G4CascadParticle Class Reference

#include <G4CascadParticle.hh>

Collaboration diagram for G4CascadParticle:

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

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

Detailed Description

Definition at line 50 of file G4CascadParticle.hh.

Constructor & Destructor Documentation

G4CascadParticle() [1/3]

G4CascadParticle::G4CascadParticle

(

)

Definition at line 47 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 57 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;
  }
}

Here is the call graph for this function:

~G4CascadParticle()

G4CascadParticle::~G4CascadParticle ( )

inline

Definition at line 60 of file G4CascadParticle.hh.

{;}         // No subclasses allowed

G4CascadParticle() [3/3]

G4CascadParticle::G4CascadParticle ( const G4CascadParticle &  cpart )

inline

Definition at line 64 of file G4CascadParticle.hh.

{ *this = cpart; }

Here is the call graph for this function:

Member Function Documentation

fill()

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

Definition at line 92 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
}

Here is the caller graph for this function:

getCurrentPath()

G4double G4CascadParticle::getCurrentPath ( ) const

inline

Definition at line 105 of file G4CascadParticle.hh.

{ return current_path; }

getCurrentZone()

G4int G4CascadParticle::getCurrentZone ( ) const

inline

Definition at line 108 of file G4CascadParticle.hh.

{ return current_zone; }

Here is the caller graph for this function:

getGeneration()

G4int G4CascadParticle::getGeneration ( ) const

inline

Definition at line 77 of file G4CascadParticle.hh.

{ return generation; }

Here is the caller graph for this function:

getHistoryId()

G4int G4CascadParticle::getHistoryId ( ) const

inline

Definition at line 80 of file G4CascadParticle.hh.

{ return historyId; }

Here is the caller graph for this function:

getMomentum()

G4LorentzVector G4CascadParticle::getMomentum ( ) const

inline

Definition at line 83 of file G4CascadParticle.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

getNumberOfReflections()

G4int G4CascadParticle::getNumberOfReflections ( ) const

inline

Definition at line 98 of file G4CascadParticle.hh.

{ return reflectionCounter; }

Here is the caller graph for this function:

getParticle() [1/2]

const G4InuclElementaryParticle& G4CascadParticle::getParticle ( ) const

inline

Definition at line 74 of file G4CascadParticle.hh.

{ return theParticle; }

Here is the caller graph for this function:

getParticle() [2/2]

G4InuclElementaryParticle& G4CascadParticle::getParticle ( )

inline

Definition at line 75 of file G4CascadParticle.hh.

{ return theParticle; }

getPathToTheNextZone()

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

Definition at line 109 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;    
}

Here is the call graph for this function:

Here is the caller graph for this function:

getPosition()

const G4ThreeVector& G4CascadParticle::getPosition ( ) const

inline

Definition at line 91 of file G4CascadParticle.hh.

{ return position; }

Here is the caller graph for this function:

incrementCurrentPath()

void G4CascadParticle::incrementCurrentPath ( G4double  npath )

inline

Definition at line 104 of file G4CascadParticle.hh.

{ current_path += npath; }

Here is the caller graph for this function:

incrementReflectionCounter()

void G4CascadParticle::incrementReflectionCounter ( )

inline

Definition at line 94 of file G4CascadParticle.hh.

                                    {
    reflectionCounter++; 
    reflected = true; 
  }

Here is the caller graph for this function:

initializePath()

void G4CascadParticle::initializePath ( G4double  npath )

inline

Definition at line 103 of file G4CascadParticle.hh.

{ current_path = npath; }

Here is the caller graph for this function:

movingInsideNuclei()

G4bool G4CascadParticle::movingInsideNuclei ( ) const

inline

Definition at line 111 of file G4CascadParticle.hh.

{ return movingIn; }

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

G4CascadParticle & G4CascadParticle::operator=

(

const G4CascadParticle &

cpart

)

Definition at line 73 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;
}

Here is the caller graph for this function:

print()

void G4CascadParticle::print

(

std::ostream &

os

)

const

Definition at line 187 of file G4CascadParticle.cc.

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

Here is the caller graph for this function:

propagateAlongThePath()

void G4CascadParticle::propagateAlongThePath

(

G4double

path

)

Definition at line 171 of file G4CascadParticle.cc.

                                                          {
  if (verboseLevel > 3) {
    G4cout << " >>> G4CascadParticle::propagateAlongThePath" << G4endl;
  }

  position += getMomentum().vect().unit()*path;
}

Here is the call graph for this function:

Here is the caller graph for this function:

reflectedNow()

G4bool G4CascadParticle::reflectedNow ( ) const

inline

Definition at line 101 of file G4CascadParticle.hh.

{ return reflected; }

Here is the caller graph for this function:

resetReflection()

void G4CascadParticle::resetReflection ( )

inline

Definition at line 100 of file G4CascadParticle.hh.

{ reflected = false; }

Here is the caller graph for this function:

setGeneration()

void G4CascadParticle::setGeneration ( G4int  gen )

inline

Definition at line 78 of file G4CascadParticle.hh.

{ generation = gen; }

Here is the caller graph for this function:

setHistoryId()

void G4CascadParticle::setHistoryId ( G4int  id )

inline

Definition at line 81 of file G4CascadParticle.hh.

{ historyId = id; }

Here is the caller graph for this function:

setMovingInsideNuclei()

void G4CascadParticle::setMovingInsideNuclei ( G4bool  isMovingIn = true )

inline

Definition at line 110 of file G4CascadParticle.hh.

{ movingIn = isMovingIn; }

Here is the caller graph for this function:

updateParticleMomentum()

void G4CascadParticle::updateParticleMomentum ( const G4LorentzVector &  mom )

inline

Definition at line 87 of file G4CascadParticle.hh.

                                                          {
    theParticle.setMomentum(mom);
  }

Here is the call graph for this function:

Here is the caller graph for this function:

updatePosition()

void G4CascadParticle::updatePosition ( const G4ThreeVector &  pos )

inline

Definition at line 92 of file G4CascadParticle.hh.

{ position = pos; }

Here is the caller graph for this function:

updateZone()

void G4CascadParticle::updateZone ( G4int  izone )

inline

Definition at line 107 of file G4CascadParticle.hh.

{ current_zone = izone; }

Here is the caller graph for this function:

young()

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

inline

Definition at line 116 of file G4CascadParticle.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

current_path

G4double G4CascadParticle::current_path

private

Definition at line 127 of file G4CascadParticle.hh.

current_zone

G4int G4CascadParticle::current_zone

private

Definition at line 126 of file G4CascadParticle.hh.

generation

G4int G4CascadParticle::generation

private

Definition at line 131 of file G4CascadParticle.hh.

historyId

G4int G4CascadParticle::historyId

private

Definition at line 132 of file G4CascadParticle.hh.

movingIn

G4bool G4CascadParticle::movingIn

private

Definition at line 128 of file G4CascadParticle.hh.

position

G4ThreeVector G4CascadParticle::position

private

Definition at line 125 of file G4CascadParticle.hh.

reflected

G4bool G4CascadParticle::reflected

private

Definition at line 130 of file G4CascadParticle.hh.

reflectionCounter

G4int G4CascadParticle::reflectionCounter

private

Definition at line 129 of file G4CascadParticle.hh.

theParticle

G4InuclElementaryParticle G4CascadParticle::theParticle

private

Definition at line 124 of file G4CascadParticle.hh.

verboseLevel

G4int G4CascadParticle::verboseLevel

private

Definition at line 123 of file G4CascadParticle.hh.

---

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

• G4CascadParticle.hh
• G4CascadParticle.cc