Esempio n. 1
0
from AthenaMonitoring.AthenaMonitoringConf import AthenaMonManager

monMan = AthenaMonManager("PhysValMonManager")
monMan.ManualDataTypeSetup = True
monMan.DataType = "monteCarlo"
monMan.Environment = "altprod"
monMan.ManualRunLBSetup = True
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()
Esempio n. 2
0
    ##################################################

    from InDetGlobalMonitoring.InDetGlobalMonitoringConf import InDetGlobalTrackMonTool
    InDetGlobalTrackMonTool = InDetGlobalTrackMonTool(
        name="InDetGlobalTrackMonTool",
        histoPathBase="/GLOBAL",
        DoIBL=InDetFlags.doIBL(),
        trackMax=75)

    TrackCollection = InDetKeys.UnslimmedTracks()

    from InDetTrackSelectionTool.InDetTrackSelectionToolConf import InDet__InDetTrackSelectionTool
    InDetTrackSelectionToolGlobalMon_Baseline = InDet__InDetTrackSelectionTool(
        name="InDetTrackSelectionToolGlobalMon_Baseline",
        UseTrkTrackTools=True,
        CutLevel="TightPrimary",
        maxNPixelHoles=1,
        minPt=5000,
        TrackSummaryTool=InDetTrackSummaryTool,
        Extrapolator=InDetExtrapolator)

    InDetTrackSelectionToolGlobalMon_TightPrimary = InDet__InDetTrackSelectionTool(
        name="InDetTrackSelectionToolGlobalMon_TightPrimary",
        UseTrkTrackTools=True,
        CutLevel="TightPrimary",
        minPt=5000,
        TrackSummaryTool=InDetTrackSummaryTool,
        Extrapolator=InDetExtrapolator)

    #ToolSvc += InDetTrackSelectionToolGlobalMon_Baseline
    #ToolSvc += InDetTrackSelectionToolGlobalMon_TightPrimary
            print "NewInDetAlignAlgSetup : setting TRT Cuts"
            minNTRTHits = 20
            maxEtaForTRTHitCuts = 1.9
            maxTRTEtaAcceptance = 0.

        from InDetTrackSelectionTool.InDetTrackSelectionToolConf import InDet__InDetTrackSelectionTool
        if not newInDetAlignAlg_Options["Cosmics"]:

            trackSelectorNew = InDet__InDetTrackSelectionTool(
                name="InDetTrackSelectionToolAlignTracks",
                UseTrkTrackTools=True,
                minPt=float(newInDetAlignAlg_Options["PtCut"]),
                maxD0=500,
                maxZ0=500,
                #minNPixelHits = 1,
                #minNPixelHitsPhysical = 1,
                #minNSiHits = 5,
                #minNSiHitsPhysical = 7,
                minNTrtHits=minNTRTHits,
                maxEtaForTrtHitCuts=maxEtaForTRTHitCuts,
                maxTrtEtaAcceptance=maxTRTEtaAcceptance,
                CutLevel="TightPrimary",
                TrackSummaryTool=InDetTrackSummaryTool,
                Extrapolator=InDetExtrapolator)
            if newInDetAlignAlg_Options["d0significanceCut"] > 0:
                trackSelectorNew.maxD0overSigmaD0 = 3
        else:
            trackSelectorNew = InDet__InDetTrackSelectionTool(
                name="InDetTrackSelectionToolAlignTracks",
                UseTrkTrackTools=True,
                minPt=float(newInDetAlignAlg_Options["PtCut"]),
                maxD0=9999.,
Esempio n. 4
0
from TrkVertexFitterUtils.TrkVertexFitterUtilsConf import Trk__ImpactPoint3dEstimator
InDetImpactPoint3dEstimator = Trk__ImpactPoint3dEstimator(
    name="InDetImpactPoint3dEstimator", Extrapolator=InDetExtrapolator)
ToolSvc += InDetImpactPoint3dEstimator

from InDetTrackSelectionTool.InDetTrackSelectionToolConf import InDet__InDetTrackSelectionTool
InDetTrackSelectorTool = InDet__InDetTrackSelectionTool(
    name="InDetDetailedTrackSelectionTool",
    CutLevel=InDetPrimaryVertexingCuts.TrackCutLevel(),
    minPt=InDetPrimaryVertexingCuts.minPT(),
    maxD0=InDetPrimaryVertexingCuts.IPd0Max(),
    maxZ0=InDetPrimaryVertexingCuts.z0Max(),
    maxZ0SinTheta=InDetPrimaryVertexingCuts.IPz0Max(),
    maxSigmaD0=InDetPrimaryVertexingCuts.sigIPd0Max(),
    maxSigmaZ0SinTheta=InDetPrimaryVertexingCuts.sigIPz0Max(),
    maxAbsEta=InDetPrimaryVertexingCuts.etaMax(),
    minNInnermostLayerHits=InDetPrimaryVertexingCuts.nHitInnermostLayer(),
    minNPixelHits=InDetPrimaryVertexingCuts.nHitPix(),
    maxNPixelHoles=InDetPrimaryVertexingCuts.nHolesPix(),
    minNSctHits=InDetPrimaryVertexingCuts.nHitSct(),
    minNTrtHits=InDetPrimaryVertexingCuts.nHitTrt(),
    minNSiHits=InDetPrimaryVertexingCuts.nHitSi(),
    TrackSummaryTool=InDetTrackSummaryTool,
    Extrapolator=InDetExtrapolator)
ToolSvc += InDetTrackSelectorTool

from TrkVertexSeedFinderUtils.TrkVertexSeedFinderUtilsConf import Trk__ImagingSeedTuningAlg
algSeq += Trk__ImagingSeedTuningAlg(
    name="SeedTuningAlg",
    useBeamConstraint=True,
    truthWindow=0.5,
Esempio n. 5
0
    PrimVtxContainerName = InDetKeys.xAODVertexContainer(),
    PassAllTracks = True, ## Uncomment this line to bypass track slection
    TrackSelectorTool    = InDetDetailedTrackSelector_Default
)




# Tracks For LoosePrimary Selection


m_TrackSelectorTool_LoosePrimary = InDet__InDetTrackSelectionTool(name         = "InDetTrackSelectorLoosePrimary",
                                                                  UseTrkTrackTools = True,
                                                                  minPt = 5000,
                                                                  maxD0 = 100000,
                                                                  maxZ0SinTheta = 150,
                                                                  minNTrtHits = 0,
                                                                  CutLevel = "LoosePrimary",
                                                                  TrackSummaryTool    = InDetTrackSummaryTool,
                                                                  Extrapolator        = InDetExtrapolator)






#ToolSvc +=  m_TrackSelectorTool_LoosePrimary


LPSelection = InDetAlignMon__TrackSelectionTool(
    name = "InDetAlignMonTracksSelectionTool_LP",
Esempio n. 6
0
def setup_eflowCaloObjectCreator(Configured, nameModifier, mlog):

    if nameModifier != "EM" and nameModifier != "LC":
        mlog.error(
            "Invalid calorimeter scale was specified : should be LC or EM, but was "
            + nameModifier)
        return False

    try:
        from eflowRec.eflowRecConf import eflowPreparation
        eflowPreparationAlgorithm = eflowPreparation("eflow" + nameModifier +
                                                     "CaloObjectBuilder")
    except:
        mlog.error("could not import eflowRec.eflowPreparation")
        print traceback.format_exc()
        return False

    Configured._eflowPreparationHandle = eflowPreparationAlgorithm

    try:
        from eflowRec.eflowRecConf import eflowTrackCaloExtensionTool
        TrackCaloExtensionTool = eflowTrackCaloExtensionTool()
    except:
        mlog.error("could not import eflowRec.eflowTrackCaloExtensionTool")
        print traceback.format_exc()
        return False

    eflowPreparationAlgorithm.TrackExtrapolatorTool = TrackCaloExtensionTool

    # sets output key of C++ algorithm equal to the python side
    eflowPreparationAlgorithm.EflowCaloObjectsOutputName = Configured.outputKey(
    )

    from eflowRec.eflowRecFlags import jobproperties

    if "EM" == nameModifier:
        eflowPreparationAlgorithm.ClustersName = "CaloTopoCluster"
        eflowPreparationAlgorithm.CalClustersName = "CaloCalTopoClusters"
    elif "LC" == nameModifier:
        eflowPreparationAlgorithm.ClustersName = "CaloCalTopoClusters"
        eflowPreparationAlgorithm.CalClustersName = ""

    if True == jobproperties.eflowRecFlags.useLeptons:
        eflowPreparationAlgorithm.useLeptons = True
        if True == jobproperties.eflowRecFlags.storeLeptonCells:
            eflowPreparationAlgorithm.storeLeptonCells = True
            if "LC" == nameModifier:
                eflowPreparationAlgorithm.eflowLeptonCellsName = "eflowRec_leptonCellContainer_LC"
    else:
        eflowPreparationAlgorithm.useLeptons = False

    if "LC" == nameModifier:
        eflowPreparationAlgorithm.EflowCaloObjectsOutputName = "eflowCaloObjects_LC"
        eflowPreparationAlgorithm.eflowElectronsName = "eflowRec_selectedElectrons_LC"

    try:
        from InDetTrackSelectionTool.InDetTrackSelectionToolConf import InDet__InDetTrackSelectionTool
        TrackSelectionTool = InDet__InDetTrackSelectionTool()
    except:
        mlog.error(
            "could not import InDetTrackSelectionTool.InDet__InDetTrackSelectionTool"
        )
        print traceback.format_exc()
        return False

    try:
        from AthenaCommon.AppMgr import ToolSvc
    except:
        mlog.error("could not import ToolSvc")
        print traceback.format_exc()
        return False

    ToolSvc += TrackSelectionTool

    TrackSelectionTool.CutLevel = "TightPrimary"
    TrackSelectionTool.minPt = 500.0

    eflowPreparationAlgorithm.TrackSelectionTool = TrackSelectionTool

    from RecExConfig.ObjKeyStore import objKeyStore
    objKeyStore.addTransient(Configured.outputType(), Configured.outputKey())

    from AthenaCommon.AlgSequence import AlgSequence
    topSequence = AlgSequence()
    topSequence += eflowPreparationAlgorithm

    return True
alg1.WhichLabel = 1  # was 3 # FLS
## FLS - to augment track info
alg1.TrackToVertexIPEstimator = ToolSvc.trkIPEstimator

algSeq += alg1

## CAN UNCOMMENT FOR DETAILED DEBUGGING
#ToolSvc.IP2DTag_Trig.OutputLevel = DEBUG
#ToolSvc.MV2c20Tag.OutputLevel = DEBUG
#ToolSvc.BTagCalibrationBrokerTool_Trig.OutputLevel = DEBUG

#ToolSvc.MV2c00Tag_Trig.MV2CalibAlias = "AntiKt4EMTopo"
#ToolSvc.MV2c00Tag_Trig.MV2CalibAlias = "AntiKt4EMTopo"
#ToolSvc.MV2c20Tag_Trig.MV2CalibAlias = "AntiKt4EMTopo"

### Additional Tools needed by the algo
from TrackVertexAssociationTool.TrackVertexAssociationToolConf import CP__TightTrackVertexAssociationTool
ToolSvc += CP__TightTrackVertexAssociationTool("TightVertexAssocTool",
                                               dzSinTheta_cut=3,
                                               doPV=True)

from InDetTrackSelectionTool.InDetTrackSelectionToolConf import InDet__InDetTrackSelectionTool
ToolSvc += InDet__InDetTrackSelectionTool("InDetTrackSelTool",
                                          CutLevel="MYINDETTRACKCUT",
                                          maxZ0SinTheta=3)

from TrigDecisionTool.TrigDecisionToolConf import Trig__TrigDecisionTool
ToolSvc += Trig__TrigDecisionTool("TrigDecisionTool")
from TrigEDMConfig.TriggerEDM import EDMLibraries
ToolSvc.TrigDecisionTool.Navigation.Dlls = EDMLibraries
Esempio n. 8
0
MUON1Seq += CfgMgr.DerivationFramework__DerivationKernel(
    "MUON1Kernel",
    AugmentationTools=MUON1AugmentTools,
    SkimmingTools=[MUON1SkimmingTool1],
    ThinningTools=MUON1ThinningTools)

#====================================================================
#       Decorate the vertices with the sum pt of their tracks
#====================================================================
alg_name = "MUONVertexDecorationAlg"
if not hasattr(DerivationFrameworkJob, alg_name):
    from DerivationFrameworkMuons.DerivationFrameworkMuonsConf import VertexDecoratorAlg
    from InDetTrackSelectionTool.InDetTrackSelectionToolConf import InDet__InDetTrackSelectionTool

    the_alg = VertexDecoratorAlg(alg_name)
    InDetTool = InDet__InDetTrackSelectionTool(
        "TrackSelectionTool_VertexDecoration")
    InDetTool.minPt = 400
    InDetTool.CutLevel = "Loose"
    ToolSvc += InDetTool
    the_alg.TrackSelectionTool = InDetTool
    DerivationFrameworkJob += the_alg

#====================================================================
# JetTagNonPromptLepton decorations
#====================================================================
import JetTagNonPromptLepton.JetTagNonPromptLeptonConfig as JetTagConfig
import LeptonTaggers.LeptonTaggersConfig as LepTagConfig
if not hasattr(MUON1Seq, "Muons_decoratePromptLepton"):
    JetTagConfig.ConfigureAntiKt4PV0TrackJets(MUON1Seq, "MUON1")
    MUON1Seq += JetTagConfig.GetDecoratePromptLeptonAlgs()
    MUON1Seq += LepTagConfig.GetDecorateImprovedPromptLeptonAlgs()
Esempio n. 9
0
def CaloExtensionBuilder(cutLevel="TightPrimary", minPT=100.0):
    try:
        from TrkExTools.AtlasExtrapolator import AtlasExtrapolator
        from TrackToCalo.TrackToCaloConf import Trk__ParticleCaloExtensionTool
    except:
        mlog.error(
            "could not import TrackToCaloConf.Trk__ParticleCaloExtensionTool")
        mlog.error("could not import TrkExTools.AtlasExtrapolator")
        mlog.error(traceback.format_exc())
    try:
        from TrackToCalo.TrackToCaloConf import Trk__CaloExtensionBuilderAlg as CaloExtensionBuilderAlg
    except:
        mlog.error(
            "could not import TrackToCaloConf.Trk__CaloExtensionBuilderAlg")
        mlog.error(traceback.format_exc())
    try:
        from InDetTrackSelectionTool.InDetTrackSelectionToolConf import InDet__InDetTrackSelectionTool
        from InDetTrackSelectorTool.InDetTrackSelectorToolConf import InDet__InDetDetailedTrackSelectorTool
    except:
        mlog.error(
            "could not import InDetTrackSelectionTool.InDet__InDetTrackSelectionTool"
        )
        mlog.error(traceback.format_exc())
    try:
        from AthenaCommon.AppMgr import ToolSvc
    except:
        mlog.error("could not import ToolSvc")
        mlog.error(traceback.format_exc())
    try:
        from AthenaCommon.AlgSequence import AlgSequence
    except:
        mlog.error("could not import AlgSequence")
        mlog.error(traceback.format_exc())

    topSequence = AlgSequence()

    theAtlasExtrapolator = AtlasExtrapolator(
        name="CaloExtensionBuilderAtlasExtrapolator")
    theAtlasExtrapolator.DoCaloDynamic = False  # this turns off dynamic

    pcExtensionTool = Trk__ParticleCaloExtensionTool(
        Extrapolator=theAtlasExtrapolator)
    ToolSvc += pcExtensionTool

    CaloExtensionBuilderTool = CaloExtensionBuilderAlg(
        LastCaloExtentionTool=pcExtensionTool)

    from AthenaCommon.BeamFlags import jobproperties

    TrackSelectionToolHC = InDet__InDetTrackSelectionTool(
        name="CaloExtensionBuilderTrackSelectionTool",
        minPt=minPT,
        CutLevel=cutLevel,
        minNSiHits=7
    )  # PixelHits + SCTHits + PixelDeadSensors + SCTDeadSensors

    CaloExtensionBuilderTool.TrkSelection = TrackSelectionToolHC

    ToolSvc += CaloExtensionBuilderTool.LastCaloExtentionTool

    topSequence += CaloExtensionBuilderTool

    return True
Esempio n. 10
0
    from ParticleJetTools.ParticleJetToolsConf import Analysis__JetPartonTruthLabel
    from ParticleJetTools.ParticleJetToolsConf import CopyTruthJetParticles
    from ParticleJetTools.ParticleJetToolsConf import ParticleJetDeltaRLabelTool

ghostScaleFactor = 1e-40

#--------------------------------------------------------------
# Track selection.
#--------------------------------------------------------------

# This is the InDet loose selection from
# https://twiki.cern.ch/twiki/bin/view/AtlasProtected/InDetTrackingPerformanceGuidelines
# October 28, 2014
#jtm += InDet__InDetDetailedTrackSelectionTool(
jtm += InDet__InDetTrackSelectionTool("trk_trackselloose",
                                      CutLevel="Loose",
                                      minPt=500)

jtm += JetTrackSelectionTool("trackselloose",
                             InputContainer=jtm.trackContainer,
                             OutputContainer="JetSelectedTracks",
                             Selector=jtm.trk_trackselloose)

jtm += InDet__InDetTrackSelectionTool("trk_trackselloose_trackjets",
                                      CutLevel="Loose",
                                      minPt=500)

jtm += JetTrackSelectionTool(
    "trackselloose_trackjets",
    InputContainer=jtm.trackContainer,
    OutputContainer="JetSelectedTracks_LooseTrackJets",
Esempio n. 11
0
#It can call the InDetDetailedTrackSelectorTool or the InDetTrackSelectionTool, that actually make the track selection

m_alignMonTrackSelectorTool = []
m_alignMonTrackSelectionTool = []

if not jobproperties.Beam.beamType() == 'cosmics':
    #
    # Collisions and single beam running
    #
    # Track selector that implement the Tracking CP supported cuts for Run2 using loosePrimary cut
    m_alignMonTrackSelectorTool.append(
        InDet__InDetTrackSelectionTool(
            name="InDetTrackSelectionTool_LoosePrimary",
            UseTrkTrackTools=True,
            minPt=5000,
            maxD0=100000,
            maxZ0SinTheta=150,
            CutLevel="LoosePrimary",
            TrackSummaryTool=InDetTrackSummaryTool,
            Extrapolator=InDetExtrapolator))

    # Track selector that implement the Tracking CP supported cuts for Run2 using tightPrimary cut
    m_alignMonTrackSelectorTool.append(
        InDet__InDetTrackSelectionTool(name="InDetTrackSelectionTool_NoTRT",
                                       UseTrkTrackTools=True,
                                       minPt=5000,
                                       maxD0=100000,
                                       maxZ0SinTheta=150,
                                       minNTrtHits=0,
                                       CutLevel="LoosePrimary",
                                       TrackSummaryTool=InDetTrackSummaryTool,