#include <G4INCLCDPP.hh>

Inheritance diagram for G4INCL::CDPP:

Public Member Functions
	CDPP ()

	~CDPP ()

G4bool	isBlocked (ParticleList const &, Nucleus const *const)

void	processOneParticle (Particle const *const p)

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

virtual	~IPauli ()

Protected Attributes
G4double	Sk

G4double	TbelowTf

NuclearPotential::INuclearPotential const *	thePotential

Detailed Description

Definition at line 46 of file G4INCLCDPP.hh.

Constructor & Destructor Documentation

G4INCL::CDPP::CDPP ( )

Definition at line 43 of file G4INCLCDPP.cc.

              :
     Sk(0.0),
     TbelowTf(0.0),
     thePotential(NULL)
   {}

G4INCL::CDPP::~CDPP ( )

Definition at line 49 of file G4INCLCDPP.cc.

49 {}

Member Function Documentation

G4bool G4INCL::CDPP::isBlocked	(	ParticleList const &	created,
		Nucleus const * const	nucleus
	)

virtual

Implements G4INCL::IPauli.

Definition at line 51 of file G4INCLCDPP.cc.

References G4INCL::Nucleus::computeSeparationEnergyBalance(), G4INCL::Nucleus::getInitialInternalEnergy(), G4INCL::Store::getParticles(), G4INCL::Nucleus::getPotential(), G4INCL::Nucleus::getStore(), processOneParticle(), Sk, TbelowTf, and thePotential.

                                                                                    {
     G4double S = nucleus->computeSeparationEnergyBalance();
 
     thePotential = nucleus->getPotential();
 
     ParticleList const &remnantParticles = nucleus->getStore()->getParticles();
 
     Sk = 0.0;
     TbelowTf = 0.0;
 
     std::for_each(remnantParticles.begin(), remnantParticles.end(), std::bind1st(std::mem_fun(&G4INCL::CDPP::processOneParticle), this));
     std::for_each(created.begin(), created.end(), std::bind1st(std::mem_fun(&G4INCL::CDPP::processOneParticle), this));
 
     const G4double Tinitial = nucleus->getInitialInternalEnergy();
     const G4double Eblock = TbelowTf - Tinitial - Sk - S;
 
     return (Eblock < 0.0);
   }

void G4INCL::CDPP::processOneParticle ( Particle const *const p )

inline

Definition at line 53 of file G4INCLCDPP.hh.

References G4INCL::NuclearPotential::INuclearPotential::getFermiEnergy(), G4INCL::Particle::getKineticEnergy(), G4INCL::Particle::getPotentialEnergy(), G4INCL::NuclearPotential::INuclearPotential::getSeparationEnergy(), G4INCL::Particle::isNucleon(), G4INCL::Particle::isPion(), G4INCL::Particle::isResonance(), Sk, TbelowTf, and thePotential.

Referenced by isBlocked().

                                                              {
       if(p->isNucleon()) {
         const G4double Tf = thePotential->getFermiEnergy(p);
         const G4double T = p->getKineticEnergy();
 
         if(T > Tf) {
           const G4double sep = thePotential->getSeparationEnergy(p);
           Sk += sep;
         } else {
           TbelowTf += T - p->getPotentialEnergy();
         }
       } else if(p->isPion() || p->isResonance()) {
         const G4double sep = thePotential->getSeparationEnergy(p);
         Sk += sep;
       }
     }