G4ReactionProduct Class Reference

#include <G4ReactionProduct.hh>

Public Member Functions
	G4ReactionProduct ()
	G4ReactionProduct (const G4ParticleDefinition *aParticleDefinition)
	~G4ReactionProduct ()
	G4ReactionProduct (const G4ReactionProduct &right)
void *	operator new (size_t)
void	operator delete (void *aReactionProduct)
G4ReactionProduct &	operator= (const G4ReactionProduct &right)
G4ReactionProduct &	operator= (const G4DynamicParticle &right)
G4ReactionProduct &	operator= (const G4HadProjectile &right)
G4bool	operator== (const G4ReactionProduct &right) const
G4bool	operator!= (const G4ReactionProduct &right) const
const G4ParticleDefinition *	GetDefinition () const
void	SetDefinition (const G4ParticleDefinition *aParticleDefinition)
void	SetDefinitionAndUpdateE (const G4ParticleDefinition *aParticleDefinition)
void	SetMomentum (const G4double x, const G4double y, const G4double z)
void	SetMomentum (const G4double x, const G4double y)
void	SetMomentum (const G4double z)
void	SetMomentum (const G4ThreeVector &mom)
G4ThreeVector	GetMomentum () const
G4double	GetTotalMomentum () const
G4double	GetTotalEnergy () const
void	SetKineticEnergy (const G4double en)
G4double	GetKineticEnergy () const
void	SetTotalEnergy (const G4double en)
void	SetMass (const G4double mas)
G4double	GetMass () const
void	SetTOF (const G4double t)
G4double	GetTOF () const
void	SetSide (const G4int sid)
G4int	GetSide () const
void	SetCreatorModel (const G4int mod)
G4int	GetCreatorModel () const
void	SetNewlyAdded (const G4bool f)
G4bool	GetNewlyAdded () const
void	SetMayBeKilled (const G4bool f)
G4bool	GetMayBeKilled () const
void	SetZero ()
void	Lorentz (const G4ReactionProduct &p1, const G4ReactionProduct &p2)
G4double	Angle (const G4ReactionProduct &p) const
void	SetPositionInNucleus (G4double x, G4double y, G4double z)
void	SetPositionInNucleus (G4ThreeVector &aPosition)
G4ThreeVector	GetPositionInNucleus () const
G4double	GetXPositionInNucleus () const
G4double	GetYPositionInNucleus () const
G4double	GetZPositionInNucleus () const
void	SetFormationTime (G4double aTime)
G4double	GetFormationTime () const
void	HasInitialStateParton (G4bool aFlag)
G4bool	HasInitialStateParton () const

Friends
G4ReactionProduct	operator+ (const G4ReactionProduct &p1, const G4ReactionProduct &p2)
G4ReactionProduct	operator- (const G4ReactionProduct &p1, const G4ReactionProduct &p2)
G4ReactionProduct	operator* (const G4double aDouble, const G4ReactionProduct &p2)

Detailed Description

Definition at line 53 of file G4ReactionProduct.hh.

Constructor & Destructor Documentation

G4ReactionProduct::G4ReactionProduct

(

)

Definition at line 36 of file G4ReactionProduct.cc.

                                      :
    theParticleDefinition(NULL),
    formationTime(0.0),
    hasInitialStateParton(false),
    mass(0.0),
    totalEnergy(0.0),
    kineticEnergy(0.0),
    timeOfFlight(0.0),
    side(0),
    theCreatorModel(-1),
    NewlyAdded(false),
    MayBeKilled(true)
  {
    SetMomentum( 0.0, 0.0, 0.0 );
    SetPositionInNucleus( 0.0, 0.0, 0.0 );
  }

Here is the call graph for this function:

G4ReactionProduct::G4ReactionProduct

(

const G4ParticleDefinition *

aParticleDefinition

)

Definition at line 53 of file G4ReactionProduct.cc.

  {
    SetMomentum( 0.0, 0.0, 0.0 );
    SetPositionInNucleus( 0.0, 0.0, 0.0 );
    formationTime = 0.0;
    hasInitialStateParton = false;
    theParticleDefinition = aParticleDefinition;
    mass = aParticleDefinition->GetPDGMass();
    totalEnergy = mass;
    kineticEnergy = 0.0;
    (aParticleDefinition->GetPDGEncoding()<0) ? timeOfFlight=-1.0 : timeOfFlight=1.0;
    side = 0;
    theCreatorModel = -1;
    NewlyAdded = false;
    MayBeKilled = true;
  }

Here is the call graph for this function:

G4ReactionProduct::~G4ReactionProduct ( )

inline

Definition at line 77 of file G4ReactionProduct.hh.

{}

G4ReactionProduct::G4ReactionProduct

(

const G4ReactionProduct &

right

)

Definition at line 71 of file G4ReactionProduct.cc.

  {
    theParticleDefinition = right.theParticleDefinition;
    positionInNucleus = right.positionInNucleus;
    formationTime = right.formationTime;
    hasInitialStateParton = right.hasInitialStateParton;
    momentum = right.momentum;
    mass = right.mass;
    totalEnergy = right.totalEnergy;
    kineticEnergy = right.kineticEnergy;
    timeOfFlight = right.timeOfFlight;
    side = right.side;
    theCreatorModel = right.theCreatorModel;
    NewlyAdded = right.NewlyAdded;
    MayBeKilled = right.MayBeKilled;
  }

Member Function Documentation

G4double G4ReactionProduct::Angle

(

const G4ReactionProduct &

p

)

const

Definition at line 225 of file G4ReactionProduct.cc.

  {
    G4ThreeVector tM = momentum;
    G4ThreeVector pM = p.momentum;
    G4double tx = tM.x(); G4double ty = tM.y(); G4double tz = tM.z();
    G4double px = pM.x(); G4double py = pM.y(); G4double pz = pM.z();
    G4double a = std::sqrt( ( px*px + py*py + pz*pz ) * ( tx*tx + ty*ty + tz*tz ) );
    if( a == 0.0 ) {
      return 0.0;
    } else {
      a = ( tx*px + ty*py + tz*pz ) / a;
      if( std::fabs(a) > 1.0 ) { a<0.0 ? a=-1.0 : a=1.0; }
      return std::acos( a );
    }
  }

Here is the call graph for this function:

G4int G4ReactionProduct::GetCreatorModel ( ) const

inline

Definition at line 168 of file G4ReactionProduct.hh.

    { return theCreatorModel; }

const G4ParticleDefinition* G4ReactionProduct::GetDefinition ( ) const

inline

Definition at line 107 of file G4ReactionProduct.hh.

    { return theParticleDefinition; }

Here is the caller graph for this function:

G4double G4ReactionProduct::GetFormationTime ( ) const

inline

Definition at line 208 of file G4ReactionProduct.hh.

{ return formationTime; }

G4double G4ReactionProduct::GetKineticEnergy ( ) const

inline

Definition at line 138 of file G4ReactionProduct.hh.

    { return kineticEnergy; }

Here is the caller graph for this function:

G4double G4ReactionProduct::GetMass ( ) const

inline

Definition at line 150 of file G4ReactionProduct.hh.

    { return mass; }

Here is the caller graph for this function:

G4bool G4ReactionProduct::GetMayBeKilled ( ) const

inline

Definition at line 180 of file G4ReactionProduct.hh.

    { return MayBeKilled; }

G4ThreeVector G4ReactionProduct::GetMomentum ( void  ) const

inline

Definition at line 123 of file G4ReactionProduct.hh.

    { return momentum; }

G4bool G4ReactionProduct::GetNewlyAdded ( ) const

inline

Definition at line 174 of file G4ReactionProduct.hh.

    { return NewlyAdded; }

G4ThreeVector G4ReactionProduct::GetPositionInNucleus ( ) const

inline

Definition at line 201 of file G4ReactionProduct.hh.

{ return positionInNucleus; }

G4int G4ReactionProduct::GetSide ( ) const

inline

Definition at line 162 of file G4ReactionProduct.hh.

    { return side; }

Here is the caller graph for this function:

G4double G4ReactionProduct::GetTOF ( ) const

inline

Definition at line 156 of file G4ReactionProduct.hh.

    { return timeOfFlight; }

Here is the caller graph for this function:

G4double G4ReactionProduct::GetTotalEnergy ( void  ) const

inline

Definition at line 129 of file G4ReactionProduct.hh.

    { return totalEnergy; }

Here is the caller graph for this function:

G4double G4ReactionProduct::GetTotalMomentum ( ) const

inline

Definition at line 126 of file G4ReactionProduct.hh.

    { return std::sqrt(std::abs(kineticEnergy*(totalEnergy+mass))); }

Here is the caller graph for this function:

G4double G4ReactionProduct::GetXPositionInNucleus ( ) const

inline

Definition at line 202 of file G4ReactionProduct.hh.

{ return positionInNucleus.x(); }

Here is the call graph for this function:

G4double G4ReactionProduct::GetYPositionInNucleus ( ) const

inline

Definition at line 203 of file G4ReactionProduct.hh.

{ return positionInNucleus.y(); }

Here is the call graph for this function:

G4double G4ReactionProduct::GetZPositionInNucleus ( ) const

inline

Definition at line 204 of file G4ReactionProduct.hh.

{ return positionInNucleus.z(); }

Here is the call graph for this function:

void G4ReactionProduct::HasInitialStateParton ( G4bool  aFlag )

inline

Definition at line 210 of file G4ReactionProduct.hh.

{ hasInitialStateParton = aFlag; }

Here is the caller graph for this function:

G4bool G4ReactionProduct::HasInitialStateParton ( ) const

inline

Definition at line 212 of file G4ReactionProduct.hh.

{ return hasInitialStateParton; }

void G4ReactionProduct::Lorentz	(	const G4ReactionProduct &	p1,
		const G4ReactionProduct &	p2
	)

Definition at line 206 of file G4ReactionProduct.cc.

  {
    G4ThreeVector p1M = p1.momentum;
    G4ThreeVector p2M = p2.momentum;
    G4double p1x = p1M.x(); G4double p1y = p1M.y(); G4double p1z = p1M.z();
    G4double p2x = p2M.x(); G4double p2y = p2M.y(); G4double p2z = p2M.z();
    G4double a = ( (p1x*p2x+p1y*p2y+p1z*p2z)/(p2.totalEnergy+p2.mass) -
                   p1.totalEnergy ) / p2.mass;
    G4double x = p1x+a*p2x;
    G4double y = p1y+a*p2y;
    G4double z = p1z+a*p2z;
    G4double p = std::sqrt(x*x+y*y+z*z);
    SetMass( p1.mass );
    SetTotalEnergy( std::sqrt( (p1.mass+p)*(p1.mass+p) - 2.*p1.mass*p ) );
    //SetTotalEnergy( std::sqrt( p1.mass*p1.mass + x*x + y*y + z*z ) );
    SetMomentum( x, y, z );
  }

Here is the call graph for this function:

Here is the caller graph for this function:

void G4ReactionProduct::operator delete ( void *  aReactionProduct )

inline

Definition at line 91 of file G4ReactionProduct.hh.

                                                        {
      aRPAllocator->FreeSingle((G4ReactionProduct*)aReactionProduct);
    }

void* G4ReactionProduct::operator new ( size_t  )

inline

Definition at line 82 of file G4ReactionProduct.hh.

                                      {
      if (!aRPAllocator) aRPAllocator = new G4Allocator<G4ReactionProduct>  ;
      return (void *)aRPAllocator->MallocSingle();
    }

G4bool G4ReactionProduct::operator!= ( const G4ReactionProduct &  right ) const

inline

Definition at line 104 of file G4ReactionProduct.hh.

    { return ( this != (G4ReactionProduct*) &right ); }

G4ReactionProduct & G4ReactionProduct::operator=

(

const G4ReactionProduct &

right

)

Definition at line 89 of file G4ReactionProduct.cc.

  {
    if( this != &right ) {
      theParticleDefinition = right.theParticleDefinition;
      positionInNucleus = right.positionInNucleus;
      formationTime = right.formationTime;
      hasInitialStateParton = right.hasInitialStateParton;
      momentum = right.momentum;
      mass = right.mass;
      totalEnergy = right.totalEnergy;
      kineticEnergy = right.kineticEnergy;
      timeOfFlight = right.timeOfFlight;
      side = right.side;
      theCreatorModel = right.theCreatorModel;
      NewlyAdded = right.NewlyAdded;
      MayBeKilled = right.MayBeKilled;
    }
    return *this;
  }

G4ReactionProduct & G4ReactionProduct::operator=

(

const G4DynamicParticle &

right

)

Definition at line 110 of file G4ReactionProduct.cc.

  {
    theParticleDefinition = right.GetDefinition();
    SetPositionInNucleus( 0.0, 0.0, 0.0 );
    formationTime = 0.0;
    hasInitialStateParton = false;
    momentum = right.GetMomentum();
    mass = right.GetDefinition()->GetPDGMass();
    totalEnergy = right.GetTotalEnergy();
    kineticEnergy = right.GetKineticEnergy();
    (right.GetDefinition()->GetPDGEncoding()<0) ? timeOfFlight=-1.0 : timeOfFlight=1.0;
    side = 0;
    theCreatorModel = -1;
    NewlyAdded = false;
    MayBeKilled = true;
    return *this;
  }

Here is the call graph for this function:

G4ReactionProduct & G4ReactionProduct::operator=

(

const G4HadProjectile &

right

)

Definition at line 129 of file G4ReactionProduct.cc.

  {
    theParticleDefinition = right.GetDefinition();
    SetPositionInNucleus( 0.0, 0.0, 0.0 );
    formationTime = 0.0;
    hasInitialStateParton = false;
    momentum = right.Get4Momentum().vect();
    mass = right.GetDefinition()->GetPDGMass();
    totalEnergy = right.Get4Momentum().e();
    kineticEnergy = right.GetKineticEnergy();
    (right.GetDefinition()->GetPDGEncoding()<0) ? timeOfFlight=-1.0 : timeOfFlight=1.0;
    side = 0;
    theCreatorModel = -1;
    NewlyAdded = false;
    MayBeKilled = true;
    return *this;
  }

Here is the call graph for this function:

G4bool G4ReactionProduct::operator== ( const G4ReactionProduct &  right ) const

inline

Definition at line 101 of file G4ReactionProduct.hh.

    { return ( this == (G4ReactionProduct*) &right ); }

void G4ReactionProduct::SetCreatorModel ( const G4int  mod )

inline

Definition at line 165 of file G4ReactionProduct.hh.

    { theCreatorModel = mod; }

void G4ReactionProduct::SetDefinition

(

const G4ParticleDefinition *

aParticleDefinition

)

Definition at line 160 of file G4ReactionProduct.cc.

  {
    theParticleDefinition = aParticleDefinition;
    mass = aParticleDefinition->GetPDGMass();
    totalEnergy = mass;
    kineticEnergy = 0.0;
    (aParticleDefinition->GetPDGEncoding()<0) ?
      timeOfFlight=-1.0 : timeOfFlight=1.0;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

void G4ReactionProduct::SetDefinitionAndUpdateE

(

const G4ParticleDefinition *

aParticleDefinition

)

Definition at line 148 of file G4ReactionProduct.cc.

  {    G4double aKineticEnergy = GetKineticEnergy();
    G4double pp = GetMomentum().mag();
    G4ThreeVector aMomentum = GetMomentum();
    SetDefinition( aParticleDefinition );
    SetKineticEnergy( aKineticEnergy );
    if( pp > DBL_MIN )
      SetMomentum( aMomentum * (std::sqrt(aKineticEnergy*aKineticEnergy +
                                    2*aKineticEnergy*GetMass())/pp) );
  }

Here is the call graph for this function:

Here is the caller graph for this function:

void G4ReactionProduct::SetFormationTime ( G4double  aTime )

inline

Definition at line 206 of file G4ReactionProduct.hh.

{ formationTime = aTime; }

Here is the caller graph for this function:

void G4ReactionProduct::SetKineticEnergy ( const G4double  en )

inline

Definition at line 132 of file G4ReactionProduct.hh.

    {
      kineticEnergy = en;
      totalEnergy = kineticEnergy + mass;
    }

Here is the caller graph for this function:

void G4ReactionProduct::SetMass ( const G4double  mas )

inline

Definition at line 147 of file G4ReactionProduct.hh.

    { mass = mas; }

Here is the caller graph for this function:

void G4ReactionProduct::SetMayBeKilled ( const G4bool  f )

inline

Definition at line 177 of file G4ReactionProduct.hh.

    { MayBeKilled = f; }

Here is the caller graph for this function:

void G4ReactionProduct::SetMomentum	(	const G4double	x,
		const G4double	y,
		const G4double	z
	)

Definition at line 171 of file G4ReactionProduct.cc.

  {
    momentum.setX( x );
    momentum.setY( y );
    momentum.setZ( z );
  }

Here is the call graph for this function:

void G4ReactionProduct::SetMomentum	(	const G4double	x,
		const G4double	y
	)

Definition at line 179 of file G4ReactionProduct.cc.

  {
    momentum.setX( x );
    momentum.setY( y );
  }

Here is the call graph for this function:

void G4ReactionProduct::SetMomentum

(

const G4double

z

)

Definition at line 186 of file G4ReactionProduct.cc.

  {
    momentum.setZ( z );
  }

Here is the call graph for this function:

void G4ReactionProduct::SetMomentum ( const G4ThreeVector &  mom )

inline

Definition at line 120 of file G4ReactionProduct.hh.

    { momentum = mom; }

void G4ReactionProduct::SetNewlyAdded ( const G4bool  f )

inline

Definition at line 171 of file G4ReactionProduct.hh.

    { NewlyAdded = f; }

Here is the caller graph for this function:

void G4ReactionProduct::SetPositionInNucleus ( G4double  x, G4double  y, G4double  z  )

inline

Definition at line 189 of file G4ReactionProduct.hh.

     {
       positionInNucleus.setX(x);
       positionInNucleus.setY(y);
       positionInNucleus.setZ(z);
     }

Here is the call graph for this function:

Here is the caller graph for this function:

void G4ReactionProduct::SetPositionInNucleus ( G4ThreeVector &  aPosition )

inline

Definition at line 196 of file G4ReactionProduct.hh.

     {
       positionInNucleus = aPosition;
     }

void G4ReactionProduct::SetSide ( const G4int  sid )

inline

Definition at line 159 of file G4ReactionProduct.hh.

    { side = sid; }

Here is the caller graph for this function:

void G4ReactionProduct::SetTOF ( const G4double  t )

inline

Definition at line 153 of file G4ReactionProduct.hh.

    { timeOfFlight = t; }

Here is the caller graph for this function:

void G4ReactionProduct::SetTotalEnergy ( const G4double  en )

inline

Definition at line 141 of file G4ReactionProduct.hh.

    {
      totalEnergy = en;
      kineticEnergy = totalEnergy - mass;
    }

Here is the caller graph for this function:

void G4ReactionProduct::SetZero

(

)

Definition at line 191 of file G4ReactionProduct.cc.

  {
    SetMomentum( 0.0, 0.0, 0.0 );
    totalEnergy = 0.0;
    kineticEnergy = 0.0;
    mass = 0.0;
    timeOfFlight = 0.0;
    side = 0;
    theCreatorModel = -1;
    NewlyAdded = false;
    SetPositionInNucleus( 0.0, 0.0, 0.0 );
    formationTime = 0.0;
    hasInitialStateParton = false;
  }

Here is the call graph for this function:

Friends And Related Function Documentation

G4ReactionProduct operator* ( const G4double  aDouble, const G4ReactionProduct &  p2  )

friend

Definition at line 61 of file G4ReactionProduct.hh.

     {
       G4ReactionProduct result;
       result.SetMomentum(aDouble*p2.GetMomentum());
       result.SetMass(p2.GetMass());
       result.SetTotalEnergy(std::sqrt(result.GetMass()*result.GetMass()+
                                  result.GetMomentum()*result.GetMomentum()));
       return result;
     }

G4ReactionProduct operator+ ( const G4ReactionProduct &  p1, const G4ReactionProduct &  p2  )

friend

Definition at line 242 of file G4ReactionProduct.cc.

  {
    G4double totEnergy = p1.totalEnergy + p2.totalEnergy;
    G4double x = p1.momentum.x() + p2.momentum.x();
    G4double y = p1.momentum.y() + p2.momentum.y();
    G4double z = p1.momentum.z() + p2.momentum.z();
    G4double newMass = totEnergy*totEnergy - ( x*x + y*y + z*z );
    if( newMass < 0.0 )
      newMass = -1. * std::sqrt( -newMass );
    else
      newMass = std::sqrt( newMass );
    G4ReactionProduct result;
    result.SetMass( newMass );
    result.SetMomentum( x, y, z );
    result.SetTotalEnergy( totEnergy );
    result.SetPositionInNucleus( 0.0, 0.0, 0.0 );
    result.SetFormationTime(0.0);
    result.HasInitialStateParton(false);
    return result;
  }

G4ReactionProduct operator- ( const G4ReactionProduct &  p1, const G4ReactionProduct &  p2  )

friend

Definition at line 264 of file G4ReactionProduct.cc.

  {
    G4double totEnergy = p1.totalEnergy - p2.totalEnergy;
    G4double x = p1.momentum.x() - p2.momentum.x();
    G4double y = p1.momentum.y() - p2.momentum.y();
    G4double z = p1.momentum.z() - p2.momentum.z();
    G4double newMass = totEnergy*totEnergy - ( x*x + y*y + z*z );
    if( newMass < 0.0 )
      newMass = -1. * std::sqrt( -newMass );
    else
      newMass = std::sqrt( newMass );
    G4ReactionProduct result;
    result.SetMass( newMass );
    result.SetMomentum( x, y, z );
    result.SetTotalEnergy( totEnergy );
    result.SetPositionInNucleus( 0.0, 0.0, 0.0 );
    result.SetFormationTime(0.0);
    result.HasInitialStateParton(false);
    return result;
  }

---

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

• geant4.10.03.p01/source/processes/hadronic/util/include/G4ReactionProduct.hh
• geant4.10.03.p01/source/processes/hadronic/util/src/G4ReactionProduct.cc