src/cpp/jail/resonanceDecay/reso_decay.cpp
/********************************************************************
*
*   Edndp3_2bodyN()
*
* transverse momentum spectrum in GeV^-2 from pions out of resonances
*********************************************************************/
Edndp3_2bodyN()
transverse momentum spectrum in GeV^-2 from pions out of resonances
Edndp3_2bodyN () transverse momentum spectrum in GeV^ - 2 from pions out of resonances
double Freeze::Edndp3_2bodyN (double y, double pt, double phi, double m1, double m2, double mr, int res_num)
{
  double mt = sqrt (pt * pt + m1 * m1);
  double norm2;         
  pblockN para;
  double res2;
  para.pt = pt;
  para.mt = mt;
  para.e = mt * cosh (y);
  para.pl = mt * sinh (y);
  para.y = y;
  para.phi = phi;
  para.m1 = m1;
  para.m2 = m2;
  para.mr = mr;
  para.res_num = res_num;
  norm2 = 1.0 / (2.0 * PI);
  res2 = norm2 * dn2ptN (m2 * m2, &para); 
  if (res2<0.) res2=0.;
  return res2;          
}
double Freeze :: Edndp3_2bodyN ( double y , double pt , double phi , double m1 , double m2 , double mr , int res_num ) { double mt = sqrt ( pt * pt + m1 * m1 ) ; double norm2 ; pblockN para ; double res2 ; para . pt = pt ; para . mt = mt ; para . e = mt * cosh ( y ) ; para . pl = mt * sinh ( y ) ; para . y = y ; para . phi = phi ; para . m1 = m1 ; para . m2 = m2 ; para . mr = mr ; para . res_num = res_num ; norm2 = 1.0 / ( 2.0 * PI ) ; res2 = norm2 * dn2ptN ( m2 * m2 , & para ) ; if ( res2 < 0. ) res2 = 0. ; return res2 ; }
Edndp3_2bodyN()
transverse momentum spectrum in GeV^-2 from pions out of resonances
Edndp3_2bodyN () 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 y;           /\* rapidity of particle 1       *\/ */ /*      double pt;          /\* transverse momentum of particle 1    *\/ */ /*      double phi;     /\* phi angle of particle 1      *\/ */ /*      double m1, m2;      /\* restmasses of decay particles in MeV *\/ */ /*      double mr;          /\* restmass of resonance MeV            *\/ */ /*      int res_num;        /\* Montecarlo number of the Resonance   */ /* 2-body normalization         */ //Calls the integration routines for 2-body /* like Ed3ndp3_2body() */