gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/deltaWidth [ Modules ]

[ Top ] [ Modules ]

NAME

module deltaWidth

PURPOSE

Implements the routines for the medium width of the Delta-Resonance P33(1232).

deltaWidth/deltawidth [ Namelists ]

[ Top ] [ deltaWidth ] [ Namelists ]

NAME

NAMELIST deltaWidth

PURPOSE

Includes the switch:

deltaWidth/deltaSwitch [ Global module-variables ]

[ Top ] [ deltaWidth ] [ Global module-variables ]

PURPOSE

Switch for different prescriptions for the delta width.

SOURCE

  integer, save :: deltaSwitch=3

NOTES

  • 1 = Oset self energies+BUU input
  • 2 = Spreading potential
  • 3 = Oset self energies
  • 4 = density dependent with BUU input

deltaWidth/delSpread [ Global module-variables ]

[ Top ] [ deltaWidth ] [ Global module-variables ]

PURPOSE

Spreading potential parameter.

SOURCE

  real, parameter :: delspread=0.04

deltaWidth/deltaOset_used [ Global module-variables ]

[ Top ] [ deltaWidth ] [ Global module-variables ]

NAME

logical function deltaOset_used

PURPOSE

return true, if deltaSwitch==3

deltaWidth/delta_nucleonPion [ Functions ]

[ Top ] [ deltaWidth ] [ Functions ]

NAME

real function delta_nucleonPion(mass,momentum,density)

PURPOSE

Returns partial width Delta->N Pi as a function of mass, momentum and density.

INPUTS

  • real :: mass -- mass of Delta
  • real :: density -- density at position
  • real, dimension(1:3) :: momentum -- momentum in LRF (only abs value used)

deltaWidth/delta_fullWidth [ Functions ]

[ Top ] [ deltaWidth ] [ Functions ]

NAME

real function delta_fullWidth(mass,momentum,density)

PURPOSE

Returns Delta width as a function of mass, momentum and density.

INPUTS

  • real :: mass -- mass of Delta
  • real :: density -- density at position
  • real, dimension(1:3) :: momentum -- momentum in LRF (only abs value used)

deltaWidth/InitDelWidth [ Subroutines ]

[ Top ] [ deltaWidth ] [ Subroutines ]

NAME

subroutine InitDelWidth

PURPOSE

Evaluates Delta width and stores it as a function of mass, momentum and density into the field medwidth.

deltaWidth/get_g1P [ Subroutines ]

[ Top ] [ deltaWidth ] [ Subroutines ]

NAME

subroutine get_g1P()

PURPOSE

Evaluates the partial width of a delta decaying into a pion nucleon systemfor a whole array and writes result to file as a function of mass, momentum and density. If it finds the information already stored in a file, it just reads in this file

NOTES

This subroutine calculates the in-medium width of a resonance in a fermi gas model.

deltaWidth/deloset [ Subroutines ]

[ Top ] [ deltaWidth ] [ Subroutines ]

NAME

subroutine deloset(mass,rho,imsig2,imsig3,imsigq)

PURPOSE

Evaluates Delta Width according to Oset and Salcedo NPA 468 (1987) 631-652

deltaWidth/deltaOset_ND_ND [ Subroutines ]

[ Top ] [ deltaWidth ] [ Subroutines ]

NAME

real function deltaOset_ND_ND(part)

PURPOSE

Evaluates ND->ND cross sections based on Delta Width according to Oset and Salcedo NPA 468 (1987) 631-652

INPUTS

  • type(particle),intent(in),dimension(1:2) :: part ! scattering particles

OUTPUT

Cross section in mb

deltaWidth/deltaOset_ND_NN [ Subroutines ]

[ Top ] [ deltaWidth ] [ Subroutines ]

NAME

real function deltaOset_ND_NN(part)

PURPOSE

Evaluates ND->NN cross sections based on Delta Width according to Oset and Salcedo NPA 468 (1987) 631-652

INPUTS

  • type(particle),intent(in),dimension(1:2) :: part ! scattering particles

OUTPUT

Cross section in mb

deltaWidth/cross [ Subroutines ]

[ Top ] [ deltaWidth ] [ Subroutines ]

NAME

function cross(part,med,switch) result(sigma)

PURPOSE

Evaluates ND->Nd (switch=1) or ND->NN (switch=2) cross sections based on Delta Width according to Oset and Salcedo NPA 468 (1987) 631-652

INPUTS

  • integer :: switch
  • type(particle),intent(in),dimension(1:2) :: part ! scattering particles

OUTPUT

Cross section in mb

deltaWidth/PrintFullWidth [ Subroutines ]

[ Top ] [ deltaWidth ] [ Subroutines ]

NAME

subroutine PrintFullWidth

PURPOSE

Print the full width in order to plot it with gnuplot as 2D plot:

  • version 1: (M,P) for different rho, file "DeltaFullWidth.gp1.dat"
  • version 2: (M,rho) for different P, file "DeltaFullWidth.gp2.dat"

In order to select the third variable, you have to use the parameter 'index' in gnuplot.

deltaWidth/GetMaxQ_Delta [ Subroutines ]

[ Top ] [ deltaWidth ] [ Subroutines ]

NAME

subroutine GetMaxQ_Delta(mass0,gamma0,rho,BinM,BinMaxQ)

PURPOSE

Calculate the maximal values of the Q weight for bins according BinM

INPUTS

  • integer :: ID -- ID of resonance
  • real :: mass0 -- pole mass
  • real :: gamma0 -- width at pole mass
  • real :: rho -- density
  • real, dimension(:) :: BinM -- array with boundaries for M binning

OUTPUT

  • real, dimension(:) :: BinMaxQ -- the maximal Q values for each bin.

NOTES

  • The size of BinMaxQ has to be at least the size of BinM minus 1.
  • It first calculates Q at the boundaries, then it iterates over the tabulated width values in order to take into account, that the Q value may be larger inbetween the boundaries.
  • if the Q value is maximal at the upper bound, we store its value as -Q.

deltaWidth/GetRhoValue [ Subroutines ]

[ Top ] [ deltaWidth ] [ Subroutines ]

NAME

subroutine GetRhoValue(rhomin,drho,nrho)

PURPOSE

return the variables concerning the density of the parametrization

INPUTS

none

OUTPUT

  • real :: rhomin, drho
  • integer :: nrho

deltaWidth/GetRhoBin [ Subroutines ]

[ Top ] [ deltaWidth ] [ Subroutines ]

NAME

subroutine GetRhoBin(rho,iRho)

PURPOSE

return the bin nummer connected to a certain density

INPUTS

  • real :: rho -- the density

OUTPUT

integer :: iRho -- the bin nummer