HeavyIon.h

Go to the documentation of this file.

/*
 *  CepGen: a central exclusive processes event generator
 *  Copyright (C) 2013-2022  Laurent Forthomme
 *
 *  This program is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#ifndef CepGen_Physics_HeavyIon_h
#define CepGen_Physics_HeavyIon_h
 
#include <ostream>
 
#include "CepGen/Physics/ParticleProperties.h"
 
namespace cepgen {
  enum class Element {
    invalid = -1,
    neutron = 0,
    H = 1,
    C = 6,
    O = 8,
    Al = 13,
    Cu = 29,
    Xe = 54,
    Au = 79,
    Pb = 82,
    U = 92
  };
  std::ostream& operator<<(std::ostream& os, const Element& elem);
 
  struct HeavyIon {
    explicit HeavyIon(unsigned short, const Element&);
 
    bool operator==(const HeavyIon& oth) const { return Z == oth.Z && A == oth.A; }
    bool operator!=(const HeavyIon& oth) const { return !(*this == oth); }
 
    double massN() const;
    double massP() const;
    double mass() const { return massN() + massP(); }
    double radius() const;
 
    static HeavyIon fromPdgId(pdgid_t);
    static bool isHI(const spdgid_t&);
    static bool isHI(const ParticleProperties&);
    static double mass(const HeavyIon& hi) { return hi.mass(); }
 
    static inline HeavyIon proton() { return HeavyIon(1, Element::H); }
    static inline HeavyIon neutron() { return HeavyIon(1, Element::neutron); }
    static inline HeavyIon Au() { return HeavyIon(197, Element::Au); }
    static inline HeavyIon Pb() { return HeavyIon(207, Element::Pb); }
 
    operator pdgid_t() const;
    friend std::ostream& operator<<(std::ostream& os, const HeavyIon& hi);
 
    Element Z{Element::invalid};
    unsigned short A{0};
  };
}  // namespace cepgen
 
#endif