G4INCL::PauliGlobal Class Reference

#include <G4INCLPauliGlobal.hh>

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

Constructor & Destructor Documentation

G4INCL::PauliGlobal::PauliGlobal

(

)

Definition at line 42 of file G4INCLPauliGlobal.cc.

{}

G4INCL::PauliGlobal::~PauliGlobal

(

)

Definition at line 43 of file G4INCLPauliGlobal.cc.

{}

Member Function Documentation

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

virtual

Implements G4INCL::IPauli.

Definition at line 45 of file G4INCLPauliGlobal.cc.

References G4INCL::NuclearPotential::INuclearPotential::getFermiMomentum(), G4INCL::Nucleus::getInitialA(), G4INCL::Nucleus::getInitialZ(), G4INCL::Store::getParticles(), G4INCL::Nucleus::getPotential(), G4INCL::Nucleus::getStore(), G4INCL::Proton, and G4INCL::Random::shoot().

                                                                               {
    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
      const ParticleList 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;

  }

---

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

• G4INCLPauliGlobal.hh
• G4INCLPauliGlobal.cc