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 50 of file G4CascadParticle.hh.

Constructor & Destructor Documentation

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::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 ( )

inline

Definition at line 60 of file G4CascadParticle.hh.

{;}         // No subclasses allowed

G4CascadParticle::G4CascadParticle ( const G4CascadParticle &  cpart )

inline

Definition at line 64 of file G4CascadParticle.hh.

{ *this = cpart; }

Member Function Documentation

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
}

G4double G4CascadParticle::getCurrentPath ( ) const

inline

Definition at line 105 of file G4CascadParticle.hh.

{ return current_path; }

G4int G4CascadParticle::getCurrentZone ( ) const

inline

Definition at line 108 of file G4CascadParticle.hh.

{ return current_zone; }

Here is the caller graph for this function:

G4int G4CascadParticle::getGeneration ( ) const

inline

Definition at line 77 of file G4CascadParticle.hh.

{ return generation; }

Here is the caller graph for this function:

G4int G4CascadParticle::getHistoryId ( ) const

inline

Definition at line 80 of file G4CascadParticle.hh.

{ return historyId; }

Here is the caller graph for this function:

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:

G4int G4CascadParticle::getNumberOfReflections ( ) const

inline

Definition at line 98 of file G4CascadParticle.hh.

{ return reflectionCounter; }

Here is the caller graph for this function:

const G4InuclElementaryParticle& G4CascadParticle::getParticle ( ) const

inline

Definition at line 74 of file G4CascadParticle.hh.

{ return theParticle; }

Here is the caller graph for this function:

G4InuclElementaryParticle& G4CascadParticle::getParticle ( )

inline

Definition at line 75 of file G4CascadParticle.hh.

{ return theParticle; }

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:

const G4ThreeVector& G4CascadParticle::getPosition ( ) const

inline

Definition at line 91 of file G4CascadParticle.hh.

{ return position; }

Here is the caller graph for this function:

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:

void G4CascadParticle::incrementReflectionCounter ( )

inline

Definition at line 94 of file G4CascadParticle.hh.

                                    {
    reflectionCounter++; 
    reflected = true; 
  }

Here is the caller graph for this function:

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:

G4bool G4CascadParticle::movingInsideNuclei ( ) const

inline

Definition at line 111 of file G4CascadParticle.hh.

{ return movingIn; }

Here is the caller graph for this function:

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

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:

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:

G4bool G4CascadParticle::reflectedNow ( ) const

inline

Definition at line 101 of file G4CascadParticle.hh.

{ return reflected; }

Here is the caller graph for this function:

void G4CascadParticle::resetReflection ( )

inline

Definition at line 100 of file G4CascadParticle.hh.

{ reflected = false; }

Here is the caller graph for this function:

void G4CascadParticle::setGeneration ( G4int  gen )

inline

Definition at line 78 of file G4CascadParticle.hh.

{ generation = gen; }

Here is the caller graph for this function:

void G4CascadParticle::setHistoryId ( G4int  id )

inline

Definition at line 81 of file G4CascadParticle.hh.

{ historyId = id; }

Here is the caller graph for this function:

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:

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:

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:

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:

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 caller graph for this function:

---

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

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