G4INCL::PauliGlobal Class Reference

#include <G4INCLPauliGlobal.hh>

Inheritance diagram for G4INCL::PauliGlobal:

Public Member Functions
G4bool	isBlocked (ParticleList const &, Nucleus const *const)
	PauliGlobal ()
	~PauliGlobal ()

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.

{}

~PauliGlobal()

G4INCL::PauliGlobal::~PauliGlobal

(

)

Definition at line 44 of file G4INCLPauliGlobal.cc.

{}

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

References CLHEP::detail::n, G4INCL::Proton, and G4INCL::Random::shoot().

---

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

• source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLPauliGlobal.hh
• source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLPauliGlobal.cc