gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


/deuterium [ Modules ]

[ Top ] [ Modules ]

NAME

module deuterium

PURPOSE

Includes subroutines to model the momentum distribution of nucleons in a Deuterium.


deuterium/waveFunction_switch [ Global module-variables ]

[ Top ] [ deuterium ] [ Global module-variables ]

SOURCE

  integer, save :: waveFunction_switch = 1

PURPOSE

Possible values are:

  • 0 -- No Wave functions! Pointlike Deuterium
  • 1 -- Wave functions according to Bonn potential
  • 2 -- Wave functions according to Argonne V18


deuterium/iParam [ Global module-variables ]

[ Top ] [ deuterium ] [ Global module-variables ]

SOURCE

  integer, save :: iParam = 1

PURPOSE

Choose parameterization of momentum distribution when using the Bonn potential. Possible values:

  • 1 -- Full Bonn (MaH87)
  • 2 -- OBEPQ (MaH87)
  • 3 -- OBEPQ-A (Mac89)
  • 4 -- OBEPQ-B (Mac89)
  • 5 -- OBEPQ-! (Mac89)
  • 6 -- OBEPR (MaH87) self-made
  • 7 -- Paris

References: MaH87: R. Machleidt et al. Phys. Rep. 149, 1 (1987) Mac89: R. Machleidt, Advances in Nucl. Phys. Vol 19


deuterium/pMax [ Global module-variables ]

[ Top ] [ deuterium ] [ Global module-variables ]

SOURCE

  real,save :: pMax = 0.5

PURPOSE

Cut-off parameter for Fermi momentum


deuterium/scaleMomentum [ Global module-variables ]

[ Top ] [ deuterium ] [ Global module-variables ]

SOURCE

  real, save :: scaleMomentum = 1.0

PURPOSE

The selected momentum is multiplied by this factor afterwards, i.e. some rescaling is done


deuterium/deuteriumFermi [ Namelists ]

[ Top ] [ deuterium ] [ Namelists ]

NAME

Namelist /deuteriumFermi/

PURPOSE

Includes the parameters:

which are used for the distribution of nucleons in deuterium.


deuterium/initDeuterium [ Subroutines ]

[ Top ] [ deuterium ] [ Subroutines ]

NAME

subroutine initDeuterium(teilchen,nuc,eventNummer,keepNumber)

PURPOSE

Represents nucleus 'nuc' in phase space by testparticles which are stored in vector 'teilchen'. The ordering in the vector teilchen is choosen to be random. The nucleus must be Deuterium.

NOTES

  • momentum distribution implemented by Thomas Falter, spatial distribution added later
  • The first nucleon is initialized at the origin, while the second nucleon gets its spatial coordinate according the squared wave function.
  • The averaged deuterium radius is <r_d>~=2.0 fm. The variable r used here is the distance between the two nucleons. Therefore: <r_d> = 1/2 * sqrt(Int_0^Infty dr r^2 psi^2)

INPUTS

  • type(nucleus) :: nuc
  • type(particle),dimension(:,:) :: teilchen
  • integer :: eventNummer -- eventNummer is given as "%event" to any initialized nucleon.
  • logical :: keepNumber -- flag whether each testparticle gets unique number or not

OUTPUT

  • type(particle),dimension(:,:) :: teilchen