HIGG1D1.py

Go to the documentation of this file.

# Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
#====================================================================
# DAOD_HIGG1D1.py
# This defines DAOD_HIGG1D1, an unskimmed DAOD format for Run 3.
# It contains the variables and objects needed for the large majority 
# of physics analyses in ATLAS.
# It requires the flag HIGG1D1 in Derivation_tf.py   
#====================================================================
 
from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
from AthenaConfiguration.ComponentFactory import CompFactory
from AthenaConfiguration.Enums import MetadataCategory
 
# Main algorithm config
def HIGG1D1KernelCfg(flags, name='HIGG1D1Kernel', **kwargs):
    """Configure the derivation framework driving algorithm (kernel) for HIGG1D1"""
    acc = ComponentAccumulator()
 
    # Common augmentations
    from DerivationFrameworkPhys.PhysCommonConfig import PhysCommonAugmentationsCfg
    acc.merge(PhysCommonAugmentationsCfg(flags, TriggerListsHelper = kwargs['TriggerListsHelper']))
 
    # Diphoton Vertex creation  
    from DerivationFrameworkHiggs.HIGG1D1CustomVertexConfig import ZeeVertexRefitterCfg, DiPhotonVertexCfg
    from DerivationFrameworkHiggs.HIGG1D1CustomJetsConfig import HIGG1D1CustomJetsCfg, HIGG1D1CustomJetsCleaningCfg
    acc.merge(ZeeVertexRefitterCfg(flags))
    acc.merge(DiPhotonVertexCfg(flags))
    
    #CustomJetsConfig
    acc.merge(HIGG1D1CustomJetsCfg(flags))
    acc.merge(HIGG1D1CustomJetsCleaningCfg(flags))
 
    # Reorder event cleaning / decoration locking to run at this
    # point in the algorithm sequence. 
    # See comments in JetCommonConfig.AddEventCleanFlagsCfg.
    seq = acc.getSequence('AthAlgSeq')
    cleanSeqs = [acc.getSequence ('EventCleanSeq'),
                 acc.getSequence ('EventCleanLockSeq')]
    seq.Members = [a for a in seq.Members if a not in cleanSeqs] + cleanSeqs
 
    from BTagging.FlavorTaggingConfig import FlavorTaggingCfg
    acc.merge(FlavorTaggingCfg(flags, 'AntiKt4EMPFlowCustomVtxJets', 'HggPrimaryVertices', trackAugmenterPrefix='btagIpHgg_'))
 
    #Custom MET
    from DerivationFrameworkJetEtMiss.METCommonConfig import METCustomVtxCfg
    acc.merge(METCustomVtxCfg(flags, 'HggPrimaryVertices', 'AntiKt4EMPFlowCustomVtx', 'CHSGCustomVtxParticleFlowObjects'))
 
    # Thinning tools...
    from DerivationFrameworkInDet.InDetToolsConfig import TrackParticleThinningCfg, MuonTrackParticleThinningCfg, TauTrackParticleThinningCfg, DiTauTrackParticleThinningCfg
    from DerivationFrameworkMCTruth.TruthDerivationToolsConfig import GenericTruthThinningCfg
    from DerivationFrameworkTools.DerivationFrameworkToolsConfig import GenericObjectThinningCfg
    
    thinningTools = []
    
    # Truth thinning
    if flags.Input.isMC:
        truth_conditions = ["((abs(TruthParticles.pdgId) >= 23) && (abs(TruthParticles.pdgId) <= 25))",  # W, Z and Higgs
                            "((abs(TruthParticles.pdgId) >= 11) && (abs(TruthParticles.pdgId) <= 16))",  # Leptons
                            "((abs(TruthParticles.pdgId) ==  6))",                                       # Top quark
                            "((abs(TruthParticles.pdgId) == 22) && (TruthParticles.pt > 1*GeV))",        # Photon
                            "(abs(TruthParticles.pdgId) >=  1000000)",                                   # BSM
                            "(TruthParticles.isGenStable)"]           # stable particles
        truth_expression = f'({" || ".join(truth_conditions)})'
 
        HIGG1D1GenericTruthThinningTool   = acc.getPrimaryAndMerge(GenericTruthThinningCfg(
            flags,
            name                          = "HIGG1D1GenericTruthThinningTool",
            StreamName                    = kwargs['StreamName'], 
            ParticleSelectionString       = truth_expression,
            PreserveDescendants           = False,
            PreserveGeneratorDescendants  = True,
            PreserveAncestors             = True))
        thinningTools.append(HIGG1D1GenericTruthThinningTool)
 
    # Inner detector group recommendations for indet tracks in analysis
    # https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/DaodRecommendations
    HIGG1D1_thinning_expression = "InDetTrackParticles.DFCommonTightPrimary && \
                                   abs(DFCommonInDetTrackZ0AtPV)*sin(InDetTrackParticles.theta) < 3.0*mm && \
                                   InDetTrackParticles.pt > 10*GeV"
    HIGG1D1TrackParticleThinningTool = acc.getPrimaryAndMerge(TrackParticleThinningCfg(
        flags,
        name                    = "HIGG1D1TrackParticleThinningTool",
        StreamName              = kwargs['StreamName'], 
        SelectionString         = HIGG1D1_thinning_expression,
        InDetTrackParticlesKey  = "InDetTrackParticles"))
    
    # Include inner detector tracks associated with muons
    HIGG1D1MuonTPThinningTool = acc.getPrimaryAndMerge(MuonTrackParticleThinningCfg(
        flags,
        name                    = "HIGG1D1MuonTPThinningTool",
        StreamName              = kwargs['StreamName'],
        MuonKey                 = "Muons",
        InDetTrackParticlesKey  = "InDetTrackParticles"))
 
    # Only keep tau tracks (and associated ID tracks) classified as charged tracks
    HIGG1D1TauTPThinningTool = acc.getPrimaryAndMerge(TauTrackParticleThinningCfg(
        flags,
        name                   = "HIGG1D1TauTPThinningTool",
        StreamName             = kwargs['StreamName'],
        TauKey                 = "TauJets",
        InDetTrackParticlesKey = "InDetTrackParticles",
        DoTauTracksThinning    = True,
        TauTracksKey           = "TauTracks"))
 
    # ID tracks associated with high-pt di-tau
    HIGG1D1DiTauTPThinningTool = acc.getPrimaryAndMerge(DiTauTrackParticleThinningCfg(
        flags,
        name                    = "HIGG1D1DiTauTPThinningTool",
        StreamName              = kwargs['StreamName'],
        DiTauKey                = "DiTauJets",
        InDetTrackParticlesKey  = "InDetTrackParticles"))
 
    
    HIGG1D1DiTauLowPtThinningTool = acc.getPrimaryAndMerge(GenericObjectThinningCfg(flags,
                                                                                 name            = "PHYSDiTauLowPtThinningTool",
                                                                                 StreamName      = kwargs['StreamName'],
                                                                                 ContainerName   = "DiTauJetsLowPt",
                                                                                 SelectionString = "DiTauJetsLowPt.nSubjets > 1"))
    
    # ID tracks associated with low-pt ditau
    HIGG1D1DiTauLowPtTPThinningTool = acc.getPrimaryAndMerge(DiTauTrackParticleThinningCfg(flags,
                                                                                        name                    = "PHYSDiTauLowPtTPThinningTool",
                                                                                        StreamName              = kwargs['StreamName'],
                                                                                        DiTauKey                = "DiTauJetsLowPt",
                                                                                        InDetTrackParticlesKey  = "InDetTrackParticles",
                                                                                        SelectionString         = "DiTauJetsLowPt.nSubjets > 1"))
    
    
    # SkimmingTool
    from DerivationFrameworkHiggs.SkimmingToolHIGG1Config import SkimmingToolHIGG1Cfg
    
    # Requires something in this list of triggers
    SkipTriggerRequirement = (flags.Input.isMC and float(flags.Beam.Energy) == 4000000.0) or not flags.Reco.EnableTrigger
    # 8 TeV MC does not have trigger information
    print("HIGG1D1.py SkipTriggerRequirement", SkipTriggerRequirement)
    TriggerExp = []
    if not SkipTriggerRequirement:
        if float(flags.Beam.Energy) == 4000000.0:
            #  8 TeV data
            TriggerExp               = ["EF_g35_loose_g25_loose"]
        if float(flags.Beam.Energy) == 6500000.0:
            # 13 TeV MC 
            TriggerExp               = ["HLT_g35_medium_g25_medium",
                                        "HLT_g35_medium_g25_medium_L12EM20VH",
                                        "HLT_g35_medium_g25_medium_L12EM15VH",
                                        "HLT_g35_loose_g25_loose",
                                        "HLT_g35_loose_g25_loose_L12EM15VH",
                                        "HLT_g120_loose",
                                        "HLT_g140_loose",
                                        "HLT_3g20_loose",
                                        "HLT_3g15_loose",
                                        "HLT_2g6_tight_icalotight_L1J50",
                                        "HLT_2g6_tight_icalotight_L1J100",
                                        "HLT_2g6_loose_L1J50",
                                        "HLT_2g6_loose_L1J100",
                                        "HLT_2g50_loose",
                                        "HLT_2g25_tight_L12EM20VH",
                                        "HLT_2g25_loose_g15_loose",
                                        "HLT_2g22_tight",
                                        "HLT_2g22_tight_L12EM15VHI",
                                        "HLT_2g22_tight_icalovloose",
                                        "HLT_2g22_tight_icalovloose_L12EM15VHI",
                                        "HLT_2g22_tight_icalotight_L12EM15VHI",
                                        "HLT_2g20_tight",
                                        "HLT_2g20_tight_icalovloose_L12EM15VHI",
                                        "HLT_2g20_tight_icalotight_L12EM15VHI",
                                        "HLT_2g20_loose",
                                        "HLT_2g20_loose_L12EM15",
                                        "HLT_2g20_loose_g15_loose",
                                        "HLT_2g50_loose_L12EM20VH"
                                       ]
        if float(flags.Beam.Energy) == 6800000.0:
            # 13.6 TeV
            TriggerExp               = ["HLT_g35_medium_g25_medium_L12EM20VH",
                                        "HLT_g35_medium_g25_medium_L12eEM24L",
                                        "HLT_g140_loose_L1EM22VHI",
                                        "HLT_g140_loose_L1eEM26M"
                                        "HLT_g120_loose_L1EM22VHI",
                                        "HLT_g120_loose_L1eEM26M",
                                        "HLT_2g9_loose_25dphiAA_invmAA80_L1DPHI-M70-2eEM9",
                                        "HLT_2g9_loose_25dphiAA_invmAA80_L12EM7",
                                        "HLT_2g50_loose_L12EM20VH",
                                        "HLT_2g50_loose_L12eEM24L",
                                        "HLT_2g25_loose_g15_loose_L12EM20VH",
                                        "HLT_2g25_loose_g15_loose_L12eEM24L",
                                        "HLT_2g22_tight_L12EM15VHI",
                                        "HLT_2g22_tight_L12eEM18M",
                                        "HLT_2g20_tight_icaloloose_L12EM15VHI",
                                        "HLT_2g20_tight_icaloloose_L12eEM18M",
                                        "HLT_2g20_loose_L12EM15VH",
                                        "HLT_2g20_loose_L12eEM18L",
                                        "HLT_2g15_tight_L1DPHI-M70-2eEM15M",
                                        "HLT_2g15_tight_25dphiAA_L1DPHI-M70-2eEM15M",
                                        "HLT_2g15_tight_25dphiAA_L12EM7",
                                        "HLT_2g15_tight_25dphiAA_invmAA80_L12EM7",
                                        "HLT_2g15_loose_25dphiAA_invmAA80_L1DPHI-M70-2eEM15M",
                                        "HLT_2g15_loose_25dphiAA_invmAA80_L12EM7",
                                       ]
    print("HIGG1D1.py Skimming Tool Triggers:", ",".join(TriggerExp))
    skimmingTool = acc.popToolsAndMerge( SkimmingToolHIGG1Cfg(flags,RequireTrigger=not SkipTriggerRequirement,Triggers=TriggerExp) )
    acc.addPublicTool(skimmingTool)
 
 
    # Finally the kernel itself
    thinningTools += [HIGG1D1TrackParticleThinningTool,
                      HIGG1D1MuonTPThinningTool,
                      HIGG1D1TauTPThinningTool,
                      HIGG1D1DiTauTPThinningTool,
                      HIGG1D1DiTauLowPtThinningTool,
                      HIGG1D1DiTauLowPtTPThinningTool]
    
    augmentationTools = [] 
 
    #====================================================================
    # Common calo decoration tools
    #====================================================================
 
    from DerivationFrameworkCalo.DerivationFrameworkCaloConfig import (
        CaloDecoratorKernelCfg, ClusterEnergyPerLayerDecoratorCfg)
    acc.merge(CaloDecoratorKernelCfg(flags))
 
    # Adding missing cluster energy decorators
    cluster_sizes = (3,5), (5,7), (7,7)
    for neta, nphi in cluster_sizes:
        cename = "ClusterEnergyPerLayerDecorator_%sx%s" % (neta, nphi)
        ClusterEnergyPerLayerDecorator = acc.popToolsAndMerge( ClusterEnergyPerLayerDecoratorCfg(flags, neta = neta, nphi=nphi, name=cename ))
        acc.addPublicTool(ClusterEnergyPerLayerDecorator)
        augmentationTools.append(ClusterEnergyPerLayerDecorator)
    
    #====================================================================
    # Truth categories decoration tool
    #====================================================================
    if flags.Input.isMC:
        from DerivationFrameworkHiggs.TruthCategoriesConfig import TruthCategoriesDecoratorCfg
        acc.merge(TruthCategoriesDecoratorCfg(flags, name="TruthCategoriesDecorator"))
    
    DerivationKernel = CompFactory.DerivationFramework.DerivationKernel
    acc.addEventAlgo(DerivationKernel(name, 
                                      SkimmingTools = [skimmingTool],
                                      ThinningTools = thinningTools,
                                      AugmentationTools = augmentationTools))     
    return acc
 
 
def HIGG1D1Cfg(flags):
 
    acc = ComponentAccumulator()
    
    from DerivationFrameworkPhys.TriggerListsHelper import TriggerListsHelper
    HIGG1D1TriggerListsHelper = TriggerListsHelper(flags)
    
    acc.merge(HIGG1D1KernelCfg(flags, name="HIGG1D1Kernel", StreamName = 'StreamDAOD_HIGG1D1', TriggerListsHelper = HIGG1D1TriggerListsHelper))
 
 
    # ============================
    # Define contents of the format
    # =============================
    from OutputStreamAthenaPool.OutputStreamConfig import OutputStreamCfg
    from xAODMetaDataCnv.InfileMetaDataConfig import SetupMetaDataForStreamCfg
    from DerivationFrameworkCore.SlimmingHelper import SlimmingHelper
    
    HIGG1D1SlimmingHelper = SlimmingHelper("HIGG1D1SlimmingHelper", NamesAndTypes = flags.Input.TypedCollections, flags = flags)
    HIGG1D1SlimmingHelper.SmartCollections = ["EventInfo",
                                              "Electrons",
                                              "Photons",
                                              "Muons",
                                              "PrimaryVertices",
                                              "InDetTrackParticles",
                                              "AntiKt4EMTopoJets",
                                              "AntiKt4EMPFlowJets",
                                              "BTagging_AntiKt4EMPFlow",
                                              "BTagging_AntiKtVR30Rmax4Rmin02Track",
                                              "MET_Baseline_AntiKt4EMTopo",
                                              "MET_Baseline_AntiKt4EMPFlow",
                                              "TauJets",
                                              "DiTauJets",
                                              "DiTauJetsLowPt",
                                              "AntiKt10LCTopoTrimmedPtFrac5SmallR20Jets",
                                              "AntiKtVR30Rmax4Rmin02PV0TrackJets"]
    
    excludedVertexAuxData = "-vxTrackAtVertex.-MvfFitInfo.-isInitialized.-VTAV"
    StaticContent = []
    StaticContent += ["xAOD::VertexContainer#SoftBVrtClusterTool_Tight_Vertices"]
    StaticContent += ["xAOD::VertexAuxContainer#SoftBVrtClusterTool_Tight_VerticesAux." + excludedVertexAuxData]
    StaticContent += ["xAOD::VertexContainer#SoftBVrtClusterTool_Medium_Vertices"]
    StaticContent += ["xAOD::VertexAuxContainer#SoftBVrtClusterTool_Medium_VerticesAux." + excludedVertexAuxData]
    StaticContent += ["xAOD::VertexContainer#SoftBVrtClusterTool_Loose_Vertices"]
    StaticContent += ["xAOD::VertexAuxContainer#SoftBVrtClusterTool_Loose_VerticesAux." + excludedVertexAuxData]
    
    HIGG1D1SlimmingHelper.StaticContent = StaticContent
    
    # Trigger content
    HIGG1D1SlimmingHelper.IncludeTriggerNavigation = False
    HIGG1D1SlimmingHelper.IncludeJetTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeMuonTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeEGammaTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeTauTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeEtMissTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeBJetTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeBPhysTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeMinBiasTriggerContent = False
 
    # Truth containers
    if flags.Input.isMC:
        HIGG1D1SlimmingHelper.AppendToDictionary = {
                    'TruthEvents':'xAOD::TruthEventContainer',
                    'TruthEventsAux':'xAOD::TruthEventAuxContainer',
                    'MET_Truth':'xAOD::MissingETContainer',
                    'MET_TruthAux':'xAOD::MissingETAuxContainer',
                    'TruthElectrons':'xAOD::TruthParticleContainer',
                    'TruthElectronsAux':'xAOD::TruthParticleAuxContainer',
                    'TruthMuons':'xAOD::TruthParticleContainer',
                    'TruthMuonsAux':'xAOD::TruthParticleAuxContainer',
                    'TruthPhotons':'xAOD::TruthParticleContainer',
                    'TruthPhotonsAux':'xAOD::TruthParticleAuxContainer',
                    'TruthTaus':'xAOD::TruthParticleContainer',
                    'TruthTausAux':'xAOD::TruthParticleAuxContainer',
                    'TruthNeutrinos':'xAOD::TruthParticleContainer',
                    'TruthNeutrinosAux':'xAOD::TruthParticleAuxContainer',
                    'TruthBSM':'xAOD::TruthParticleContainer',
                    'TruthBSMAux':'xAOD::TruthParticleAuxContainer',
                    'TruthBoson':'xAOD::TruthParticleContainer',
                    'TruthBosonAux':'xAOD::TruthParticleAuxContainer',
                    'TruthTop':'xAOD::TruthParticleContainer',
                    'TruthTopAux':'xAOD::TruthParticleAuxContainer',
                    'TruthForwardProtons':'xAOD::TruthParticleContainer',
                    'TruthForwardProtonsAux':'xAOD::TruthParticleAuxContainer',
                    'BornLeptons':'xAOD::TruthParticleContainer',
                    'BornLeptonsAux':'xAOD::TruthParticleAuxContainer',
                    'TruthBosonsWithDecayParticles':'xAOD::TruthParticleContainer',
                    'TruthBosonsWithDecayParticlesAux':'xAOD::TruthParticleAuxContainer',
                    'TruthBosonsWithDecayVertices':'xAOD::TruthVertexContainer',
                    'TruthBosonsWithDecayVerticesAux':'xAOD::TruthVertexAuxContainer',
                    'TruthBSMWithDecayParticles':'xAOD::TruthParticleContainer',
                    'TruthBSMWithDecayParticlesAux':'xAOD::TruthParticleAuxContainer',
                    'TruthBSMWithDecayVertices':'xAOD::TruthVertexContainer',
                    'TruthBSMWithDecayVerticesAux':'xAOD::TruthVertexAuxContainer',
                    'TruthHFWithDecayParticles':'xAOD::TruthParticleContainer',
                    'TruthHFWithDecayParticlesAux':'xAOD::TruthParticleAuxContainer',
                    'TruthHFWithDecayVertices':'xAOD::TruthVertexContainer',
                    'TruthHFWithDecayVerticesAux':'xAOD::TruthVertexAuxContainer',
                    'TruthCharm':'xAOD::TruthParticleContainer',
                    'TruthCharmAux':'xAOD::TruthParticleAuxContainer',
                    'TruthPrimaryVertices':'xAOD::TruthVertexContainer',
                    'TruthPrimaryVerticesAux':'xAOD::TruthVertexAuxContainer'
                }
 
        from DerivationFrameworkMCTruth.MCTruthCommonConfig import addTruth3ContentToSlimmerTool
        addTruth3ContentToSlimmerTool(HIGG1D1SlimmingHelper)
        HIGG1D1SlimmingHelper.AllVariables += ["TruthHFWithDecayParticles",
                                               "TruthHFWithDecayVertices",
                                               "TruthCharm",
                                               "TruthPileupParticles",
                                               "InTimeAntiKt4TruthJets",
                                               "OutOfTimeAntiKt4TruthJets",
                                               "AntiKt4TruthDressedWZJets",
                                               "AntiKt4TruthWZJets",
                                               "TruthEvents",
                                               "TruthPrimaryVertices",
                                               "TruthVertices",
                                               "TruthParticles"]
 
        HIGG1D1SlimmingHelper.ExtraVariables += ["Electrons.TruthLink",
                                                 "Muons.TruthLink",
                                                 "Photons.TruthLink",
                                                 "AntiKt4EMTopoJets.DFCommonJets_QGTagger_truthjet_nCharged.DFCommonJets_QGTagger_truthjet_pt.DFCommonJets_QGTagger_truthjet_eta.DFCommonJets_QGTagger_NTracks.DFCommonJets_QGTagger_TracksWidth.DFCommonJets_QGTagger_TracksC1.ConeExclBHadronsFinal.ConeExclCHadronsFinal.GhostBHadronsFinal.GhostCHadronsFinal.GhostBHadronsFinalCount.GhostBHadronsFinalPt.GhostCHadronsFinalCount.GhostCHadronsFinalPt",
                                                 "AntiKt4EMPFlowJets.QGTransformer_ConstScore.DFCommonJets_QGTagger_truthjet_nCharged.DFCommonJets_QGTagger_truthjet_pt.DFCommonJets_QGTagger_truthjet_eta.DFCommonJets_QGTagger_NTracks.DFCommonJets_QGTagger_TracksWidth.DFCommonJets_QGTagger_TracksC1.ConeExclBHadronsFinal.ConeExclCHadronsFinal.GhostBHadronsFinal.GhostCHadronsFinal.GhostBHadronsFinalCount.GhostBHadronsFinalPt.GhostCHadronsFinalCount.GhostCHadronsFinalPt",
                                                 "TruthPrimaryVertices.t.x.y.z",
                                                 "EventInfo.DFCommonJetsCustomVtx_eventClean_LooseBad.DFCommonJetsCustomVtx_eventClean_TightBad.hardScatterVertexLink.timeStampNSOffset",
                                                 "TauJets.dRmax.etOverPtLeadTrk"]
        if flags.Reco.EnableTrigger:
            HIGG1D1SlimmingHelper.ExtraVariables += [
                "HLT_xAOD__TrigMissingETContainer_TrigEFMissingET.ex.ey",
                "HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_mht.ex.ey"
            ]
 
    #  Additional content for HIGG1D1   
    HIGG1D1SlimmingHelper.AppendToDictionary.update({
                "AntiKt4EMPFlowCustomVtxJets": "xAOD::JetContainer",
                "AntiKt4EMPFlowCustomVtxJetsAux":"xAOD::JetAuxContainer",
                "METAssoc_AntiKt4EMPFlowCustomVtx": "xAOD::MissingETAssociationMap",
                "METAssoc_AntiKt4EMPFlowCustomVtxAux":"xAOD::MissingETAuxAssociationMap",
                "MET_Core_AntiKt4EMPFlowCustomVtx": "xAOD::MissingETContainer",
                "MET_Core_AntiKt4EMPFlowCustomVtxAux":"xAOD::MissingETAuxContainer",
                "HggPrimaryVertices":"xAOD::VertexContainer",
                "HggPrimaryVerticesAux":"xAOD::ShallowAuxContainer",
                "Kt4EMPFlowCustomVtxEventShape":"xAOD::EventShape",
                "Kt4EMPFlowCustomVtxEventShapeAux":"xAOD::EventShapeAuxInfo",
                "Kt4EMPFlowEventShape":"xAOD::EventShape",
                "Kt4EMPFlowEventShapeAux":"xAOD::EventShapeAuxInfo",
                "ZeeRefittedPrimaryVertices":"xAOD::VertexContainer",
                "ZeeRefittedPrimaryVerticesAux":"xAOD::VertexAuxContainer",
                "AFPSiHitContainer":"xAOD::AFPSiHitContainer",
                "AFPSiHitContainerAux":"xAOD::AFPSiHitAuxContainer",
                "AFPToFHitContainer":"xAOD::AFPToFHitContainer",
                "AFPToFHitContainerAux":"xAOD::AFPToFHitAuxContainer",
                "AFPVertexContainer":"xAOD::AFPVertexContainer",
                "AFPVertexContainerAux":"xAOD::AFPVertexAuxContainer",
                "AFPToFTrackContainer":"xAOD::AFPToFTrackContainer",
                "AFPToFTrackContainerAux":"xAOD::AFPToFTrackAuxContainer",
                "BTagging_AntiKt4EMPFlowCustomVtx":"xAOD::BTaggingContainer",
                "BTagging_AntiKt4EMPFlowCustomVtxAux":"xAOD::BTaggingAuxContainer"
             })
 
    HIGG1D1SlimmingHelper.AllVariables += [
            "HggPrimaryVertices",
            "ZeeRefittedPrimaryVertices",
            "AntiKt4EMPFlowCustomVtxJets",
            "Kt4EMPFlowCustomVtxEventShape",
            "Kt4EMPFlowEventShape",
            "METAssoc_AntiKt4EMPFlowCustomVtx",
            "MET_Core_AntiKt4EMPFlowCustomVtx"
        ]
    
    # Add AFP information
    HIGG1D1SlimmingHelper.AllVariables += ["AFPSiHitContainer",
                                           "AFPToFHitContainer",
                                           "AFPVertexContainer",
                                           "AFPToFTrackContainer"]
    # Add Btagging information
    from DerivationFrameworkFlavourTag.BTaggingContent import BTaggingStandardContent
    HIGG1D1SlimmingHelper.ExtraVariables += BTaggingStandardContent("AntiKt4EMPFlowCustomVtxJets", flags)
    HIGG1D1SlimmingHelper.ExtraVariables += BTaggingStandardContent("AntiKt4EMPFlowJets", flags)
 
    # is this really needed given Photons are in the AllVariables list ?
    from DerivationFrameworkEGamma.PhotonsCPDetailedContent import PhotonsCPDetailedContent
    HIGG1D1SlimmingHelper.ExtraVariables += PhotonsCPDetailedContent
 
 
    # Add the variables for Gain and Cluster energy
    from DerivationFrameworkCalo.DerivationFrameworkCaloConfig import (
        getGainDecorations, getClusterEnergyPerLayerDecorations )
    gainDecorations = getGainDecorations(acc, flags, 'HIGG1D1Kernel')
    HIGG1D1SlimmingHelper.ExtraVariables.extend(gainDecorations)
    clusterEnergyDecorations = getClusterEnergyPerLayerDecorations(
        acc, 'HIGG1D1Kernel' )
    HIGG1D1SlimmingHelper.ExtraVariables.extend(clusterEnergyDecorations)
        
    # Add HTXS variables
    HIGG1D1SlimmingHelper.ExtraVariables.extend(["EventInfo.HTXS_prodMode",
                                                 "EventInfo.HTXS_errorCode",
                                                 "EventInfo.HTXS_Stage0_Category",
                                                 "EventInfo.HTXS_Stage1_Category_pTjet25",
                                                 "EventInfo.HTXS_Stage1_Category_pTjet30",
                                                 "EventInfo.HTXS_Stage1_FineIndex_pTjet30",
                                                 "EventInfo.HTXS_Stage1_FineIndex_pTjet25",
                                                 "EventInfo.HTXS_Stage1_2_Category_pTjet25",
                                                 "EventInfo.HTXS_Stage1_2_Category_pTjet30",
                                                 "EventInfo.HTXS_Stage1_2_FineIndex_pTjet30",
                                                 "EventInfo.HTXS_Stage1_2_FineIndex_pTjet25",
                                                 "EventInfo.HTXS_Stage1_2_Fine_Category_pTjet25",
                                                 "EventInfo.HTXS_Stage1_2_Fine_Category_pTjet30",
                                                 "EventInfo.HTXS_Stage1_2_Fine_FineIndex_pTjet30",
                                                 "EventInfo.HTXS_Stage1_2_Fine_FineIndex_pTjet25",
                                                 "EventInfo.HTXS_Njets_pTjet25",
                                                 "EventInfo.HTXS_Njets_pTjet30",
                                                 "EventInfo.HTXS_isZ2vvDecay",
                                                 "EventInfo.HTXS_Higgs_eta",
                                                 "EventInfo.HTXS_Higgs_m",
                                                 "EventInfo.HTXS_Higgs_phi",
                                                 "EventInfo.HTXS_Higgs_pt",
                                                 "EventInfo.HTXS_V_jets30_eta",
                                                 "EventInfo.HTXS_V_jets30_m",
                                                 "EventInfo.HTXS_V_jets30_phi",
                                                 "EventInfo.HTXS_V_jets30_pt",
                                                 "EventInfo.HTXS_V_pt"])
    
    # Add variables from MaxCell decorator and PhotonPointing decorator
    HIGG1D1SlimmingHelper.ExtraVariables.extend(["Electrons.maxEcell_time",
                                                 "Electrons.maxEcell_energy",
                                                 "Electrons.maxEcell_gain",
                                                 "Electrons.maxEcell_onlId",
                                                 "Photons.maxEcell_time",
                                                 "Photons.maxEcell_energy",
                                                 "Photons.maxEcell_gain",
                                                 "Photons.maxEcell_onlId",
                                                 "Photons.zvertex"])
    # Add TTVA variables
    HIGG1D1SlimmingHelper.ExtraVariables.extend(["InDetTrackParticles.TTVA_AMVFVertices.TTVA_AMVFWeights.TTVA_AMVFVertices_forReco.TTVA_AMVFWeights_forReco.TTVA_AMVFVertices_forHiggs.TTVA_AMVFWeights_forHiggs.eProbabilityHT.numberOfTRTHits.numberOfTRTOutliers"])
 
    # Trigger content
    HIGG1D1SlimmingHelper.IncludeTriggerNavigation = False
    HIGG1D1SlimmingHelper.IncludeJetTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeMuonTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeEGammaTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeTauTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeEtMissTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeBJetTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeBPhysTriggerContent = False
    HIGG1D1SlimmingHelper.IncludeMinBiasTriggerContent = False
 
    # Trigger matching
    # Run 2
    if flags.Trigger.EDMVersion == 2:
        from DerivationFrameworkPhys.TriggerMatchingCommonConfig import AddRun2TriggerMatchingToSlimmingHelper
        AddRun2TriggerMatchingToSlimmingHelper(SlimmingHelper = HIGG1D1SlimmingHelper, 
                                               OutputContainerPrefix = "TrigMatch_", 
                                               TriggerList = HIGG1D1TriggerListsHelper.Run2TriggerNamesTau)
        AddRun2TriggerMatchingToSlimmingHelper(SlimmingHelper = HIGG1D1SlimmingHelper, 
                                               OutputContainerPrefix = "TrigMatch_",
                                               TriggerList = HIGG1D1TriggerListsHelper.Run2TriggerNamesNoTau)
    # Run 3, or Run 2 with navigation conversion
    if flags.Trigger.EDMVersion == 3 or (flags.Trigger.EDMVersion == 2 and flags.Trigger.doEDMVersionConversion):
        from TrigNavSlimmingMT.TrigNavSlimmingMTConfig import AddRun3TrigNavSlimmingCollectionsToSlimmingHelper
        AddRun3TrigNavSlimmingCollectionsToSlimmingHelper(HIGG1D1SlimmingHelper)
    
    # Output stream    
    HIGG1D1ItemList = HIGG1D1SlimmingHelper.GetItemList()
    acc.merge(OutputStreamCfg(flags, "DAOD_HIGG1D1", ItemList=HIGG1D1ItemList, AcceptAlgs=["HIGG1D1Kernel"]))
    acc.merge(SetupMetaDataForStreamCfg(flags, "DAOD_HIGG1D1", AcceptAlgs=["HIGG1D1Kernel"], createMetadata=[MetadataCategory.CutFlowMetaData, MetadataCategory.TruthMetaData]))
 
    return acc