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
#-------------------------------------------------------------
from InDetPhysValMonitoring.InDetPhysValMonitoringConf import InDetPhysValMonitoringTool
InDetPhysValMonitoringTool = InDetPhysValMonitoringTool("InDetPhysValMonitoringTool")
InDetPhysValMonitoringTool.useTrackSelection = True
InDetPhysValMonitoringTool.TrackSelectionTool = InDetTrackSelectorTool
#InDetPhysValMonitoringTool.TruthSelectionTool = TrackTruthSelectionTool
InDetPhysValMonitoringTool.TruthParticleContainerName = "TruthParticles"
#InDetPhysValMonitoringTool.PileupSwitch = "HardScatter"
InDetPhysValMonitoringTool.OutputLevel = INFO
ToolSvc += InDetPhysValMonitoringTool
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
ToolSvc += tool1
monMan.AthenaMonTools += [tool1]
'''
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()
tool1.TruthSelectionTool = AthTruthSelectionTool
tool1.useTrackSelection = usingTrackSelection
tool1.TrackSelectionTool = InDetTrackSelectorTool
tool1.FillTrackInJetPlots = True
print tool1