#include <G4INCLPauliGlobal.hh>

Inheritance diagram for G4INCL::PauliGlobal:

Collaboration diagram for G4INCL::PauliGlobal:

Public Member Functions
	PauliGlobal ()

	~PauliGlobal ()

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

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

virtual	~IPauli ()

Detailed Description

Definition at line 44 of file G4INCLPauliGlobal.hh.

Constructor & Destructor Documentation

◆ PauliGlobal()

G4INCL::PauliGlobal::PauliGlobal ( )

Definition at line 43 of file G4INCLPauliGlobal.cc.

43 {}

◆ ~PauliGlobal()

G4INCL::PauliGlobal::~PauliGlobal ( )

Definition at line 44 of file G4INCLPauliGlobal.cc.

44 {}

Member Function Documentation

◆ isBlocked()

G4bool G4INCL::PauliGlobal::isBlocked	(	ParticleList const &	pL,
		Nucleus const * const	n
	)

virtual

Implements G4INCL::IPauli.

Definition at line 46 of file G4INCLPauliGlobal.cc.

                                                                                {
     for(ParticleIter p=pL.begin(), e=pL.end(); p!=e; ++p) {
       // Pauli blocking only applies to nucleons
       if(!(*p)->isNucleon()) continue;
 
       // If the particle is above T_F, it is never blocked
       const ParticleType t = (*p)->getType();
       const G4double pFermi = n->getPotential()->getFermiMomentum(t);
       const G4double pFermiSquared = pFermi*pFermi;
       if((*p)->getMomentum().mag2() > pFermiSquared) continue;
 
       // Count particles of the same type as p below the Fermi sea
       ParticleList const &particles = n->getStore()->getParticles();
       G4int nSea = 0;
       for(ParticleIter i=particles.begin(), end=particles.end(); i!=end; ++i) {
         if((*i)->getType() != t) continue;
         const G4double pmod2 = (*i)->getMomentum().mag2();
         if(pmod2<pFermiSquared) nSea++;
       }
 
       // Compute the blocking probability
       G4double probBlocking;
       if(t==Proton)
         probBlocking = ((G4double) nSea)/((G4double) n->getInitialZ());
       else
         probBlocking = ((G4double) nSea)/((G4double) (n->getInitialA() - n->getInitialZ()));
 
       // The avatar is blocked if any particle is blocked
       if(Random::shoot() < probBlocking) return true;
 
     }
 
     // Not blocked
     return false;
 
   }

Here is the call graph for this function: