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] ''' from InDetTrackHoleSearch.InDetTrackHoleSearchConf import InDet__InDetTrackHoleSearchTool InDetHoleSearchTool = InDet__InDetTrackHoleSearchTool(name = "InDetHoleSearchTool", Extrapolator = InDetExtrapolator, usePixel = True, useSCT= True, CountDeadModulesAfterLastHit = True) ToolSvc += InDetHoleSearchTool print InDetHoleSearchTool ''' from GaudiSvc.GaudiSvcConf import THistSvc
FileKey="PhysValMon",
ManualDataTypeSetup=True,
DataType="userDefined", # use this for collision data for now
Environment="user",
ManualRunLBSetup=True,
Run=1,
LumiBlock=1)
topSequence += InDetPhysValMonManager
if (InDetFlags.doPrintConfigurables()):
print(InDetPhysValMonManager)
from InDetPhysValMonitoring.InDetPhysValMonitoringConf import InDetPhysValMonitoringTool
if InDetFlags.doDBMstandalone():
InDetPhysValMonTool = InDetPhysValMonitoringTool(
useTrackSelection=False,
TrackParticleContainerName="InDetDBMTrackParticles",
TrackSelectionTool=None,
TruthSelectionTool=None,
jetContainerName="")
else:
InDetPhysValMonTool = InDetPhysValMonitoringTool(
useTrackSelection=True, TrackSelectionTool=InDetTrackSelectorTool)
# TrackParticleContainerName = InDetKeys.Tracks())
# InputTrackCollectionTruth = InDetKeys.TracksTruth()
# if InDetFlags.doDBM():
# InDetPhysValMonToolDBM = InDetPhysValMonitoringTool (useTrackSelection = False,
# TrackParticleContainerName = "InDetDBMTrackParticles",
# TruthParticleContainerName = "DBMTracksTruth",
# TrackSelectionTool = None,
# TruthSelectionTool = None,
# jetContainerName = "")
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
#ToolSvc += tool1
monMan.AthenaMonTools += [tool1]
from InDetBoundaryCheckTool.InDetBoundaryCheckToolConf import InDet__InDetBoundaryCheckTool
InDetBoundaryCheckTool = InDet__InDetBoundaryCheckTool(
name="InDetBoundaryCheckTool")
ToolSvc += InDetBoundaryCheckTool
from InDetTrackHoleSearch.InDetTrackHoleSearchConf import InDet__InDetTrackHoleSearchTool
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
monMan.AthenaMonTools += [InDetPhysValMonitoringTool]
# set up output file
from GaudiSvc.GaudiSvcConf import THistSvc
ServiceMgr += THistSvc()
svcMgr.THistSvc.Output += ["MyPhysVal DATAFILE='MyPhysVal.root' OPT='RECREATE'"]
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") tool2.TruthSelectionTool = AthTruthSelectionTool tool2.useTrackSelection = True tool2.TrackSelectionTool = InDetTrackSelectorTool tool2.FillTrackInJetPlots = False tool2.DirName = "TightPrimary/" ToolSvc += tool2