met::ORMETMakerAlg Class Reference

#include <ORMETMakerAlg.h>

Inheritance diagram for met::ORMETMakerAlg:

Collaboration diagram for met::ORMETMakerAlg:

Public Member Functions
	ORMETMakerAlg (const std::string &name, ISvcLocator *pSvcLocator)
	Constructor with parameters: More...
virtual	~ORMETMakerAlg ()
	Destructor: More...
virtual StatusCode	initialize () override
	Athena algorithm's Hooks. More...
virtual StatusCode	execute () override
virtual StatusCode	finalize () override
virtual StatusCode	sysInitialize () override
	Override sysInitialize. More...
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies. More...
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
virtual StatusCode	sysStart () override
	Handle START transition. More...
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
MsgStream &	msg (const MSG::Level lvl) const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Protected Attributes
std::string	m_softclname
std::string	m_softtrkname
SG::ReadHandleKey< xAOD::ElectronContainer >	m_ElectronContainerKey
SG::ReadHandleKey< xAOD::PhotonContainer >	m_PhotonContainerKey
SG::ReadHandleKey< xAOD::TauJetContainer >	m_TauJetContainerKey
SG::ReadHandleKey< xAOD::MuonContainer >	m_MuonContainerKey
SG::ReadHandleKey< xAOD::JetContainer >	m_JetContainerKey
SG::ReadHandleKey< xAOD::MissingETContainer >	m_CoreMetKey
SG::WriteHandleKey< xAOD::MissingETContainer >	m_metKey
SG::ReadHandleKey< xAOD::MissingETAssociationMap >	m_metMapKey
bool	m_doTruthLep
ToolHandle< IMETMaker >	m_metmaker
	Athena configured tools. More...
ToolHandle< CP::IMuonSelectionTool >	m_muonSelTool
ToolHandle< IAsgElectronLikelihoodTool >	m_elecSelLHTool
ToolHandle< IAsgPhotonIsEMSelector >	m_photonSelIsEMTool
ToolHandle< TauAnalysisTools::ITauSelectionTool >	m_tauSelTool

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
	ORMETMakerAlg ()
	Default constructor: More...
virtual bool	accept (const xAOD::Electron *el) final
virtual bool	accept (const xAOD::Photon *ph) final
virtual bool	accept (const xAOD::TauJet *tau) final
virtual bool	accept (const xAOD::Muon *muon) final
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
std::string	m_soft
SG::ReadHandleKey< xAOD::MissingETAssociationMap >	m_ORMetMapKey
SG::WriteHandleKey< xAOD::PFOContainer >	m_chargedPFOContainerWriteHandleKey {this,"PFOChargedOutputName","OverlapRemovedCHSChargedParticleFlowObjects","WriteHandleKey for charged PFO"}
SG::WriteHandleKey< xAOD::PFOContainer >	m_neutralPFOContainerWriteHandleKey {this,"PFONeutralOutputName","OverlapRemovedCHSNeutralParticleFlowObjects","WriteHandleKey for charged PFO"}
SG::WriteHandleKey< xAOD::PFOContainer >	m_PFOContainerWriteHandleKey {this,"PFOOutputName","OverlapRemovedCHSParticleFlowObjects","WriteHandleKey for PFO"}
SG::ReadHandleKey< xAOD::PFOContainer >	m_inPFOKey {this, "InPFOKey", "", "ReadHandleKey for modified PFlow Objects"}
SG::WriteHandleKey< xAOD::PFOContainer >	m_outPFOKey {this, "OutPFOKey", "", "WriteHandleKey for modified PFlow Objects"}
bool	m_doRetrieveORconstit
bool	m_retainMuonConstit
bool	m_doORMet
double	m_electronPT
double	m_muonPT
double	m_photonPT
double	m_tauPT
double	m_electronETA
double	m_muonETA
double	m_photonETA
double	m_tauETA
bool	m_useUnmatched
bool	m_doJVT
DataObjIDColl	m_extendedExtraObjects
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 43 of file ORMETMakerAlg.h.

Member Typedef Documentation

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

ORMETMakerAlg() [1/2]

met::ORMETMakerAlg::ORMETMakerAlg	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Constructor with parameters:

Definition at line 34 of file ORMETMakerAlg.cxx.

    : METMakerAlg(name, pSvcLocator) 
 {
    declareProperty("ORMETMapName",    m_ORMetMapKey = "ORMETAssoc" );
    declareProperty( "DoRetrieveORconstit", m_doRetrieveORconstit = false              );
    declareProperty( "RetainMuonConstit", m_retainMuonConstit = false              );
    declareProperty( "DoORMet", m_doORMet = false              );
    //pT thresholds
    declareProperty("ElectronPT", m_electronPT = 10000, "pT for electrons");
    declareProperty("MuonPT", m_muonPT = 10000, "pT for muons");
    declareProperty("PhotonPT", m_photonPT = 10000, "pT for photons");
    declareProperty("TauPT", m_tauPT = 10000, "pT for photons");
    //eta thresholds
    declareProperty("ElectronETA", m_electronETA = 2.47, "eta for electrons");
    declareProperty("MuonETA", m_muonETA = 2.7, "eta for muons");
    declareProperty("PhotonETA", m_photonETA = 2.47, "eta for photons");
    declareProperty("TauETA", m_tauETA = 2.47, "eta for photons");
 
    declareProperty("UseUnmatched", m_useUnmatched = true, "Include or reject egamma with unmatched clusters");
    declareProperty("DoJVT", m_doJVT = false);
    declareProperty("Soft", m_soft = "Clus");
 
 
  }

~ORMETMakerAlg()

met::ORMETMakerAlg::~ORMETMakerAlg ( )

virtualdefault

Destructor:

ORMETMakerAlg() [2/2]

met::ORMETMakerAlg::ORMETMakerAlg ( )

private

Default constructor:

Member Function Documentation

accept() [1/4]

bool met::ORMETMakerAlg::accept ( const xAOD::Electron *  el )

finalprivatevirtual

Reimplemented from met::METMakerAlg.

Definition at line 351 of file ORMETMakerAlg.cxx.

                                                  {
    if( fabs(el->eta())>m_electronETA || el->pt()<m_electronPT) return false;
    return static_cast<bool> (m_elecSelLHTool->accept(el));
  }

accept() [2/4]

bool met::ORMETMakerAlg::accept ( const xAOD::Muon *  muon )

finalprivatevirtual

Reimplemented from met::METMakerAlg.

Definition at line 346 of file ORMETMakerAlg.cxx.

                                              {
        if( mu->pt()<m_muonPT && fabs(mu->eta())>m_muonETA)  return false;
    return static_cast<bool>(m_muonSelTool->accept(*mu));
  }

accept() [3/4]

bool met::ORMETMakerAlg::accept ( const xAOD::Photon *  ph )

finalprivatevirtual

Reimplemented from met::METMakerAlg.

Definition at line 356 of file ORMETMakerAlg.cxx.

                                                {
    if( !(ph->author()&20) || fabs(ph->eta())>m_photonETA || ph->pt()<m_photonPT) return false;
    return static_cast<bool> (m_photonSelIsEMTool->accept(ph));
  }

accept() [4/4]

bool met::ORMETMakerAlg::accept ( const xAOD::TauJet *  tau )

finalprivatevirtual

Reimplemented from met::METMakerAlg.

Definition at line 361 of file ORMETMakerAlg.cxx.

                                                 { 
    return static_cast<bool>(m_tauSelTool->accept( *tau ));
  }

declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl, const SG::VarHandleKeyArrayType &    )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl, const SG::VarHandleKeyType &    )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl, const SG::VarHandleType &    )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
  }

declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  t, const SG::NotHandleType &    )

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(t);
  }

declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleBase &  hndl, const std::string &  doc, const SG::VarHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

  {
    this->declare(hndl.vhKey());
    hndl.vhKey().setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleKey &  hndl, const std::string &  doc, const SG::VarHandleKeyType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

  {
    this->declare(hndl);
    hndl.setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleKeyArray &  hndArr, const std::string &  doc, const SG::VarHandleKeyArrayType &    )

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

  {
 
    // std::ostringstream ost;
    // ost << Algorithm::name() << " VHKA declareProp: " << name 
    //     << " size: " << hndArr.keys().size() 
    //     << " mode: " << hndArr.mode() 
    //     << "  vhka size: " << m_vhka.size()
    //     << "\n";
    // debug() << ost.str() << endmsg;
 
    hndArr.setOwner(this);
    m_vhka.push_back(&hndArr);
 
    Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
    if (p != 0) {
      p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
    } else {
      ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                    << name);
    }
 
    return p;
 
  }

declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc, const SG::NotHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(name, property, doc);
  }

declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc = "none"  )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

  {
    typedef typename SG::HandleClassifier<T>::type htype;
    return declareProperty (name, property, doc, htype());
  }

declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( Gaudi::Property< T > &  t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                 {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( ) const

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode met::ORMETMakerAlg::execute ( )

overridevirtual

MET

Jets

Electrons

Photons

Taus

Muons

PFOs for OR constituents

Reimplemented from met::METMakerAlg.

Definition at line 89 of file ORMETMakerAlg.cxx.

                                    {
    ATH_MSG_VERBOSE("Executing " << name() << "...");
    // Create a MissingETContainer with its aux store
    auto ctx = getContext();
    auto metHandle= SG::makeHandle (m_metKey,ctx);
    ATH_CHECK( metHandle.record (std::make_unique<xAOD::MissingETContainer>(),                      std::make_unique<xAOD::MissingETAuxContainer>()) );
    xAOD::MissingETContainer* newMet=metHandle.ptr();
 
 
    SG::ReadHandle<xAOD::MissingETAssociationMap> metMap(m_metMapKey);
    if (!metMap.isValid()) {
      ATH_MSG_WARNING("Unable to retrieve MissingETAssociationMap: " << m_metMapKey.key());
      return StatusCode::FAILURE;
    }
 
    MissingETAssociationHelper metHelper(&(*metMap));
 
    SG::ReadHandle<xAOD::MissingETAssociationMap> ORMetMap(m_ORMetMapKey);
    if (m_doORMet && !ORMetMap.isValid()) {
      ATH_MSG_WARNING("Unable to retrieve MissingETAssociationMap: " << m_ORMetMapKey.key());
      return StatusCode::FAILURE;
    }
 
    
    MissingETAssociationHelper ORMetHelper((m_doORMet ? &(*ORMetMap) : nullptr));
 
    // Retrieve containers ***********************************************
 
    SG::ReadHandle<xAOD::MissingETContainer> coreMet(m_CoreMetKey);
    if (!coreMet.isValid()) {
      ATH_MSG_WARNING("Unable to retrieve MissingETContainer: " << m_CoreMetKey.key());
      return StatusCode::FAILURE;
    }
 
    SG::ReadHandle<xAOD::JetContainer> Jets(m_JetContainerKey);
    if (!Jets.isValid()) {
      ATH_MSG_WARNING("Unable to retrieve JetContainer: " << Jets.key());
      return StatusCode::FAILURE;
    }
 
    SG::ReadHandle<xAOD::ElectronContainer> Electrons(m_ElectronContainerKey);
    if (!Electrons.isValid()) {
      ATH_MSG_WARNING("Unable to retrieve ElectronContainer: " << Electrons.key());
      return StatusCode::FAILURE;
    }
 
    SG::ReadHandle<xAOD::PhotonContainer> Gamma(m_PhotonContainerKey);
    if (!Gamma.isValid()) {
      ATH_MSG_WARNING("Unable to retrieve GammaContainer: " << Gamma.key());
      return StatusCode::FAILURE;
    }
 
    SG::ReadHandle<xAOD::TauJetContainer> TauJets(m_TauJetContainerKey);
    if (!TauJets.isValid()) {
      ATH_MSG_WARNING("Unable to retrieve TauJetContainer: " << TauJets.key());
      return StatusCode::FAILURE;
    }
 
    SG::ReadHandle<xAOD::MuonContainer> Muons(m_MuonContainerKey);
    if (!Muons.isValid()) {
      ATH_MSG_WARNING("Unable to retrieve MuonContainer: "  << Muons.key());
      return StatusCode::FAILURE;
    }
 
    SG::ReadHandle<xAOD::PFOContainer> PFOs("CHSParticleFlowObjects");
    if (!PFOs.isValid() && m_doRetrieveORconstit ) {
      ATH_MSG_WARNING("Unable to retrieve PFOContainer: " << PFOs.key());
      return StatusCode::FAILURE;
    }
 
    // Select and flag objects for final MET building ***************************
 
    MissingETBase::UsageHandler::Policy objScale = MissingETBase::UsageHandler::PhysicsObject;
    if(m_doTruthLep) objScale = MissingETBase::UsageHandler::TruthParticle;
    // Electrons
    if(!m_ElectronContainerKey.empty()) {
      ConstDataVector<ElectronContainer> metElectrons(SG::VIEW_ELEMENTS);
      for(const auto *const el : *Electrons) {
        if(accept(el)) {
          metElectrons.push_back(el);
 
        }
      }
      if( m_metmaker->rebuildMET("RefEle", xAOD::Type::Electron, newMet,
                     metElectrons.asDataVector(),
                     metHelper, objScale).isFailure() ) {
        ATH_MSG_WARNING("Failed to build electron term.");
      }
      ATH_MSG_DEBUG("Selected " << metElectrons.size() << " MET electrons. "
                << acc_constitObjLinks(*(*newMet)["RefEle"]).size() << " are non-overlapping.");
    }
 
    // Photons
    if(!m_PhotonContainerKey.empty()) {
      ConstDataVector<PhotonContainer> metPhotons(SG::VIEW_ELEMENTS);
      for(const auto *const ph : *Gamma) {
        if(accept(ph)) {
          metPhotons.push_back(ph);
        }
      }
      if( m_metmaker->rebuildMET("RefGamma", xAOD::Type::Photon, newMet,
                     metPhotons.asDataVector(),
                     metHelper, objScale).isFailure() ) {
        ATH_MSG_WARNING("Failed to build photon term.");
      }
      ATH_MSG_DEBUG("Selected " << metPhotons.size() << " MET photons. "
                << acc_constitObjLinks(*(*newMet)["RefGamma"]).size() << " are non-overlapping.");
    }
 
    // Taus
    if(!m_TauJetContainerKey.empty()) {
      ConstDataVector<TauJetContainer> metTaus(SG::VIEW_ELEMENTS);
      for(const auto *const tau : *TauJets) {
        if(accept(tau)) {
          metTaus.push_back(tau);
        }
      }
      if( m_metmaker->rebuildMET("RefTau", xAOD::Type::Tau, newMet,
                     metTaus.asDataVector(),
                     metHelper, objScale).isFailure() ){
        ATH_MSG_WARNING("Failed to build tau term.");
      }
      ATH_MSG_DEBUG("Selected " << metTaus.size() << " MET taus. "
                << acc_constitObjLinks(*(*newMet)["RefTau"]).size() << " are non-overlapping.");
    }
 
    // Muons
    ConstDataVector<MuonContainer> metMuons(SG::VIEW_ELEMENTS);
    if(!m_MuonContainerKey.empty()) {
      for(const auto *const mu : *Muons) {
        if(accept(mu)) {
          metMuons.push_back(mu);
        }
      }
      
      if(m_doTruthLep) objScale = MissingETBase::UsageHandler::OnlyTrack;
 
      if (m_doORMet && m_retainMuonConstit) {
        if( m_metmaker->rebuildMET("Muons", xAOD::Type::Muon, newMet,
                     metMuons.asDataVector(),
                     ORMetHelper, objScale).isFailure() ) {
          ATH_MSG_WARNING("Failed to build muon term.");
        }
      } else {
        if( m_metmaker->rebuildMET("Muons", xAOD::Type::Muon, newMet,
                     metMuons.asDataVector(),
                     metHelper, objScale).isFailure() ) {
          ATH_MSG_WARNING("Failed to build muon term.");
        }
      }
      ATH_MSG_DEBUG("Selected " << metMuons.size() << " MET muons. "
                << acc_constitObjLinks(*(*newMet)["Muons"]).size() << " are non-overlapping.");
    }
 
    if (m_doRetrieveORconstit) {
 
      /*** Fill overlapRemovedPFlowObjects ***/
      SG::WriteHandle<xAOD::PFOContainer> PFOContainerWriteHandle(m_PFOContainerWriteHandleKey,ctx); 
      ATH_CHECK(PFOContainerWriteHandle.record(std::make_unique<xAOD::PFOContainer>(SG::VIEW_ELEMENTS)));
 
      ATH_MSG_DEBUG("Retrieve OR constituents");
      ConstDataVector<PFOContainer> met_PFO(SG::VIEW_ELEMENTS);
      for(const auto *const pfo : *PFOs) {met_PFO.push_back(pfo);}
 
      const xAOD::PFOContainer *OR_pfos; 
   
      if (!metMuons.empty() && m_retainMuonConstit) {ATH_CHECK(m_metmaker->retrieveOverlapRemovedConstituents(met_PFO.asDataVector(),metHelper, &OR_pfos, true, metMuons.asDataVector()));}
      else {ATH_CHECK(m_metmaker->retrieveOverlapRemovedConstituents(met_PFO.asDataVector(),metHelper, &OR_pfos, false));}
 
      *PFOContainerWriteHandle=*OR_pfos;
 
     /**** If want to write CHScharged/neutralOverlapRemovedPFlowObjects container to AOD ****
      SG::WriteHandle<xAOD::PFOContainer> PFOContainerWriteHandle(m_PFOContainerWriteHandleKey,ctx); 
      ATH_CHECK(PFOContainerWriteHandle.record(std::make_unique<xAOD::PFOContainer>(SG::VIEW_ELEMENTS)));
      //cPFO
      SG::ReadHandle<xAOD::PFOContainer> cPFOs("CHSChargedParticleFlowObjects");
      if (!cPFOs.isValid()) {
        ATH_MSG_WARNING("Unable to retrieve PFOContainer: " << cPFOs.key());
        return StatusCode::FAILURE;
      }
      //nPFO
      SG::ReadHandle<xAOD::PFOContainer> nPFOs("CHSNeutralParticleFlowObjects");
      if (!nPFOs.isValid()) {
        ATH_MSG_WARNING("Unable to retrieve PFOContainer: " << nPFOs.key());
         return StatusCode::FAILURE;
      }
      SG::WriteHandle<xAOD::PFOContainer> chargedPFOContainerWriteHandle(m_chargedPFOContainerWriteHandleKey,ctx); 
      SG::WriteHandle<xAOD::PFOContainer> neutralPFOContainerWriteHandle(m_neutralPFOContainerWriteHandleKey,ctx); 
      ATH_CHECK(chargedPFOContainerWriteHandle.record(std::make_unique<xAOD::PFOContainer>(),std::make_unique<xAOD::PFOAuxContainer>())); 
      ATH_CHECK(neutralPFOContainerWriteHandle.record(std::make_unique<xAOD::PFOContainer>(),std::make_unique<xAOD::PFOAuxContainer>())); 
      
 
      ConstDataVector<PFOContainer> met_cPFO(SG::VIEW_ELEMENTS);
      for(const auto& pfo : *cPFOs) {met_cPFO.push_back(pfo);}
      ConstDataVector<PFOContainer> met_nPFO(SG::VIEW_ELEMENTS);
      for(const auto& pfo : *nPFOs) {met_nPFO.push_back(pfo);}
 
      xAOD::PFOContainer* or_cPFO = chargedPFOContainerWriteHandle.ptr();
      xAOD::PFOContainer* or_nPFO = neutralPFOContainerWriteHandle.ptr();
 
      if (metMuons.size()>0 && m_retainMuonConstit) {ATH_CHECK(m_metmaker->retrieveOverlapRemovedConstituents(met_cPFO.asDataVector(), met_nPFO.asDataVector(),metHelper,or_cPFO,or_nPFO,true, metMuons.asDataVector()));}
      else {ATH_CHECK(m_metmaker->retrieveOverlapRemovedConstituents(met_cPFO.asDataVector(), met_nPFO.asDataVector(),metHelper,or_cPFO,or_nPFO,false));}
 
    (*PFOContainerWriteHandle).assign((*neutralPFOContainerWriteHandle).begin(), (*neutralPFOContainerWriteHandle).end());
    (*PFOContainerWriteHandle).insert((*PFOContainerWriteHandle).end(),
          (*chargedPFOContainerWriteHandle).begin(), 
          (*chargedPFOContainerWriteHandle).end());***/
 
     for (auto tmp_constit : *PFOContainerWriteHandle){ATH_MSG_VERBOSE("ORMETMaker constit with index " << tmp_constit->index() << ", charge " << tmp_constit->charge()<< " pT " << tmp_constit->pt() << " in OverlapRemovedCHSParticleFlowObjects");}
 
    }
 
    m_doJVT= (m_soft != "Clus");
    if (m_doORMet) {
      if( m_metmaker->rebuildJetMET("RefJet", (m_soft=="Clus" ? m_softclname : m_softtrkname), newMet,
                  Jets.cptr(), coreMet.cptr(), ORMetHelper, m_doJVT ).isFailure() ) {
        ATH_MSG_WARNING("Failed to build jet and soft terms.");
      } 
      ATH_MSG_DEBUG("Of " << Jets.cptr()->size()  << " jets, "
          << acc_constitObjLinks(*(*newMet)["RefJet"]).size() << " are non-overlapping, "
          << acc_constitObjLinks(*(*newMet)[(m_soft=="Clus" ? m_softclname : m_softtrkname)]).size() << " are soft");
    } else {
      if( m_metmaker->rebuildJetMET("RefJet", (m_soft=="Clus" ? m_softclname : m_softtrkname), newMet,
                  Jets.cptr(), coreMet.cptr(), metHelper, m_doJVT ).isFailure() ) {
        ATH_MSG_WARNING("Failed to build jet and soft terms.");
      }
      ATH_MSG_DEBUG("Of " << Jets.cptr()->size()  << " jets, "
          << acc_constitObjLinks(*(*newMet)["RefJet"]).size() << " are non-overlapping, "
          << acc_constitObjLinks(*(*newMet)[(m_soft=="Clus" ? m_softclname : m_softtrkname)]).size() << " are soft");
    }
 
 
    ATH_CHECK( buildMETSum("Final"+m_soft, 
                newMet,
                (m_soft=="Clus")? (*newMet)["SoftClus"]->source() : (*newMet)["PVSoftTrk"]->source()) 
         ); 
 
    xAOD::MissingET* metSum = (*newMet)["Final"+m_soft];
    ATH_MSG_DEBUG("MET sum = " << 0.001 * metSum->met() << " GeV");
    ATH_MSG_DEBUG("MET jet = " << 0.001 * (*newMet)["RefJet"]->met() << " GeV");
    ATH_MSG_DEBUG("MET mu = " << 0.001 * (*newMet)["Muons"]->met() << " GeV");
    ATH_MSG_DEBUG("MET el = " << 0.001 * (*newMet)["RefEle"]->met() << " GeV");
    ATH_MSG_DEBUG("MET ph = " << 0.001 * (*newMet)["RefGamma"]->met() << " GeV");
    ATH_MSG_DEBUG("MET tau = " << 0.001 * (*newMet)["RefTau"]->met() << " GeV");
    ATH_MSG_DEBUG("MET soft = " << 0.001 * (*newMet)["Soft"+m_soft]->met() << " GeV");
 
    return StatusCode::SUCCESS;
  }

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return Algorithm::extraOutputDeps();
}

finalize()

StatusCode met::ORMETMakerAlg::finalize ( )

overridevirtual

Reimplemented from met::METMakerAlg.

Definition at line 82 of file ORMETMakerAlg.cxx.

                                     {
    ATH_MSG_INFO ("Finalizing " << name() << "...");
    return StatusCode::SUCCESS;
  }

initialize()

StatusCode met::ORMETMakerAlg::initialize ( )

overridevirtual

Athena algorithm's Hooks.

Reimplemented from met::METMakerAlg.

Definition at line 66 of file ORMETMakerAlg.cxx.

                                       {
 
    ATH_CHECK( METMakerAlg::initialize() );
    ATH_MSG_INFO("Initializing " << name() << "...");
    ATH_MSG_INFO("Retrieving tools...");
 
    ATH_CHECK(m_ORMetMapKey.initialize());
    ATH_CHECK(m_chargedPFOContainerWriteHandleKey.initialize()); 
    ATH_CHECK(m_neutralPFOContainerWriteHandleKey.initialize()); 
    ATH_CHECK(m_PFOContainerWriteHandleKey.initialize()); 
 
    return StatusCode::SUCCESS;
  }

inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

msg() [1/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msg() [2/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                  {
    return this->msgStream(lvl);
  }

msgLvl()

bool AthCommonMsg< Algorithm >::msgLvl ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< Algorithm > >::renounce ( T &  h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::renounceArray ( SG::VarHandleKeyArray &  handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, PyAthena::Alg, and AthHistogramAlgorithm.

Definition at line 66 of file AthAlgorithm.cxx.

                                       {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;
}

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

updateVHKA()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase &  )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_chargedPFOContainerWriteHandleKey

SG::WriteHandleKey<xAOD::PFOContainer> met::ORMETMakerAlg::m_chargedPFOContainerWriteHandleKey {this,"PFOChargedOutputName","OverlapRemovedCHSChargedParticleFlowObjects","WriteHandleKey for charged PFO"}

private

Definition at line 73 of file ORMETMakerAlg.h.

m_CoreMetKey

SG::ReadHandleKey<xAOD::MissingETContainer> met::METMakerAlg::m_CoreMetKey

protectedinherited

Definition at line 77 of file METMakerAlg.h.

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_doJVT

bool met::ORMETMakerAlg::m_doJVT

private

Definition at line 98 of file ORMETMakerAlg.h.

m_doORMet

bool met::ORMETMakerAlg::m_doORMet

private

Definition at line 83 of file ORMETMakerAlg.h.

m_doRetrieveORconstit

bool met::ORMETMakerAlg::m_doRetrieveORconstit

private

Definition at line 81 of file ORMETMakerAlg.h.

m_doTruthLep

bool met::METMakerAlg::m_doTruthLep

protectedinherited

Definition at line 83 of file METMakerAlg.h.

m_elecSelLHTool

ToolHandle<IAsgElectronLikelihoodTool> met::METMakerAlg::m_elecSelLHTool

protectedinherited

Definition at line 89 of file METMakerAlg.h.

m_ElectronContainerKey

SG::ReadHandleKey<xAOD::ElectronContainer> met::METMakerAlg::m_ElectronContainerKey

protectedinherited

Definition at line 71 of file METMakerAlg.h.

m_electronETA

double met::ORMETMakerAlg::m_electronETA

private

Definition at line 91 of file ORMETMakerAlg.h.

m_electronPT

double met::ORMETMakerAlg::m_electronPT

private

Definition at line 86 of file ORMETMakerAlg.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects

privateinherited

Definition at line 79 of file AthAlgorithm.h.

m_inPFOKey

SG::ReadHandleKey<xAOD::PFOContainer> met::ORMETMakerAlg::m_inPFOKey {this, "InPFOKey", "", "ReadHandleKey for modified PFlow Objects"}

private

Definition at line 78 of file ORMETMakerAlg.h.

m_JetContainerKey

SG::ReadHandleKey<xAOD::JetContainer> met::METMakerAlg::m_JetContainerKey

protectedinherited

Definition at line 75 of file METMakerAlg.h.

m_metKey

SG::WriteHandleKey<xAOD::MissingETContainer> met::METMakerAlg::m_metKey

protectedinherited

Definition at line 79 of file METMakerAlg.h.

m_metmaker

ToolHandle<IMETMaker> met::METMakerAlg::m_metmaker

protectedinherited

Athena configured tools.

Definition at line 86 of file METMakerAlg.h.

m_metMapKey

SG::ReadHandleKey<xAOD::MissingETAssociationMap> met::METMakerAlg::m_metMapKey

protectedinherited

Definition at line 80 of file METMakerAlg.h.

m_MuonContainerKey

SG::ReadHandleKey<xAOD::MuonContainer> met::METMakerAlg::m_MuonContainerKey

protectedinherited

Definition at line 74 of file METMakerAlg.h.

m_muonETA

double met::ORMETMakerAlg::m_muonETA

private

Definition at line 92 of file ORMETMakerAlg.h.

m_muonPT

double met::ORMETMakerAlg::m_muonPT

private

Definition at line 87 of file ORMETMakerAlg.h.

m_muonSelTool

ToolHandle<CP::IMuonSelectionTool> met::METMakerAlg::m_muonSelTool

protectedinherited

Definition at line 88 of file METMakerAlg.h.

m_neutralPFOContainerWriteHandleKey

SG::WriteHandleKey<xAOD::PFOContainer> met::ORMETMakerAlg::m_neutralPFOContainerWriteHandleKey {this,"PFONeutralOutputName","OverlapRemovedCHSNeutralParticleFlowObjects","WriteHandleKey for charged PFO"}

private

Definition at line 74 of file ORMETMakerAlg.h.

m_ORMetMapKey

SG::ReadHandleKey<xAOD::MissingETAssociationMap> met::ORMETMakerAlg::m_ORMetMapKey

private

Definition at line 71 of file ORMETMakerAlg.h.

m_outPFOKey

SG::WriteHandleKey<xAOD::PFOContainer> met::ORMETMakerAlg::m_outPFOKey {this, "OutPFOKey", "", "WriteHandleKey for modified PFlow Objects"}

private

Definition at line 79 of file ORMETMakerAlg.h.

m_PFOContainerWriteHandleKey

SG::WriteHandleKey<xAOD::PFOContainer> met::ORMETMakerAlg::m_PFOContainerWriteHandleKey {this,"PFOOutputName","OverlapRemovedCHSParticleFlowObjects","WriteHandleKey for PFO"}

private

Definition at line 75 of file ORMETMakerAlg.h.

m_PhotonContainerKey

SG::ReadHandleKey<xAOD::PhotonContainer> met::METMakerAlg::m_PhotonContainerKey

protectedinherited

Definition at line 72 of file METMakerAlg.h.

m_photonETA

double met::ORMETMakerAlg::m_photonETA

private

Definition at line 93 of file ORMETMakerAlg.h.

m_photonPT

double met::ORMETMakerAlg::m_photonPT

private

Definition at line 88 of file ORMETMakerAlg.h.

m_photonSelIsEMTool

ToolHandle<IAsgPhotonIsEMSelector> met::METMakerAlg::m_photonSelIsEMTool

protectedinherited

Definition at line 90 of file METMakerAlg.h.

m_retainMuonConstit

bool met::ORMETMakerAlg::m_retainMuonConstit

private

Definition at line 82 of file ORMETMakerAlg.h.

m_soft

std::string met::ORMETMakerAlg::m_soft

private

Definition at line 69 of file ORMETMakerAlg.h.

m_softclname

std::string met::METMakerAlg::m_softclname

protectedinherited

Definition at line 67 of file METMakerAlg.h.

m_softtrkname

std::string met::METMakerAlg::m_softtrkname

protectedinherited

Definition at line 68 of file METMakerAlg.h.

m_tauETA

double met::ORMETMakerAlg::m_tauETA

private

Definition at line 94 of file ORMETMakerAlg.h.

m_TauJetContainerKey

SG::ReadHandleKey<xAOD::TauJetContainer> met::METMakerAlg::m_TauJetContainerKey

protectedinherited

Definition at line 73 of file METMakerAlg.h.

m_tauPT

double met::ORMETMakerAlg::m_tauPT

private

Definition at line 89 of file ORMETMakerAlg.h.

m_tauSelTool

ToolHandle<TauAnalysisTools::ITauSelectionTool> met::METMakerAlg::m_tauSelTool

protectedinherited

Definition at line 91 of file METMakerAlg.h.

m_useUnmatched

bool met::ORMETMakerAlg::m_useUnmatched

private

Definition at line 97 of file ORMETMakerAlg.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

The documentation for this class was generated from the following files:

• ORMETMakerAlg.h
• ORMETMakerAlg.cxx