Beispiel #1
0
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)
Beispiel #2
0
 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()]
Beispiel #3
0
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()]
Beispiel #4
0
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()]
Beispiel #5
0
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
Beispiel #6
0
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
Beispiel #7
0
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