G4CascadParticle Class Reference

#include <G4CascadParticle.hh>

Public Member Functions
	G4CascadParticle ()
	G4CascadParticle (const G4InuclElementaryParticle &particle, const G4ThreeVector &pos, G4int izone, G4double cpath, G4int gen)
	~G4CascadParticle ()
	G4CascadParticle (const G4CascadParticle &cpart)
G4CascadParticle &	operator= (const G4CascadParticle &cpart)
void	fill (const G4InuclElementaryParticle &particle, const G4ThreeVector &pos, G4int izone, G4double cpath, G4int gen)
const G4InuclElementaryParticle &	getParticle () const
G4InuclElementaryParticle &	getParticle ()
G4int	getGeneration () const
void	setGeneration (G4int gen)
G4int	getHistoryId () const
void	setHistoryId (G4int id)
G4LorentzVector	getMomentum () const
void	updateParticleMomentum (const G4LorentzVector &mom)
const G4ThreeVector &	getPosition () const
void	updatePosition (const G4ThreeVector &pos)
void	incrementReflectionCounter ()
G4int	getNumberOfReflections () const
void	resetReflection ()
G4bool	reflectedNow () const
void	initializePath (G4double npath)
void	incrementCurrentPath (G4double npath)
G4double	getCurrentPath () const
void	updateZone (G4int izone)
G4int	getCurrentZone () const
void	setMovingInsideNuclei (G4bool isMovingIn=true)
G4bool	movingInsideNuclei () const
G4double	getPathToTheNextZone (G4double rz_in, G4double rz_out)
void	propagateAlongThePath (G4double path)
G4bool	young (G4double young_path_cut, G4double cpath) const
void	print (std::ostream &os) const

Detailed Description

Definition at line 49 of file G4CascadParticle.hh.

Constructor & Destructor Documentation

G4CascadParticle() [1/3]

G4CascadParticle::G4CascadParticle

(

)

Definition at line 46 of file G4CascadParticle.cc.

  : verboseLevel(G4CascadeParameters::verbose()),
    current_zone(-1), current_path(-1.), movingIn(false),
    reflectionCounter(0), reflected(false), generation(-1), historyId(-1) {
  if (verboseLevel > 3) {
    G4cout << " >>> G4CascadParticle::G4CascadParticle" << G4endl;
  }
}

G4CascadParticle() [2/3]

G4CascadParticle::G4CascadParticle	(	const G4InuclElementaryParticle &	particle,
		const G4ThreeVector &	pos,
		G4int	izone,
		G4double	cpath,
		G4int	gen )

Definition at line 55 of file G4CascadParticle.cc.

    : verboseLevel(G4CascadeParameters::verbose()),
      theParticle(particle), position(pos), 
      current_zone(izone), current_path(cpath), movingIn(true),
      reflectionCounter(0), reflected(false), generation(gen), historyId(-1) {
  if (verboseLevel > 3) {
    G4cout << " >>> G4CascadParticle::G4CascadParticle "
       << particle.getDefinition()->GetParticleName() 
       << " @ " << pos << G4endl;
  }
}

~G4CascadParticle()

G4CascadParticle::~G4CascadParticle ( )

inline

Definition at line 59 of file G4CascadParticle.hh.

{;}         // No subclasses allowed

G4CascadParticle() [3/3]

G4CascadParticle::G4CascadParticle ( const G4CascadParticle & cpart )

inline

Definition at line 63 of file G4CascadParticle.hh.

{ *this = cpart; }

Member Function Documentation

fill()

void G4CascadParticle::fill	(	const G4InuclElementaryParticle &	particle,
		const G4ThreeVector &	pos,
		G4int	izone,
		G4double	cpath,
		G4int	gen )

Definition at line 91 of file G4CascadParticle.cc.

                                           {
  if (verboseLevel > 3) G4cout << " >>> G4CascadParticle::fill" << G4endl;
 
  theParticle = particle;
  position = pos;
  current_zone = izone;
  current_path = cpath;
  movingIn = true;
  reflectionCounter = 0;
  reflected = false;
  generation = gen;
  historyId = -1;   // New particle, no history entry yet
}

getCurrentPath()

G4double G4CascadParticle::getCurrentPath ( ) const

inline

Definition at line 104 of file G4CascadParticle.hh.

{ return current_path; }

getCurrentZone()

G4int G4CascadParticle::getCurrentZone ( ) const

inline

Definition at line 107 of file G4CascadParticle.hh.

{ return current_zone; }

Referenced by G4NucleiModel::boundaryTransition(), G4NucleiModel::choosePointAlongTraj(), G4IntraNucleiCascader::decayTrappedParticle(), G4NucleiModel::generateInteractionPartners(), G4NucleiModel::generateParticleFate(), G4NucleiModel::inverseMeanFreePath(), G4CascadeHistory::PrintEntry(), G4NucleiModel::stillInside(), and G4NucleiModel::worthToPropagate().

getGeneration()

G4int G4CascadParticle::getGeneration ( ) const

inline

Definition at line 76 of file G4CascadParticle.hh.

{ return generation; }

Referenced by G4IntraNucleiCascader::decayTrappedParticle(), G4NucleiModel::generateParticleFate(), G4NucleiModel::isProjectile(), and G4CascadeHistory::PrintEntry().

getHistoryId()

G4int G4CascadParticle::getHistoryId ( ) const

inline

Definition at line 79 of file G4CascadParticle.hh.

{ return historyId; }

Referenced by G4CascadeHistory::AddEntry(), G4CascadeHistory::AssignHistoryID(), and G4CascadeHistory::DropEntry().

getMomentum()

G4LorentzVector G4CascadParticle::getMomentum ( ) const

inline

Definition at line 82 of file G4CascadParticle.hh.

                                      {     // Can't return ref; temporary
    return theParticle.getMomentum(); 
  }

Referenced by G4NucleiModel::boundaryTransition(), G4NucleiModel::choosePointAlongTraj(), G4NucleiModel::generateInteractionPartners(), getPathToTheNextZone(), G4CascadeHistory::PrintEntry(), and propagateAlongThePath().

getNumberOfReflections()

G4int G4CascadParticle::getNumberOfReflections ( ) const

inline

Definition at line 97 of file G4CascadParticle.hh.

{ return reflectionCounter; }

Referenced by G4IntraNucleiCascader::generateCascade().

getParticle() [1/2]

G4InuclElementaryParticle & G4CascadParticle::getParticle ( )

inline

Definition at line 74 of file G4CascadParticle.hh.

{ return theParticle; }

getParticle() [2/2]

const G4InuclElementaryParticle & G4CascadParticle::getParticle ( ) const

inline

Definition at line 73 of file G4CascadParticle.hh.

{ return theParticle; }

Referenced by G4CollisionOutput::addOutgoingParticle(), G4NucleiModel::boundaryTransition(), G4IntraNucleiCascader::decayTrappedParticle(), G4NucleiModel::forceFirst(), G4IntraNucleiCascader::generateCascade(), G4NucleiModel::generateInteractionPartners(), G4NucleiModel::generateParticleFate(), G4CascadeHistory::GuessTarget(), G4NucleiModel::inverseMeanFreePath(), G4ParticleLargerBeta::operator()(), G4ParticleLargerEkin::operator()(), G4IntraNucleiCascader::particleCanInteract(), G4CascadeHistory::PrintEntry(), G4IntraNucleiCascader::processSecondary(), G4IntraNucleiCascader::processTrappedParticle(), and G4NucleiModel::worthToPropagate().

getPathToTheNextZone()

G4double G4CascadParticle::getPathToTheNextZone	(	G4double	rz_in,
		G4double	rz_out )

Definition at line 108 of file G4CascadParticle.cc.

                                         {
  if (verboseLevel > 3) {
    G4cout << " >>> G4CascadParticle::getPathToTheNextZone rz_in " << rz_in
       << " rz_out " << rz_out << G4endl;
  }
 
  const G4LorentzVector& mom = getMomentum();
 
  G4double path = -1.0;
  G4double rp = mom.vect().dot(position);
  G4double rr = position.mag2();
  G4double pp = mom.vect().mag2();
 
  if (std::abs(pp) < 1e-9) {    // Cut-off for "at rest" is 1 eV momentum
    if (verboseLevel > 3) G4cout << " at rest; path length is zero" << G4endl;
 
    if (current_zone == 0) movingIn = false;    // Allow 'stuck' to escape
    return 0.;
  }
 
  G4double ra = rr - rp * rp / pp;
  pp = std::sqrt(pp);
  G4double ds;
  G4double d2;
 
  if (verboseLevel > 3) {
    G4cout << " current_zone " << current_zone << " rr " << rr
       << " rp " << rp << " pp " << pp << " ra " << ra << G4endl;
  }
 
  if (current_zone == 0 || rp > 0.0) {
    d2 = rz_out * rz_out - ra;
    if (d2 > 0.0) {
      ds = 1.0;
      movingIn = false;
    } else {
      d2 = rz_in * rz_in - ra;
      ds = -1.0;
      movingIn = true;
    }
  } else { 
    d2 = rz_in * rz_in - ra;
    if (d2 > 0.0) {
      ds = -1.0;
      movingIn = true;
    } else {
      d2 = rz_out * rz_out - ra;
      ds = 1.0;
      movingIn = false;
    }
  }
 
  if (verboseLevel > 3) G4cout << " ds " << ds << " d2 " << d2 << G4endl;
 
  if (d2 < 0.0 && d2 > -1e-6) d2 = 0.;      // Account for round-off
 
  if (d2 > 0.0) path = ds * std::sqrt(d2) - rp / pp;    // Avoid FPE failure
 
  return path;    
}

Referenced by G4NucleiModel::generateInteractionPartners().

getPosition()

const G4ThreeVector & G4CascadParticle::getPosition ( ) const

inline

Definition at line 90 of file G4CascadParticle.hh.

{ return position; }

Referenced by G4NucleiModel::boundaryTransition(), G4NucleiModel::choosePointAlongTraj(), G4IntraNucleiCascader::decayTrappedParticle(), G4NucleiModel::generateParticleFate(), and G4CascadeHistory::PrintEntry().

incrementCurrentPath()

void G4CascadParticle::incrementCurrentPath ( G4double npath )

inline

Definition at line 103 of file G4CascadParticle.hh.

{ current_path += npath; }

Referenced by G4NucleiModel::generateParticleFate().

incrementReflectionCounter()

void G4CascadParticle::incrementReflectionCounter ( )

inline

Definition at line 93 of file G4CascadParticle.hh.

                                    {
    reflectionCounter++; 
    reflected = true; 
  }

Referenced by G4NucleiModel::boundaryTransition().

initializePath()

void G4CascadParticle::initializePath ( G4double npath )

inline

Definition at line 102 of file G4CascadParticle.hh.

{ current_path = npath; }

Referenced by G4IntraNucleiCascader::processSecondary().

movingInsideNuclei()

G4bool G4CascadParticle::movingInsideNuclei ( ) const

inline

Definition at line 110 of file G4CascadParticle.hh.

{ return movingIn; }

Referenced by G4NucleiModel::boundaryTransition().

operator=()

G4CascadParticle & G4CascadParticle::operator=

(

const G4CascadParticle &

cpart

)

Definition at line 72 of file G4CascadParticle.cc.

                                                                           {
  if (&cpart != this) {             // Avoid unnecessary work
    verboseLevel = cpart.verboseLevel;
    theParticle = cpart.theParticle;
    position = cpart.position;
    current_zone = cpart.current_zone;
    current_path = cpart.current_path;
    movingIn = cpart.movingIn;
    reflectionCounter = cpart.reflectionCounter;
    reflected = cpart.reflected;
    generation = cpart.generation;
    historyId = cpart.historyId;
  }
 
  return *this;
}

print()

void G4CascadParticle::print

(

std::ostream &

os

)

const

Definition at line 186 of file G4CascadParticle.cc.

                                                 {
  os << " pos " << position << " zone " << current_zone
     << " current_path " << current_path
     << " reflectionCounter " << reflectionCounter << G4endl
     << theParticle << G4endl;     
}

Referenced by operator<<().

propagateAlongThePath()

void G4CascadParticle::propagateAlongThePath

(

G4double

path

)

Definition at line 170 of file G4CascadParticle.cc.

                                                          {
  if (verboseLevel > 3) {
    G4cout << " >>> G4CascadParticle::propagateAlongThePath" << G4endl;
  }
 
  position += getMomentum().vect().unit()*path;
}

Referenced by G4NucleiModel::generateParticleFate().

reflectedNow()

G4bool G4CascadParticle::reflectedNow ( ) const

inline

Definition at line 100 of file G4CascadParticle.hh.

{ return reflected; }

Referenced by G4NucleiModel::worthToPropagate().

resetReflection()

void G4CascadParticle::resetReflection ( )

inline

Definition at line 99 of file G4CascadParticle.hh.

{ reflected = false; }

Referenced by G4NucleiModel::boundaryTransition().

setGeneration()

void G4CascadParticle::setGeneration ( G4int gen )

inline

Definition at line 77 of file G4CascadParticle.hh.

{ generation = gen; }

Referenced by G4IntraNucleiCascader::processSecondary().

setHistoryId()

void G4CascadParticle::setHistoryId ( G4int id )

inline

Definition at line 80 of file G4CascadParticle.hh.

{ historyId = id; }

Referenced by G4CascadeHistory::AssignHistoryID().

setMovingInsideNuclei()

void G4CascadParticle::setMovingInsideNuclei ( G4bool isMovingIn = true )

inline

Definition at line 109 of file G4CascadParticle.hh.

{ movingIn = isMovingIn; }

Referenced by G4IntraNucleiCascader::processSecondary().

updateParticleMomentum()

void G4CascadParticle::updateParticleMomentum ( const G4LorentzVector & mom )

inline

Definition at line 86 of file G4CascadParticle.hh.

                                                          {
    theParticle.setMomentum(mom);
  }

Referenced by G4NucleiModel::boundaryTransition().

updatePosition()

void G4CascadParticle::updatePosition ( const G4ThreeVector & pos )

inline

Definition at line 91 of file G4CascadParticle.hh.

{ position = pos; }

Referenced by G4NucleiModel::choosePointAlongTraj(), G4NucleiModel::generateParticleFate(), and G4IntraNucleiCascader::processSecondary().

updateZone()

void G4CascadParticle::updateZone ( G4int izone )

inline

Definition at line 106 of file G4CascadParticle.hh.

{ current_zone = izone; }

Referenced by G4NucleiModel::boundaryTransition(), G4NucleiModel::choosePointAlongTraj(), and G4IntraNucleiCascader::processSecondary().

young()

G4bool G4CascadParticle::young ( G4double young_path_cut, G4double cpath ) const

inline

Definition at line 115 of file G4CascadParticle.hh.

                                                              { 
    return ((current_path < 1000.) && (cpath < young_path_cut));
  }

Referenced by G4NucleiModel::generateInteractionLength().

---

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

• G4CascadParticle.hh
• G4CascadParticle.cc