minNInnermostLayerHits = InDetPrimaryVertexingCuts.nHitInnermostLayer(),
minNPixelHits = InDetPrimaryVertexingCuts.nHitPix(),
maxNPixelHoles = InDetPrimaryVertexingCuts.nHolesPix(),
minNSctHits = InDetPrimaryVertexingCuts.nHitSct(),
minNTrtHits = InDetPrimaryVertexingCuts.nHitTrt(),
minNSiHits = InDetPrimaryVertexingCuts.nHitSi(),
TrackSummaryTool = InDetTrackSummaryTool,
Extrapolator = InDetExtrapolator)
ToolSvc += InDetTrackSelectorTool
'''
#This section should control TTST 7-12-16
mode = "FWD" #Set this to "Back" for backtracking
from InDetPhysValMonitoring.InDetPhysValMonitoringConf import AthTruthSelectionTool
truthSelection = AthTruthSelectionTool()
if mode == "Back":
# max prod. vertex radius for secondaries [mm]
# < 0 corresponds to : do not require decay before pixel
truthSelection.maxProdVertRadius = -999.9
truthSelection.maxBarcode = -1
ToolSvc += truthSelection
from InDetPhysValMonitoring.InDetPhysValMonitoringConf import InDetPhysValMonitoringTool
tool1 = InDetPhysValMonitoringTool()
tool1.TruthSelectionTool = truthSelection
tool1.useTrackSelection = False
#tool1.TrackSelectionTool=InDetTrackSelectorTool
tool1.FillTrackInJetPlots = False
minNTrtHits = InDetPrimaryVertexingCuts.nHitTrt(),
minNSiHits = InDetPrimaryVertexingCuts.nHitSi(),
TrackSummaryTool = InDetTrackSummaryTool,
Extrapolator = InDetExtrapolator)
ToolSvc += InDetTrackSelectorTool
'''
from InDetTrackSelectionTool.InDetTrackSelectionToolConf import InDet__InDetTrackSelectionTool
InDetTrackSelectorTool = InDet__InDetTrackSelectionTool()
ToolSvc += InDetTrackSelectorTool
InDetTrackSelectorTool.CutLevel = "TightPrimary"
#This section should control TTST 7-12-16
from InDetPhysValMonitoring.InDetPhysValMonitoringConf import AthTruthSelectionTool
AthTruthSelectionTool = AthTruthSelectionTool()
if mode == "Back":
# max prod. vertex radius for secondaries [mm]
AthTruthSelectionTool.minPt = 5000
AthTruthSelectionTool.maxProdVertRadius = 4000
AthTruthSelectionTool.maxBarcode = -1
AthTruthSelectionTool.hasNoGrandparent = True
AthTruthSelectionTool.poselectronfromgamma = True
os.environ["BACKTRACKDEBUG"] = "1"
print AthTruthSelectionTool
#ToolSvc += AthTruthSelectionTool
from InDetPhysValMonitoring.InDetPhysValMonitoringConf import InDetPhysValMonitoringTool
tool1 = InDetPhysValMonitoringTool()
InDetTrackSelectorTool.minPt = 400 # Mev InDetTrackSelectorTool.maxD0 = 1 # mm InDetTrackSelectorTool.maxZ0 = 150 # mm InDetTrackSelectorTool.minNSiHits = 9 # Pixel + SCT InDetTrackSelectorTool.maxNPixelHoles = 2 # Pixel only #eta dependant hit cut below #InDetTrackSelectorTool.vecEtaCutoffsForSiHitsCut = [0,1.0,1.2,1.8,2.2] #InDetTrackSelectorTool.vecMinNSiHitsAboveEta = [11,11,11,13,10] ToolSvc += InDetTrackSelectorTool print "Set Up InDetTrackSelectorTool" #------------------------------------------------------------- # Set up truth selection tool #------------------------------------------------------------- from InDetPhysValMonitoring.InDetPhysValMonitoringConf import AthTruthSelectionTool TrackTruthSelectionTool = AthTruthSelectionTool() # todo: manually adapt this acc. to the fiducial acceptance of your detector TrackTruthSelectionTool.maxEta = 4.0 TrackTruthSelectionTool.maxPt = -1 TrackTruthSelectionTool.minPt = 900 # default 400 MeV TrackTruthSelectionTool.maxBarcode = int(200e3) TrackTruthSelectionTool.pdgId = -1 TrackTruthSelectionTool.requireCharged = True TrackTruthSelectionTool.requireStatus1 = True TrackTruthSelectionTool.maxProdVertRadius = 260. #max prod. vertex radius of secondaries [mm] TrackTruthSelectionTool.OutputLevel = INFO ToolSvc += TrackTruthSelectionTool #------------------------------------------------------------- # Set up InDetPhysValMonitoring tool #-------------------------------------------------------------
monMan.Run = 1 monMan.LumiBlock = 1 monMan.FileKey = "M_output" topSequence += monMan from InDetTrackSelectionTool.InDetTrackSelectionToolConf import InDet__InDetTrackSelectionTool InDetTrackSelectorTool = InDet__InDetTrackSelectionTool() ToolSvc += InDetTrackSelectorTool InDetTrackSelectorTool.CutLevel = "TightPrimary" #This section should control TTST 7-12-16 from InDetPhysValMonitoring.InDetPhysValMonitoringConf import AthTruthSelectionTool AthTruthSelectionTool = AthTruthSelectionTool() print AthTruthSelectionTool ToolSvc += AthTruthSelectionTool from InDetPhysValMonitoring.InDetPhysValMonitoringConf import InDetPhysValMonitoringTool tool1 = InDetPhysValMonitoringTool() tool1.TruthSelectionTool = AthTruthSelectionTool tool1.useTrackSelection = False #tool1.TrackSelectionTool=InDetTrackSelectorTool tool1.FillTrackInJetPlots = False print tool1 ToolSvc += tool1 # tool2 = InDetPhysValMonitoringTool(name="TightPrimary")