double pyIn = 0., double pzIn = 0., double eIn = 0.,
double mIn = 0., double scaleIn = 0., ParticleData* ptr = 0)
: Particle(idIn, statusIn, mother1In, mother2In, daughter1In, daughter2In,
- colIn, acolIn, pxIn, pyIn, pzIn, eIn, mIn, scaleIn) {
+ colIn, acolIn, pxIn, pyIn, pzIn, eIn, mIn, scaleIn), idx(0), direction(1) {
if (ptr) { setPDTPtr(ptr); setPDEPtr(); }
rho = vector< vector<complex> >(spinStates(),
vector<complex>(spinStates(), 0));
int daughter1In, int daughter2In, int colIn, int acolIn, Vec4 pIn,
double mIn = 0., double scaleIn = 0., ParticleData* ptr = 0)
: Particle(idIn, statusIn, mother1In, mother2In, daughter1In, daughter2In,
- colIn, acolIn, pIn, mIn, scaleIn) {
+ colIn, acolIn, pIn, mIn, scaleIn), idx(0), direction(1) {
if (ptr) { setPDTPtr(ptr); setPDEPtr();}
rho = vector< vector<complex> >(spinStates(),
vector<complex>(spinStates(), 0));
for (int i = 0; i < spinStates(); i++) { rho[i][i] = 0.5; D[i][i] = 1;}
direction = 1; }
HelicityParticle(const Particle& ptIn, ParticleData* ptr = 0)
- : Particle(ptIn) {
+ : Particle(ptIn), idx(0), direction(1) {
if (ptr) { setPDTPtr(ptr); setPDEPtr();}
rho = vector< vector<complex> >(spinStates(),
vector<complex>(spinStates(), 0));