Fixed bug where electron masses were subtracted twice in kinematic correction

This commit is contained in:
Gordon McCann 2022-01-06 14:35:05 -05:00
parent 220dae97b8
commit be2ff0470c

View File

@ -64,10 +64,10 @@ namespace EventBuilder {
B /= 10000; //convert to tesla B /= 10000; //convert to tesla
angle *= DEGTORAD; angle *= DEGTORAD;
MT = MASS.FindMass(ZT, AT) - ZT*RESTMASS_ELECTRON*UTOMEV; MT = MASS.FindMass(ZT, AT);
MP = MASS.FindMass(ZP, AP) - ZP*RESTMASS_ELECTRON*UTOMEV; MP = MASS.FindMass(ZP, AP);
ME = MASS.FindMass(ZE, AE) - ZE*RESTMASS_ELECTRON*UTOMEV; ME = MASS.FindMass(ZE, AE);
MR = MASS.FindMass(ZR, AR) - ZR*RESTMASS_ELECTRON*UTOMEV; MR = MASS.FindMass(ZR, AR);
if (MT*MP*ME*MR == 0) if (MT*MP*ME*MR == 0)
{ {