#include <G4INCLPauliGlobal.hh>

Inheritance diagram for G4INCL::PauliGlobal:

Public Member Functions
	PauliGlobal ()

	~PauliGlobal ()

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 45 of file G4INCLPauliGlobal.hh.

Constructor & Destructor Documentation

◆ PauliGlobal()

G4INCL::PauliGlobal::PauliGlobal ( )

Definition at line 44 of file G4INCLPauliGlobal.cc.

44{}

◆ ~PauliGlobal()

G4INCL::PauliGlobal::~PauliGlobal ( )

Definition at line 45 of file G4INCLPauliGlobal.cc.

45{}

Member Function Documentation

◆ isBlocked()

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

virtual

Implements G4INCL::IPauli.

Definition at line 47 of file G4INCLPauliGlobal.cc.

                                                                                    {
    for(ParticleIter p=pL.begin(); p!=pL.end(); ++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
      const ParticleList particles = n->getStore()->getParticles();
      G4int nSea = 0;
      for(ParticleIter i=particles.begin(); i!=particles.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;
 
  }