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.

References G4cout, and G4endl.

  : 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.

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

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

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.

References G4cout, and G4endl.

                                           {
  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.

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

{ return current_zone; }

G4int G4CascadParticle::getGeneration ( ) const

inline

Definition at line 77 of file G4CascadParticle.hh.

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

{ return generation; }

G4int G4CascadParticle::getHistoryId ( ) const

inline

Definition at line 80 of file G4CascadParticle.hh.

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

{ return historyId; }

G4LorentzVector G4CascadParticle::getMomentum ( ) const

inline

Definition at line 83 of file G4CascadParticle.hh.

References G4InuclParticle::getMomentum().

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

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

G4int G4CascadParticle::getNumberOfReflections ( ) const

inline

Definition at line 98 of file G4CascadParticle.hh.

Referenced by G4IntraNucleiCascader::generateCascade().

{ return reflectionCounter; }

const G4InuclElementaryParticle& G4CascadParticle::getParticle ( ) const

inline

Definition at line 74 of file G4CascadParticle.hh.

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

{ return theParticle; }

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.

References CLHEP::Hep3Vector::dot(), G4cout, G4endl, getMomentum(), CLHEP::Hep3Vector::mag2(), G4InuclParticleNames::pp, and CLHEP::HepLorentzVector::vect().

Referenced by G4NucleiModel::generateInteractionPartners().

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

const G4ThreeVector& G4CascadParticle::getPosition ( ) const

inline

Definition at line 91 of file G4CascadParticle.hh.

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

{ return position; }

void G4CascadParticle::incrementCurrentPath ( G4double  npath )

inline

Definition at line 104 of file G4CascadParticle.hh.

Referenced by G4NucleiModel::generateParticleFate().

{ current_path += npath; }

void G4CascadParticle::incrementReflectionCounter ( )

inline

Definition at line 94 of file G4CascadParticle.hh.

Referenced by G4NucleiModel::boundaryTransition().

                                    {
    reflectionCounter++; 
    reflected = true; 
  }

void G4CascadParticle::initializePath ( G4double  npath )

inline

Definition at line 103 of file G4CascadParticle.hh.

Referenced by G4IntraNucleiCascader::processSecondary().

{ current_path = npath; }

G4bool G4CascadParticle::movingInsideNuclei ( ) const

inline

Definition at line 111 of file G4CascadParticle.hh.

Referenced by G4NucleiModel::boundaryTransition().

{ return movingIn; }

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 185 of file G4CascadParticle.cc.

References G4endl.

Referenced by operator<<().

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

void G4CascadParticle::propagateAlongThePath

(

G4double

path

)

Definition at line 169 of file G4CascadParticle.cc.

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

Referenced by G4NucleiModel::generateParticleFate().

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

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

G4bool G4CascadParticle::reflectedNow ( ) const

inline

Definition at line 101 of file G4CascadParticle.hh.

Referenced by G4NucleiModel::worthToPropagate().

{ return reflected; }

void G4CascadParticle::resetReflection ( )

inline

Definition at line 100 of file G4CascadParticle.hh.

Referenced by G4NucleiModel::boundaryTransition().

{ reflected = false; }

void G4CascadParticle::setGeneration ( G4int  gen )

inline

Definition at line 78 of file G4CascadParticle.hh.

Referenced by G4IntraNucleiCascader::processSecondary().

{ generation = gen; }

void G4CascadParticle::setHistoryId ( G4int  id )

inline

Definition at line 81 of file G4CascadParticle.hh.

Referenced by G4CascadeHistory::AssignHistoryID().

{ historyId = id; }

void G4CascadParticle::setMovingInsideNuclei ( G4bool  isMovingIn = true )

inline

Definition at line 110 of file G4CascadParticle.hh.

Referenced by G4IntraNucleiCascader::processSecondary().

{ movingIn = isMovingIn; }

void G4CascadParticle::updateParticleMomentum ( const G4LorentzVector &  mom )

inline

Definition at line 87 of file G4CascadParticle.hh.

References G4InuclParticle::setMomentum().

Referenced by G4NucleiModel::boundaryTransition().

                                                          {
    theParticle.setMomentum(mom);
  }

void G4CascadParticle::updatePosition ( const G4ThreeVector &  pos )

inline

Definition at line 92 of file G4CascadParticle.hh.

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

{ position = pos; }

void G4CascadParticle::updateZone ( G4int  izone )

inline

Definition at line 107 of file G4CascadParticle.hh.

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

{ current_zone = izone; }

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

inline

Definition at line 116 of file G4CascadParticle.hh.

Referenced by G4NucleiModel::generateInteractionLength().

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

---

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

• G4CascadParticle.hh
• G4CascadParticle.cc