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

Constructor & Destructor Documentation

◆ CDPP()

G4INCL::CDPP::CDPP ( )

Definition at line 44 of file G4INCLCDPP.cc.

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

◆ ~CDPP()

G4INCL::CDPP::~CDPP ( )

Definition at line 50 of file G4INCLCDPP.cc.

50{}

Member Function Documentation

◆ isBlocked()

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(),
                  [this](Particle const * const p){ this->G4INCL::CDPP::processOneParticle(p); }
                  );
    std::for_each(created.begin(), created.end(),
                  [this](Particle const * const p){ this->G4INCL::CDPP::processOneParticle(p); }
                  );
 
    const G4double Tinitial = nucleus->getInitialInternalEnergy();
    const G4double Eblock = TbelowTf - Tinitial - Sk - S;
 
    return (Eblock < 0.0);
  }

◆ processOneParticle()

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

inline

Definition at line 54 of file G4INCLCDPP.hh.

                                                             {
      if(p->isBaryon()) {
        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->isMeson() || p->isResonance()) {
        const G4double sep = thePotential->getSeparationEnergy(p);
        Sk += sep;
      }
    }