NucleusProperties.cc

Go to the documentation of this file.

#include <utl/NucleusProperties.h>
#include <utl/Particle.h>
#include <utl/PhysicalConstants.h>
#include <utl/Math.h>
#include <boost/format.hpp>

using boost::format;
using boost::io::str;


namespace {

  const char* gNuclearNames[] = {
    "neutron",
    "H",
    "He",
    "Li",
    "Be",
    "B",
    "C",
    "N",
    "O",
    "F",
    "Ne",
    "Na",
    "Mg",
    "Al",
    "Si",
    "P",
    "S",
    "Cl",
    "Ar",
    "K",
    "Ca",
    "Sc",
    "Ti",
    "V",
    "Cr",
    "Mn",
    "Fe"
  };

  unsigned int kMaxNuclearNameIndex = utl::Length(gNuclearNames) - 1;

}


namespace utl {

  std::string
  NucleusProperties::GetName()
    const
  {
    const unsigned int z = (fType / kChargeFactor) % kChargeMask;
    const unsigned int a = (fType / kAtomicNumberFactor) % kAtomicNumberMask;

    return (z <= kMaxNuclearNameIndex) ?
      str(format("%1% %2%") % gNuclearNames[z] % a) :
      str(format("N(%1%) %2%") % z % a);
  }


  double
  NucleusProperties::GetAtomicNumber()
    const
  {
    return (fType / kAtomicNumberFactor) % kAtomicNumberMask;
  }


  double
  NucleusProperties::GetMass()
    const
  {
    return GetAtomicNumber() * kProtonMass;
  }


  int
  NucleusProperties::TypeCode(const unsigned int charge,
                              const unsigned int atomicNumber)
  {
    return
      kNucleusBase
      + kChargeFactor       * charge
      + kAtomicNumberFactor * atomicNumber;
  }


  bool
  NucleusProperties::IsNucleus(const int type)
  {
    return kNucleusBase < type && type < kNucleusMax;
  }

}