---

ChargeService.cxx

Go to the documentation of this file.

//=================================================
//
//
// find particle charge
//
//=================================================
#include "AtlfastUtils/ChargeService.h"
namespace Atlfast {
  using std::abs;

  ChargeService::ChargeService(){
    HepMC::IO_PDG_ParticleDataTable pdg_io("PDGTABLE");
    m_particleDataTable = *pdg_io.read_particle_data_table();
    m_particleDataTable.make_antiparticles_from_particles();
  }
  ChargeService::ChargeService(const ChargeService & other):
    m_particleDataTable(other.m_particleDataTable){}
  
  double ChargeService::operator()(const HepMC::GenParticle* p)const {
    
    //does particle exist in table only positive ids for neutrals !?
    HepMC::ParticleData* ap = m_particleDataTable.find( abs( p->pdg_id() ) );
    
    if(!ap){ 
      std::cout << "ChargeService WARNING: abs(id) " 
           << abs(p->pdg_id())
           << " is not in particle data table"<<std::endl;
      return 0.;
    }
    
    //if in table, and the charge is 0, return 0. 
    double c = ap->charge();
    if(c==0) return 0.;
    
    
    //if in table, and the charge is not 0, use signed id. 
    
    HepMC::ParticleData* pp  = m_particleDataTable.find( p->pdg_id() );
    
    if(!pp){ 
      std::cout << "ChargeService WARNING: id " 
           << p->pdg_id()<< " is not in particle data table"<<std::endl;
      return 0.;
    }
    
    return pp->charge();
    
  }
 
  double ChargeService::operator()(const ReconstructedParticle* p)const {
    return (*this)(p->truth());
  }
 
}

---