diff --git a/NPLib/Physics/NPReaction.h b/NPLib/Physics/NPReaction.h
index 71388278fdde2cbc1408432c84c3071436e447df..7200a8edffdec255810ff1ab9350cb0414e56237 100644
--- a/NPLib/Physics/NPReaction.h
+++ b/NPLib/Physics/NPReaction.h
@@ -243,7 +243,7 @@ namespace NPL {
 
     // Return Excitation Energy
     double ReconstructRelativistic(double EnergyLab, double ThetaLab, double PhiLab = 0);
-    
+
     // Return Lorentz Vector of Heavy Recoil after reaction
     TLorentzVector LorentzAfterReaction(double EnergyLab, double ThetaLab, double PhiLab = 0);
 
@@ -261,6 +261,11 @@ namespace NPL {
 
     void SetParticle3(double EnergyLab, double ThetaLab);
 
+    void SetParticle1(Beam p) { fParticle1 = p; };
+    void SetParticle2(Particle p) { fParticle2 = p; };
+    void SetParticle3(Particle p) { fParticle3 = p; };
+    void SetParticle4(Particle p) { fParticle4 = p; };
+
     TGraph* GetKinematicLine3(double AngleStep_CM = 1);
     TGraph* GetKinematicLine4(double AngleStep_CM = 1);
     TGraph* GetBrhoLine3(double AngleStep_CM = 1);