genie::GHepParticle Class Reference

STDHEP-like event record entry that can fit a particle or a nucleus. More...

#include <GHepParticle.h>

Inheritance diagram for genie::GHepParticle:

Public Member Functions
	GHepParticle ()
	GHepParticle (const GHepParticle &particle)
	GHepParticle (int pdg, GHepStatus_t status, int mother1, int mother2, int daughter1, int daughter2, const TLorentzVector &p, const TLorentzVector &v)
	GHepParticle (int pdg, GHepStatus_t status, int mother1, int mother2, int daughter1, int daughter2, double px, double py, double pz, double E, double x, double y, double z, double t)
	GHepParticle (TRootIOCtor *)
	~GHepParticle ()
int	Pdg (void) const
GHepStatus_t	Status (void) const
int	RescatterCode (void) const
int	FirstMother (void) const
int	LastMother (void) const
int	FirstDaughter (void) const
int	LastDaughter (void) const
bool	HasDaughters (void) const
bool	IsBound (void) const
string	Name (void) const
	Name that corresponds to the PDG code. More...
double	Mass (void) const
	Mass that corresponds to the PDG code. More...
double	Charge (void) const
	Chrg that corresponds to the PDG code. More...
const TLorentzVector *	P4 (void) const
const TLorentzVector *	X4 (void) const
TLorentzVector *	P4 (void)
TLorentzVector *	X4 (void)
TLorentzVector *	GetP4 (void) const
TLorentzVector *	GetX4 (void) const
double	Px (void) const
	Get Px. More...
double	Py (void) const
	Get Py. More...
double	Pz (void) const
	Get Pz. More...
double	E (void) const
	Get energy. More...
double	Energy (void) const
	Get energy. More...
double	KinE (bool mass_from_pdg=false) const
	Get kinetic energy. More...
double	Vx (void) const
	Get production x. More...
double	Vy (void) const
	Get production y. More...
double	Vz (void) const
	Get production z. More...
double	Vt (void) const
	Get production time. More...
double	RemovalEnergy (void) const
	Get removal energy. More...
bool	Compare (const GHepParticle *p) const
bool	ComparePdgCodes (const GHepParticle *p) const
bool	CompareStatusCodes (const GHepParticle *p) const
bool	CompareFamily (const GHepParticle *p) const
bool	CompareMomentum (const GHepParticle *p) const
bool	IsOnMassShell (void) const
bool	IsOffMassShell (void) const
int	Z (void) const
int	A (void) const
double	PolzPolarAngle (void) const
double	PolzAzimuthAngle (void) const
bool	PolzIsSet (void) const
void	GetPolarization (TVector3 &polz)
void	SetPdgCode (int c)
void	SetStatus (GHepStatus_t s)
void	SetRescatterCode (int code)
void	SetFirstMother (int m)
void	SetLastMother (int m)
void	SetFirstDaughter (int d)
void	SetLastDaughter (int d)
void	SetMomentum (const TLorentzVector &p4)
void	SetPosition (const TLorentzVector &v4)
void	SetMomentum (double px, double py, double pz, double E)
void	SetPosition (double x, double y, double z, double t)
void	SetEnergy (double E)
void	SetPx (double px)
void	SetPy (double py)
void	SetPz (double pz)
void	SetPolarization (double theta, double phi)
void	SetPolarization (const TVector3 &polz)
void	SetBound (bool bound)
void	SetRemovalEnergy (double Erm)
void	CleanUp (void)
void	Reset (void)
void	Clear (Option_t *option)
void	Copy (const GHepParticle &particle)
void	Print (ostream &stream) const
void	Print (Option_t *opt) const
bool	operator== (const GHepParticle &p) const
GHepParticle &	operator= (const GHepParticle &p)

Private Member Functions
void	Init (void)
void	AssertIsKnownParticle (void) const

Private Attributes
int	fPdgCode
	particle PDG code More...
GHepStatus_t	fStatus
	particle status More...
int	fRescatterCode
	rescattering code More...
int	fFirstMother
	first mother idx More...
int	fLastMother
	last mother idx More...
int	fFirstDaughter
	first daughter idx More...
int	fLastDaughter
	last daughter idx More...
TLorentzVector *	fP4
	momentum 4-vector (GeV) More...
TLorentzVector *	fX4
	position 4-vector (in the target nucleus coordinate system / x,y,z in fm / t from the moment of the primary interaction in ys(yocto second = 10^-24 s) More...
double	fPolzTheta
	polar polarization angle (rad) More...
double	fPolzPhi
	azimuthal polarization angle (rad) More...
double	fRemovalEnergy
	removal energy for bound nucleons (GeV) More...
bool	fIsBound
	'is it a bound particle?' flag More...

Friends
ostream &	operator<< (ostream &stream, const GHepParticle &p)

Detailed Description

STDHEP-like event record entry that can fit a particle or a nucleus.

Author

Costas Andreopoulos <constantinos.andreopoulos cern.ch> University of Liverpool & STFC Rutherford Appleton Laboratory

November 18, 2004

Copyright (c) 2003-2020, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 39 of file GHepParticle.h.

Constructor & Destructor Documentation

GHepParticle::GHepParticle

(

)

Definition at line 50 of file GHepParticle.cxx.

                           :
TObject()
{
  this->Init();
}

GHepParticle::GHepParticle

(

const GHepParticle &

particle

)

Definition at line 104 of file GHepParticle.cxx.

                                                        :
TObject()
{
  this->Init();
  this->Copy(particle);
}

GHepParticle::GHepParticle	(	int	pdg,
		GHepStatus_t	status,
		int	mother1,
		int	mother2,
		int	daughter1,
		int	daughter2,
		const TLorentzVector &	p,
		const TLorentzVector &	v
	)

Definition at line 57 of file GHepParticle.cxx.

                                                            :
TObject(),
fStatus(status),
fFirstMother(mother1),
fLastMother(mother2),
fFirstDaughter(daughter1),
fLastDaughter(daughter2)
{
  this->SetPdgCode(pdg);

  fP4 = new TLorentzVector(p);
  fX4 = new TLorentzVector(v);

  fRescatterCode  = -1;
  fPolzTheta      = -999;
  fPolzPhi        = -999;
  fIsBound        = false;
  fRemovalEnergy  = 0.;
}

GHepParticle::GHepParticle	(	int	pdg,
		GHepStatus_t	status,
		int	mother1,
		int	mother2,
		int	daughter1,
		int	daughter2,
		double	px,
		double	py,
		double	pz,
		double	E,
		double	x,
		double	y,
		double	z,
		double	t
	)

Definition at line 80 of file GHepParticle.cxx.

                                                :
TObject(),
fStatus(status),
fFirstMother(mother1),
fLastMother(mother2),
fFirstDaughter(daughter1),
fLastDaughter(daughter2)
{
  this->SetPdgCode(pdg);

  fP4 = new TLorentzVector(px,py,pz,En);
  fX4 = new TLorentzVector(x,y,z,t);

  fRescatterCode  = -1;
  fPolzTheta      = -999;
  fPolzPhi        = -999;
  fIsBound        = false;
  fRemovalEnergy  = 0.;
}

GHepParticle::GHepParticle

(

TRootIOCtor *

)

Definition at line 111 of file GHepParticle.cxx.

                                       :
TObject(),
fPdgCode(0),
fStatus(kIStUndefined),
fRescatterCode(-1),
fFirstMother(-1),
fLastMother(-1),
fFirstDaughter(-1),
fLastDaughter(-1),
fP4(0),
fX4(0),
fPolzTheta(-999.),
fPolzPhi(-999.),
fRemovalEnergy(0),
fIsBound(false)
{

}

GHepParticle::~GHepParticle

(

)

Definition at line 130 of file GHepParticle.cxx.

{
  this->CleanUp();
}

Member Function Documentation

int GHepParticle::A

(

void

)

const

Definition at line 184 of file GHepParticle.cxx.

{
// Decoding A from the PDG code

  if(!pdg::IsIon(fPdgCode))
    return -1;
  else
    return pdg::IonPdgCodeToA(fPdgCode);
}

void GHepParticle::AssertIsKnownParticle ( void  ) const

private

Definition at line 535 of file GHepParticle.cxx.

{
  TParticlePDG * p = PDGLibrary::Instance()->Find(fPdgCode);
  if(!p) {
    LOG("GHepParticle", pFATAL)
      << "\n** You are attempting to insert particle with PDG code = "
      << fPdgCode << " into the event record."
      << "\n** This particle can not be found in "
      << "$GENIE/data/evgen/catalogues/pdg/genie_pdg_table.txt";
    gAbortingInErr = true;
    exit(1);
  }
}

double GHepParticle::Charge

(

void

)

const

Chrg that corresponds to the PDG code.

Definition at line 151 of file GHepParticle.cxx.

{
  this->AssertIsKnownParticle();

  TParticlePDG * p = PDGLibrary::Instance()->Find(fPdgCode);
  return p->Charge();
}

void GHepParticle::CleanUp

(

void

)

Definition at line 405 of file GHepParticle.cxx.

{
// deallocate memory

  if(fP4) delete fP4;
  if(fX4) delete fX4;
  fP4 = 0;
  fX4 = 0;
}

void GHepParticle::Clear

(

Option_t *

option

)

Definition at line 423 of file GHepParticle.cxx.

{
// implement the Clear(Option_t *) method so that the GHepParticle when is a
// member of a GHepRecord, gets deleted properly when calling TClonesArray's
// Clear("C")

  this->CleanUp();
}

bool GHepParticle::Compare

(

const GHepParticle *

p

)

const

Definition at line 463 of file GHepParticle.cxx.

{
// Do the comparisons in steps & put the ones that cost most
// in the inner-most {}

  bool same_particle = (this->fPdgCode == p->fPdgCode);
  bool same_status   = (this->fStatus  == p->fStatus );

  if( !same_particle || !same_status )  return false;
  else {
     if ( ! this->CompareFamily(p) )    return false;
     else {
       if( ! this->CompareMomentum(p) ) return false;
       else                             return true;
     }
  }
}

bool GHepParticle::CompareFamily

(

const GHepParticle *

p

)

const

Definition at line 491 of file GHepParticle.cxx.

{
  bool same_family  = (
          this->fFirstMother == p->fFirstMother &&
                  this->fLastMother  == p->fLastMother &&
                          this->fFirstDaughter == p->fFirstDaughter &&
                                   this->fLastDaughter  == p->fLastDaughter
  );
  return same_family;
}

bool GHepParticle::CompareMomentum

(

const GHepParticle *

p

)

const

Definition at line 502 of file GHepParticle.cxx.

{
  double dE  = TMath::Abs( this->E()  - p->E()  );
  double dPx = TMath::Abs( this->Px() - p->Px() );
  double dPy = TMath::Abs( this->Py() - p->Py() );
  double dPz = TMath::Abs( this->Pz() - p->Pz() );

  bool same_momentum =
       (dE < kPCutOff && dPx < kPCutOff && dPy < kPCutOff && dPz < kPCutOff);

  return same_momentum;
}

bool GHepParticle::ComparePdgCodes

(

const GHepParticle *

p

)

const

Definition at line 481 of file GHepParticle.cxx.

{
  return (this->fPdgCode == p->fPdgCode);
}

bool GHepParticle::CompareStatusCodes

(

const GHepParticle *

p

)

const

Definition at line 486 of file GHepParticle.cxx.

{
  return (this->fStatus == p->fStatus);
}

void GHepParticle::Copy

(

const GHepParticle &

particle

)

Definition at line 515 of file GHepParticle.cxx.

{
  this->SetStatus           (particle.Status()          );
  this->SetPdgCode          (particle.Pdg()             );
  this->SetRescatterCode    (particle.RescatterCode()   );
  this->SetFirstMother      (particle.FirstMother()     );
  this->SetLastMother       (particle.LastMother()      );
  this->SetFirstDaughter    (particle.FirstDaughter()   );
  this->SetLastDaughter     (particle.LastDaughter()    );

  this->SetMomentum (*particle.P4());
  this->SetPosition (*particle.X4());

  this->fPolzTheta = particle.fPolzTheta;
  this->fPolzPhi   = particle.fPolzPhi;

  this->fIsBound       = particle.fIsBound;
  this->fRemovalEnergy = particle.fRemovalEnergy;
}

double genie::GHepParticle::E ( void  ) const

inline

Get energy.

Definition at line 91 of file GHepParticle.h.

double genie::GHepParticle::Energy ( void  ) const

inline

Get energy.

Definition at line 92 of file GHepParticle.h.

int genie::GHepParticle::FirstDaughter ( void  ) const

inline

Definition at line 68 of file GHepParticle.h.

{ return  fFirstDaughter;      }

int genie::GHepParticle::FirstMother ( void  ) const

inline

Definition at line 66 of file GHepParticle.h.

{ return  fFirstMother;        }

TLorentzVector * GHepParticle::GetP4

(

void

)

const

Definition at line 194 of file GHepParticle.cxx.

{
// see GHepParticle::P4() for a method that does not create a new object and
// transfers its ownership

  if(fP4) {
     TLorentzVector * p4 = new TLorentzVector(*fP4);
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
     LOG("GHepParticle", pDEBUG)
          << "Return vp = " << utils::print::P4AsShortString(p4);
#endif
     return p4;
  } else {
    LOG("GHepParticle", pWARN) << "NULL 4-momentum TLorentzVector";
    return 0;
  }
}

void GHepParticle::GetPolarization

(

TVector3 &

polz

)

Definition at line 315 of file GHepParticle.cxx.

{
// gets the polarization vector

  if(! this->PolzIsSet() ) {
      polz.SetXYZ(0.,0.,0.);
      return;
  }
  polz.SetX( TMath::Sin(fPolzTheta) * TMath::Cos(fPolzPhi) );
  polz.SetY( TMath::Sin(fPolzTheta) * TMath::Sin(fPolzPhi) );
  polz.SetZ( TMath::Cos(fPolzTheta) );
}

TLorentzVector * GHepParticle::GetX4

(

void

)

const

Definition at line 212 of file GHepParticle.cxx.

{
// see GHepParticle::X4() for a method that does not create a new object and
// transfers its ownership

  if(fX4) {
     TLorentzVector * x4 = new TLorentzVector(*fX4);
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
     LOG("GHepParticle", pDEBUG)
         << "Return x4 = " << utils::print::X4AsString(x4);
#endif
     return x4;
  } else {
    LOG("GHepParticle", pWARN) << "NULL 4-position TLorentzVector";
    return 0;
  }
}

bool genie::GHepParticle::HasDaughters ( void  ) const

inline

Definition at line 70 of file GHepParticle.h.

{ return (fFirstDaughter!=-1); }

void GHepParticle::Init ( void  )

private

Definition at line 388 of file GHepParticle.cxx.

{
  fPdgCode       = 0;
  fStatus        = kIStUndefined;
  fRescatterCode = -1;
  fFirstMother   = -1;
  fLastMother    = -1;
  fFirstDaughter = -1;
  fLastDaughter  = -1;
  fPolzTheta     = -999;
  fPolzPhi       = -999;
  fIsBound       = false;
  fRemovalEnergy = 0.;
  fP4            = new TLorentzVector(0,0,0,0);
  fX4            = new TLorentzVector(0,0,0,0);
}

bool genie::GHepParticle::IsBound ( void  ) const

inline

Definition at line 71 of file GHepParticle.h.

{ return  fIsBound;            }

bool GHepParticle::IsOffMassShell

(

void

)

const

Definition at line 303 of file GHepParticle.cxx.

{
  return (! this->IsOnMassShell());
}

bool GHepParticle::IsOnMassShell

(

void

)

const

Definition at line 289 of file GHepParticle.cxx.

{
  this->AssertIsKnownParticle();

  TParticlePDG * p = PDGLibrary::Instance()->Find(fPdgCode);

  double Mpdg = p->Mass();
  double M4p  = (fP4) ? fP4->M() : 0.;

//  return utils::math::AreEqual(Mpdg, M4p);

  return (TMath::Abs(M4p-Mpdg) < kOffShellDm);
}

double GHepParticle::KinE

(

bool

mass_from_pdg = false

)

const

Get kinetic energy.

Definition at line 159 of file GHepParticle.cxx.

{
  if(!fP4) {
    LOG("GHepParticle", pWARN) << "4-momentum not yet set!";
    return 0;
  }

  double En = fP4->Energy();
  double M = ( (mass_from_pdg) ? this->Mass() : fP4->M() );
  double K = En - M;

  K = TMath::Max(K,0.);
  return K;
}

int genie::GHepParticle::LastDaughter ( void  ) const

inline

Definition at line 69 of file GHepParticle.h.

{ return  fLastDaughter;       }

int genie::GHepParticle::LastMother ( void  ) const

inline

Definition at line 67 of file GHepParticle.h.

{ return  fLastMother;         }

double GHepParticle::Mass

(

void

)

const

Mass that corresponds to the PDG code.

Definition at line 143 of file GHepParticle.cxx.

{
  this->AssertIsKnownParticle();

  TParticlePDG * p = PDGLibrary::Instance()->Find(fPdgCode);
  return p->Mass();
}

string GHepParticle::Name

(

void

)

const

Name that corresponds to the PDG code.

Definition at line 135 of file GHepParticle.cxx.

{
  this->AssertIsKnownParticle();

  TParticlePDG * p = PDGLibrary::Instance()->Find(fPdgCode);
  return p->GetName();
}

GHepParticle & GHepParticle::operator=

(

const GHepParticle &

p

)

Definition at line 554 of file GHepParticle.cxx.

{
  this->Copy(p);
  return (*this);
}

bool GHepParticle::operator==

(

const GHepParticle &

p

)

const

Definition at line 549 of file GHepParticle.cxx.

{
  return (this->Compare(&p));
}

const TLorentzVector* genie::GHepParticle::P4 ( void  ) const

inline

Definition at line 78 of file GHepParticle.h.

{ return fP4; }

TLorentzVector* genie::GHepParticle::P4 ( void  )

inline

Definition at line 80 of file GHepParticle.h.

{ return fP4; }

int genie::GHepParticle::Pdg ( void  ) const

inline

Definition at line 63 of file GHepParticle.h.

{ return  fPdgCode;            }

double genie::GHepParticle::PolzAzimuthAngle ( void  ) const

inline

Definition at line 120 of file GHepParticle.h.

{ return fPolzPhi; }

bool GHepParticle::PolzIsSet

(

void

)

const

Definition at line 308 of file GHepParticle.cxx.

{
// checks whether the polarization angles have been set

  return (fPolzTheta > -999 && fPolzPhi > -999);
}

double genie::GHepParticle::PolzPolarAngle ( void  ) const

inline

Definition at line 119 of file GHepParticle.h.

{ return fPolzTheta; }

void GHepParticle::Print

(

ostream &

stream

)

const

Definition at line 432 of file GHepParticle.cxx.

{
  stream << "\n |";
  stream << setfill(' ') << setw(14) << this->Name()           << " | ";
  stream << setfill(' ') << setw(14) << this->Pdg()            << " | ";
  stream << setfill(' ') << setw(6)  << this->Status()         << " | ";
  stream << setfill(' ') << setw(3)  << this->FirstMother()    << " | ";
  stream << setfill(' ') << setw(3)  << this->LastMother()     << " | ";
  stream << setfill(' ') << setw(3)  << this->FirstDaughter()  << " | ";
  stream << setfill(' ') << setw(3)  << this->LastDaughter()   << " | ";
  stream << std::fixed  << setprecision(3);
  stream << setfill(' ') << setw(6)  << this->Px()             << " | ";
  stream << setfill(' ') << setw(6)  << this->Py()             << " | ";
  stream << setfill(' ') << setw(6)  << this->Pz()             << " | ";
  stream << setfill(' ') << setw(6)  << this->E()              << " | ";
  stream << setfill(' ') << setw(6)  << this->Mass()           << " | ";

  int rescat_code = this->RescatterCode();
  if( rescat_code != -1 )
  {
     stream << setfill(' ') << setw(5)  << rescat_code << " | ";
  }
}

void GHepParticle::Print

(

Option_t *

opt

)

const

Definition at line 456 of file GHepParticle.cxx.

{
// implement the TObject's Print(Option_t *) method

  this->Print(cout);
}

double genie::GHepParticle::Px ( void  ) const

inline

Get Px.

Definition at line 88 of file GHepParticle.h.

double genie::GHepParticle::Py ( void  ) const

inline

Get Py.

Definition at line 89 of file GHepParticle.h.

double genie::GHepParticle::Pz ( void  ) const

inline

Get Pz.

Definition at line 90 of file GHepParticle.h.

double genie::GHepParticle::RemovalEnergy ( void  ) const

inline

Get removal energy.

Definition at line 100 of file GHepParticle.h.

int genie::GHepParticle::RescatterCode ( void  ) const

inline

Definition at line 65 of file GHepParticle.h.

{ return  fRescatterCode;      }

void GHepParticle::Reset

(

void

)

Definition at line 415 of file GHepParticle.cxx.

{
// deallocate memory + initialize

  this->CleanUp();
  this->Init();
}

void GHepParticle::SetBound

(

bool

bound

)

Definition at line 361 of file GHepParticle.cxx.

{
  // only set it for p or n
  bool is_nucleon = pdg::IsNeutronOrProton(fPdgCode);
  if(!is_nucleon && bound) {
    LOG("GHepParticle", pERROR)
       << "Refusing to set the bound flag for particles other than nucleons";
    LOG("GHepParticle", pERROR)
       << "(Requested for pdg = " << fPdgCode << ")";
    return;
  }
  // if the particles isn't bound then make sure that its removal energy = 0
  if(!bound) {
   fRemovalEnergy = 0;
  }
  // set the flag
  fIsBound = bound;
}

void GHepParticle::SetEnergy

(

double

E

)

Definition at line 269 of file GHepParticle.cxx.

{
  this->SetMomentum(this->Px(), this->Py(), this->Pz(), En);
}

void genie::GHepParticle::SetFirstDaughter ( int  d )

inline

Definition at line 134 of file GHepParticle.h.

{ fFirstDaughter = d; }

void genie::GHepParticle::SetFirstMother ( int  m )

inline

Definition at line 132 of file GHepParticle.h.

{ fFirstMother   = m; }

void genie::GHepParticle::SetLastDaughter ( int  d )

inline

Definition at line 135 of file GHepParticle.h.

{ fLastDaughter  = d; }

void genie::GHepParticle::SetLastMother ( int  m )

inline

Definition at line 133 of file GHepParticle.h.

{ fLastMother    = m; }

void GHepParticle::SetMomentum

(

const TLorentzVector &

p4

)

Definition at line 236 of file GHepParticle.cxx.

{
  if(fP4)
      fP4->SetPxPyPzE( p4.Px(), p4.Py(), p4.Pz(), p4.Energy() );
  else
      fP4 = new TLorentzVector(p4);
}

void GHepParticle::SetMomentum	(	double	px,
		double	py,
		double	pz,
		double	E
	)

Definition at line 244 of file GHepParticle.cxx.

{
  if(fP4)
      fP4->SetPxPyPzE(px, py, pz, En);
  else
      fP4 = new TLorentzVector(px, py, pz, En);
}

void GHepParticle::SetPdgCode

(

int

c

)

Definition at line 230 of file GHepParticle.cxx.

{
  fPdgCode = code;
  this->AssertIsKnownParticle();
}

void GHepParticle::SetPolarization	(	double	theta,
		double	phi
	)

Definition at line 328 of file GHepParticle.cxx.

{
// sets the polarization angles

  if(theta>=0 && theta<=kPi && phi>=0 && phi<2*kPi)
  {
    fPolzTheta = theta;
    fPolzPhi   = phi;

  } else {
    LOG("GHepParticle", pERROR)
      << "Invalid polarization angles (polar = " << theta
      << ", azimuthal = " << phi << ")";
  }
}

void GHepParticle::SetPolarization

(

const TVector3 &

polz

)

Definition at line 344 of file GHepParticle.cxx.

{
// sets the polarization angles

  double p = polz.Mag();
  if(! (p>0) ) {
    LOG("GHepParticle", pERROR)
           << "Input polarization vector has non-positive norm! Ignoring it";
    return;
  }

  double theta = TMath::ACos(polz.z()/p);
  double phi   = kPi + TMath::ATan2(-polz.y(), -polz.x());

  this->SetPolarization(theta,phi);
}

void GHepParticle::SetPosition

(

const TLorentzVector &

v4

)

Definition at line 252 of file GHepParticle.cxx.

{
  this->SetPosition(v4.X(), v4.Y(), v4.Z(), v4.T());
}

void GHepParticle::SetPosition	(	double	x,
		double	y,
		double	z,
		double	t
	)

Definition at line 257 of file GHepParticle.cxx.

{
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GHepParticle", pDEBUG)
            << "Setting position to (x = " << x << ", y = "
                               << y << ", z = " << z << ", t = " << t << ")";
#endif

  if(fX4) fX4->SetXYZT(x,y,z,t);
  else    fX4 = new TLorentzVector(x,y,z,t);
}

void GHepParticle::SetPx

(

double

px

)

Definition at line 274 of file GHepParticle.cxx.

{
  this->SetMomentum(px, this->Py(), this->Pz(), this->E());
}

void GHepParticle::SetPy

(

double

py

)

Definition at line 279 of file GHepParticle.cxx.

{
  this->SetMomentum(this->Px(), py, this->Pz(), this->E());
}

void GHepParticle::SetPz

(

double

pz

)

Definition at line 284 of file GHepParticle.cxx.

{
  this->SetMomentum(this->Px(), this->Py(), pz, this->E());
}

void GHepParticle::SetRemovalEnergy

(

double

Erm

)

Definition at line 380 of file GHepParticle.cxx.

{
  fRemovalEnergy = TMath::Max(Erm, 0.); // non-negative

  // if a value was set, make sure that the IsBound flag is turned on
  if(fRemovalEnergy>0) this->SetBound(true);
}

void genie::GHepParticle::SetRescatterCode ( int  code )

inline

Definition at line 129 of file GHepParticle.h.

{ fRescatterCode = code; }

void genie::GHepParticle::SetStatus ( GHepStatus_t  s )

inline

Definition at line 126 of file GHepParticle.h.

{ fStatus = s; }

GHepStatus_t genie::GHepParticle::Status ( void  ) const

inline

Definition at line 64 of file GHepParticle.h.

{ return  fStatus;             }

double genie::GHepParticle::Vt ( void  ) const

inline

Get production time.

Definition at line 97 of file GHepParticle.h.

double genie::GHepParticle::Vx ( void  ) const

inline

Get production x.

Definition at line 94 of file GHepParticle.h.

double genie::GHepParticle::Vy ( void  ) const

inline

Get production y.

Definition at line 95 of file GHepParticle.h.

double genie::GHepParticle::Vz ( void  ) const

inline

Get production z.

Definition at line 96 of file GHepParticle.h.

const TLorentzVector* genie::GHepParticle::X4 ( void  ) const

inline

Definition at line 79 of file GHepParticle.h.

{ return fX4; }

TLorentzVector* genie::GHepParticle::X4 ( void  )

inline

Definition at line 81 of file GHepParticle.h.

{ return fX4; }

int GHepParticle::Z

(

void

)

const

Definition at line 174 of file GHepParticle.cxx.

{
// Decoding Z from the PDG code

  if(!pdg::IsIon(fPdgCode))
    return -1;
  else
    return pdg::IonPdgCodeToZ(fPdgCode);
}

Friends And Related Function Documentation

ostream& operator<< ( ostream &  stream, const GHepParticle &  p  )

friend

Member Data Documentation

int genie::GHepParticle::fFirstDaughter

private

first daughter idx

Definition at line 179 of file GHepParticle.h.

int genie::GHepParticle::fFirstMother

private

first mother idx

Definition at line 177 of file GHepParticle.h.

bool genie::GHepParticle::fIsBound

private

'is it a bound particle?' flag

Definition at line 186 of file GHepParticle.h.

int genie::GHepParticle::fLastDaughter

private

last daughter idx

Definition at line 180 of file GHepParticle.h.

int genie::GHepParticle::fLastMother

private

last mother idx

Definition at line 178 of file GHepParticle.h.

TLorentzVector* genie::GHepParticle::fP4

private

momentum 4-vector (GeV)

Definition at line 181 of file GHepParticle.h.

int genie::GHepParticle::fPdgCode

private

particle PDG code

Definition at line 174 of file GHepParticle.h.

double genie::GHepParticle::fPolzPhi

private

azimuthal polarization angle (rad)

Definition at line 184 of file GHepParticle.h.

double genie::GHepParticle::fPolzTheta

private

polar polarization angle (rad)

Definition at line 183 of file GHepParticle.h.

double genie::GHepParticle::fRemovalEnergy

private

removal energy for bound nucleons (GeV)

Definition at line 185 of file GHepParticle.h.

int genie::GHepParticle::fRescatterCode

private

rescattering code

Definition at line 176 of file GHepParticle.h.

GHepStatus_t genie::GHepParticle::fStatus

private

particle status

Definition at line 175 of file GHepParticle.h.

TLorentzVector* genie::GHepParticle::fX4

private

position 4-vector (in the target nucleus coordinate system / x,y,z in fm / t from the moment of the primary interaction in ys(yocto second = 10^-24 s)

Definition at line 182 of file GHepParticle.h.

---

The documentation for this class was generated from the following files:

• Generator/src/Framework/GHEP/GHepParticle.h
• Generator/src/Framework/GHEP/GHepParticle.cxx