#include "PHASIC++/Channels/Decay_Dalitz.H" #include "PHASIC++/Channels/Channel_Elements.H" using namespace PHASIC; using namespace ATOOLS; using namespace std; Decay_Dalitz::Decay_Dalitz( const ATOOLS::Flavour * fl, const double& mass, const double& width, size_t dir, size_t p1, size_t p2, const ATOOLS::Mass_Selector* masssel) : Single_Channel(1,3,fl), m_decvec(Vec4D(fl[0].HadMass(),0.,0.,0.)), m_pmass(mass), m_pwidth(width), m_sexp(.5), m_p1(p1), m_p2(p2), m_dir(dir), m_mode(0), p_masssel(masssel) { for (short int i=0;iMass2(fl[i]); m_smin = ATOOLS::sqr(p_masssel->Mass(fl[m_p1])+p_masssel->Mass(fl[m_p2])); m_smax = ATOOLS::sqr(p_masssel->Mass(fl[0])-p_masssel->Mass(fl[m_dir])); if (sqrt(m_smin)