From edf25adcab10ecfeadc25d3ce74f1aaa531eb5c8 Mon Sep 17 00:00:00 2001
From: Gordon McCann <43148247+gwm17@users.noreply.github.com>
Date: Fri, 28 Jan 2022 14:09:07 -0500
Subject: [PATCH] Fixed double subtraction of electron masses in FP_kinematics

---
 src/evb/FP_kinematics.cpp | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/evb/FP_kinematics.cpp b/src/evb/FP_kinematics.cpp
index 9a558c9..cdb0f5b 100644
--- a/src/evb/FP_kinematics.cpp
+++ b/src/evb/FP_kinematics.cpp
@@ -61,10 +61,10 @@ double Delta_Z(int ZT, int AT, int ZP, int AP, int ZE, int AE,
   B /= 10000; //convert to tesla
   angle *= DEGTORAD;
 
-  MT = MASS.FindMass(ZT, AT) - ZT*RESTMASS_ELECTRON*UTOMEV;
-  MP = MASS.FindMass(ZP, AP) - ZP*RESTMASS_ELECTRON*UTOMEV;
-  ME = MASS.FindMass(ZE, AE) - ZE*RESTMASS_ELECTRON*UTOMEV;
-  MR = MASS.FindMass(ZR, AR) - ZR*RESTMASS_ELECTRON*UTOMEV;
+  MT = MASS.FindMass(ZT, AT);
+  MP = MASS.FindMass(ZP, AP);
+  ME = MASS.FindMass(ZE, AE);
+  MR = MASS.FindMass(ZR, AR);
   
   if (MT*MP*ME*MR == 0) {
     std::cerr << "***WARNING: error loading one or more masses; returning 0\n";