diff --git a/Projects/Inelastic/RunToTreat.txt b/Projects/Inelastic/RunToTreat.txt
index 2a165eaafa002675ea5d8bfb5f9f9b6ca95de8dc..1ad4fdee63c263faa0fb252bd2361ed0c04b33ae 100644
--- a/Projects/Inelastic/RunToTreat.txt
+++ b/Projects/Inelastic/RunToTreat.txt
@@ -1,5 +1,4 @@
 TTreeName
   SimulatedTree
 RootFileName
-  ../../Outputs/Simulation/vendeta_el_*.root 
-  ../../Outputs/Simulation/vendeta_1st_inel_*.root 
+  ./root/simulation/neutron_56Fe_gs.root 
diff --git a/Projects/Inelastic/Vendeta_inelastic.detector b/Projects/Inelastic/Vendeta_inelastic.detector
index a28b4cb0c31fd35cfaf40d48dfee78b0ce984a6e..1cf7a74c1b207d03ecf0a8b3b6a94a00aecdb252 100644
--- a/Projects/Inelastic/Vendeta_inelastic.detector
+++ b/Projects/Inelastic/Vendeta_inelastic.detector
@@ -2,458 +2,83 @@
 Target 
   THICKNESS= 10 micrometer
   RADIUS= 20 mm
-  MATERIAL= 238U 
+  MATERIAL= Fe 
   ANGLE= 0 deg 
   X= 0 mm 
   Y= 0 mm 
   Z= 0 mm 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 20 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 30 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 40 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 50 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 60 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 70 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 80 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 90 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 100 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 110 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 120 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 130 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 140 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 150 deg
   PHI= 0 deg
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Vendeta
-  R= 2500 mm
+  R= 2000 mm
   THETA= 160 deg
   PHI= 0 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 20 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 30 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 40 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 50 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 60 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 70 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 80 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 90 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 100 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 110 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 120 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 130 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 140 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 150 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 160 deg
-  PHI= 20 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 20 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 30 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 40 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 50 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 60 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 70 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 80 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 90 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 100 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 110 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 120 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 130 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 140 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 150 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 160 deg
-  PHI= 40 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 20 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 30 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 40 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 50 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 60 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 70 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 80 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 90 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 100 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 110 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 120 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 130 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 140 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 150 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 160 deg
-  PHI= 140 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 20 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 30 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 40 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 50 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 60 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 70 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 80 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 90 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 100 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 110 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 120 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 130 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 140 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 150 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 160 deg
-  PHI= 160 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 20 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 30 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 40 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 50 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 60 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 70 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 80 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 90 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 100 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 110 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 120 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 130 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 140 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 150 deg
-  PHI= 180 deg
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-Vendeta
-  R= 2500 mm
-  THETA= 160 deg
-  PHI= 180 deg
diff --git a/Projects/Inelastic/macro/CheckKinematics.C b/Projects/Inelastic/macro/CheckKinematics.C
new file mode 100644
index 0000000000000000000000000000000000000000..09c5612cf441eb497615c9d574d734f0c2edf42b
--- /dev/null
+++ b/Projects/Inelastic/macro/CheckKinematics.C
@@ -0,0 +1,41 @@
+TRandom3 ra;
+
+void CheckKinematics(){
+
+  TFile* ofile = new TFile("kine.root","recreate");
+
+  NPL::Reaction* r1 = new NPL::Reaction("7Li(p,1n)7Be@13.95");
+
+  double beam_energy_init = r1->GetBeamEnergy();
+  double eloss = 0.070;
+  
+  double beam_energy;
+  double ThetaCM;
+  double E3, Theta3, E4, Theta4;
+  TTree* tree = new TTree("KineTree","KineTree");
+  tree->Branch("beam_energy",&beam_energy,"beam_energy/D");
+  tree->Branch("ThetaCM",&ThetaCM,"ThetaCM/D");
+  tree->Branch("E3",&E3,"E3/D");
+  tree->Branch("Theta3",&Theta3,"Theta3/D");
+  tree->Branch("E4",&E4,"E4/D");
+  tree->Branch("Theta4",&Theta4,"Theta4/D");
+
+
+  for(int i=0; i<1e6; i++){
+    beam_energy = ra.Uniform(beam_energy_init - eloss, beam_energy_init);
+
+    ThetaCM = ra.Uniform(0,180);
+    r1->SetThetaCM(ThetaCM*deg);
+    r1->SetBeamEnergy(beam_energy);
+    r1->KineRelativistic(Theta3,E3,Theta4,E4);
+
+    Theta3 = Theta3/deg;
+    Theta4 = Theta4/deg;
+    
+    tree->Fill();
+  }
+
+  tree->Write();
+
+  ofile->Close();
+}
diff --git a/Projects/Inelastic/macro/DirectKinematics.C b/Projects/Inelastic/macro/DirectKinematics.C
new file mode 100644
index 0000000000000000000000000000000000000000..a8d33489ffcdede6b2ccc68079f8e6b8c04b532f
--- /dev/null
+++ b/Projects/Inelastic/macro/DirectKinematics.C
@@ -0,0 +1,52 @@
+TRandom3 ra;
+double density_7Li = 0.534; // g/cm3
+///////////////////////////////////////////////////
+void DirectKinematics(double thickness_mg_cm2 = 0.3){
+
+  NPL::EnergyLoss* p7Li_loss;
+  p7Li_loss = new NPL::EnergyLoss("../EnergyLossTable/proton_7Li.G4table","G4Table",100);
+
+  TFile* ofile = new TFile("direct_kine.root","recreate");
+
+  NPL::Reaction* r1 = new NPL::Reaction("p(7Li,1n)7Be@2.8");
+
+  double thickness_micrometer;
+  thickness_micrometer = thickness_mg_cm2/(density_7Li*1e3)*1e4;
+  cout << thickness_mg_cm2 << " mg/cm2 -> " << thickness_micrometer << " um" << endl;
+
+  double beam_energy_init = r1->GetBeamEnergy();
+  double eloss;
+  eloss = p7Li_loss->EnergyLossCalculation(beam_energy_init,thickness_micrometer*micrometer,0);
+
+  cout << "Energy Loss = " << eloss << " MeV" << endl;
+
+  double beam_energy;
+  double ThetaCM;
+  double E3, Theta3, E4, Theta4;
+  TTree* tree = new TTree("KineTree","KineTree");
+  tree->Branch("beam_energy",&beam_energy,"beam_energy/D");
+  tree->Branch("ThetaCM",&ThetaCM,"ThetaCM/D");
+  tree->Branch("E3",&E3,"E3/D");
+  tree->Branch("Theta3",&Theta3,"Theta3/D");
+  tree->Branch("E4",&E4,"E4/D");
+  tree->Branch("Theta4",&Theta4,"Theta4/D");
+
+
+  for(int i=0; i<1e6; i++){
+    beam_energy = ra.Uniform(beam_energy_init - eloss, beam_energy_init);
+
+    ThetaCM = ra.Uniform(0,180);
+    r1->SetThetaCM(ThetaCM*deg);
+    r1->SetBeamEnergy(beam_energy);
+    r1->KineRelativistic(Theta3,E3,Theta4,E4);
+
+    Theta3 = Theta3/deg;
+    Theta4 = Theta4/deg;
+
+    tree->Fill();
+  }
+
+  tree->Write();
+
+  ofile->Close();
+}
diff --git a/Projects/Inelastic/neutron_inelastic.reaction b/Projects/Inelastic/neutron_inelastic.reaction
new file mode 100644
index 0000000000000000000000000000000000000000..66ecd1606323e1b52941d20874eca7b3a2da806f
--- /dev/null
+++ b/Projects/Inelastic/neutron_inelastic.reaction
@@ -0,0 +1,28 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%%%% Reaction file for 11Li(d,3He)10He reaction %%%%%%%%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+Beam
+ Particle= 1n
+ Energy= 1.1 MeV
+ SigmaEnergy= 0.016 MeV
+ SigmaThetaX= 0. deg
+ SigmaPhiY= 0. deg
+ SigmaX= 0 mm
+ SigmaY= 0 mm
+ MeanThetaX= 0 deg
+ MeanPhiY= 0 deg
+ MeanX= 0 mm
+ MeanY= 0 mm
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+TwoBodyReaction
+ Beam= 1n
+ Target= 56Fe
+ Light= 1n
+ Heavy= 56Fe
+ ExcitationEnergy3= 0.0 MeV
+ ExcitationEnergy4= 0.0 MeV
+ CrossSectionPath= flat.txt CS10He
+ ShootLight= 1
+ ShootHeavy= 1
+  
diff --git a/Projects/Inelastic/project.config b/Projects/Inelastic/project.config
new file mode 100644
index 0000000000000000000000000000000000000000..c19d2a29e602f4c6a6e78f6e87a02866c785ff81
--- /dev/null
+++ b/Projects/Inelastic/project.config
@@ -0,0 +1,5 @@
+Project e850
+ AnalysisOutput= ./root/analysis 
+ SimulationOutput= ./root/simulation
+ EnergyLoss= ./EnergyLossTable
+