gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


/pauliBlockingModule [ Modules ]

[ Top ] [ Modules ]

NAME

module PauliBlockingModule

PURPOSE

Contains all information and routines which are necessary for the Pauli blocking of the test-particles.


pauliBlockingModule/pauliSwitch [ Global module-variables ]

[ Top ] [ pauliBlockingModule ] [ Global module-variables ]

SOURCE

  integer, save :: pauliSwitch=1

PURPOSE

  • 0 : No Pauli blocking
  • 1 : dynamic Pauli blocking (use actual phase space densities)
  • 2 : analytic Pauli blocking (use ground state assumption) (not possible for Heavy Ions!)
  • 3 : dynamic Pauli blocking in a box


pauliBlockingModule/densDepMomCutFlag [ Global module-variables ]

[ Top ] [ pauliBlockingModule ] [ Global module-variables ]

SOURCE

  logical, save :: densDepMomCutFlag=.false.

PURPOSE

if .true. - the radius in momentum space for selecting nucleons around given nucleon will depend on local Fermi momentum

NOTES

Used only for dynamic pauli blocking.


pauliBlockingModule/Gauss [ Global module-variables ]

[ Top ] [ pauliBlockingModule ] [ Global module-variables ]

SOURCE

  real, save :: Gauss=1.0

PURPOSE

Smearing for dynamic pauli blocking


pauliBlockingModule/cutGauss [ Global module-variables ]

[ Top ] [ pauliBlockingModule ] [ Global module-variables ]

SOURCE

  real, save :: cutGauss=2.2

PURPOSE

Cutoff for gauss Smearing


pauliBlockingModule/cutMom [ Global module-variables ]

[ Top ] [ pauliBlockingModule ] [ Global module-variables ]

SOURCE

  real, save :: cutMom=0.08

PURPOSE


pauliBlockingModule/cutPos [ Global module-variables ]

[ Top ] [ pauliBlockingModule ] [ Global module-variables ]

SOURCE

  real, save :: cutPos=1.86

PURPOSE

Radius of phase space box in position space


pauliBlockingModule/nGridPos [ Global module-variables ]

[ Top ] [ pauliBlockingModule ] [ Global module-variables ]

SOURCE

  integer, save :: nGridPos=30

PURPOSE

number of points in position space to save weights on


pauliBlockingModule/ensembleJump [ Global module-variables ]

[ Top ] [ pauliBlockingModule ] [ Global module-variables ]

SOURCE

  integer, save :: ensembleJump=5

PURPOSE

Parameter for speedup. Only every "ensemblejump"th ensemble is considered to evaluate the probability for pauli blocking.


pauliBlockingModule/DoHistogram [ Global module-variables ]

[ Top ] [ pauliBlockingModule ] [ Global module-variables ]

SOURCE

  logical, save :: DoHistogram=.false.

PURPOSE

if .true., a historgram is filled representing the blocking probability as function of the fermi momentum. You have to call 'WriteBlockMom' explicitely for writing the histogram


pauliBlockingModule/pauliBlocking [ Functions ]

[ Top ] [ pauliBlockingModule ] [ Functions ]

NAME

logical function pauliBlocking(momentum,position,Nukcharge,Teilchen,probabilityOut,weight)

PURPOSE

  • evaluates Pauli blocking for neutrons or protons
  • returns .true. if position and momentum are blocked by Pauli principle
  • returns .false. if no Pauli blocking

INPUTS

  • real, dimension(0:3) :: momentum
  • real, dimension(1:3) :: position
  • integer :: Nukcharge -- charge of nucleon
  • type(particle), dimension(:,:), optional :: Teilchen -- full particle vector
  • real, optional :: weight -- particle weight (for documentation purposes)

OUTPUT

  • real, optional :: probabilityOut -- probability of blocking, i.e. occupation number


pauliBlockingModule/WriteBlockMom [ Subroutines ]

[ Top ] [ pauliBlockingModule ] [ Subroutines ]

NAME

subroutine WriteBlockMom(mul)

PURPOSE

Write Histogram 'PauliBlocking.BlockMom.dat'


pauliBlockingModule/init [ Subroutines ]

[ Top ] [ pauliBlockingModule ] [ Subroutines ]

NAME

subroutine init

PURPOSE

Reads data out of namelist 'initPauli'.


pauliBlockingModule/initPauli [ Namelists ]

[ Top ] [ pauliBlockingModule ] [ Namelists ]

NAME

NAMELIST initPauli

PURPOSE

Includes the input switches and variables:


pauliBlockingModule/checkPauli [ Functions ]

[ Top ] [ pauliBlockingModule ] [ Functions ]

NAME

logical function checkPauli(teilchen,realParticles)

PURPOSE

Administrates the Pauli-blocking decision.

INPUTS

OUTPUT

  • .true. = momentum of particle "teilchen" IS NOT pauli blocked, .false.= momentum of particle "teilchen" IS pauli blocked

NOTE The position and momentum of the outgoing particles must be set!