GeoSwimParticle.cxx

Go to the documentation of this file.

//
// GeoSwimParticle
//
// Adapted from example particle class provided with geant4_vmc package. 
// January 20, 2006  M.Ishitsuka First version for GeoSwimmer
//                   ref. SwimPartile.cxx
#include <iostream>
using namespace std;

#include <TParticle.h>
#include <TObjArray.h>
#include <TMCProcess.h>

#include "GeoSwimmer/GeoSwimParticle.h"

ClassImp(GeoSwimParticle)

//__________________________________________________________
GeoSwimParticle::GeoSwimParticle(Int_t pdg, TParticle* particle) :
  fPDG(pdg),fS(0),fRange(0),fParticle(particle)
{ 
}

//__________________________________________________________
GeoSwimParticle::GeoSwimParticle(const TVector3 position,
                                       const TVector3 momentum,
                                       Double_t       charge,
                                       Int_t          pdg,
                                       Double_t       mass) :
  fInitPosition(position),fPosition(position),fMomentum(momentum),
  fMass(mass),fCharge(charge),fPDG(pdg),fS(0),fRange(0),fParticle(0)
{ 
}

//_________________________________________________________
GeoSwimParticle::GeoSwimParticle() : 
  fPDG(-1),fS(0),fRange(0),fParticle(0)
{
}

//__________________________________________________________
GeoSwimParticle::~GeoSwimParticle() {
// Destructor
}

//_________________________________________________________

const TVector3 GeoSwimParticle::GetInitPosition() const
{
  return fInitPosition;
} 

//_________________________________________________________

const TVector3 GeoSwimParticle::GetPosition() const
{
  return fPosition;
} 

//_________________________________________________________

const TVector3 GeoSwimParticle::GetMomentum() const
{
  return fMomentum;
} 

//_________________________________________________________

TVector3 GeoSwimParticle::GetDirection() const
{
  double modulus = GeoSwimParticle::GetMomentumModulus();  
  if (modulus!=0.0) return fMomentum.Unit();
  else return TVector3(0.0,0.0,0.0);
}

//_________________________________________________________

double GeoSwimParticle::GetMomentumModulus() const
{
  return fMomentum.Mag();
} 

//_________________________________________________________

double GeoSwimParticle::GetEnergy() const
{
  return sqrt(fMass*fMass + fMomentum.Mag2());
}

//_________________________________________________________

double GeoSwimParticle::GetMass() const 
{ 
  return fMass; 
} 

//_________________________________________________________

double GeoSwimParticle::GetCharge() const 
{ 
  return fCharge; 
} 

//__________________________________________________________
Int_t GeoSwimParticle:: GetID() const {
  return fPDG;
}  

//__________________________________________________________

double GeoSwimParticle::GetS() const 
{ 
  return fS; 
} 

//__________________________________________________________

double GeoSwimParticle::GetRange() const 
{ 
  return fRange; 
} 

//__________________________________________________________
TParticle*  GeoSwimParticle::GetParticle() const {
  return fParticle;
}  

//......................................................................

void GeoSwimParticle::SetPosition(const TVector3 position) 
{ 
  fPosition = position;
}

//......................................................................

void GeoSwimParticle::SetMomentum(const TVector3 momentum)
{
  fMomentum = momentum;
}

//......................................................................

void GeoSwimParticle::SetMass(double mass)
{
  fMass = mass;
}

//......................................................................

void GeoSwimParticle::SetCharge(double charge)
{
  fCharge = charge;
}

//......................................................................

void GeoSwimParticle::SetID(int pdg)
{
  fPDG = pdg;
}

//......................................................................

void GeoSwimParticle::SetS(double s)
{
  fS = s;
}

//......................................................................

void GeoSwimParticle::SetRange(double range)
{
  fRange = range;
}

---