InDetFlags.doV0VtxNtuple = False
InDetFlags.doRefit = doRefitTracks
InDetFlags.doLowBetaFinder = False
InDetFlags.doPrintConfigurables.set_Value_and_Lock(True)
InDetFlags.doStatistics = True
##### TEST FOR USING PRIMARY VERTEXING TO FIND DVs
##InDetFlags.primaryVertexCutSetup="StartUp"
##InDetFlags.useBeamConstraint.set_Value_and_Lock(False)
### END TEST
from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
#from ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
#print "SiSpacePointsSeedMaker",SiSpacePointsSeedMaker
InDetNewTrackingCutsLowPt = ConfiguredNewTrackingCuts("LowPt")
#InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minPT= 100. * Units.GeV
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minPT = 500. * Units.MeV
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxPT = 1. * Units.TeV
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxEta = 5
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxPrimaryImpact = 300.0 * Units.mm
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxZImpact = 1500.0 * Units.mm
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxSecondaryImpact = 300.0 * Units.mm
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minSecondaryPt = 500.0 * Units.MeV
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minClusters = 7
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minSiNotShared = 5
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxShared = 2
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxTracksPerSharedPRD = 2
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__seedFilterLevel = 1
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minPixel = 0
rec.doMuonCombined.set_Value_and_Lock(False)
rec.doEgamma.set_Value_and_Lock(False)
rec.doTau.set_Value_and_Lock(False)
from RecExConfig.RecAlgsFlags import recAlgs
recAlgs.doEgammaBremReco.set_Value_and_Lock(False)
from TrackD3PDMaker.VertexD3PDAnalysisFlags import VertexD3PDAnalysisFlags
VertexD3PDAnalysisFlags.useBeamspot.set_Value_and_Lock(True)
#Change tracking reconstruction cuts (Increase maximum z of tracking accpetance)
from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
InDetNewTrackingCuts = ConfiguredNewTrackingCuts("Offline")
InDetNewTrackingCuts._ConfiguredNewTrackingCuts__maxZImpact = 1700.
InDetNewTrackingCuts.printInfo()
InDetNewTrackingCutsLowPt = ConfiguredNewTrackingCuts("LowPt")
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxZImpact = 1700.
InDetNewTrackingCutsLowPt.printInfo()
#Change vertex reconstruction cuts (Increase maximum z of tracking accpetance)
from InDetRecExample.ConfiguredVertexingCuts import ConfiguredVertexingCuts
InDetPrimaryVertexingCuts = ConfiguredVertexingCuts("LowPt")
InDetPrimaryVertexingCuts._ConfiguredVertexingCuts__z0Max = 1700.
InDetPrimaryVertexingCuts._ConfiguredVertexingCuts__IPz0Max = 1700.
InDetPrimaryVertexingCuts._ConfiguredVertexingCuts__sigIPd0Max = 5.0
InDetPrimaryVertexingCuts.printInfo()
#
# ----------- import, lock and print InDetKeys
#
from InDetRecExample.InDetKeys import InDetKeys
InDetKeys.lockAllExceptAlias()
print "Printing InDetKeys"
InDetKeys.print_JobProperties()
#
# --- setup of cut values for NewTracking
#
if (not 'InDetNewTrackingCuts' in dir()):
print "InDetRec_jobOptions: InDetNewTrackingCuts not set before - import them now"
from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
if InDetFlags.doCosmics():
InDetNewTrackingCuts = ConfiguredNewTrackingCuts("Cosmics")
elif InDetFlags.doHeavyIon():
InDetNewTrackingCuts = ConfiguredNewTrackingCuts("HeavyIon")
elif InDetFlags.doSLHC():
InDetNewTrackingCuts = ConfiguredNewTrackingCuts("SLHC")
elif InDetFlags.doMinBias():
InDetNewTrackingCuts = ConfiguredNewTrackingCuts("MinBias")
else:
InDetNewTrackingCuts = ConfiguredNewTrackingCuts("Offline")
InDetNewTrackingCuts.printInfo()
# --- backward compatible
InDetCutValues = InDetNewTrackingCuts
#
# --- setup of cut values for Primary vertexing
#
SpacePointsOverlapName=InDetKeys.OverlapSpacePoints(),
ProcessPixels=DetFlags.haveRIO.pixel_on(),
ProcessSCTs=DetFlags.haveRIO.SCT_on(),
ProcessOverlaps=DetFlags.haveRIO.pixel_on() and DetFlags.haveRIO.SCT_on(),
OverrideBeamSpot=not doBeamSpot)
if numThreads >= 2:
InDetSiTrackerSpacePointFinder.Cardinality = numThreads
topSequence += InDetSiTrackerSpacePointFinder
# Set up ConfiguredNewTrackingCuts
# Taken from InDetRecExample/share/InDetRec_jobOptions.py
from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
NewTrackingCuts = None
if doPixel and doSCT:
NewTrackingCuts = ConfiguredNewTrackingCuts("Offline")
elif doPixel:
NewTrackingCuts = ConfiguredNewTrackingCuts("Pixel")
elif doSCT:
NewTrackingCuts = ConfiguredNewTrackingCuts("SCT")
# Set up InDet__InDetPRD_AssociationToolGangedPixels (public)
if (NewTrackingCuts.mode() == "LowPt" or NewTrackingCuts.mode() == "VeryLowPt"
or NewTrackingCuts.mode() == "LargeD0"
or NewTrackingCuts.mode() == "LowPtLargeD0"
or NewTrackingCuts.mode() == "BeamGas"
or NewTrackingCuts.mode() == "ForwardTracks"
or NewTrackingCuts.mode() == "ForwardSLHCTracks"
or NewTrackingCuts.mode() == "Disappearing"
or NewTrackingCuts.mode() == "VeryForwardSLHCTracks"
or NewTrackingCuts.mode() == "SLHCConversionFinding"):
#InDetFlags.doiPatRec = False
#
# ----------- lock and printout the setup
#
InDetFlags.init()
InDetFlags.printInfo()
InDetFlags.print_JobProperties()
if not 'InDetNewTrackingCuts' in dir():
#
# --- setup physics cuts with default values
#
print "InDetRec_jobOptions: InDetCutValues not set - setting to defaults"
from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
InDetNewTrackingCuts = ConfiguredNewTrackingCuts("Offline")
from InDetRecExample.ConfiguredVertexingCuts import ConfiguredVertexingCuts
InDetPrimaryVertexingCuts = ConfiguredVertexingCuts("NoLumi")
#
# ----------- printout the setup
#
InDetNewTrackingCuts.printInfo()
if not 'InDetKeys' in dir():
#
# --- setup StoreGate keys (JobProperties!)
#
print "InDetRec_jobOptions: InDetKeys not set - setting to defaults"
from InDetRecExample.InDetKeys import InDetKeys
InDetKeys.lock_JobProperties()
IDDerivationSequence += xAOD_TRT_PrepDataToxAOD
if (printIdTrkDxAODConf):
print(xAOD_TRT_PrepDataToxAOD)
print(xAOD_TRT_PrepDataToxAOD.properties())
# to store dEdx info
from TRT_ElectronPidTools.TRT_ElectronPidToolsConf import TRT_ToT_dEdx
TRT_dEdx_Tool = TRT_ToT_dEdx(name="TRT_ToT_dEdx")
from InDetRecExample.TrackingCommon import getInDetTRT_LocalOccupancy
TRT_dEdx_Tool.TRT_LocalOccupancyTool = getInDetTRT_LocalOccupancy()
ToolSvc += TRT_dEdx_Tool
TrackCollectionKeys = []
from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
InDetNewTrackingCutsPixel = ConfiguredNewTrackingCuts("Pixel")
print(InDetNewTrackingCutsPixel)
if not hasattr(topSequence, 'InDetTrackCollectionMerger'):
from InDetRecExample import TrackingCommon
from InDetTrackPRD_Association.InDetTrackPRD_AssociationConf import InDet__InDetTrackPRD_Association
InDetPRD_Association = InDet__InDetTrackPRD_Association(
name='InDetPRD_Association' +
InDetNewTrackingCutsPixel.extension(),
AssociationTool=TrackingCommon.
getInDetPRDtoTrackMapToolGangedPixels(),
TracksName=['Tracks'],
AssociationMapName="PRDtoTrackMap" + InDetKeys.UnslimmedTracks(),
OutputLevel=INFO)
topSequence += InDetPRD_Association
print(InDetPRD_Association)
InDetFlags.doSctTrkNtuple = False
InDetFlags.doTrtTrkNtuple = False
InDetFlags.doPixelClusterNtuple = False
InDetFlags.doSctClusterNtuple = False
InDetFlags.doTrtDriftCircleNtuple = False
InDetFlags.doVtxNtuple = False
InDetFlags.doConvVtxNtuple = False
InDetFlags.doV0VtxNtuple = False
InDetFlags.doRefit = doRefitTracks
InDetFlags.doLowBetaFinder = False
InDetFlags.doPrintConfigurables.set_Value_and_Lock(True)
InDetFlags.doStatistics = True
from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
InDetNewTrackingCutsHighD0 = ConfiguredNewTrackingCuts("HighD0")
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__minPT = 500. * Units.MeV
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxPT = 1. * Units.TeV
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxEta = 5
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxPrimaryImpact = 300.0 * Units.mm
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxZImpact = 1500.0 * Units.mm
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxSecondaryImpact = 300.0 * Units.mm
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__minSecondaryPt = 500.0 * Units.MeV
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__minClusters = 7
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__minSiNotShared = 5
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxShared = 2
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxTracksPerSharedPRD = 2
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__seedFilterLevel = 1
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__minPixel = 0
InDetNewTrackingCutsHighD0.printInfo()
#from tauRec.tauRecFlags import jobproperties
#jobproperties.tauRecFlags.doTauRec = False
# AO (this turns off all MissingET and EFlow)
rec.doJetMissingETTag = False
recAlgs.doMissingET = False
recAlgs.doMissingETSig = False
recAlgs.doEFlow = False
recAlgs.doEFlowJet = False
recAlgs.doMuGirl = False
recAlgs.doMuidLowPt = False
recAlgs.doCaloTrkMuId = False
from MuonCombinedRecExample.MuonCombinedRecFlags import muonCombinedRecFlags
muonCombinedRecFlags.doMuGirlLowBeta = False
ptCutVal = 0.5
if jobproperties.HIRecExampleFlags.ptCutOn():
ptCutVal = jobproperties.HIRecExampleFlags.ptCut()
print("ptCutOn: Modify ptCutVal to ", ptCutVal, " GeV")
if (rec.triggerStream() == "HardProbes") or (globalflags.DataSource
== "geant4"):
print("Modify ConfiguredNewTrackingCuts__minPT to ", ptCutVal,
" GeV")
from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
InDetNewTrackingCuts = ConfiguredNewTrackingCuts("HeavyIon")
InDetNewTrackingCuts._ConfiguredNewTrackingCuts__minPT = ptCutVal * Units.GeV
InDetNewTrackingCuts.printInfo()
############################################################################