InDetFlags.doSctTrkNtuple     = False 
InDetFlags.doTrtTrkNtuple     = False 
InDetFlags.doPixelClusterNtuple = False
InDetFlags.doSctClusterNtuple   = False
InDetFlags.doTrtDriftCircleNtuple = False
InDetFlags.doVtxNtuple        = False 
InDetFlags.doConvVtxNtuple    = False 
InDetFlags.doV0VtxNtuple      = False 
InDetFlags.doRefit            = doRefitTracks
InDetFlags.doLowBetaFinder    = False
InDetFlags.doPrintConfigurables.set_Value_and_Lock(True)
InDetFlags.doStatistics = True

from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts

InDetNewTrackingCutsHighD0 = ConfiguredNewTrackingCuts("HighD0")
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__minPT= 500. * Units.MeV
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxPT= 1. * Units.TeV
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxEta= 5
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxPrimaryImpact=300.0 * Units.mm
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxZImpact=1500.0 * Units.mm
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxSecondaryImpact=300.0 * Units.mm
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__minSecondaryPt=500.0 * Units.MeV
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__minClusters=7
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__minSiNotShared=5
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxShared=2
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__maxTracksPerSharedPRD=2
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__seedFilterLevel=1
InDetNewTrackingCutsHighD0._ConfiguredNewTrackingCuts__minPixel=0

InDetNewTrackingCutsHighD0.printInfo()
        #
        # ----------- import, lock and print InDetKeys
        #
        from InDetRecExample.InDetKeys import InDetKeys
        InDetKeys.lockAllExceptAlias()
        print "Printing InDetKeys"
        InDetKeys.print_JobProperties()

        #
        # --- setup of cut values for NewTracking
        #
        if (not 'InDetNewTrackingCuts' in dir()):
            print "InDetRec_jobOptions: InDetNewTrackingCuts not set before - import them now"
            from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
            if InDetFlags.doCosmics():
                InDetNewTrackingCuts = ConfiguredNewTrackingCuts("Cosmics")
            elif InDetFlags.doHeavyIon():
                InDetNewTrackingCuts = ConfiguredNewTrackingCuts("HeavyIon")
            elif InDetFlags.doSLHC():
                InDetNewTrackingCuts = ConfiguredNewTrackingCuts("SLHC")
            elif InDetFlags.doMinBias():
                InDetNewTrackingCuts = ConfiguredNewTrackingCuts("MinBias")
            else:
                InDetNewTrackingCuts = ConfiguredNewTrackingCuts("Offline")
        InDetNewTrackingCuts.printInfo()
        # --- backward compatible
        InDetCutValues = InDetNewTrackingCuts

        #
        # --- setup of cut values for  Primary vertexing
        #
예제 #3
0
InDetFlags.doV0VtxNtuple = False
InDetFlags.doRefit = doRefitTracks
InDetFlags.doLowBetaFinder = False
InDetFlags.doPrintConfigurables.set_Value_and_Lock(True)
InDetFlags.doStatistics = True

##### TEST FOR USING PRIMARY VERTEXING TO FIND DVs
##InDetFlags.primaryVertexCutSetup="StartUp"
##InDetFlags.useBeamConstraint.set_Value_and_Lock(False)
### END TEST

from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
#from ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
#print "SiSpacePointsSeedMaker",SiSpacePointsSeedMaker

InDetNewTrackingCutsLowPt = ConfiguredNewTrackingCuts("LowPt")
#InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minPT= 100. * Units.GeV
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minPT = 500. * Units.MeV
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxPT = 1. * Units.TeV
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxEta = 5
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxPrimaryImpact = 300.0 * Units.mm
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxZImpact = 1500.0 * Units.mm
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxSecondaryImpact = 300.0 * Units.mm
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minSecondaryPt = 500.0 * Units.MeV
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minClusters = 7
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minSiNotShared = 5
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxShared = 2
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxTracksPerSharedPRD = 2
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__seedFilterLevel = 1
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__minPixel = 0
예제 #4
0
rec.doMuonCombined.set_Value_and_Lock(False)
rec.doEgamma.set_Value_and_Lock(False)
rec.doTau.set_Value_and_Lock(False)

from RecExConfig.RecAlgsFlags import recAlgs

recAlgs.doEgammaBremReco.set_Value_and_Lock(False)

from TrackD3PDMaker.VertexD3PDAnalysisFlags import VertexD3PDAnalysisFlags

VertexD3PDAnalysisFlags.useBeamspot.set_Value_and_Lock(True)

#Change tracking reconstruction cuts (Increase maximum z of tracking accpetance)
from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts

InDetNewTrackingCuts = ConfiguredNewTrackingCuts("Offline")
InDetNewTrackingCuts._ConfiguredNewTrackingCuts__maxZImpact = 1700.
InDetNewTrackingCuts.printInfo()

InDetNewTrackingCutsLowPt = ConfiguredNewTrackingCuts("LowPt")
InDetNewTrackingCutsLowPt._ConfiguredNewTrackingCuts__maxZImpact = 1700.
InDetNewTrackingCutsLowPt.printInfo()

#Change vertex reconstruction cuts (Increase maximum z of tracking accpetance)
from InDetRecExample.ConfiguredVertexingCuts import ConfiguredVertexingCuts

InDetPrimaryVertexingCuts = ConfiguredVertexingCuts("LowPt")
InDetPrimaryVertexingCuts._ConfiguredVertexingCuts__z0Max = 1700.
InDetPrimaryVertexingCuts._ConfiguredVertexingCuts__IPz0Max = 1700.
InDetPrimaryVertexingCuts._ConfiguredVertexingCuts__sigIPd0Max = 5.0
InDetPrimaryVertexingCuts.printInfo()
예제 #5
0
    SpacePointsOverlapName=InDetKeys.OverlapSpacePoints(),
    ProcessPixels=DetFlags.haveRIO.pixel_on(),
    ProcessSCTs=DetFlags.haveRIO.SCT_on(),
    ProcessOverlaps=DetFlags.haveRIO.pixel_on() and DetFlags.haveRIO.SCT_on(),
    OverrideBeamSpot=not doBeamSpot)
if numThreads >= 2:
    InDetSiTrackerSpacePointFinder.Cardinality = numThreads
topSequence += InDetSiTrackerSpacePointFinder

# Set up ConfiguredNewTrackingCuts
# Taken from InDetRecExample/share/InDetRec_jobOptions.py
from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts

NewTrackingCuts = None
if doPixel and doSCT:
    NewTrackingCuts = ConfiguredNewTrackingCuts("Offline")
elif doPixel:
    NewTrackingCuts = ConfiguredNewTrackingCuts("Pixel")
elif doSCT:
    NewTrackingCuts = ConfiguredNewTrackingCuts("SCT")

# Set up InDet__InDetPRD_AssociationToolGangedPixels (public)
if (NewTrackingCuts.mode() == "LowPt" or NewTrackingCuts.mode() == "VeryLowPt"
        or NewTrackingCuts.mode() == "LargeD0"
        or NewTrackingCuts.mode() == "LowPtLargeD0"
        or NewTrackingCuts.mode() == "BeamGas"
        or NewTrackingCuts.mode() == "ForwardTracks"
        or NewTrackingCuts.mode() == "ForwardSLHCTracks"
        or NewTrackingCuts.mode() == "Disappearing"
        or NewTrackingCuts.mode() == "VeryForwardSLHCTracks"
        or NewTrackingCuts.mode() == "SLHCConversionFinding"):
예제 #6
0
#InDetFlags.doiPatRec     = False

#
# ----------- lock and printout the setup
#
InDetFlags.init()
InDetFlags.printInfo()
InDetFlags.print_JobProperties()

if not 'InDetNewTrackingCuts' in dir():
    #
    # --- setup physics cuts with default values
    #
    print "InDetRec_jobOptions: InDetCutValues not set - setting to defaults"
    from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
    InDetNewTrackingCuts = ConfiguredNewTrackingCuts("Offline")
    from InDetRecExample.ConfiguredVertexingCuts import ConfiguredVertexingCuts
    InDetPrimaryVertexingCuts = ConfiguredVertexingCuts("NoLumi")

#
# ----------- printout the setup
#
InDetNewTrackingCuts.printInfo()

if not 'InDetKeys' in dir():
    #
    # --- setup StoreGate keys (JobProperties!)
    #
    print "InDetRec_jobOptions: InDetKeys not set - setting to defaults"
    from InDetRecExample.InDetKeys import InDetKeys
    InDetKeys.lock_JobProperties()
예제 #7
0
    IDDerivationSequence += xAOD_TRT_PrepDataToxAOD
    if (printIdTrkDxAODConf):
        print(xAOD_TRT_PrepDataToxAOD)
        print(xAOD_TRT_PrepDataToxAOD.properties())

    # to store dEdx info
    from TRT_ElectronPidTools.TRT_ElectronPidToolsConf import TRT_ToT_dEdx
    TRT_dEdx_Tool = TRT_ToT_dEdx(name="TRT_ToT_dEdx")
    from InDetRecExample.TrackingCommon import getInDetTRT_LocalOccupancy
    TRT_dEdx_Tool.TRT_LocalOccupancyTool = getInDetTRT_LocalOccupancy()
    ToolSvc += TRT_dEdx_Tool

    TrackCollectionKeys = []

    from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
    InDetNewTrackingCutsPixel = ConfiguredNewTrackingCuts("Pixel")
    print(InDetNewTrackingCutsPixel)

    if not hasattr(topSequence, 'InDetTrackCollectionMerger'):
        from InDetRecExample import TrackingCommon
        from InDetTrackPRD_Association.InDetTrackPRD_AssociationConf import InDet__InDetTrackPRD_Association
        InDetPRD_Association = InDet__InDetTrackPRD_Association(
            name='InDetPRD_Association' +
            InDetNewTrackingCutsPixel.extension(),
            AssociationTool=TrackingCommon.
            getInDetPRDtoTrackMapToolGangedPixels(),
            TracksName=['Tracks'],
            AssociationMapName="PRDtoTrackMap" + InDetKeys.UnslimmedTracks(),
            OutputLevel=INFO)
        topSequence += InDetPRD_Association
        print(InDetPRD_Association)
예제 #8
0
        #from tauRec.tauRecFlags import jobproperties
        #jobproperties.tauRecFlags.doTauRec = False

        # AO (this turns off all MissingET and EFlow)
        rec.doJetMissingETTag = False
        recAlgs.doMissingET = False
        recAlgs.doMissingETSig = False
        recAlgs.doEFlow = False
        recAlgs.doEFlowJet = False

        recAlgs.doMuGirl = False
        recAlgs.doMuidLowPt = False
        recAlgs.doCaloTrkMuId = False
        from MuonCombinedRecExample.MuonCombinedRecFlags import muonCombinedRecFlags
        muonCombinedRecFlags.doMuGirlLowBeta = False

        ptCutVal = 0.5
        if jobproperties.HIRecExampleFlags.ptCutOn():
            ptCutVal = jobproperties.HIRecExampleFlags.ptCut()
            print("ptCutOn: Modify ptCutVal to ", ptCutVal, " GeV")
            if (rec.triggerStream() == "HardProbes") or (globalflags.DataSource
                                                         == "geant4"):
                print("Modify ConfiguredNewTrackingCuts__minPT to ", ptCutVal,
                      " GeV")
                from InDetRecExample.ConfiguredNewTrackingCuts import ConfiguredNewTrackingCuts
                InDetNewTrackingCuts = ConfiguredNewTrackingCuts("HeavyIon")
                InDetNewTrackingCuts._ConfiguredNewTrackingCuts__minPT = ptCutVal * Units.GeV
                InDetNewTrackingCuts.printInfo()

############################################################################