def combatRun():
from Configurables import TbCombatBuilder
TbCombatBuilder().CombatInputFile0 = "/afs/cern.ch/user/v/vifranco/public/Run31100003_SAMPLE.txt"
TbCombatBuilder().ReadoutFormat = "RelaxD"
# TbCombatBuilder().CombatInputFile1 = "/afs/cern.ch/user/v/vifranco/cmtuser/KEPLER/KEPLER_HEAD/240614_30V_TH500_100ms_merge.dat"
seq = GaudiSequencer("Telescope")
seq.Members = [TbCombatBuilder(), TbClustering(), TbTracking()]
TbCombatBuilder().PrintFreq = 5000
TbCombatBuilder().NumberArms = 2
TbTracking().MinNClusters = 6
TbTracking().SearchRadius = 2
TbTracking().VolumeAngle = 0.0015
TbTracking().SearchPlanes = [3,5]
TbTracking().SearchVolume = "diabolo"
TbTracking().MaxChi2 = 2000000.
TbTracking().ViewerEventNum = 100
TbTracking().ViewerTimeLow = 0
TbTracking().ViewerTimeUp = 2000000
#TbTracking().ViewerOutput = True
TbTracking().CombatRun = True
seq = GaudiSequencer("Monitoring")
from Configurables import TbHitMonitor, TbClusterPlots, TbTrackPlots
seq.Members = [TbHitMonitor(), TbClusterPlots(),TbTrackPlots()]
TbTrackPlots().ParametersResidualsXY = ("", -1.0, 1.0, 200)
TbTrackPlots().ParametersXY = ("", -2, 16, 200)
TbClusterPlots().ParametersXY = ("", -2, 16, 200)
def configureMonitoring(self, seq):
from Configurables import TbHitMonitor, TbClusterPlots, TbTrackPlots
if self.getProp("Tracking") == True:
seq.Members += [TbHitMonitor(), TbClusterPlots(), TbTrackPlots()]
from Configurables import TbTriggerMonitor
seq.Members += [TbTriggerMonitor()]
from Configurables import TbDUTMonitor
seq.Members += [TbDUTMonitor()]
else:
seq.Members += [TbHitMonitor(), TbClusterPlots()]
def egRun():
from Configurables import TbEventBuilder, TbTrackPlots, TbCalibration, TbVisualiserOutput
from Configurables import TbClustering, TbClusterPlots, TbSimpleTracking, TbHitMonitor
from Configurables import TbClusterAssociator, TbEfficiency, TbDUTMonitor
TbEventBuilder().PrintFreq = 10
TbSimpleTracking().TimeWindow = 6
TbSimpleTracking().MaxDistance = 0
TbSimpleTracking().MinPlanes = 7
TbSimpleTracking().MaskedPlanes = [DUT_id, 4]
TbSimpleTracking().MaxOpeningAngle = 0.005
TbSimpleTracking().RecheckTrack = True
TbSimpleTracking().ChargeCutLow = 0
TbSimpleTracking().DoOccupancyCut = True
TbSimpleTracking().MaxClusterSize = 20
TbSimpleTracking().MaxOccupancy = 7 # not inclusive
TbTrackPlots().MaskedPlanes = [DUT_id, 4]
TbClusterAssociator().DUTs = [DUT_id]
TbDUTMonitor().DUTs = [DUT_id]
TbClusterAssociator().TimeWindow = 10
TbClusterAssociator().XWindow = allowance
TbEfficiency().CheckHitDUT = True
TbEfficiency().CheckHitAlivePixel = True
TbEfficiency().DUT = DUT_id
TbEfficiency().TakeDeadPixelsFromFile = True
TbEfficiency().nTotalTracks =500000
TbEfficiency().MaxChi = 2
seq = GaudiSequencer("Telescope")
seq.Members = [TbEventBuilder(), TbClustering(), TbSimpleTracking(), TbClusterAssociator(), TbEfficiency()]
#seq.Members = [TbEventBuilder()]
seq = GaudiSequencer("Monitoring")
from Configurables import TbHitMonitor, TbClusterPlots, TbTrackPlots
#seq.Members = [TbHitMonitor(), TbClusterPlots(), TbTrackPlots(), TbDUTMonitor(), TbVisualiserOutput()]
seq.Members = [TbDUTMonitor()]
def sim():
seq = GaudiSequencer("Telescope")
# TestMC options.
from Configurables import TbTestMC
TbTestMC().doMisAlign = False
TbTestMC().NTracks = 200
TbTestMC().RunDuration = 10000
TbTestMC().NNoiseClusters = 0
TbTestMC().HitTimeJitter = 20
TbTestMC().ClusterPosnError = 0.001
TbTestMC().ClusterADCMPV = 300.0 # Ish.
TbTestMC().ChargeSharingWidth = 0.033 # Best not to change much.
TbTestMC().ThresholdCut = 20
TbTestMC().ForceEfficiency = True
#TbTestMC().InitAlignment = "Alignment_perfect.dat"
#TbTestMC().InitAlignment = Alignment_raw.dat
# TbClustering options.
from Configurables import TbClustering
TbClustering().TimeWindow = 200
# TbTracking options.
from Configurables import TbTracking
TbTracking().TimeWindow = 300
TbTracking().SearchRadius = 1
TbTracking().MinNClusters = 7
TbTracking().Monitoring = True
TbTracking().SearchVolume = "cylinder" # {cylinder, diabolo}
# TbTracking speed options.
TbTracking().nComboCut = 1000 # O(10) for speed.
TbTracking(
).SearchVolumeFillAlgorithm = "adap_seq" # {seq, adap_seq}. Adap faster.
seq.Members = [TbTestMC(), TbClustering(), TbTracking()]
seq = GaudiSequencer("Monitoring")
from Configurables import TbClusterPlots, TbTrackPlots
seq.Members = [TbClusterPlots(), TbTrackPlots()]
from Gaudi.Configuration import * trackFitToolName = "TbKalmanTrackFit" from Configurables import TbTracking TbTracking().TrackFitTool = trackFitToolName from Configurables import TbTrackPlots TbTrackPlots().TrackFitTool = trackFitToolName from Configurables import TbSimpleTracking #TbSimpleTracking().TrackFitTool = trackFitToolName from Configurables import TbVertexTracking #TbVertexTracking().TrackFitTool = trackFitToolName from Configurables import TbAlignment # TbAlignment().TrackFitTool = trackFitToolName
from Configurables import Kepler
from Configurables import TbClusterAssociator
TbClusterAssociator().DUTs = [4]
TbClusterAssociator().XWindow = 0.1
TbClusterAssociator().TimeWindow = 35
from Configurables import TbEventBuilder
TbEventBuilder().PrintHeader = True
TbEventBuilder().MinPlanesWithHits = 8
TbEventBuilder().PrintFreq = 50
TbEventBuilder().Monitoring = True
TbEventBuilder().ForceCaching = False
from Configurables import TbClustering
TbClustering().SearchDist = 2
from Configurables import TbSimpleTracking
Kepler().UseSimpleTracking = True
TbSimpleTracking().MinPlanes = 8
TbSimpleTracking().MaskedPlanes = [4]
from Configurables import TbTrackPlots
TbTrackPlots().ParametersSlope = ("", -0.01, 0.01, 1000)
TbTrackPlots().ParametersResidualsXY = ("", -0.055, 0.055, 500)
TbTrackPlots().ParametersXY = ("", 0, 15, 500)
from Configurables import TbTriggerMonitor
TbTriggerMonitor().OutputLevel = ERROR
from Gaudi.Configuration import *
from Configurables import Kepler
Kepler().InputFiles = ["eos/lhcb/testbeam/velo/timepix3/July2014/RawData/Run1062/"]
Kepler().PixelConfigFile = "eos/lhcb/testbeam/velo/timepix3/July2014/RootFiles/Run1062/Conditions/PixelMask.dat"
Kepler().EvtMax = 100
Kepler().WriteTuples = False
Kepler().Monitoring = False
Kepler().Alignment = True
Kepler().AlignmentFile = "Alignment.dat"
from Configurables import TbTrackPlots
TbTrackPlots().ParametersResidualsXY = ("",-0.1,0.1,2000)
from Configurables import TbAlignment, TbTracking
TbAlignment().PrintConfiguration = True
TbAlignment().FitStrategy = 1
TbAlignment().ReferencePlane = 3
TbAlignment().MaskedPlanes = TbTracking().MaskedPlanes
TbAlignment().OutputLevel = INFO