Particle.cc

Go to the documentation of this file.

#include <cmath>

#include <boost/tuple/tuple_io.hpp>

#include <utl/Particle.h>
#include <utl/ParticlePropertiesFactory.h>
#include <utl/NucleusProperties.h>

#include <utl/CoordinateSystem.h>
#include <utl/ErrorLogger.h>
#include <utl/Math.h>

using namespace utl;


Particle::Particle(const int type, const Source& source,
                   const Point& position, const Vector& direction,
                   const TimeInterval& time,
                   const double weight, const double kineticEnergy) :
  fProperties(ParticlePropertiesFactory::Create(type)),
  fSource(source),
  fPosition(position),
  fDirection(direction),
  fTime(time),
  fWeight(weight),
  fKineticEnergy(kineticEnergy)
{ }


Particle::Particle(const int type, const Source& source,
                   const Point& position, const Vector& momentum,
                   const TimeInterval& time,
                   const double weight) :
  fProperties(ParticlePropertiesFactory::Create(type)),
  fSource(source),
  fPosition(position),
  //fDirection
  fTime(time),
  fWeight(weight)
{
  SetMomentum(momentum);
}


utl::Vector
Particle::GetMomentum()
  const
{
  const double m = GetMass();
  const double p = std::sqrt(Sqr(fKineticEnergy + m) - Sqr(m));
  return p * fDirection;
}


void
Particle::SetMomentum(const utl::Vector& momentum)
{
  fDirection = Normalized(momentum);
  const double m = GetMass();
  fKineticEnergy = std::sqrt(momentum.GetMag2() + Sqr(m)) - m;
}


int
Particle::NucleusCode(const int charge, const int atomicNumber)
{
  return NucleusProperties::TypeCode(charge, atomicNumber);
}