def _preparePedestalRun(self):
# dummy aligment just to run Kepler
Kepler().PixelConfigFile = ["eos/lhcb/testbeam/velo/timepix3/July2014/RootFiles/Run1236/Conditions/PixelConfig.dat"]
Kepler().AlignmentFile = "eos/lhcb/testbeam/velo/timepix3/Oct2014/RootFiles/Run2176/Conditions/Alignment2176mille.dat"
Kepler().InputFiles = ['eos/lhcb/testbeam/velo/timepix3/Oct2014/RawData/Run2176/']
Kepler().EvtMax = self.eventMax
Kepler().HistogramFile="MambaPedestal.root" # should be set more correctly
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
from Gaudi.Configuration import * from Configurables import Kepler from Configurables import TbEventBuilder TbEventBuilder().MinPlanesWithHits = 2 TbEventBuilder().PrintFreq = 100 TbEventBuilder().Monitoring = True Kepler().EvtMax = 5000 Kepler( ).AlignmentFile = "eos/lhcb/testbeam/velo/timepix3/Dec2014/RootFiles/Run4012/Conditions/Alignment4012.dat" from Configurables import TbTracking TbTracking().PrintConfiguration = True TbTracking().MinNClusters = 3 TbTracking().SearchRadius = 0.5 TbTracking().VolumeAngle = 0.2 TbTracking().TimeWindow = 75 TbTracking().Monitoring = False TbTracking().SearchPlanes = [2, 3] TbTracking().MaskedPlanes = [] from Configurables import TbTrackPlots
from Gaudi.Configuration import *
from Configurables import Kepler
# Add TbAlignment to the Telescope sequence.
Kepler().Alignment = True
# Set the number of events to process.
Kepler().EvtMax = 100
from Configurables import TbEventBuilder
# Skip noise events.
TbEventBuilder().MinPlanesWithHits = 5
from Configurables import TbClusterPlots
# Set the reference plane.
TbClusterPlots().ReferencePlane = 3
# Widen the range of the difference histograms if needed.
TbClusterPlots().ParametersDifferenceXY = ('', -10., 10., 200)
from Configurables import TbAlignment
TbAlignment().AlignmentTechnique = "survey"
from Gaudi.Configuration import *
from Configurables import Kepler
Kepler().Alignment = True
from Configurables import TbAlignment
from Configurables import TbMillepede, TbAlignmentMinuit2, TbAlignmentMinuit1
Kepler().addAlignment(
TbAlignmentMinuit1(DOFs=[1, 1, 0, 0, 0, 1],
MaxChi2=9999999999999999999,
ReferencePlane=3))
Kepler().addAlignment(
TbMillepede(DOFs=[1, 1, 0, 1, 1, 1],
MaxChi2=200,
ResCutInit=0.5,
ResCut=0.1))
Kepler().addAlignment(
TbMillepede(DOFs=[1, 1, 1, 1, 1, 1],
MaxChi2=50,
ResCutInit=1.0,
ResCut=0.2))
Kepler().addAlignment(
TbMillepede(Monitoring=True,
DOFs=[1, 1, 1, 1, 1, 1],
MaxChi2=15,
ResCutInit=1.0,
ResCut=0.2))
TbAlignment().PrintConfiguration = True
TbAlignment().NTracks = 8000
from Gaudi.Configuration import *
from Configurables import Kepler
Kepler().Alignment = True
from Configurables import TbAlignment
from Configurables import TbMillepede, TbAlignmentMinuit2, TbAlignmentMinuit1
Kepler().addAlignment(
TbAlignmentMinuit1(MaskedPlanes=[4],
DOFs=[1, 1, 0, 0, 0, 1],
MaxChi2=9999999999999999999,
ReferencePlane=3))
Kepler().addAlignment(
TbMillepede(MaskedPlanes=[4],
DOFs=[1, 1, 0, 1, 1, 1],
MaxChi2=2000,
ResCutInit=0.5,
ResCut=0.1))
Kepler().addAlignment(
TbMillepede(MaskedPlanes=[4],
DOFs=[1, 1, 0, 1, 1, 1],
MaxChi2=15,
ResCutInit=1.0,
ResCut=0.2))
Kepler().addAlignment(
TbMillepede(Monitoring=True,
MaskedPlanes=[4],
DOFs=[1, 1, 1, 1, 1, 1],
MaxChi2=15,
ResCutInit=1.0,
from Gaudi.Configuration import *
from Configurables import Kepler
# Set the path to the directory/files to be processed
Kepler().InputFiles = [
"/afs/cern.ch/work/t/tevans/public/Kepler_data_1MHz_120s_trigger/"
]
Kepler().EvtMax = 4000
def sim():
telseq = GaudiSequencer("Telescope")
from Configurables import TbTrackFitter
telseq.Members = [TbTrackFitter()]
# scattering estimate (theta0 squared?)
TbTrackFitter().scat2 = 1.e-9
#TbTrackFitter().scat2 = 0
TbTrackFitter().direction = -1
# divergence of true track
TbTrackFitter().theta0 = 1.e-4
#TbTrackFitter().hiterror2 = 9.0e-6
#TbTrackFitter().theta0 = 0
monseq = GaudiSequencer("Monitoring")
monseq.Members = []
appendPostConfigAction(sim)