I started playing with the phase space generator TGenPhaseSpace inside ROOT and wanted to compute the 3-particle phase space volume for a muon decaying into 3 massless particles with the simple attached macro. The phase space volume should depend on the total CM energy like:
E^2/(4(4pi)^3)
as derived for example in
http://www.ippp.dur.ac.uk/~krauss/Lectures/QuarksLeptons/Basics/PS_2.html
However, when adding up the weights with the macro one always gets a value around ~0.25 regardless of how big the CM energy really is.
Being the code already used in a number of applications I guess something is missing in how the phase space generator should be used.
Does somebody know what am I doing wrong?
Many thanks in advance!
Gionata
PSvolume.cc (1.42 KB)