### 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