#include "METOOLS/SpinCorrelations/Amplitude2_Matrix.H" #include "ATOOLS/Org/Exception.H" #include "ATOOLS/Phys/Particle.H" using namespace METOOLS; using namespace ATOOLS; Amplitude2_Matrix::Amplitude2_Matrix(const ATOOLS::Particle* p) : p_part(p) { m_nhel=p_part->RefFlav().IntSpin()+1; if (m_nhel==3 && IsZero(p_part->RefFlav().Mass())) m_nhel=2; resize(m_nhel*m_nhel, Complex(0.0,0.0)); } Amplitude2_Matrix::~Amplitude2_Matrix() { } Complex Amplitude2_Matrix::Trace() const { Complex sum(0.,0.); for (size_t i(0); iRefFlav()<<":"<