TABLE OF CONTENTS
/threeBodyPhaseSpace [ Modules ]
NAME
module threeBodyPhaseSpace
NOTES
Includes routines which are necessary for the three-body phase-space. The subroutine momenta_in_3BodyPS is put into the module nBodyPhaseSpace.
threeBodyPhaseSpace/Integrate_3bodyPS [ Functions ]
[ Top ] [ threeBodyPhaseSpace ] [ Functions ]
NAME
function Integrate_3bodyPS (srts, mass1, mass2, mass3) result (ps)
PURPOSE
Evaluates Integral over three body phase space in vacuum
ps=d\Phi_3 *16*(2*pi)**7 (d\Phi_3 from PDG)
INPUTS
- real :: srts -- sqrt(s) in the problem
- real :: mass1,mass2,mass3 -- masses of the three particles
OUTPUT
- real :: ps = Integral of phase space
NOTES
Formerly known as "bops3".Is now faster than bops3. The code is so hard to read, since I wanted to make it faster. This routine is called often, therefore the optimization became important!
threeBodyPhaseSpace/Integrate_3bodyPS_Resonance [ Functions ]
[ Top ] [ threeBodyPhaseSpace ] [ Functions ]
NAME
function Integrate_3bodyPS_Resonance (srts, mass1, mass2, resonanceID, scalarPotential) result (ps)
PURPOSE
Evaluates Integral over three body phase space in vacuum with a resonance
among the three particles. Therefore one has to Integrate over the
mass of the resonance as well. ps=Integral d(massResonance) d\Phi_3 *16*(2*pi)**7 (d\Phi_3 from PDG)
INPUTS
- real :: mass1,mass2 = masses of the three particles
- integer :: idRes = Id of the resonance
- real,optional :: scalarPotential = scalarPotential of the resonance
- srts : sqrt(s) in the problem
OUTPUT
- real,dimension(1:2) :: ps = Integral of phase space
- ps(1) : Full width in the nominator of the spectral function
- ps(2) : Only (nucleon kaonBar) width in nominator of spectral function
NOTES
Formerly known as "massInt"
Be careful : since there is a mass cut off on the masses, there is no normalization off the spectral functions any more