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