G4INCL::CDPP Class Reference

#include <G4INCLCDPP.hh>

Inheritance diagram for G4INCL::CDPP:

Collaboration 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 47 of file G4INCLCDPP.hh.

Constructor & Destructor Documentation

G4INCL::CDPP::CDPP

(

)

Definition at line 44 of file G4INCLCDPP.cc.

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

G4INCL::CDPP::~CDPP

(

)

Definition at line 50 of file G4INCLCDPP.cc.

{}

Member Function Documentation

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

virtual

Implements G4INCL::IPauli.

Definition at line 52 of file G4INCLCDPP.cc.

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

Here is the call graph for this function:

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

inline

Definition at line 54 of file G4INCLCDPP.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

G4double G4INCL::CDPP::Sk

protected

Definition at line 72 of file G4INCLCDPP.hh.

G4double G4INCL::CDPP::TbelowTf

protected

Definition at line 73 of file G4INCLCDPP.hh.

NuclearPotential::INuclearPotential const* G4INCL::CDPP::thePotential

protected

Definition at line 74 of file G4INCLCDPP.hh.

---

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

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