Re: phase space decays

Hi Rene,

> TGenPhaseSpace::SetDecay is overestimating the maximum weigth
> for your process. This slows down the rejection technique used in
> TGenPhaseSpace::Generate, but the results should be correct.

ok, I see. Well, since there is no bug and TGenPhaseSpace::Generate is quite fast I can afford to do a quick search for the maximum weight myself.

Doing that, producing 'unweighted' events became almost as fast as producing 'weighted' events and gives identical results as before. I'll stick with that.

> Note that you can slightly speed up your script by running it
> with ACLIC.

Oh yes, we don't run the generator in interpreted mode from CINT. That was just an example.

Thank you,

        Costas

-- 
Dr. Constantinos Andreopoulos, CCLRC - Rutherford Lab.
http://hepunx.rl.ac.uk/~candreop/