gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/baryonWidthMedium [ Modules ]

[ Top ] [ Modules ]

NAME

module baryonWidthMedium

PURPOSE

Basically this module calls the routines of baryonWidth and adds medium modifications to it.

USES

module baryonWidth

baryonWidthMedium/mediumSwitch [ Global module-variables ]

[ Top ] [ baryonWidthMedium ] [ Global module-variables ]

SOURCE

  logical, save :: mediumSwitch=.false.

PURPOSE

Switch on and off the in-medium width of all baryons at once. If .false., the vacuum width are used.

baryonWidthMedium/mediumSwitch_Delta [ Global module-variables ]

[ Top ] [ baryonWidthMedium ] [ Global module-variables ]

SOURCE

  logical, save :: mediumSwitch_Delta=.false.

PURPOSE

Only meaningful if mediumSwitch=.true.: Switch on and off the in-medium width of the Delta. (.false.=off)

Note that in that case the Delta is treated specially: what is used for the in-medium width is determined by the flag in deltaWidth. This switch is not consistent with mediumSwitch_coll!

baryonWidthMedium/mediumSwitch_coll [ Global module-variables ]

[ Top ] [ baryonWidthMedium ] [ Global module-variables ]

SOURCE

  logical, save :: mediumSwitch_coll=.false.

PURPOSE

Only meaningful if mediumSwitch=.true.: Use in-medium width according to collision term.

NOTES

if set to TRUE, then also UseOffShellPotentialBaryons (see module offShellPotential) must be .true.

baryonWidthMedium/mediumSwitch_proton_neutron [ Global module-variables ]

[ Top ] [ baryonWidthMedium ] [ Global module-variables ]

SOURCE

  logical, save :: mediumSwitch_proton_neutron=.false.

PURPOSE

Only meaningful if mediumSwitch=.true.: Switch on and off the in-medium width of the proton and the neutron. (.false.=off)

Note that in that case the nucleons are treated specially. This switch is not consistent with mediumSwitch_coll!

baryonWidthMedium/verboseInit [ Global module-variables ]

[ Top ] [ baryonWidthMedium ] [ Global module-variables ]

SOURCE

  logical, save :: verboseInit = .false.

PURPOSE

switch on/off informational messages during initialization

baryonWidthMedium/verboseInitStop [ Global module-variables ]

[ Top ] [ baryonWidthMedium ] [ Global module-variables ]

SOURCE

  logical, save :: verboseInitStop = .false.

PURPOSE

Stop after informational messages during initialization or not.

baryonWidthMedium/get_MediumSwitch [ Functions ]

[ Top ] [ baryonWidthMedium ] [ Functions ]

NAME

logical function get_MediumSwitch()

PURPOSE

Returns the value of mediumSwitch

baryonWidthMedium/get_MediumSwitch_coll [ Functions ]

[ Top ] [ baryonWidthMedium ] [ Functions ]

NAME

logical function get_MediumSwitch_coll()

PURPOSE

Returns the value of mediumSwitch_coll

baryonWidthMedium/get_MediumSwitch_Delta [ Functions ]

[ Top ] [ baryonWidthMedium ] [ Functions ]

NAME

logical function get_MediumSwitch_Delta()

PURPOSE

Returns the value of mediumSwitch_Delta

baryonWidthMedium/get_MediumSwitch_proton_neutron [ Functions ]

[ Top ] [ baryonWidthMedium ] [ Functions ]

NAME

logical function get_MediumSwitch_proton_neutron()

PURPOSE

Returns the value of mediumSwitch_proton_neutron

baryonWidthMedium/readInput [ Subroutines ]

[ Top ] [ baryonWidthMedium ] [ Subroutines ]

NAME

subroutine readInput

PURPOSE

Reads input in jobcard out of namelist "width_baryon".

baryonWidthMedium/width_Baryon [ Namelists ]

[ Top ] [ baryonWidthMedium ] [ Namelists ]

NAME

NAMELIST width_Baryon

PURPOSE

Includes the switches:

baryonWidthMedium/partialWidthBaryonMedium [ Functions ]

[ Top ] [ baryonWidthMedium ] [ Functions ]

NAME

real function partialWidthBaryonMedium(particleID,mass,inWidth,mesonID,baryonID,momentumLRF,mediumATposition,baryonMass,mesonMass)

PURPOSE

This function calculates the partial width with energy dependence and medium modifications. If the medium modification is not yet implemented for a specific case, partialWidthBaryon (vacuum) is called!

INPUTS

  • integer :: particleID -- id of resonance
  • real :: mass -- bare mass of the resonance (offshell), not including the potentials!
  • logical :: inWidth -- .true.=> in-width (only important for channels with unstable particles); .false.=> out-width
  • integer :: mesonID, baryonID -- ID's of the decay products which one is interested in
  • type(medium) :: mediumATposition -- Medium information
  • real, dimension(0:3) :: momentumLRF -- Momentum im LRF of the resonance
  • real, OPTIONAL :: baryonMass,mesonMass -- Possibility to define the masses of the incoming baryon and meson, needed in the case of the In-Width if one of them is off-shell. Otherwise not relevant.

baryonWidthMedium/decayWidthBaryonMedium [ Subroutines ]

[ Top ] [ baryonWidthMedium ] [ Subroutines ]

NAME

subroutine decayWidthBaryonMedium(particleID,mass,momentumLRF,mediumATposition, decayWidth, pauliFlag)

PURPOSE

This function returns the partial out-widths with energy dependence and medium modifications. If the out-width has already Pauli blocking effects included, then pauliFlag is set to .true. .

INPUTS

  • integer :: particleID -- id of resonance
  • real :: mass -- bare mass of the resonance (offshell), not including potentials!
  • real, dimension(1:3) :: position -- position vector of the resonance

OUTPUT

  • real, dimension(:) :: decayWidth -- array of decay widths for all decay channels

baryonWidthMedium/WidthBaryonMedium [ Functions ]

[ Top ] [ baryonWidthMedium ] [ Functions ]

NAME

real function WidthBaryonMedium(particleID,mass,momentumLRF,mediumATposition,outOfBounds)

PURPOSE

This function calculates the full width with energy dependence and medium modifications. If for a specific resonance, medium modifications are not yet implemented, then WidthBaryon is called and the vacuum width is used.

INPUTS

  • integer :: particleID -- id of resonance
  • real :: mass -- bare mass of the resonance (offshell), not including potentials
  • type(medium) :: mediumATposition -- Medium at decay vertex
  • real,dimension(0:3):: momentumLRF -- momentum of resonance in LRF

OUTPUT

  • logical, OPTIONAL :: outOfBounds

NOTES

For the Delta, the in-medium width 'delta_fullWidth' does not converge to the vacuum values 'FullWidthBaryon(2)' for large M in the limit rho->0 ! This is due to upper bound of the tabulation in both cases (for the inMedium cases, the width is kept constant for M>2GeV, while inVacuum the upper mass is larger.)

In order to get a smooth behaviour, we do not respect "medium%useMedium" anymore, if "mediumSwitch" is set: Also very tiny densities are treated as 'in-medium' and not as vacuum.

baryonWidthMedium/GetMassAssInfo_Baryon [ Subroutines ]

[ Top ] [ baryonWidthMedium ] [ Subroutines ]

NAME

subroutine GetMassAssInfo_Baryon(MassAssInfo, ID,momLRF,mediumAtPos)

PURPOSE

Set the information necessary for improved MassAss

INPUTS

  • integer :: ID -- ID of particle
  • real, dimension(0:3) :: momLRF -- momentum of resonance in LRF
  • type(medium) :: mediumAtPos -- medium at position

OUTPUT

NOTES

  • The type 'MassAssInfo' is given as 'intent(inout)' since we have to bother with the allocation of its components

baryonWidthMedium/InitMassAssInfo_Baryon [ Subroutines ]

[ Top ] [ baryonWidthMedium ] [ Subroutines ]

NAME

subroutine InitMassAssInfo_Baryon

PURPOSE

This routine initializes the array of information used by the improved massass routines.