/cvmfs/atlas.cern.ch/repo/sw/ASG/AnalysisBase/2.4.31/TruthUtils/TruthUtils/HepMCHelpers.h

Go to the documentation of this file.

#pragma once


#include "TruthUtils/PIDHelpers.h"
#include "TruthUtils/TruthParticleHelpers.h"
#include "MCUtils/HepMCUtils.h"
// #include "MCUtils/Clustering.h"

// Common imports from external namespaces
using HepMC::GenEvent;
using HepMC::GenParticle;
using HepMC::GenVertex;
//using fastjet::PseudoJet;
//using fastjet::ClusterSequence;

// Alias BOOST_FOREACH as foreach in this semi-safe way
#include <boost/foreach.hpp>
#ifndef foreach
namespace boost { namespace BOOST_FOREACH = foreach; }
#define foreach BOOST_FOREACH
#endif


namespace MC {

  // Use the MCUtils and HEPUtils functions as if they were defined in the ATLAS MC and MC::PID namespaces
  using namespace MCUtils;
  using namespace HEPUtils;




  inline bool isGenStable(const HepMC::GenParticle* p) {
    return isGenStable(p->status(), p->barcode());
  }




  inline bool isSimStable(const HepMC::GenParticle* p) {
    if (p->status() != 1) return false;
    if (isGenStable(p)) return p->end_vertex() == NULL;
    return true;
  }

  inline bool isGenSimulStable(const HepMC::GenParticle* p) {
    return isSimStable(p);
  }


  inline bool isNonInteracting(const HepMC::GenParticle* p) {
    return MC::isNonInteracting(p->pdg_id()); //< From TruthUtils/PIDHelpers.h
  }

  // inline bool isSimInteracting(int pid) {
  //   if (! MC::isGenStable(pid)) return false; //skip particles which the simulation would not see
  //   return !MC::isNonInteracting(pid);
  // }
  inline bool isSimInteracting(const HepMC::GenParticle* p) {
    if (! MC::isGenStable(p)) return false; //skip particles which the simulation would not see
    return !MC::isNonInteracting(p);
  }

  inline bool isGenInteracting(const HepMC::GenParticle* p) {
    return isSimInteracting(p);
  }



}

---