Program Listing for File ParticleID.cpp

↰ Return to documentation for file (src/ParticleID.cpp)

#include "crpropa/ParticleID.h"

#include "HepPID/ParticleIDMethods.hh"
#include "HepPID/ParticleName.hh"
#include "kiss/convert.h"

#include <string>

namespace crpropa {

int nucleusId(int a, int z) {
        if (z < 0)
                throw std::runtime_error(
                                "crpropa::Nucleus: no nucleus with Z < 0, A=" + kiss::str(a) + " Z="
                                                + kiss::str(z));
        if (a < 1)
                throw std::runtime_error(
                                "crpropa::Nucleus: no nucleus with A < 1, A=" + kiss::str(a) + " Z="
                                                + kiss::str(z));
        if (a < z)
                throw std::runtime_error(
                                "crpropa::Nucleus: no nucleus with A < Z, A=" + kiss::str(a) + " Z="
                                                + kiss::str(z));
        return 1000000000 + z * 10000 + a * 10;
}

int chargeNumber(int id) {
        return HepPID::Z(id);
}

int massNumber(int id) {
        if (id == 2112)
                return 1;
        return HepPID::A(id);
}

bool isNucleus(int id) {
        if (id == 2112)
                return true; // consider neutron as nucleus
        return HepPID::isNucleus(id);
}

std::string convertIdToName(int id) {
        // handle a few extra cases that HepPID doesn't like
        if (id == 1000000010) // neutron
                id = 2112;
        if (id == -1000000010) // anti-neutron
                id = -2112;
        if (id == -1000010010) // anti-proton
                id = -2212;
        return HepPID::particleName(id);
}

}