MetMakerAlg.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/

 
 
 
#ifndef MET_ANALYSIS_ALGORITHMS__MET_MAKER_ALG_H
#define MET_ANALYSIS_ALGORITHMS__MET_MAKER_ALG_H
 
#include <AnaAlgorithm/AnaAlgorithm.h>
#include <SystematicsHandles/SysReadHandle.h>
#include <SystematicsHandles/SysHandleArray.h>
#include <SystematicsHandles/SysWriteHandle.h>
#include <SystematicsHandles/SysListHandle.h>
#include <SelectionHelpers/SysReadSelectionHandle.h>
#include <METInterface/IMETMaker.h>
#include <METInterface/IMETSystematicsTool.h>
#include <xAODBase/IParticleContainer.h>
#include <xAODMissingET/MissingETContainer.h>
#include <xAODMissingET/MissingETAuxContainer.h>
#include <xAODMissingET/MissingETAssociationMap.h>
#include <AsgDataHandles/ReadHandleKey.h>
#include <AsgTools/PropertyWrapper.h>
 
namespace CP
{
 
  class MetMakerAlg final : public EL::AnaAlgorithm
  {
  public:
    using EL::AnaAlgorithm::AnaAlgorithm;
    StatusCode initialize () override;
    StatusCode execute (const EventContext& ctx) override;
 
 

  private:
    ToolHandle<IMETMaker> m_makerTool {this, "makerTool", "METMaker", "the METMaker tool we apply"};

  private:
    ToolHandle<IMETSystematicsTool> m_systematicsTool {this, "systematicsTool", "", "the systematics tool we apply"};

  private:
    SG::ReadHandleKey<xAOD::MissingETContainer> m_metCoreKey {this, "metCore", "",  "the name of the core MissingETContainer"};

  private:
    SG::ReadHandleKey<xAOD::MissingETAssociationMap> m_metAssociationKey {this, "metAssociation", "", "the name of the MissingETAssociationMap"};

  private:
    SysListHandle m_systematicsList {this};

  private:
    SysReadHandle<xAOD::IParticleContainer> m_electronsHandle {
      this, "electrons", "", "the electron container to use" };

  private:
    SysReadSelectionHandle m_electronsSelection {
      this, "electronsSelection", "", "the selection on the input electrons" };

  private:
    Gaudi::Property<std::string> m_electronsKey {this, "electronsKey", "RefEle", "the key for the electrons"};

  private:
    SysReadHandle<xAOD::IParticleContainer> m_photonsHandle {
      this, "photons", "", "the photon container to use" };

  private:
    SysReadSelectionHandle m_photonsSelection {
      this, "photonsSelection", "", "the selection on the input photons" };

  private:
    Gaudi::Property<std::string> m_photonsKey {this, "photonsKey", "RefGamma", "the key for the photons"};

  private:
    SysReadHandle<xAOD::IParticleContainer> m_muonsHandle {
      this, "muons", "", "the muon container to use" };

  private:
    SysReadSelectionHandle m_muonsSelection {
      this, "muonsSelection", "", "the selection on the input muons" };

  private:
    Gaudi::Property<std::string> m_muonsKey {this, "muonsKey", "Muons", "the key for the muons"};

  private:
    SysReadHandle<xAOD::IParticleContainer> m_tausHandle {
      this, "taus", "", "the tau container to use" };

  private:
    SysReadSelectionHandle m_tausSelection {
      this, "tausSelection", "", "the selection on the input taus" };

  private:
    Gaudi::Property<std::string> m_tausKey {this, "tausKey", "RefTau", "the key for the taus"};

  private:
    SysReadHandle<xAOD::JetContainer> m_jetsHandle {
      this, "jets", "", "the jet collection we use"};

  private:
    SysHandleArray<SysReadHandle<xAOD::IParticleContainer>> m_invisHandles {
      this, "invisible", {}, "Any particles to treat as invisible."};

  private:
    Gaudi::Property<std::vector<std::string>> m_invisSelectionKeys {
      this, "invisibleSelection", {}, "the selection on the particles to be treated as invisible." };
  private:
    std::vector<SysReadSelectionHandle> m_invisSelections;

  private:
    Gaudi::Property<std::string> m_jetsKey {this, "jetsKey", "RefJet", "the key for the jets"};

  private:
    Gaudi::Property<std::string> m_softTermKey {this, "softTermKey", "PVSoftTrk", "the soft term key"};

  private:
    Gaudi::Property<bool> m_doTrackMet {this, "doTrackMet", false, "whether to use track-met instead of jet-met"};

  private:
    Gaudi::Property<bool> m_doJetJVT {this, "doJetJVT", true, "whether to do jet JVT"};

  private:
    SysWriteHandle<xAOD::MissingETContainer,xAOD::MissingETAuxContainer> m_metHandle {
      this, "met", "MissingET_%SYS%", "the met collection we produce"};

  private:
    Gaudi::Property<bool> m_switchTauMuonOrder {this, "switchTauMu", false, "whether to switch order of taus and muons"};

  private:
    Gaudi::Property<bool> m_evaluateNNMET {this, "evaluateNNMET", false,
      "run NN MET inference after rebuildJetMET; requires a maker tool that "
      "overrides evaluateNNMET (e.g. met::METNet). Skips the need for MetBuilderAlg."};

  private:
    Gaudi::Property<std::string> m_finalKey {this, "finalKey", "Final",
      "name of the total MET term written by evaluateNNMET"};
  };
}
 
#endif