[ROOT] TGenPhaseSpace bug?

From: F. Sabatie (sabatie@hep.saclay.cea.fr)
Date: Wed Mar 07 2001 - 11:34:19 MET


Hello,

has anyone used the TGenPhaseSpace class by any chance? I currently use
ROOT_v2.25 on RH6.1. My problem comes when I try to Generate() the phase
space, I get a core dump.

The interesting part of the program is really simple:

  TGenPhaseSpace ps;
  TLorentzVector cm_mom;
  Double_t masses;

  masses[0]=0.;
  masses[1]=mproton;

  cm_mom.SetPxPyPzE(0,0,-qmod,nu+mproton); // qmod, nu and mproton are
Double_t

  if ( ps.SetDecay(cm_mom,2,masses,"normal") ) cout << "SetDecay Success
!" << endl;
// This does produce a success, so SetDecay works

  ps.Generate();
// This core dumps !!!!!!!!



Anyone has an idea what's going on?

Regards,
Franck



This archive was generated by hypermail 2b29 : Tue Jan 01 2002 - 17:50:39 MET