00001 #include <vector> 00002 00003 #include "TruthHelper/IsGenNonInteracting.h" 00004 #include "HepMC/GenParticle.h" 00005 using std::vector; 00006 using std::abs; 00007 00008 //**************************************************************** 00009 //* IsGenNonInteracting * 00010 //**************************************************************** 00011 IsGenNonInteracting::IsGenNonInteracting() {} 00012 00013 bool IsGenNonInteracting::operator()( const Particle* const p ) const { 00014 return (abs(p->pdg_id())==12 || abs(p->pdg_id())==14 00015 || abs(p->pdg_id())==16 || 00016 (abs(p->pdg_id())==1000039 && p->status()%1000==1 ) || 00017 (abs(p->pdg_id())==1000022 && p->status()%1000==1 ) ) ? true:false; 00018 } 00019 00020 bool IsGenNonInteracting::operator() ( const Particle& p ) const { 00021 return this->operator()(&p); 00022 } 00023 GenIMCselector* IsGenNonInteracting::create() const {return new IsGenNonInteracting(*this);} 00024