TABLE OF CONTENTS
- 1. /mesonWidthVacuum
- 1.1. mesonWidthVacuum/interactionRadius
- 1.2. mesonWidthVacuum/srts_srt_switch
- 1.3. mesonWidthVacuum/omega_width
- 1.4. mesonWidthVacuum/readInput
- 1.5. mesonWidthVacuum/MesonWidthVacuum
- 1.6. mesonWidthVacuum/vacuumWidth
- 1.7. mesonWidthVacuum/stableFinalState
- 1.8. mesonWidthVacuum/semistableFinalState
- 1.9. mesonWidthVacuum/rho_AB_Integrand
- 1.10. mesonWidthVacuum/omegaVacuum
- 1.11. mesonWidthVacuum/dileptonWidth
/mesonWidthVacuum [ Modules ]
NAME
module mesonWidthVacuum
NOTES
Module which calculates the partial and full widths of the meson resonances in dependence of their mass. Their pole mass is given by their mass in 'particleProperties' and the widths at this pole mass as well. Everything corresponds to the vacuum situation. The resulting width is therefore always the width in vacuum! As mass of the resonance we use the four-vector definition: p_mu p^mu= mass**2 Prescription according to Manley et al. Phys. Rev. D45 (1992) 4002. The In-Width is assumed to be the outwidth.
mesonWidthVacuum/interactionRadius [ Global module-variables ]
[ Top ] [ mesonWidthVacuum ] [ Global module-variables ]
SOURCE
real, parameter, public :: interactionRadius = 1./hbarc
PURPOSE
interaction radius for Blatt-Weisskopf functions: r = 1fm
mesonWidthVacuum/srts_srt_switch [ Global module-variables ]
[ Top ] [ mesonWidthVacuum ] [ Global module-variables ]
SOURCE
logical, save, public :: srts_srt_switch=.false.
PURPOSE
Modifies the width according to S. Leupold's definition of the width, one especially has to exchange s against sqrt(s) in the denominator of Formula 2.76 of Effenbergers Phd
NOTES
The default value is .false. and the power of the mass resp. sqrt(s) is 1. If the flag is .true., the power is 2.
mesonWidthVacuum/omega_width [ Global module-variables ]
[ Top ] [ mesonWidthVacuum ] [ Global module-variables ]
SOURCE
integer, save :: omega_width = 1
PURPOSE
Select a parametrization for the omega vacuum width:
mesonWidthVacuum/readInput [ Subroutines ]
[ Top ] [ mesonWidthVacuum ] [ Subroutines ]
NAME
subroutine readInput
PURPOSE
Reads input in jobcard out of namelist "MesonWidthVacuum".
mesonWidthVacuum/MesonWidthVacuum [ Namelists ]
[ Top ] [ mesonWidthVacuum ] [ Namelists ]
NAME
NAMELIST /MesonWidthVacuum/
PURPOSE
Includes the input switches:
mesonWidthVacuum/vacuumWidth [ Subroutines ]
[ Top ] [ mesonWidthVacuum ] [ Subroutines ]
NAME
function vacuumWidth(mass, partID, ratio) result (gammaTotal)
NOTES
This routine calculates the total vacuum decay width of a mesonic resonance, i.e. the sum of all partial decay widths, and the branching ratios for each channel as a function of the offshell mass (m^2 = p_mu p^mu). Parameters taken from Manley et al. Phys. Rev. D45 (1992) 4002 and PDG. The mass dependence of the resonances is treated according to Manley.
INPUTS
- real, intent(in) :: mass -- offshell mass of Resonance in GeV
- integer, intent(in) :: partID -- ID of resonance
OUTPUT
- real, dimension(1:nDecays), optional :: ratio -- branching ratio, i.e. partial width / full width for all decay channels of the mesons
- real :: gammaTotal -- full width (mass) in GeV at the mass of "mass"
mesonWidthVacuum/stableFinalState [ Functions ]
[ Top ] [ mesonWidthVacuum ] [ Functions ]
NAME
real function stableFinalState(mass, poleMass, mass1, mass2, L, partialWidth_pole)
INPUTS
- integer :: L --- angular momentum
- real :: polemass --- pole mass of mother resonance
- real :: mass --- mass of mother resonance
- real :: mass1 --- mass of first decay product
- real :: mass2 --- mass of second decay product
- real :: partialWidth_pole --- partial widht of decayChannel at pole mass
NOTES
Resonance decays into stable particles. Decay only allowed if:
Manley et al. Phys. Rev. D45 (1992) 4002.
OUTPUT
mesonWidthVacuum/semistableFinalState [ Functions ]
[ Top ] [ mesonWidthVacuum ] [ Functions ]
NAME
real function semistableFinalState(mass,polemass,ID,L,partialWidth_pole)
INPUTS
- real :: mass --- mass of mother resonance
- real :: polemass --- polemass of mother resonance
- integer :: ID(1:2) --- outgoing IDs: 1=stable, 2=unstable
- integer :: L --- angular momentum of final state
- real :: partialWidth_pole --- partial width at pole of mother resonance
NOTES
Resonance decays into one stable and one unstable particle. Calculates the partial width dependend on mass of a meson resonance decaying into one stable and one unstable decay product. According to Manley and Effenberger' Dr. Thesis equation 2.76.
OUTPUT
mesonWidthVacuum/rho_AB_Integrand [ Functions ]
[ Top ] [ mesonWidthVacuum ] [ Functions ]
NAME
real function rho_AB_Integrand(mu)
NOTES
This is the integrand of rho_ab(mu) in Effenbergers Dr.-Thesis, page 27, formula 2.76.
Only in the case of only one unstable particle, currently only used for decays into "rho pi"!
mesonWidthVacuum/omegaVacuum [ Functions ]
[ Top ] [ mesonWidthVacuum ] [ Functions ]
NAME
real function omegaVacuum(s, r)
INPUTS
real, intent(in) :: s ! m^2 of the omega meson integer, intent(in) :: r ! select channel
NOTES
Returns the partial omega widths in the vacuum:
- r=0 : total width
- r=1 : omega -> pi pi
- r=2 : omega -> pi0 gamma
- r=3 : omega -> rho pi (as an approximation for omega -> 3pi)
See: Pascal Muehlich, PhD thesis, chapter 8.2.1.
mesonWidthVacuum/dileptonWidth [ Subroutines ]
[ Top ] [ mesonWidthVacuum ] [ Subroutines ]
NAME
real function dileptonWidth(ID, mass)
NOTES
Calculate the dilepton decay width of the vector mesons rho, omega and phi (V -> e+e-) according to strict vector-meson dominance (VMD). References:
- http://arxiv.org/abs/1203.3557v2, eq. (11) and table 3.
- Effenberger PhD, chapter 2.7.1
- Corrections due to finite electron mass see, for example, in https://arxiv.org/abs/nucl-th/0005041, eq. (26)
INPUTS
- integer :: ID --- vector meson ID (should be 103, 105 or 107)
- real :: mass --- vector meson mass in GeV
OUTPUT
- width in GeV