TABLE OF CONTENTS
- 1. /deltaWidth
- 1.1. deltaWidth/deltawidth
- 1.2. deltaWidth/deltaSwitch
- 1.3. deltaWidth/delSpread
- 1.4. deltaWidth/osetDelta_used
- 1.5. deltaWidth/delta_nucleonPion
- 1.6. deltaWidth/delta_fullWidth
- 1.7. deltaWidth/InitDelWidth
- 1.8. deltaWidth/get_g1P
- 1.9. deltaWidth/deloset
- 1.10. deltaWidth/deltaOset_ND_ND
- 1.11. deltaWidth/deltaOset_ND_NN
- 1.12. deltaWidth/cross
- 1.13. deltaWidth/PrintFullWidth
- 1.14. deltaWidth/GetMaxQ_Delta
- 1.15. deltaWidth/GetRhoValue
- 1.16. deltaWidth/GetRhoBin
/deltaWidth [ 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 = use Oset self energies+BUU input
- 2 = Spreading potential
- 3 = use 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/osetDelta_used [ Global module-variables ]
[ Top ] [ deltaWidth ] [ Global module-variables ]
NAME
logical function osetDelta_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