Ticket #125: decays-flux.patch

PHASIC++/Main/Phase_Space_Handler.C

@@ -208 +208 @@
 {
   if (m_nin == 1) {
     m_E = m_m[0];
-    m_flux = 1./(2.*m_E);
     m_s = m_E*m_E;
     p[0] = Vec4D(m_E,0.,0.,0.);
     return 1;

PHASIC++/Main/Phase_Space_Handler.H

@@ -97 +97 @@
     long int m_trials, m_maxtrials, m_sumtrials, m_events;
 
     double m_E, m_m[2], m_m2[2], m_s, m_smin, m_psweight;
-    double m_result_1, m_result_2, m_flux, m_weight, m_error; 
+    double m_result_1, m_result_2, m_weight, m_error; 
     
     std::vector<std::vector<double> > m_stats;
 

PHASIC++/Process/Single_Process.C

@@ -71 +71 @@
 {
   if (!m_use_biweight) return 1.;
   double wgt(1.0);
-  if (m_nin!=2) return wgt;
+  if (m_nin!=2) return 0.5/p_int->Momenta()[0].Mass();
   if (p_int->ISR()) {
     wgt*=p_int->ISR()->Weight
       (mode,p_int->Momenta()[0],p_int->Momenta()[1],