#LHCbApp().DDDBtag = "head-20090330"
#LHCbApp().CondDBtag = "head-20090402"
LHCbApp().DDDBtag = "head-20100119"
LHCbApp().CondDBtag = "head-20091112"
# Latest cosmic run, with CALO, OT and (!!) RICH2 (35569 events)
from Configurables import CountingPrescaler
CountingPrescaler("EscherPrescaler").OutputLevel = 5
# do not do the TES check because these data do not have the required lists
GaudiSequencer("InitReprocSeq").Enable = False
Escher().DataType = "2009"
Escher().Simulation = False
Escher().SpecialData += ["fieldOff", "veloOpen", "earlyData"]
Escher().InputType = "DST"
Escher().EvtMax = 100
Escher().DatasetName = "63596_1"
TAlignment().TrackLocation = "Event/Alignment/OTTracks"
ATrackSelector().UseWeights = True
# add 2009 collision data files
EventSelector().Input = [
"DATAFILE='file:/lhcb/data/2009/DST/00005843/0000/00005843_00000001_1.dst' TYP='POOL_ROOTTREE' OPT='READ'"
# "DATAFILE='file:/data/data/2009/DST/00005714/0000/00005714_00000001_1.dst' SVC='LHCb::MDFSelector'",
# "DATAFILE='file:/media/disk-2/data/2009/DST/00005738/0000/00005738_00000001_1.dst' SVC='LHCb::MDFSelector'"
]
from Configurables import (ATrackSelector) ATrackSelector().MinPCut = 0.0 ATrackSelector().MinPtCut = 0.0 ATrackSelector().MinChi2Cut = 0.0 ATrackSelector().MaxChi2Cut = 15.0 ATrackSelector().Charge = 0 ATrackSelector().MaxPCut = -1 ATrackSelector().MaxPtCut = -1 ATrackSelector().MinITHitCut = 0 ATrackSelector().MinOTHitCut = 15 ATrackSelector().MinTTHitCut = 0 ATrackSelector().MinEnergyCut = 0.0 ATrackSelector().OutputLevel = 1
#otCalib = CondDBAccessSvc( 'OTCalib' ) #otCalib.ConnectionString ='sqlite_file:/data/work/DataBases/Collision09_OT_QuarterT0s_071209.db/LHCBCOND' #CondDB().addLayer( otCalib ) # #from Configurables import UpdateManagerSvc #UpdateManagerSvc().ConditionsOverride += ["Conditions/Online/LHCb/Magnet/Measured := double Current = 5849.9936523438 ; int Polarity = -1;"] # #aligndb = '/data/work/DataBases/TTLayerITLayerOTModules_CollisionsDecember2009.db' #from Configurables import ( CondDB, CondDBAccessSvc ) #alignCond = CondDBAccessSvc( 'AlignCond' ) #alignCond.ConnectionString = 'sqlite_file:' + aligndb + '/LHCBCOND' #CondDB().addLayer( alignCond ) from Configurables import (ATrackSelector, Escher) ATrackSelector().MinPCut = 0.0 # value is in MeV ATrackSelector().MinPtCut = 0.0 ATrackSelector().MinChi2Cut = 0.0 ATrackSelector().MaxChi2Cut = 15.0 # chi2perdof ATrackSelector().Charge = 0 ATrackSelector().MaxPCut = -1 # value = -1 -> no p info available (Boff) ATrackSelector().MaxPtCut = -1 ATrackSelector().MinITHitCut = 0 ATrackSelector().MinOTHitCut = 15 ATrackSelector().MinTTHitCut = 0 ATrackSelector().MinEnergyCut = 0.0 Escher().Detectors = ["OT"] Escher().EvtMax = 100 Escher().AlignmentLevel = "layers" Escher().Millepede = True
TAConfig().IT_ladder = False TAConfig().MUON_system = False TAConfig().MUON_stations = False TAConfig().MUON_chambers = False TAConfig().TT_system = False TAConfig().TT_station = False TAConfig().TT_layer = False TAConfig().IT_alignment_condition_location = "" TAConfig().TT_alignment_condition_location = "" TAConfig().OT_alignment_condition_location = "" WriteAlignmentConditions().topElement = "" WriteAlignmentConditions().footer = "</DDDB>" WriteAlignmentConditions().startTag = "<condition" WriteAlignmentConditions().outputFile = "alignment.xml" WriteAlignmentConditions().depths = [] WriteAlignmentConditions().precision = 16 ATrackSelector().MinPCut = 0.0 ATrackSelector().MinPtCut = 0.0 ATrackSelector().MinEnergyCut = 0.0 ATrackSelector().MinChi2Cut = 0.0 ATrackSelector().Charge = 0 ATrackSelector().MaxPCut = [] ATrackSelector().MaxPtCut = [] ATrackSelector().MaxChi2Cut = [] ATrackSelector().MaxHitCut = []
Escher().PrintFreq = 1000 Escher().AlignmentLevel = "layers" Escher().Millepede = True Escher().Kalman = False Escher().Incident = "GlobalMPedeFit" Escher().TrackContainer = "Rec/Track/Best" Escher().skipBigCluster = True Escher().OutputLevel = 3 AlignTrTools().nTrackModelParameters = 4 AlignTrTools().Chi2Scale = 10 AlignTrTools().minChi2 = 5 AlignTrTools().Outlier = 3 AlignTrTools().Sim = True # simulated data or real data ATrackSelector().UseWeights = True ATrackSelector().UniformCutOff = 6 ATrackSelector().MinEnergyCut = 0.0 ATrackSelector().MinPCut = 0.0 ATrackSelector().MinPtCut = 0.0 ATrackSelector().MinChi2Cut = 0.0 ATrackSelector().MaxChi2Cut = 15.0 ATrackSelector().Charge = 0 ATrackSelector().MaxPCut = -1 ATrackSelector().MaxPtCut = -1 ATrackSelector().MinITHitCut = 0 ATrackSelector().MinOTHitCut = 15 ATrackSelector().MinTTHitCut = 0 #ATrackSelector().YCutMin_at_T1 = -2500.0 #ATrackSelector().YCutMax_at_T1 = -830.0 ATrackSelector().OutputLevel = 3
GaudiSequencer("InitReprocSeq").Enable = False
Escher().DataType = "2009"
Escher().Simulation = False
Escher().SpecialData += ["fieldOff", "cosmics"]
Escher().skipBigCluster = True
Escher().InputType = 'DST'
Escher().Detectors = ["OT"]
Escher().EvtMax = 1000
Escher().AlignmentLevel = "layers"
Escher().Millepede = True
Escher().Kalman = False
Escher().Incident = "GlobalMPedeFit"
Escher().DatasetName = "Cosmics"
ATrackSelector().UseWeights = False
ApplicationMgr().ExtSvc += ["NTupleSvc"]
tupleFile = "2009-Cosmic.root"
tuple = "FILE1 DATAFILE='" + tupleFile + "' TYP='ROOT' OPT='NEW'"
NTupleSvc().Output = [tuple]
data = [
'58739', '58744', '59051', '59052', '59054', '59075', '59077', '59081',
'59088', '59092', '59093', '59096', '59097', '59195', '59198', '59306',
'59307', '59565', '59568'
]
EventSelector().Input = []
for d in data:
name = "DATAFILE='" + d + "' TYP='POOL_ROOTTREE' OPT='READ'"
def TrackSelection(self):
from Configurables import ATrackSelector
ATrackSelector().MinPCut = 0.0
ATrackSelector().MinPtCut = 0.0
ATrackSelector().MinChi2Cut = 0.0
ATrackSelector().MaxChi2Cut = 15.0
ATrackSelector().Charge = 0
ATrackSelector().MaxPCut = -1
ATrackSelector().MaxPtCut = -1
ATrackSelector().MinITHitCut = 0
ATrackSelector().MinOTHitCut = 15
ATrackSelector().MinTTHitCut = 0
ATrackSelector().MinEnergyCut = 0.0
ATrackSelector().OutputLevel = 1
ATrackSelector().YCutMax_at_T1 = -3000.0
ATrackSelector().YCutMin_at_T1 = 3000.0