egammaGenParticleAssociationTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
// $Id$
 
 
#include "egammaGenParticleAssociationTool.h"
#include "xAODEgamma/Egamma.h"
#include "xAODEgamma/Electron.h"
#include "xAODEgamma/Photon.h"
#include "xAODTruth/TruthParticle.h"
#include "FourMomUtils/P4Helpers.h"
#include "AthenaKernel/errorcheck.h"
#include "AtlasHepMC/GenEvent.h"
#include <cstdlib>
 
 
namespace D3PD {
 

egammaGenParticleAssociationTool::egammaGenParticleAssociationTool
  (const std::string& type,
   const std::string& name,
   const IInterface* parent)
    : Base (type, name, parent),
      m_classifier ("MCTruthClassifier")
{
  declareProperty ("Classifier", m_classifier, "Classifier tool instance.");
  declareProperty ("DRVar", m_drvar,
                   "If not empty, the variable name to use for DR.");
 
  m_dr = 0;
}
 

StatusCode egammaGenParticleAssociationTool::initialize()
{
  CHECK( Base::initialize() );
  CHECK( m_classifier.retrieve() );
  return StatusCode::SUCCESS;
}
 

StatusCode egammaGenParticleAssociationTool::book ()
{
  if (!m_drvar.empty())
    CHECK( addVariable (m_drvar, m_dr) );
  return StatusCode::SUCCESS;
}
 

const xAOD::TruthParticle*
egammaGenParticleAssociationTool::get (const xAOD::Egamma& p)
{
  const xAOD::TruthParticle* out = nullptr;
  if (const xAOD::Electron* q =
      dynamic_cast<const xAOD::Electron*> (&p))
  {
    std::tie(std::ignore, std::ignore, out, std::ignore) = m_classifier->particleTruthClassifier_full (q);
  }
  else if (const xAOD::Photon* q =
           dynamic_cast<const xAOD::Photon*> (&p))
  {
    std::tie(std::ignore, std::ignore, out, std::ignore) = m_classifier->particleTruthClassifier_full (q);
  }
  else
    std::abort();
 
  if (!m_drvar.empty()) {
    if (out)
      *m_dr = p.p4().DeltaR (out->p4());
    else
      *m_dr = -999;
  }
 
  return out;
}
 
 
} // namespace D3PD