gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/nucleus [ Modules ]

[ Top ] [ Modules ]

NAME

module nucleus

PURPOSE

Provides access to the projectile and target nuclei.

INPUTS

Namelists "target" and "projectile".

nucleus/getTarget [ Functions ]

[ Top ] [ nucleus ] [ Functions ]

NAME

function getTarget()

PURPOSE

Returns an initialized target nucleus resting at 0. with velocity=0.

INPUTS

  • NONE

OUTPUT

nucleus/getProjectile [ Functions ]

[ Top ] [ nucleus ] [ Functions ]

NAME

function getProjectile()

PURPOSE

Returns an initialized projectile nucleus resting at 0. with velocity=0.

INPUTS

  • NONE

OUTPUT

nucleus/initProjectile [ Subroutines ]

[ Top ] [ nucleus ] [ Subroutines ]

NAME

subroutine initProjectile

PURPOSE

Initializes the projectile nucleus resting in the origin according to Information in namelist "projectile" in jobcard.

initProjectile/A [ Global module-variables ]

[ Top ] [ initProjectile ] [ Global module-variables ]

SOURCE

    integer, save :: A = 0

PURPOSE

Mass A of projectile nucleus ( = number of nucleons). If zero, a default isotope is chosen for the given Z.

initProjectile/Z [ Global module-variables ]

[ Top ] [ initProjectile ] [ Global module-variables ]

SOURCE

    integer, save :: Z = 20

PURPOSE

Charge Z of projectile nucleus ( = number of protons).

initProjectile/Projectile_A [ Global module-variables ]

[ Top ] [ initProjectile ] [ Global module-variables ]

SOURCE

    integer, save :: Projectile_A = -99

PURPOSE

deprecated, use 'A' instead

initProjectile/Projectile_Z [ Global module-variables ]

[ Top ] [ initProjectile ] [ Global module-variables ]

SOURCE

    integer, save :: Projectile_Z = -99

PURPOSE

deprecated, use 'Z' instead

initProjectile/fermiMotion [ Global module-variables ]

[ Top ] [ initProjectile ] [ Global module-variables ]

SOURCE

    logical,save :: fermiMotion=.true.

PURPOSE

Determines whether particles in projectile nucleus feel Fermi motion.

initProjectile/densitySwitch_static [ Global module-variables ]

[ Top ] [ initProjectile ] [ Global module-variables ]

SOURCE

    integer,save :: densitySwitch_static=3

PURPOSE

This switch is important, because it decides, which static density is used to set up the testparticles in the nuclei before the first time-step.

Possible values:

  • 0 : density=0.0
  • 1 : Static density using Woods-Saxon according to H. Lenske
  • 2 : Static density according to Nieves93
  • 3 : Static density using Woods-Saxon according to H. Lenske, different radii for neutrons and protons
  • 4 : Static density according oscillator shell model
  • 5 : Density distribution is a sphere with density according to the input value of "fermiMomentum_input".
  • 6 : Static Density based on LDA, implemented by Birger Steinmueller
  • 7 : Static Density based on LDA + Welke potential
  • 8 : Static Density prescription according Relativistic Thomas-Fermi (valid only in RMF-mode)

Possible nuclei for the different prescriptions:

  • 1 : A > 2 (only A > 16 makes sense)
  • 2 : Be(9), C(12), O(16,18), Al(27), Ar(40), Ca(40,44), Fe(56), Cu(63), As(75), Ce(142), Sn(112,116,120,124), Ta(181), Au(197), Pb(208) see densityStatic.f90, subroutine denspar for more info
  • 3 : 6->C(12), 8->O(16),O(18), 13->Al(27), 20->Ca(40),Ca(44), 79->Au(197) 82->Pb(208)
  • 4: 2->He(4), 4->Be(9), 5->B(11), 6->C(12), 8->O(16)

Refs:

  • J.Nieves, E.Oset, C.Garcia-Recio, Nucl.Phys.A 554 (1993) 509

initProjectile/fermiMomentum_input [ Global module-variables ]

[ Top ] [ initProjectile ] [ Global module-variables ]

SOURCE

    real,save :: fermiMomentum_input=0.225

PURPOSE

Input value of the fermi momentum for densitySwitch_static=5 (in GeV).

initProjectile/anti [ Global module-variables ]

[ Top ] [ initProjectile ] [ Global module-variables ]

SOURCE

    logical,save :: anti=.FALSE.

PURPOSE

Indicate, whether it is a anti-nucleus

initProjectile/projectile [ Namelists ]

[ Top ] [ initProjectile ] [ Namelists ]

NAME

NAMELIST /projectile/

PURPOSE

Includes the input parameters for the projectile nucleus:

nucleus/DensTab_projectile.dat [ Output files ]

[ Top ] [ nucleus ] [ Output files ]

NAME

file DensTab_projectile.dat

PURPOSE

Density tabulation of the projectile nucleus at initialization.

nucleus/initTarget [ Subroutines ]

[ Top ] [ nucleus ] [ Subroutines ]

NAME

subroutine initTarget

PURPOSE

Initializes the target nucleus resting in the origin according to Information in namelist "target" in jobcard.

initTarget/A [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    integer, save :: A = 0

PURPOSE

Mass A of target nucleus ( = number of nucleons). If zero, a default isotope is chosen for the given Z.

initTarget/Z [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    integer, save :: Z = 20

PURPOSE

Charge Z of target nucleus ( = number of protons).

initTarget/Target_A [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    integer, save :: Target_A = -99

PURPOSE

deprecated, use 'A' instead

initTarget/Target_Z [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    integer, save :: Target_Z = -99

PURPOSE

deprecated, use 'Z' instead

initTarget/fermiMotion [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    logical,save :: fermiMotion=.true.

PURPOSE

Determines whether particles in target nucleus feel Fermi motion.

initTarget/densitySwitch_static [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    integer,save :: densitySwitch_static=3

PURPOSE

This switch is important, because it decides, which static density is used to set up the testparticles in the nuclei before the first time-step.

Possible values:

  • 0 : density=0.0
  • 1 : Static density using Woods-Saxon according to H. Lenske
  • 2 : Static density according to Nieves93
  • 3 : Static density using Woods-Saxon according to H. Lenske, different radii for neutrons and protons
  • 4 : Static density according oscillator shell model
  • 5 : Density distribution is a sphere with density according to the input value of "fermiMomentum_input".
  • 6 : Static density based on LDA, implemented by Birger Steinmueller
  • 7 : Static density based on LDA + Welke potential
  • 8 : Static density prescription according Relativistic Thomas-Fermi (valid only in RMF-mode)

Possible nuclei for the different prescriptions:

  • 1 : A > 2 (only A > 16 makes sense)
  • 2 : Be(9), C(12), O(16,18), Al(27), Ar(40), Ca(40,44), Fe(56), Cu(63), As(75), Ce(142), Sn(112,116,120,124), Ta(181), Au(197), Pb(208) see densityStatic.f90, subroutine denspar for more info
  • 3 : 6->C(12), 8->O(16),O(18), 13->Al(27), 20->Ca(40),Ca(44), 79->Au(197) 82->Pb(208)
  • 4: 2->He(4), 4->Be(9), 5->B(11), 6->C(12), 8->O(16)

Refs:

  • J.Nieves, E.Oset, C.Garcia-Recio, Nucl.Phys.A 554 (1993) 509

initTarget/fermiMomentum_input [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    real,save :: fermiMomentum_input=0.225

PURPOSE

Input value of the fermi momentum for densitySwitch_static=5 (in GeV).

initTarget/radius_input [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    real,dimension(1:2),save :: radius_input = -99.9

PURPOSE

Input value(s) of the radius of the nucleus in a Woods-Saxon parametrization.

If densitySwitch_static=...:

  • 1: one value is used both for protons and neutrons.
  • 3: two values: 1) proton, 2) neutron

Overrides default values only if given and positive.

initTarget/surface_input [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    real,dimension(1:2),save :: surface_input = -99.9

PURPOSE

Input value(s) of the surface of the nucleus in a Woods-Saxon parametrization

If densitySwitch_static=...:

  • 1: one value is used both for protons and neutrons.
  • 3: two values: 1) proton, 2) neutron

Overrides default values only if given and positive.

initTarget/ReAdjustForConstBinding [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    logical, save :: ReAdjustForConstBinding = .true.

PURPOSE

If this flag is set to true, we use the selected density distribution only for a preliminary step, where we calculate the baryonic potential as function of r (which depends on the density distribution). From the condition, that the binding energy has to be constant, we deduce the distribution of the fermi momentum and thus the 'new' density distribution.

The tabulated density distribution is replaced via the 'new' one and all behaviour is as usual.

Only valid in Skyrme mode (and also not for HeavyIon init).

initTarget/ConstBinding [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    real, save :: ConstBinding = -0.008

PURPOSE

if 'ReAdjustForConstBinding' equals true, we a trying to readjust the fermi momentum and the density such, we quarantee this value for the binding energy (in GeV).

Only valid in Skyrme mode (and also not for HeavyIon init)

initTarget/doPrintGlauber [ Global module-variables ]

[ Top ] [ initTarget ] [ Global module-variables ]

SOURCE

    logical, save :: doPrintGlauber = .false.

PURPOSE

switch on/off printing of simple Glauber calculations after init returning sigma_hA as function of sigma_hN (cf. Falter PhD, eq. (5.15))

initTarget/target [ Namelists ]

[ Top ] [ initTarget ] [ Namelists ]

NAME

NAMELIST /target/

PURPOSE

Includes the input parameters for the target nucleus:

nucleus/initNucleus [ Subroutines ]

[ Top ] [ nucleus ] [ Subroutines ]

NAME

subroutine initNucleus(Nuc, fermiMomentum_input)

PURPOSE

Initializes a nucleus resting in the frame of calculation at r=0.

INPUTS

NOTES

In 'Nuc', mass and charge have to be set as input variables.

we have to set nuc%radius, since this is referred to in further init routines (e.g. for calculating the distance)

nucleus/adjustTimesteps [ Subroutines ]

[ Top ] [ nucleus ] [ Subroutines ]

NAME

subroutine adjustTimesteps(Nuc)

PURPOSE

Adjust the overall running time of the simulation if the target is a helium target