* Fixing Kinematic and Example1

        - If statement was commented out for inverse kinematic case

Examples/Example1/RunToTreat.txt

 TTreeName 
 	SimulatedTree
 RootFileName 
-	../../Outputs/Simulation/Example1.root
+	./root/simulation/Example1.root
 

Examples/Example1/ShowResults.C

@@ -6,6 +6,7 @@ TChain* chain=NULL ;
 
 ////////////////////////////////////////////////////////////////////////////////
 void LoadCuts(){
+    std::cout << "Loading Cuts" << std::endl;
     TFile* File_ETOF = new TFile("cuts/ETOF.root","READ");
     ETOF = (TCutG*) File_ETOF->FindObjectAny("ETOF");
     
@@ -15,6 +16,8 @@ void LoadCuts(){
 
 ////////////////////////////////////////////////////////////////////////////////
 void LoadChain(){
+
+    std::cout << "Loading Chain" << std::endl;
     chain = new TChain("PhysicsTree");
     chain->Add("root/analysis/Example1.root");
 }

Examples/Example1/configs/ConfigMust2.dat

@@ -14,4 +14,4 @@ ConfigMust2
 	DISABLE_ALL MM6
 	DISABLE_ALL MM7
 	DISABLE_ALL MM8
-  CSI_SIZE 256
+  CSI_SIZE 32

Inputs/DetectorConfiguration/Example1.detector

@@ -94,8 +94,6 @@ M2Telescope
  CSI=    0
  VIS=    all
 
-%%%%%%%%%%%%%%%%%%%%%     
-SSSDArray
 %%%%%%%%%% Det 1 %%%%%%%%     
 SSSD
  A= 17.61 9.85  104.11 mm

NPLib/Physics/NPReaction.cxx

@@ -261,10 +261,10 @@ void Reaction::KineRelativistic(double& ThetaLab3, double& KineticEnergyLab3, do
   // case of inverse kinematics
 
   double theta = fThetaCM;
-  /*if (m1 > m2){
+  if (m1 > m2){
     theta = M_PI - fThetaCM;
     //fThetaCM = M_PI - fThetaCM;
-  }*/
+  }
   fEnergyImpulsionCM_3 = TLorentzVector(pCM_3 * sin(theta), 0, pCM_3 * cos(theta), ECM_3);
   fEnergyImpulsionCM_4 = fTotalEnergyImpulsionCM - fEnergyImpulsionCM_3;