EXP1803

Проблему понял, TGetPhaseSpace действительно ломается при этом условии, что описано в его методах: 

 fTeCmTm = P.Mag();           // total energy in C.M. minus the sum of the masses
   for (n=0;nfNt;n++) {
      fMass[n]  = mass[n];
      fTeCmTm  -= mass[n];
   }

   if (fTeCmTmreturn kFALSE;    // not enough energy for this decay

Можно решить проблему по аналогии с тем, как это было устроено в SIMONE: при условии, что продуктов реакции 2 кинематика рассчитывается одним способои, при условии что 3/5... другим. 

В связи с этим, скорее всего, нужно немного переделать механизм разыгрывания кинематики (создания 4-векторов) продуктов реакции. 

      Double_t mass5H = (fIs5HUserMassSet) ? f5HMass : f5H->Mass();
      LOG(DEBUG)

      Double_t reactMasses[2];
      reactMasses[0] = f3He->Mass();
      reactMasses[1] = mass5H;

      fReactionPhaseSpace->SetDecay(lvReaction, 2, reactMasses);
      fReactionPhaseSpace->Generate();

      TLorentzVector *lv3He = fReactionPhaseSpace->GetDecay(0);
      TLorentzVector *lv5H  = fReactionPhaseSpace->GetDecay(1);

эти строки нужно переделать: 

С точки зрения физики это правильно: "Метод TGenPhaseSpace::SetDecay() возвращает false, если суммарная энергия He6+H2 в ЦМ меньше суммы масс H5 и He3." Вас, Иван, попросли объяснить с т.з. физики, а вы о поправке кода. В чем физический смыс ваших изменений, и почему вы считаете, что указанное условие неверно?