G4INCL::ReflectionChannel Class Reference

#include <G4INCLReflectionChannel.hh>

Inheritance diagram for G4INCL::ReflectionChannel:

Collaboration diagram for G4INCL::ReflectionChannel:

Public Member Functions
	ReflectionChannel (Nucleus *n, Particle *p)
virtual	~ReflectionChannel ()
void	fillFinalState (FinalState *fs)
Public Member Functions inherited from G4INCL::IChannel
	IChannel ()
virtual	~IChannel ()
FinalState *	getFinalState ()

Detailed Description

Definition at line 48 of file G4INCLReflectionChannel.hh.

Constructor & Destructor Documentation

G4INCL::ReflectionChannel::ReflectionChannel	(	Nucleus *	n,
		Particle *	p
	)

Definition at line 49 of file G4INCLReflectionChannel.cc.

    :theNucleus(n),theParticle(p)
  {
  }

G4INCL::ReflectionChannel::~ReflectionChannel ( )

virtual

Definition at line 54 of file G4INCLReflectionChannel.cc.

  {
  }

Member Function Documentation

void G4INCL::ReflectionChannel::fillFinalState ( FinalState *  fs )

virtual

Implements G4INCL::IChannel.

Definition at line 58 of file G4INCLReflectionChannel.cc.

                                                       {
    fs->setTotalEnergyBeforeInteraction(theParticle->getEnergy() - theParticle->getPotentialEnergy());

    const ThreeVector &oldMomentum = theParticle->getMomentum();
    const ThreeVector thePosition = theParticle->getPosition();
    G4double pspr = thePosition.dot(oldMomentum);
    if(pspr>=0) { // This means that the particle is trying to leave; perform a reflection
      const G4double x2cour = thePosition.mag2();
      const ThreeVector newMomentum = oldMomentum - (thePosition * (2.0 * pspr/x2cour));
      const G4double deltaP2 = (newMomentum-oldMomentum).mag2();
      theParticle->setMomentum(newMomentum);
      const G4double minDeltaP2 = sinMinReflectionAngleSquaredOverFour * newMomentum.mag2();
      if(deltaP2 < minDeltaP2) { // Avoid extremely small reflection angles
        theParticle->setPosition(thePosition * positionScalingFactor);
        INCL_DEBUG("Reflection angle for particle " << theParticle->getID() << " was too tangential: " << '\n'
            << "  " << deltaP2 << "=deltaP2<minDeltaP2=" << minDeltaP2 << '\n'
            << "  Resetting the particle position to ("
            << thePosition.getX() << ", "
            << thePosition.getY() << ", "
            << thePosition.getZ() << ")" << '\n');
      }
      theNucleus->updatePotentialEnergy(theParticle);
    } else { // The particle momentum is already directed towards the inside of the nucleus; do nothing
      // ...but make sure this only happened because of the frozen propagation
// assert(theParticle->getPosition().dot(theParticle->getPropagationVelocity())>0.);
    }

    theParticle->thawPropagation();
    fs->addModifiedParticle(theParticle);
  }

Here is the call graph for this function:

---

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

• geant4.10.03.p01/source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLReflectionChannel.hh
• geant4.10.03.p01/source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLReflectionChannel.cc