gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/baryonWidthMedium_tables [ Modules ]

[ Top ] [ Modules ]

NAME

module baryonWidthMedium_tables

PURPOSE

Implements the routines for the medium width of the baryons

NOTES

The corresponding tables are created by inMediumWidth.f90

baryonWidthMedium_tables/debugFlag [ Global module-variables ]

[ Top ] [ baryonWidthMedium_tables ] [ Global module-variables ]

SOURCE

  logical, parameter :: debugFlag=.false.

PURPOSE

Switch for debug information

baryonWidthMedium_tables/maxRes [ Global module-variables ]

[ Top ] [ baryonWidthMedium_tables ] [ Global module-variables ]

SOURCE

  integer,save  :: maxRes=1000

PURPOSE

Read the data table up to a maximum resonance ID. ONLY FOR TESTING!!!

baryonWidthMedium_tables/minRes [ Global module-variables ]

[ Top ] [ baryonWidthMedium_tables ] [ Global module-variables ]

SOURCE

  integer,save  :: minRes=-1000

PURPOSE

Read the data table starting at this minimal resonance ID. ONLY FOR TESTING!!!

baryonWidthMedium_tables/deltaOset [ Global module-variables ]

[ Top ] [ baryonWidthMedium_tables ] [ Global module-variables ]

SOURCE

  logical,save  :: deltaOset=.false.

PURPOSE

Use delta width according to Oset et al. NPA 468 (1987)

baryonWidthMedium_tables/minimalWidth [ Global module-variables ]

[ Top ] [ baryonWidthMedium_tables ] [ Global module-variables ]

SOURCE

  real,parameter  :: minimalWidth=0.001

PURPOSE

The minimal width a particle can get. Must be finite for numerics.

baryonWidthMedium_tables/inMediumParameterset [ Global module-variables ]

[ Top ] [ baryonWidthMedium_tables ] [ Global module-variables ]

SOURCE

  integer,save  :: inMediumParameterset=2

PURPOSE

chooses the parameters for the inMediumWidth (1 electron, 2 neutrino)

baryonWidthMedium_tables/onlyNucleon [ Global module-variables ]

[ Top ] [ baryonWidthMedium_tables ] [ Global module-variables ]

SOURCE

  logical, save :: onlyNucleon=.false.

PURPOSE

Only for debugging: only nucleon width is read in.

baryonWidthMedium_tables/extrapolateAbsP [ Global module-variables ]

[ Top ] [ baryonWidthMedium_tables ] [ Global module-variables ]

SOURCE

  logical,save  :: extrapolateAbsP=.false.

PURPOSE

if(true) then set absP to maxAbsP if absP is larger

baryonWidthMedium_tables/FileNameMask [ Global module-variables ]

[ Top ] [ baryonWidthMedium_tables ] [ Global module-variables ]

PURPOSE

The absolute filename of the file containing the in-medium tables.

This is only the basename, the final name is given by 

   '[FileNameMask].iii_nnn.dat.bz2'

where iii is the ID of the particle and nnn=1 or 2

possible values:

NOTE if you want to use the files 'XXX.iii_nnn.dat.bz2' in the actual directory, give this mask as './XXX'

SOURCE

  character(1000), save :: FileNameMask = ''

baryonWidthMedium_tables/readInput [ Subroutines ]

[ Top ] [ baryonWidthMedium_tables ] [ Subroutines ]

NAME

subroutine readInput

PURPOSE

Reads input in jobcard from namelist "BaryonWidthMedium_tables".

baryonWidthMedium_tables/BaryonWidthMedium_tables [ Namelists ]

[ Top ] [ baryonWidthMedium_tables ] [ Namelists ]

NAME

NAMELIST /BaryonWidthMedium_tables/

PURPOSE

Includes the input switches:

baryonWidthMedium_tables/get_inMediumWidth [ Functions ]

[ Top ] [ baryonWidthMedium_tables ] [ Functions ]

NAME

function get_inMediumWidth(particleID,momentum,baremass,rhoN,rhoP,switch,outOfBounds) result(width)

PURPOSE

  • Returns the in-medium width of baryons according to Gamma=Gamma_free*Pauli_Blocking+sigma*rho*v *Pauli_Blocking
  • An average over the Fermi see is performed
  • Sigma is given by the actual full collision term, with the incoming particles in the vacuum
  • Simplification: An average over the charge of the baryon is performed.

INPUTS

  • integer, intent(in) :: particleID -- ID of baryon
  • real, dimension (0:3), intent(in) :: Momentum -- Momentum (in LRF)
  • real, intent(in) :: baremass -- bare mass of baryon
  • real, intent(in) :: rhoN,rhoP -- proton and neutron density in fm^-3
  • integer, intent(in) :: switch -- see OUTPUT for consequences of this switch
  • logical, optional, intent(in) :: doMassCorrect_in -- only for debugging: switch on/off mass correction

OUTPUT

  • real :: width
  • logical, optional :: outOfBounds -- .true. if the input variables are out of bounds

If switch equals...:

  • 1=only collisional width
  • 2=only pauli blocked free width
  • else= sum of collisional width and pauli blocked free width

get_inMediumWidth/readTable [ Subroutines ]

[ Top ] [ get_inMediumWidth ] [ Subroutines ]

NAME

subroutine readTable

PURPOSE

Reads tabulated arrays from files:

  • buuinput/inMediumWidth/InMediumWidth."particleID"_001.dat.bz2 and
  • buuinput/inMediumWidth/InMediumWidth."particleID"_002.dat.bz2 .

OUTPUT

Initialized arrays widthTable, min_mass, max_mass, delta_mass

NOTES

The file name may be given via the jobcard, see fileNameMask