Ejemplo n.º 1
0
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
Ejemplo n.º 2
0
        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 = "")
Ejemplo n.º 3
0
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'"]
Ejemplo n.º 5
0
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