G4INCL::ReflectionChannel Class Reference

#include <G4INCLReflectionChannel.hh>

Inheritance diagram for G4INCL::ReflectionChannel:

Public Member Functions
	ReflectionChannel (Nucleus *n, Particle *p)
virtual	~ReflectionChannel ()
FinalState *	getFinalState ()

---

Detailed Description

Definition at line 46 of file G4INCLReflectionChannel.hh.

---

Constructor & Destructor Documentation

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

Definition at line 48 of file G4INCLReflectionChannel.cc.

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

G4INCL::ReflectionChannel::~ReflectionChannel

(

)

[virtual]

Definition at line 53 of file G4INCLReflectionChannel.cc.

  {
  }

---

Member Function Documentation

FinalState * G4INCL::ReflectionChannel::getFinalState

(

)

[virtual]

Implements G4INCL::IChannel.

Definition at line 57 of file G4INCLReflectionChannel.cc.

References G4INCL::FinalState::addModifiedParticle(), DEBUG, G4INCL::ThreeVector::dot(), G4INCL::Particle::getEnergy(), G4INCL::Particle::getID(), G4INCL::Particle::getMomentum(), G4INCL::Particle::getPosition(), G4INCL::Particle::getPotentialEnergy(), G4INCL::ThreeVector::getX(), G4INCL::ThreeVector::getY(), G4INCL::ThreeVector::getZ(), G4INCL::ThreeVector::mag2(), G4INCL::Particle::setMomentum(), G4INCL::Particle::setPosition(), G4INCL::FinalState::setTotalEnergyBeforeInteraction(), G4INCL::Particle::thawPropagation(), and G4INCL::Nucleus::updatePotentialEnergy().

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

---

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

• G4INCLReflectionChannel.hh
• G4INCLReflectionChannel.cc

---