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