G4CascadParticle Class Reference

#include <G4CascadParticle.hh>

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

Private Attributes
G4double	current_path
G4int	current_zone
G4int	generation
G4int	historyId
G4bool	movingIn
G4ThreeVector	position
G4bool	reflected
G4int	reflectionCounter
G4InuclElementaryParticle	theParticle
G4int	verboseLevel

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;
  }
}

References G4cout, G4endl, and verboseLevel.

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;
  }
}

References G4cout, G4endl, G4InuclParticle::getDefinition(), G4ParticleDefinition::GetParticleName(), pos, and verboseLevel.

~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
}

References current_path, current_zone, G4cout, G4endl, generation, historyId, movingIn, pos, reflected, reflectionCounter, theParticle, and verboseLevel.

getCurrentPath()

G4double G4CascadParticle::getCurrentPath ( ) const

inline

Definition at line 104 of file G4CascadParticle.hh.

{ return current_path; }

References current_path.

getCurrentZone()

G4int G4CascadParticle::getCurrentZone ( ) const

inline

Definition at line 107 of file G4CascadParticle.hh.

{ return current_zone; }

References 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; }

References 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; }

References 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(); 
  }

References G4InuclParticle::getMomentum(), and theParticle.

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; }

References reflectionCounter.

Referenced by G4IntraNucleiCascader::generateCascade().

getParticle() [1/2]

G4InuclElementaryParticle & G4CascadParticle::getParticle ( )

inline

Definition at line 74 of file G4CascadParticle.hh.

{ return theParticle; }

References theParticle.

getParticle() [2/2]

const G4InuclElementaryParticle & G4CascadParticle::getParticle ( ) const

inline

Definition at line 73 of file G4CascadParticle.hh.

{ return theParticle; }

References 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;    
}

References current_zone, d2, CLHEP::Hep3Vector::dot(), anonymous_namespace{G4QuasiElRatios.cc}::ds, G4cout, G4endl, getMomentum(), CLHEP::Hep3Vector::mag2(), movingIn, G4InuclParticleNames::pp, CLHEP::HepLorentzVector::vect(), and verboseLevel.

Referenced by G4NucleiModel::generateInteractionPartners().

getPosition()

const G4ThreeVector & G4CascadParticle::getPosition ( ) const

inline

Definition at line 90 of file G4CascadParticle.hh.

{ return position; }

References 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; }

References current_path.

Referenced by G4NucleiModel::generateParticleFate().

incrementReflectionCounter()

void G4CascadParticle::incrementReflectionCounter ( )

inline

Definition at line 93 of file G4CascadParticle.hh.

                                    {
    reflectionCounter++; 
    reflected = true; 
  }

References reflected, and reflectionCounter.

Referenced by G4NucleiModel::boundaryTransition().

initializePath()

void G4CascadParticle::initializePath ( G4double  npath )

inline

Definition at line 102 of file G4CascadParticle.hh.

{ current_path = npath; }

References current_path.

Referenced by G4IntraNucleiCascader::processSecondary().

movingInsideNuclei()

G4bool G4CascadParticle::movingInsideNuclei ( ) const

inline

Definition at line 110 of file G4CascadParticle.hh.

{ return movingIn; }

References 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;
}

References current_path, current_zone, generation, historyId, movingIn, position, reflected, reflectionCounter, theParticle, and verboseLevel.

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;     
}

References current_path, current_zone, G4endl, reflectionCounter, and theParticle.

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;
}

References G4cout, G4endl, getMomentum(), CLHEP::Hep3Vector::unit(), CLHEP::HepLorentzVector::vect(), and verboseLevel.

Referenced by G4NucleiModel::generateParticleFate().

reflectedNow()

G4bool G4CascadParticle::reflectedNow ( ) const

inline

Definition at line 100 of file G4CascadParticle.hh.

{ return reflected; }

References reflected.

Referenced by G4NucleiModel::worthToPropagate().

resetReflection()

void G4CascadParticle::resetReflection ( )

inline

Definition at line 99 of file G4CascadParticle.hh.

{ reflected = false; }

References reflected.

Referenced by G4NucleiModel::boundaryTransition().

setGeneration()

void G4CascadParticle::setGeneration ( G4int  gen )

inline

Definition at line 77 of file G4CascadParticle.hh.

{ generation = gen; }

References generation.

Referenced by G4IntraNucleiCascader::processSecondary().

setHistoryId()

void G4CascadParticle::setHistoryId ( G4int  id )

inline

Definition at line 80 of file G4CascadParticle.hh.

{ historyId = id; }

References historyId.

Referenced by G4CascadeHistory::AssignHistoryID().

setMovingInsideNuclei()

void G4CascadParticle::setMovingInsideNuclei ( G4bool  isMovingIn = true )

inline

Definition at line 109 of file G4CascadParticle.hh.

{ movingIn = isMovingIn; }

References movingIn.

Referenced by G4IntraNucleiCascader::processSecondary().

updateParticleMomentum()

void G4CascadParticle::updateParticleMomentum ( const G4LorentzVector &  mom )

inline

Definition at line 86 of file G4CascadParticle.hh.

                                                          {
    theParticle.setMomentum(mom);
  }

References G4InuclParticle::setMomentum(), and theParticle.

Referenced by G4NucleiModel::boundaryTransition().

updatePosition()

void G4CascadParticle::updatePosition ( const G4ThreeVector &  pos )

inline

Definition at line 91 of file G4CascadParticle.hh.

{ position = pos; }

References 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; }

References current_zone.

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));
  }

References current_path.

Referenced by G4NucleiModel::generateInteractionLength().

Field Documentation

current_path

G4double G4CascadParticle::current_path

private

Definition at line 126 of file G4CascadParticle.hh.

Referenced by fill(), getCurrentPath(), incrementCurrentPath(), initializePath(), operator=(), print(), and young().

current_zone

G4int G4CascadParticle::current_zone

private

Definition at line 125 of file G4CascadParticle.hh.

Referenced by fill(), getCurrentZone(), getPathToTheNextZone(), operator=(), print(), and updateZone().

generation

G4int G4CascadParticle::generation

private

Definition at line 130 of file G4CascadParticle.hh.

Referenced by fill(), getGeneration(), operator=(), and setGeneration().

historyId

G4int G4CascadParticle::historyId

private

Definition at line 131 of file G4CascadParticle.hh.

Referenced by fill(), getHistoryId(), operator=(), and setHistoryId().

movingIn

G4bool G4CascadParticle::movingIn

private

Definition at line 127 of file G4CascadParticle.hh.

Referenced by fill(), getPathToTheNextZone(), movingInsideNuclei(), operator=(), and setMovingInsideNuclei().

position

G4ThreeVector G4CascadParticle::position

private

Definition at line 124 of file G4CascadParticle.hh.

Referenced by getPosition(), and operator=().

reflected

G4bool G4CascadParticle::reflected

private

Definition at line 129 of file G4CascadParticle.hh.

Referenced by fill(), incrementReflectionCounter(), operator=(), reflectedNow(), and resetReflection().

reflectionCounter

G4int G4CascadParticle::reflectionCounter

private

Definition at line 128 of file G4CascadParticle.hh.

Referenced by fill(), getNumberOfReflections(), incrementReflectionCounter(), operator=(), and print().

theParticle

G4InuclElementaryParticle G4CascadParticle::theParticle

private

Definition at line 123 of file G4CascadParticle.hh.

Referenced by fill(), getMomentum(), getParticle(), operator=(), print(), and updateParticleMomentum().

verboseLevel

G4int G4CascadParticle::verboseLevel

private

Definition at line 122 of file G4CascadParticle.hh.

Referenced by fill(), G4CascadParticle(), getPathToTheNextZone(), operator=(), and propagateAlongThePath().

---

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

• source/processes/hadronic/models/cascade/cascade/include/G4CascadParticle.hh
• source/processes/hadronic/models/cascade/cascade/src/G4CascadParticle.cc