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: charges accessible with electron
  • 2: charges accessible with neutrino (= all)

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_n_s.dat.bz2'

where iii is the ID of the particle, n=1,2, and s=1,2

possible values:

NOTE if you want to use the files 'XXX.iii_n_s.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, where the possible charges were chosen according inMediumParameterset during tabulation

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"_1_?.dat.bz2 and
  • buuinput/inMediumWidth/InMediumWidth."particleID"_2_?.dat.bz2 .

OUTPUT

Initialized arrays widthTable, min_mass, max_mass, delta_mass

NOTES

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

baryonWidthMedium_tables/dumpTable [ Subroutines ]

[ Top ] [ baryonWidthMedium_tables ] [ Subroutines ]

NAME

subroutine dumpTable(iFile,particleID,iAbsP,iMass,iRhoN,iRhoP, iColumn)

PURPOSE

write table in special format for plotting

INPUTS

  • integer :: iFile -- output file
  • integer :: particleID -- ID of baryon
  • integer :: iAbsP -- bin of absP to use
  • integer :: iMass -- bin of mass to use
  • integer :: iRhoN -- bin of rhoN to use
  • integer :: iRhoP -- bin of rhoP to use
  • integer :: iColumn -- select the output, see below

OUTPUT

according the value of iColumn, the running variable is * 1: absP * 2: mass * 12,21: 2D-Array, absP, mass

NOTES

  • if iColumn=1, then the value iMass=-1 indicates, that the average over all masses is printed (this is useful for the nucleon only, where the width does not depend on the mass)