src/cpp/jail/resonanceDecay/reso_decay.cpp
/********************************************************************
*
*   Edndp3_3bodyN()
*
* transverse momentum spectrum in GeV^-2 from pions out of resonances
*********************************************************************/
Edndp3_3bodyN()
transverse momentum spectrum in GeV^-2 from pions out of resonances
Edndp3_3bodyN () transverse momentum spectrum in GeV^ - 2 from pions out of resonances
double Freeze::Edndp3_3bodyN (double y, double pt, double phi, double m1, double m2,
              double m3, double mr, double norm3, int res_num)
{
  double mt = sqrt (pt * pt + m1 * m1);
  pblockN para;
  double wmin, wmax;
  double res3;
  para.pt = pt;
  para.mt = mt;
  para.y = y;
  para.e = mt * cosh (y);
  para.pl = mt * sinh (y);
  para.phi = phi;
  para.m1 = m1;
  para.m2 = m2;
  para.m3 = m3;
  para.mr = mr;
  para.res_num = res_num;
  wmin = (m2 + m3) * (m2 + m3); 
  wmax = (mr - m1) * (mr - m1);
  res3 = 2.0 * norm3 * gauss (PTS4, &Freeze::dn3ptN, wmin, wmax, &para) / mr;  
  if (res3<0.) res3=0.; 
  return res3;
}
double Freeze :: Edndp3_3bodyN ( double y , double pt , double phi , double m1 , double m2 , double m3 , double mr , double norm3 , int res_num ) { double mt = sqrt ( pt * pt + m1 * m1 ) ; pblockN para ; double wmin , wmax ; double res3 ; para . pt = pt ; para . mt = mt ; para . y = y ; para . e = mt * cosh ( y ) ; para . pl = mt * sinh ( y ) ; para . phi = phi ; para . m1 = m1 ; para . m2 = m2 ; para . m3 = m3 ; para . mr = mr ; para . res_num = res_num ; wmin = ( m2 + m3 ) * ( m2 + m3 ) ; wmax = ( mr - m1 ) * ( mr - m1 ) ; res3 = 2.0 * norm3 * gauss ( PTS4 , & Freeze :: dn3ptN , wmin , wmax , & para ) / mr ; if ( res3 < 0. ) res3 = 0. ; return res3 ; }
Edndp3_3bodyN()
transverse momentum spectrum in GeV^-2 from pions out of resonances
Edndp3_3bodyN () transverse momentum spectrum in GeV^ - 2 from pions out of resonances
/* in units of GeV^-2,includes phasespace and volume,
           does not include degeneracy factors  */ //   double slope;          /* slope of resonance for high mt */ //   int pn; //   pn = partid[MHALF + res_num]; //Integrates "W" using gaussian 