gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/initHadron [ Modules ]

[ Top ] [ Modules ]

NAME

module initHadron

PURPOSE

Initialises the hadron projectile for a hadron-nucleus collision. The eventtype 'hadron' always uses real particles.

NOTES

The target nucleus has to be initialised separately (by using initNucPhaseSpace routine) before calling the initHadronInduced routine.

initHadron/impactParameter [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  real, save :: impactParameter=0.

PURPOSE

smaller 0: Impact parameter will be chosen randomly in the interval [0;abs(impactParameter)] (see subroutine setGeometry). It is recommended to take very large negative value of impactParameter in order to have good automatic random choice, e.g. impactParameter=-100.

initHadron/bRaiseFlag [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  logical, save :: bRaiseFlag=.false.

PURPOSE

if .true.: actual impact parameter will be raised by deltaB after nRunsPerB subsequent runs. Starting value is given by the impactParameter variable.

initHadron/deltaB [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  real, save :: deltaB=0.

PURPOSE

impact parameter step (relevant if bRaiseFlag=.true.)

initHadron/nRunsPerB [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  integer, save :: nRunsPerB=1

PURPOSE

number of subsequent runs per impact parameter (relevant if bRaiseFlag=.true.)

initHadron/perturbative [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  logical, save :: perturbative=.false.

PURPOSE

if .true. then the hadron is a perturbative particle

initHadron/numberParticles [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  integer,save   :: numberParticles = 200

PURPOSE

Number of projectile testparticles per ensemble in the case of a perturbative treatment

initHadron/particleId [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  integer, save :: particleId=1

PURPOSE

Identity of the projectile hadron.

initHadron/antiParticle [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  logical, save :: antiParticle=.false.

PURPOSE

if .true. then the hadron is an antiparticle

initHadron/particleCharge [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  integer, save :: particleCharge=0

PURPOSE

Charge of the hadron

initHadron/ekin_lab [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  real,    save :: ekin_lab=1.

PURPOSE

kinetic energy of the hadron in the rest frame of the target nucleus (GeV)

NOTES

If ekin_lab < 0. --- initialization according to the binding energy

initHadron/E_bind [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  real,    save :: E_bind=0.

PURPOSE

binding energy of initialized hadron (GeV)

NOTES

Active for iniType= 0,2 if ekin_lab < 0. is set.

initHadron/iniType [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  integer, save :: iniType=0

PURPOSE

initHadron/zChoice [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  integer, save :: zChoice=1

PURPOSE

Relevant for iniType=0 or iniType=2.

initHadron/delta [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  real, save :: delta=0.5

PURPOSE

Relevant for iniType=0 or iniType=2.

initHadron/deltaZ [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  real, save :: deltaZ=5.

PURPOSE

z = -deltaZ - R_nucleus, where z is z-coordinate of the hadron Relevant for iniType=0,2 and zChoice=2.

initHadron/width [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  real, save ::    width=1.

PURPOSE

Width of a gaussian density profile [fm]. Only relevant for iniType= 1 and 2.

initHadron/debug [ Global module-variables ]

[ Top ] [ initHadron ] [ Global module-variables ]

SOURCE

  logical, parameter :: debug=.false.

PURPOSE

if .true. then additional printouts are done for debugging

initHadron/initHadronInduced [ Subroutines ]

[ Top ] [ initHadron ] [ Subroutines ]

NAME

subroutine initHadronInduced (teilchenReal, teilchenPert)

PURPOSE

Provides initial conditions for a hadron.

INPUTS

  • type(particle),dimension(:,:),intent(inout) :: teilchenReal --- array of real particles to store the hadron if it is real
  • type(particle),dimension(:,:),intent(inout) :: teilchenPert --- array of perturbative particles to store the hadron if it is perturbative

NOTES

  • The user has to provide the namelist 'hadron', which contains the impact parameter, lab. energy, id, charge and antiparticle- parameters of the hadron. If the input impact parameter is negative, then the actual impact parameter is chosen by Monte-Carlo between 0 and abs(input impact parameter).
  • The target nucleus has to be initialised before calling this routine.
  • The program initialises the real (not perturbative) particle.

initHadronInduced/initInput [ Subroutines ]

[ Top ] [ initHadronInduced ] [ Subroutines ]

NAME

subroutine initInput

PURPOSE

Reads input out of jobcard. Namelist 'hadron'.

initHadron/hadron [ Namelists ]

[ Top ] [ initHadron ] [ Namelists ]

NAME

NAMELIST hadron

PURPOSE

Includes parameters for initialization of a hadron in the case of the hadron-nucleus collision:

initHadronInduced/setKinematics [ Subroutines ]

[ Top ] [ initHadronInduced ] [ Subroutines ]

NAME

subroutine setKinematics

PURPOSE

Sets basic kinematics of a hadron colliding with a nucleus.

initHadronInduced/setPosition [ Subroutines ]

[ Top ] [ initHadronInduced ] [ Subroutines ]

NAME

subroutine setPosition

PURPOSE

Sets position of the hadron.

NOTES

In the case of usual initialisation (initype=0), the actual impact parameter b and the coordinate z has to be defined before by calling the subroutine setGeometry.

initHadronInduced/setGeometry [ Subroutines ]

[ Top ] [ initHadronInduced ] [ Subroutines ]

NAME

subroutine setGeometry (flagPrint)

PURPOSE

Sets the impact parameter b and coordinate z of the hadron.

INPUTS

  • logical, intent(in) :: flagPrint -- if .true., then print output parameters to standard output file

NOTES

If the input impactParameter is less than zero, than the actual impact parameter is choosen by a Monte-Carlo-decision.