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")