HepMC Reference Documentation HepMC

ParticleData.cc

Go to the documentation of this file.

// Matt.Dobbs@Cern.CH, September 1999
//
// particle properties common to all particles of a given PDG id

#include "HepMC/ParticleData.h"
#include <cstdio>       // needed for formatted output using sprintf 

namespace HepMC {

    ParticleData::ParticleData( std::string name, int id, double charge, 
                                double mass, double clifetime, double spin ) :
        m_name(name), m_pdg_id(id), m_mass(mass), m_clifetime(clifetime)
    {
        //
        set_charge(charge);
        set_spin(spin);
        s_counter++;
    }

    ParticleData::ParticleData( const char* name, int id, double charge, 
                                double mass, double clifetime, double spin ):
        m_name(name), m_pdg_id(id), m_mass(mass), m_clifetime(clifetime)
    {
        //
        set_charge(charge);
        set_spin(spin); 
        s_counter++;
    }

    ParticleData::~ParticleData() {
        s_counter--;
    }

    void ParticleData::print( std::ostream& ostr ) const {
        ostr << "ParticleData: " << name() << "\t"
             << " ID[pdg]:" << pdg_id()
             << " Charge[e+]:" << charge()
             << " Mass:" << mass()
             << " Tau:" << clifetime()
             << " J:" << spin() << std::endl;
    }

    // access methods //

    int ParticleData::model_independent_pdg_id_() const { 
        int id = m_pdg_id;
        if ( id/1000000 >=1 && id/1000000 <= 4 ) id %= 1000000;
        return id;
    }

    double ParticleData::width() const { 
        double width;
        if ( m_clifetime > 0 ) { 
            width = HepMC_hbarc/m_clifetime;
        } else if ( m_clifetime == 0 ) { 
            width = -1;
        } else { 
            width = 0; 
        }
        return width;
    }

    // Static  //
    unsigned int ParticleData::counter() { return s_counter; }
    unsigned int ParticleData::s_counter = 0; 

    // Friends //

    std::ostream& operator<<( std::ostream& ostr, const ParticleData& pdata ) {
        char outline[80];
        sprintf( outline,"%+9d%21s%+6.2f%19.11e%19.11e%5.1f",
                 pdata.pdg_id(),
                 pdata.name().substr(0,21).c_str(),
                 pdata.charge(),
                 pdata.mass(),
                 pdata.clifetime(),
                 pdata.spin() );
        return ostr << outline;
    }

    // Related Methods  //

    double clifetime_from_width( double width ) {
        // i.e. new ParticleData(
        if ( width > 0 ) return HepMC_hbarc/width;
        if ( width == 0. ) return -1.;
        return 0.;
    }
    
} // HepMC

---