#include <G4INCLCDPP.hh>

Inheritance diagram for G4INCL::CDPP:

Public Member Functions
	CDPP ()

	~CDPP ()

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

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

virtual	~IPauli ()

virtual G4bool	isBlocked (ParticleList const, Nucleus const *const) const =0

Detailed Description

Definition at line 47 of file G4INCLCDPP.hh.

Constructor & Destructor Documentation

◆ CDPP()

G4INCL::CDPP::CDPP ( )

Definition at line 43 of file G4INCLCDPP.cc.

43{}

◆ ~CDPP()

G4INCL::CDPP::~CDPP ( )

Definition at line 45 of file G4INCLCDPP.cc.

45{}

Member Function Documentation

◆ isBlocked()

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

virtual

Implements G4INCL::IPauli.

Definition at line 47 of file G4INCLCDPP.cc.

                                                                                        {
    G4double S = nucleus->computeSeparationEnergyBalance();
 
    const G4double Sp = nucleus->getPotential()->getSeparationEnergy(Proton);
    const G4double Sn = nucleus->getPotential()->getSeparationEnergy(Neutron);
 
    ParticleList remnantParticles = nucleus->getStore()->getParticles();
    remnantParticles.insert(remnantParticles.end(), created.begin(), created.end());
 
    G4double Sk = 0.0;
    G4double TbelowTf = 0.0;
    for(ParticleIter i = remnantParticles.begin(); i != remnantParticles.end();
        ++i) {
 
      if((*i)->isNucleon()) {
        const G4double Tf = nucleus->getPotential()->getFermiEnergy(*i);
        const G4double T = (*i)->getKineticEnergy();
 
        if(T > Tf) {
          const G4double sep = nucleus->getPotential()->getSeparationEnergy(*i);
          Sk += sep;
        } else {
          TbelowTf += T - (*i)->getPotentialEnergy();
        }
      } else if((*i)->isResonance()) {
        const G4double Tf = nucleus->getPotential()->getFermiEnergy(*i);
        const G4double T = (*i)->getKineticEnergy();
 
        if(T > Tf) {
          const G4double sep = nucleus->getPotential()->getSeparationEnergy(*i);
          Sk += sep;
        } else { // Ugly! We should use total energies everywhere!
          TbelowTf += (*i)->getEnergy() - ParticleTable::getINCLMass(Proton) - (*i)->getPotentialEnergy();
        }
      } else if((*i)->getType() == PiPlus)
        Sk += Sp - Sn;
      else if((*i)->getType() == PiMinus)
        Sk += Sn - Sp;
 
    }
    G4double Tinitial = nucleus->getInitialInternalEnergy();
    G4double Eblock = TbelowTf - Tinitial - Sk - S;
 
    return (Eblock < 0.0);
  }