TVector3 momTot(Px(),Py(),Pz());
Double_t lQlNeg = momNeg.Dot(momTot)/momTot.Mag();
- Double_t lQlPos = momNeg.Dot(momTot)/momTot.Mag();
+ Double_t lQlPos = momPos.Dot(momTot)/momTot.Mag();
return 1.-2./(1.+lQlNeg/lQlPos);
}
//
// calculate effective mass
//
- const Float_t kpmass[5] = {5.10000000000000037e-04,1.05660000000000004e-01,1.39570000000000000e-01,
- 4.93599999999999983e-01, 9.38270000000000048e-01};
+ const Float_t kpmass[5] = {TDatabasePDG::Instance()->GetParticle(kElectron)->Mass(),
+ TDatabasePDG::Instance()->GetParticle(kMuonMinus)->Mass(),
+ TDatabasePDG::Instance()->GetParticle(kPiPlus)->Mass(),
+ TDatabasePDG::Instance()->GetParticle(kKPlus)->Mass(),
+ TDatabasePDG::Instance()->GetParticle(kProton)->Mass()};
if (p1>4) return -1;
if (p2>4) return -1;
Float_t mass1 = kpmass[p1];