### TABLE OF CONTENTS

## /twoBodyPhaseSpace [ Modules ]

NAME

module **twoBodyPhaseSpace**

NOTES

Includes all the routines which are necessary for the two-body phase-space.

## twoBodyPhaseSpace/setMaxSqrts [ Subroutines ]

[ Top ] [ twoBodyPhaseSpace ] [ Subroutines ]

NAME

subroutine **setMaxSqrts**(srtsMax)

PURPOSE

set the upper bound of the tabulation

## twoBodyPhaseSpace/Integrate_2bodyPS_resonance [ Functions ]

[ Top ] [ twoBodyPhaseSpace ] [ Functions ]

NAME

function **Integrate_2bodyPS_resonance**(resID, srts, massStable, scalarPotential) result (ps)

PURPOSE

Returns Integral over the CM-momentum with a resonance among the two particles. Therefore one has to Integrate over the mass of the resonance as well:

ps=Integral { p_final(m_R) * Spectral function(m_R} d(m_R)

with

m_R = massResonance

It uses an internal routine "Calculate" which is either called directly or we store its values to a field.

This routine is not giving the real Phase space !!!

INPUTS

- integer, intent(in) :: resID -- Id of the resonance
- real, intent(in) :: srts -- sqrt(s) in the problem
- real, intent(in) :: massStable -- mass of the stable particle
- real, intent(in) :: scalarPotential -- scalarPotential of the resonance

OUTPUT

- real,dimension(1:5) :: ps -- Integral as given above

Here the different components are:

- ps(1): With full width evaluated in the nominator of spectral function
- ps(2): With partial width (pion N) evaluated in the nominator of spectral function
- ps(3): With partial width (eta N) evaluated in the nominator of spectral function
- ps(4): With partial width (rho N) evaluated in the nominator of spectral function
- ps(5): With partial width (omega N) evaluated in the nominator of spectral function

NOTES

Formerly known as "massInt2"

## twoBodyPhaseSpace/nnRR [ Subroutines ]

[ Top ] [ twoBodyPhaseSpace ] [ Subroutines ]

NAME

real function **nnRR**(srts, ID)

PURPOSE

Procedure for calculation of integrals of Resonance Resonance' CM-momentum

Evaluates Integral over two body phase space in vacuum for two baryon resonances. Therefore one has to Integrate over the mass of the two resonances as well.

ps=Integral d(mass_A ) d(mass_B ) p_AB * Spectralfunction_A Spectralfunction_B

INPUTS

- integer :: ID(1:2) -- Ids of the resonances
- real :: srts -- sqrt(s) of the reaction

OUTPUT

- Integral as given above

NOTES

To increase speed, we tabulate the possible output for Delta+R.