MissingET_v1.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
 
#include "xAODMissingET/versions/MissingET_v1.h"
#include "xAODMissingET/versions/MissingETBase.h"
 
using namespace xAOD;
 
// Static accessors
const SG::AuxElement::Accessor<float> MissingET_v1::m_acc_mpx("mpx");
const SG::AuxElement::Accessor<float> MissingET_v1::m_acc_mpy("mpy");
const SG::AuxElement::Accessor<float> MissingET_v1::m_acc_sumet("sumet");
namespace {
   // register attributes' types early (before I/O) for schema evolution
   static const SG::AuxElement::Accessor<MissingETBase::Types::bitmask_t> dummy("source");
}
 
MissingET_v1::MissingET_v1( bool createStore )
  : SG::AuxElement() {
 
   if( createStore ) {
      createPrivateStore();
      setName(this->name());
   }
}
 
MissingET_v1::MissingET_v1( const std::string& name,
                            MissingETBase::Types::bitmask_t src )
   : SG::AuxElement() {
 
   createPrivateStore();
   setName( name );
   setSource( src );
 
}
 
MissingET_v1::MissingET_v1( float mpx, float mpy, float sumet,
                            const std::string& name,
                            MissingETBase::Types::bitmask_t src )
   : SG::AuxElement() {
 
   createPrivateStore();
   setName( name );
   setSource( src );
   this->f_mpx() = mpx;
   this->f_mpy() = mpy;
   this->f_sumet() = sumet;
}
 
MissingET_v1::MissingET_v1( const IParticle* particle, const std::string& name,
                            MissingETBase::Types::bitmask_t src )
   : SG::AuxElement() {
 
   createPrivateStore();
   setName( name );
   setSource( src );
   add( particle->p4().Px(), particle->p4().Py(), particle->pt() );
}
 
// CHECK use makePrivateStore to control that behaviour ??
MissingET_v1::MissingET_v1(const MissingET_v1& met)
  : SG::AuxElement(met) {
  
  this->makePrivateStore(&met);
  updateHash();
}
 
MissingET_v1::~MissingET_v1()
{ }
 
// Note: the add(...) methods are virtual and should therefore not be inlined.
 
void MissingET_v1::add(const IParticle* particle)
{
  // retrieve kinematics
  EXTRACT_PX( float, px, (*particle) );
  EXTRACT_PY( float, py, (*particle) );
  // add
  this->add(px,py,particle->pt());
}
 
void MissingET_v1::add(const IParticle* particle,float scale)
{
  // retrieve kinematics
  EXTRACT_SCALED_PX( float, px, (*particle), scale);
  EXTRACT_SCALED_PY( float ,py, (*particle), scale);
  float pt(particle->pt()*scale);
  // add
  this->add(px,py,pt);
}
 
void xAOD::MissingET_v1::add(float px,float py, float pt)
{ this->f_mpx() -= px; this->f_mpy() -= py; this->f_sumet() += pt; }
 
MissingET_v1& MissingET_v1::operator=(const MissingET_v1& met)
{
  setName( met.name() );
  setSource( met.source() );
  this->f_mpx() = met.mpx();
  this->f_mpy() = met.mpy();
  this->f_sumet() = met.sumet();
  return *this;
}
 
MissingET_v1& MissingET_v1::operator-=(const IParticle* particle)
{
  EXTRACT_PX( float, px, (*particle) );
  EXTRACT_PY( float, py, (*particle) );
  this->add(-px,-py,-particle->pt());
  return *this;
}
 
void MissingET_v1::createPrivateStore() {
 
   makePrivateStore();
   setName( MissingETBase::Strings::unknownString() );
   setSource( MissingETBase::Source::unknown() );
   clear();
}