G4INCL::TransmissionChannel Class Reference

#include <G4INCLTransmissionChannel.hh>

Inheritance diagram for G4INCL::TransmissionChannel:

Public Member Functions
void	fillFinalState (FinalState *fs)
FinalState *	getFinalState ()
	TransmissionChannel (Nucleus *n, Particle *p)
	TransmissionChannel (Nucleus *n, Particle *p, const G4double kOut, const G4double cosR)
	TransmissionChannel (Nucleus *n, Particle *p, const G4double TOut)
virtual	~TransmissionChannel ()

Private Member Functions
G4double	initializeKineticEnergyOutside ()
	Kinetic energy of the transmitted particle. More...
void	particleLeaves ()
	Modify particle that leaves the nucleus. More...

Private Attributes
const G4double	cosRefractionAngle
	Cosine of the refraction angle. More...
const G4double	kineticEnergyOutside
	Kinetic energy of the particle outside the nucleus. More...
const G4double	pOutMag
	Momentum of the particle outside the nucleus. More...
const G4bool	refraction
	True if refraction should be applied. More...
Nucleus *const	theNucleus
Particle *const	theParticle

Detailed Description

Definition at line 48 of file G4INCLTransmissionChannel.hh.

Constructor & Destructor Documentation

TransmissionChannel() [1/3]

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

Definition at line 42 of file G4INCLTransmissionChannel.cc.

    : theNucleus(nucleus), theParticle(particle),
    refraction(false),
    pOutMag(0.),
    kineticEnergyOutside(initializeKineticEnergyOutside()),
    cosRefractionAngle(1.)
  {}

TransmissionChannel() [2/3]

G4INCL::TransmissionChannel::TransmissionChannel	(	Nucleus *	n,
		Particle *	p,
		const G4double	TOut
	)

Definition at line 50 of file G4INCLTransmissionChannel.cc.

    : theNucleus(nucleus), theParticle(particle),
    refraction(false),
    pOutMag(0.),
    kineticEnergyOutside(TOut),
    cosRefractionAngle(1.)
  {}

TransmissionChannel() [3/3]

G4INCL::TransmissionChannel::TransmissionChannel	(	Nucleus *	n,
		Particle *	p,
		const G4double	kOut,
		const G4double	cosR
	)

Definition at line 58 of file G4INCLTransmissionChannel.cc.

    : theNucleus(nucleus), theParticle(particle),
    refraction(true),
    pOutMag(kOut),
    kineticEnergyOutside(initializeKineticEnergyOutside()),
    cosRefractionAngle(cosR)
  {}

~TransmissionChannel()

G4INCL::TransmissionChannel::~TransmissionChannel ( )

virtual

Definition at line 66 of file G4INCLTransmissionChannel.cc.

{}

Member Function Documentation

fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 113 of file G4INCLTransmissionChannel.cc.

                                                         {
    G4double initialEnergy = 0.0;
    initialEnergy = theParticle->getEnergy() - theParticle->getPotentialEnergy();
    // Correction for real masses
    const G4int AParent = theNucleus->getA();
    const G4int ZParent = theNucleus->getZ();
    const G4int SParent = theNucleus->getS();
    initialEnergy += theParticle->getTableMass() - theParticle->getMass()
      + theParticle->getEmissionQValueCorrection(AParent,ZParent,SParent);
    particleLeaves();
    fs->setTotalEnergyBeforeInteraction(initialEnergy);
    fs->addOutgoingParticle(theParticle); // We write the particle down as outgoing
  }

References G4INCL::FinalState::addOutgoingParticle(), G4INCL::Particle::getA(), G4INCL::Particle::getEmissionQValueCorrection(), G4INCL::Particle::getEnergy(), G4INCL::Particle::getMass(), G4INCL::Particle::getPotentialEnergy(), G4INCL::Particle::getS(), G4INCL::Particle::getTableMass(), G4INCL::Particle::getZ(), particleLeaves(), G4INCL::FinalState::setTotalEnergyBeforeInteraction(), theNucleus, and theParticle.

getFinalState()

FinalState * G4INCL::IChannel::getFinalState ( )

inherited

Definition at line 50 of file G4INCLIChannel.cc.

                                      {
    FinalState *fs = new FinalState;
    fillFinalState(fs);
    return fs;
  }

References G4INCL::IChannel::fillFinalState().

initializeKineticEnergyOutside()

G4double G4INCL::TransmissionChannel::initializeKineticEnergyOutside ( )

private

Kinetic energy of the transmitted particle.

Calculate the kinetic energy of the particle outside the nucleus, if the value has not been provided as a pre-calculated argument to the constructor.

Definition at line 68 of file G4INCLTransmissionChannel.cc.

                                                               {
    // The particle energy outside the nucleus. Subtract the nuclear
    // potential from the kinetic energy when leaving the nucleus
    G4double TOut = theParticle->getEnergy()
      - theParticle->getPotentialEnergy()
      - theParticle->getMass();
 
    // Correction for real masses
    const G4int AParent = theNucleus->getA();
    const G4int ZParent = theNucleus->getZ();
    const G4int SParent = theNucleus->getS();
    const G4double theQValueCorrection = theParticle->getEmissionQValueCorrection(AParent,ZParent,SParent);
    TOut += theQValueCorrection;
    return TOut;
  }

References G4INCL::Particle::getA(), G4INCL::Particle::getEmissionQValueCorrection(), G4INCL::Particle::getEnergy(), G4INCL::Particle::getMass(), G4INCL::Particle::getPotentialEnergy(), G4INCL::Particle::getS(), G4INCL::Particle::getZ(), theNucleus, and theParticle.

particleLeaves()

void G4INCL::TransmissionChannel::particleLeaves ( )

private

Modify particle that leaves the nucleus.

Modify the particle momentum and/or position when the particle leaves the nucleus.

Definition at line 84 of file G4INCLTransmissionChannel.cc.

                                           {
 
    // Use the table mass in the outside world
    theParticle->setTableMass();
    theParticle->setPotentialEnergy(0.);
 
    if(refraction) {
      // Change the momentum direction
      // The magnitude of the particle momentum outside the nucleus will be
      // fixed by the kineticEnergyOutside variable. This is done in order to
      // avoid numerical inaccuracies.
      const ThreeVector &position = theParticle->getPosition();
      const G4double r2 = position.mag2();
      ThreeVector normal;
      if(r2>0.)
        normal = position / std::sqrt(r2);
 
      const ThreeVector &momentum = theParticle->getMomentum();
 
      const ThreeVector pOut = normal * (pOutMag * cosRefractionAngle) + momentum - normal * normal.dot(momentum);
// assert(std::fabs(pOut.mag()-pOutMag)<1.e-5);
 
      theParticle->setMomentum(pOut);
    }
    // Scaling factor for the particle momentum
    theParticle->setEnergy(kineticEnergyOutside + theParticle->getMass());
    theParticle->adjustMomentumFromEnergy();
  }

References G4INCL::Particle::adjustMomentumFromEnergy(), cosRefractionAngle, G4INCL::Particle::getMass(), G4INCL::Particle::getMomentum(), G4INCL::Particle::getPosition(), kineticEnergyOutside, CLHEP::normal(), pOutMag, refraction, G4INCL::Particle::setEnergy(), G4INCL::Particle::setMomentum(), G4INCL::Particle::setPotentialEnergy(), G4INCL::Particle::setTableMass(), and theParticle.

Referenced by fillFinalState().

Field Documentation

cosRefractionAngle

const G4double G4INCL::TransmissionChannel::cosRefractionAngle

private

Cosine of the refraction angle.

Definition at line 86 of file G4INCLTransmissionChannel.hh.

Referenced by particleLeaves().

kineticEnergyOutside

const G4double G4INCL::TransmissionChannel::kineticEnergyOutside

private

Kinetic energy of the particle outside the nucleus.

Definition at line 83 of file G4INCLTransmissionChannel.hh.

Referenced by particleLeaves().

pOutMag

const G4double G4INCL::TransmissionChannel::pOutMag

private

Momentum of the particle outside the nucleus.

Definition at line 80 of file G4INCLTransmissionChannel.hh.

Referenced by particleLeaves().

refraction

const G4bool G4INCL::TransmissionChannel::refraction

private

True if refraction should be applied.

Definition at line 77 of file G4INCLTransmissionChannel.hh.

Referenced by particleLeaves().

theNucleus

Nucleus* const G4INCL::TransmissionChannel::theNucleus

private

Definition at line 73 of file G4INCLTransmissionChannel.hh.

Referenced by fillFinalState(), and initializeKineticEnergyOutside().

theParticle

Particle* const G4INCL::TransmissionChannel::theParticle

private

Definition at line 74 of file G4INCLTransmissionChannel.hh.

Referenced by fillFinalState(), initializeKineticEnergyOutside(), and particleLeaves().

---

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

• source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLTransmissionChannel.hh
• source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLTransmissionChannel.cc