#include <G4INCLReflectionChannel.hh>

Inheritance diagram for G4INCL::ReflectionChannel:

Public Member Functions
	ReflectionChannel (Nucleus n, Particle p)

virtual	~ReflectionChannel ()

FinalState *	getFinalState ()

Public Member Functions inherited from G4INCL::IChannel
	IChannel ()

virtual	~IChannel ()

virtual G4INCL::FinalState *	getFinalState ()=0

Detailed Description

Definition at line 48 of file G4INCLReflectionChannel.hh.

Constructor & Destructor Documentation

◆ ReflectionChannel()

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

Definition at line 50 of file G4INCLReflectionChannel.cc.

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

◆ ~ReflectionChannel()

G4INCL::ReflectionChannel::~ReflectionChannel ( )

virtual

Definition at line 55 of file G4INCLReflectionChannel.cc.

56 {

57 }

Member Function Documentation

◆ getFinalState()

FinalState * G4INCL::ReflectionChannel::getFinalState ( )

virtual

Implements G4INCL::IChannel.

Definition at line 59 of file G4INCLReflectionChannel.cc.

  {
    FinalState *fs = new FinalState(); // Create final state for the output
    fs->setTotalEnergyBeforeInteraction(theParticle->getEnergy() - theParticle->getPotentialEnergy());
 
    const ThreeVector &oldMomentum = theParticle->getMomentum();
    G4double pspr = theParticle->getPosition().dot(oldMomentum);
    if(pspr>=0) { // This means that the particle is trying to leave; perform a reflection
      const G4double x2cour = theParticle->getPosition().mag2();
      const ThreeVector newMomentum = oldMomentum - (theParticle->getPosition() * (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(theParticle->getPosition() * positionScalingFactor);
        DEBUG("Reflection angle for particle " << theParticle->getID() << " was too tangential: " << std::endl
            << "  " << deltaP2 << "=deltaP2<minDeltaP2=" << minDeltaP2 << std::endl
            << "  Resetting the particle position to ("
            << theParticle->getPosition().getX() << ", "
            << theParticle->getPosition().getY() << ", "
            << theParticle->getPosition().getZ() << ")" << std::endl);
      }
      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);
    return fs;
  }