gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/inMediumWidth [ Modules ]

[ Top ] [ Modules ]

NAME

module inMediumWidth

PURPOSE

  • Implements the routines for the medium width of the baryons
  • Generates tabulated input files
  • see also code/inMediumWidth/testRun/tabulateImagPart.f90 which creates an executable using this routine to tabulate the width

inMediumWidth/debugFlag [ Global module-variables ]

[ Top ] [ inMediumWidth ] [ Global module-variables ]

SOURCE

  logical,save  :: debugFlag=.false.

PURPOSE

  • Switch for debug information

inMediumWidth/writeLocal [ Global module-variables ]

[ Top ] [ inMediumWidth ] [ Global module-variables ]

SOURCE

  logical,save  :: writeLocal=.true.

PURPOSE

  • If .true. then output files are written to ./"filename", else to buuinput.

inMediumWidth/maxRes [ Global module-variables ]

[ Top ] [ inMediumWidth ] [ Global module-variables ]

SOURCE

  integer,save  :: maxRes=1000

PURPOSE

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

inMediumWidth/minRes [ Global module-variables ]

[ Top ] [ inMediumWidth ] [ Global module-variables ]

SOURCE

  integer,save  :: minRes=-1000

PURPOSE

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

inMediumWidth/num_MonteCarlo_Points [ Global module-variables ]

[ Top ] [ inMediumWidth ] [ Global module-variables ]

SOURCE

  integer, save :: num_MonteCarlo_Points = 2000

PURPOSE

number of monte-carlo points for the sampling of the phase space

usual values: 2000 for baryons, 250 for mesons

inMediumWidth/readInput [ Subroutines ]

[ Top ] [ inMediumWidth ] [ Subroutines ]

NAME

subroutine readInput

PURPOSE

Reads input in jobcard out of namelist "inMediumWidth".

inMediumWidth/tabulate_inMediumWidth_baryons [ Subroutines ]

[ Top ] [ inMediumWidth ] [ Subroutines ]

NAME

subroutine tabulate_inMediumWidth_baryons()

PURPOSE

  • Tabulates the in-medium width of baryons according to Gamma=Gamma_free*Pauli_Blocking+sigma*rho*v *Pauli_Blocking
  • An average over the Fermi sea 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

  • NONE

OUTPUT

inMediumWidth/tabulate_inMediumWidth_mesons [ Subroutines ]

[ Top ] [ inMediumWidth ] [ Subroutines ]

NAME

subroutine tabulate_inMediumWidth_mesons()

PURPOSE

  • Tabulates the in-medium width of mesons according to Gamma=Gamma_free+sigma*rho*v
  • An average over the Fermi sea 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 meson is performed.

INPUTS

  • NONE

OUTPUT

inMediumWidth/calcPauliBlockedDecayWidth [ Functions ]

[ Top ] [ inMediumWidth ] [ Functions ]

NAME

function calcPauliBlockedDecayWidth(particleID,momentumLRF,mass,rhoN,rhoP) result(gammaDecay)

PURPOSE

  • Returns the in-medium width of baryons according to Gamma=Gamma_free*Pauli_Blocking
  • An average over the Fermi sea is performed
  • Simplification: An average over the charge of the baryon is performed.

INPUTS

  • integer, intent(in) :: particleID -- ID of baryon
  • real, intent(in) :: absP -- absolute Momentum
  • real, intent(in) :: mass -- Mass of baryon
  • real, intent (in) :: rhoN,rhoP -- proton, neutron density in fm^-3

OUTPUT

  • real :: gammaDecay

inMediumWidth/calcCollBroadening_Bar [ Functions ]

[ Top ] [ inMediumWidth ] [ Functions ]

NAME

function calcCollBroadening_Bar(particleID,momentumLRF,mass,rhoN,rhoP) RESULT(gcoll)

PURPOSE

  • Returns the in-medium width of baryons according to Gamma=sigma*rho*v *Pauli_Blocking
  • An average over the Fermi sea 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, intent(in) :: absP -- absolute Momentum
  • real, intent(in) :: mass -- Mass of baryon
  • real, intent (in) :: rhoN,rhoP -- proton and neutron density in GeV^3
  • integer, optional, intent (in) :: monte_in -- if this input is given, then the number of Monte Carlo points is chosen according to this input; otherwise it is equal to "num_MonteCarlo_Points".

OUTPUT

  • real :: gcoll
  • real, optional, intent(out) :: gcoll_elastic_out -- elastic collisional width

inMediumWidth/calcCollBroadening_Mes [ Functions ]

[ Top ] [ inMediumWidth ] [ Functions ]

NAME

function calcCollBroadening_Mes(particleID,momentumLRF,mass,rhoN,rhoP) RESULT(gcoll)

PURPOSE

  • Returns the in-medium width of mesons according to Gamma=sigma*rho*v
  • An average over the Fermi sea 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 meson is performed.

INPUTS

  • integer, intent(in) :: particleID -- ID of particle
  • real, intent(in) :: absP -- absolute Momentum
  • real, intent(in) :: mass -- Mass of meson
  • real, intent(in) :: rhoN,rhoP -- proton and neutron density in GeV^3

OUTPUT

  • real :: gcoll

inMediumWidth/fermiMom [ Functions ]

[ Top ] [ inMediumWidth ] [ Functions ]

NAME

real function fermiMom(rho)

PURPOSE

  • Returns the fermi momentum.

INPUTS

  • real, intent(in) :: rho -- proton or neutron density in GeV^-3

OUTPUT