9 #ifndef ThePEG_RhoDMatrix_H
10 #define ThePEG_RhoDMatrix_H
13 #include "ThePEG/PDT/PDT.h"
14 #include "ThePEG/Helicity/HelicityDefinitions.h"
46 for(
size_t ix=0; ix<
_ispin; ++ix)
47 for(
size_t iy=0; iy<
_ispin; ++iy)
48 _matrix[ix][iy] = (average && ix==iy) ? 1./_ispin : 0.;
79 static const double epsa=1e-30, epsb=1e-10;
82 for(
size_t ix=0; ix<
_ispin; ++ix)
84 assert(norm.real() > epsa);
85 assert(norm.imag()/norm.real() < epsb);
86 double invnorm = 1./norm.real();
87 for(
size_t ix=0; ix<
_ispin; ++ix)
88 for(
size_t iy=0; iy<
_ispin; ++iy)
122 enum { MAXSPIN = 5 };
135 for (
size_t ix = 0; ix < rd.
_ispin; ++ix) {
136 for (
size_t iy = 0; iy < rd.
_ispin; ++iy)
137 os << rd.
_matrix[ix][iy] <<
" ";