diff --git a/NPSimulation/Core/NPIonIonInelasticPhysic.cc b/NPSimulation/Core/NPIonIonInelasticPhysic.cc
index 71a72955a8a02ba5e2b4087c75f13c742a8d8ae2..7993915244bce612f939fd1a892e7a589fe4bdda 100755
--- a/NPSimulation/Core/NPIonIonInelasticPhysic.cc
+++ b/NPSimulation/Core/NPIonIonInelasticPhysic.cc
@@ -64,6 +64,7 @@
 #include "G4IonInelasticProcess.hh"
 #include "G4GeneralSpaceNNCrossSection.hh"
 
+using namespace std;
 /////////////////////////////////////////////////////////////////////////////
 NPIonIonInelasticPhysic::NPIonIonInelasticPhysic(const G4String& name):
 G4VPhysicsConstructor(name)
@@ -99,7 +100,7 @@ void NPIonIonInelasticPhysic::ConstructProcess()
     G4IonsShenCrossSection* ShenCrossSections = new G4IonsShenCrossSection;
     G4ComponentGGHadronNucleusXsc* GlauberGribovCrossSection = new G4ComponentGGHadronNucleusXsc;
     
-    G4CrossSectionInelastic* GlauberGribovDataSet = new G4CrossSectionInelastic(GlauberGribovCrossSection);
+    //G4CrossSectionInelastic* GlauberGribovDataSet = new G4CrossSectionInelastic(GlauberGribovCrossSection);
     
     // ******************
     // **** Elastic ****
diff --git a/NPSimulation/Detectors/NeutronWall/NeutronWall.cc b/NPSimulation/Detectors/NeutronWall/NeutronWall.cc
index 7c770cabc05616760a8a8ac139f319a83a1c1a8e..aab029a161b1d281fde81e732fb8a13a8fbf4052 100644
--- a/NPSimulation/Detectors/NeutronWall/NeutronWall.cc
+++ b/NPSimulation/Detectors/NeutronWall/NeutronWall.cc
@@ -60,9 +60,9 @@ using namespace CLHEP;
 namespace NeutronWall_NS{
     // Energy and time Resolution
     const double EnergyThreshold = 0.1*MeV;
-    const double ResoTime = 0*ns ;
-    const double ResoEnergy = 0*MeV ;
-    const double ResoPosition = 0*cm;
+    //const double ResoTime = 0*ns ;
+    //const double ResoEnergy = 0*MeV ;
+    //const double ResoPosition = 0*cm;
     //The size of NS should depend on the distance between NeutronWall and plastic Bar right now
     double NS_X = 2020.0*mm;
     double NS_Y = 2020.0*mm;
diff --git a/NPSimulation/Detectors/Tiara/Tiara.cc b/NPSimulation/Detectors/Tiara/Tiara.cc
index 813fe4a2ebf5d7260621e5ad1c9265fc813ecb4d..89dd7a7c04b11ae13460bfa0e474dc470c7de923 100644
--- a/NPSimulation/Detectors/Tiara/Tiara.cc
+++ b/NPSimulation/Detectors/Tiara/Tiara.cc
@@ -283,7 +283,7 @@ void Tiara::ReadSensitive(const G4Event* event){
     double EF = RandGauss::shoot(Info[0],ResoEnergyHyball);
     if(EF>EnergyThreshold){
       int RingNumber=Info[8];
-      RingNumber=abs(RingNumber-17);                                                                                 //
+      RingNumber=abs(RingNumber-17);    
       Info[8]=RingNumber; 
       m_EventHyball->SetRingE(Info[7],Info[8],EF);
       m_EventHyball->SetRingT(Info[7],Info[8],Info[1]); 
diff --git a/NPSimulation/Detectors/Tiara/Tiara.hh b/NPSimulation/Detectors/Tiara/Tiara.hh
index 58b838ac0429f4a63c337460d4765ef3c51fb4c0..d24a277030a10b5368cfa9e9243feaa93b0f99f3 100644
--- a/NPSimulation/Detectors/Tiara/Tiara.hh
+++ b/NPSimulation/Detectors/Tiara/Tiara.hh
@@ -48,9 +48,9 @@ using namespace CLHEP;
 namespace TIARA{
   // Energy and time Resolution
   const G4double ResoTime    = 0      ;
-  const G4double ResoEnergyInnerBarrel  = 0.058*MeV ;// = 136keV of Resolution   //   Unit is MeV/2.35
-  const G4double ResoEnergyOuterBarrel  = 0.058*MeV ;// = 136keV of Resolution   //   Unit is MeV/2.35
-  const G4double ResoEnergyHyball       = 0.029*MeV ;// =  70keV of Resolution   //   Unit is MeV/2.35
+  const G4double ResoEnergyInnerBarrel  = 0.058*MeV ;// = 136keV FWHM 
+  const G4double ResoEnergyOuterBarrel  = 0.058*MeV ;// = 136keV FWHM
+  const G4double ResoEnergyHyball       = 0.029*MeV ;// =  70keV FWHM
 
   const G4double EnergyThreshold = 200*keV;
 
diff --git a/Projects/T40/22Nadp.reaction b/Projects/T40/22Nadp.reaction
index 09c0acb3555ba507d3cebba4852ecd849fa5f0d3..099142ed7d6973786ae3614d8292e558cd794e70 100644
--- a/Projects/T40/22Nadp.reaction
+++ b/Projects/T40/22Nadp.reaction
@@ -5,7 +5,7 @@ Beam
   Particle= 22Na
   ExcitationEnergy= 0
   Energy= 220
-  SigmaEnergy= 1
+  SigmaEnergy= 0
   SigmaThetaX= 0.01 
   SigmaPhiY= 0.01
   SigmaX= 0.5
diff --git a/Projects/T40/Analysis.cxx b/Projects/T40/Analysis.cxx
index 9245f7e61968ae86ce6fa027f0a821ae0594628a..e412f92306bc6069d57f491f69588dfaf80d8866 100644
--- a/Projects/T40/Analysis.cxx
+++ b/Projects/T40/Analysis.cxx
@@ -119,7 +119,8 @@ void Analysis::TreatEvent(){
     Energy = Si_E_TH;
 
     // Evaluate energy using the thickness 
-    ELab = LightAl.EvaluateInitialEnergy( Energy ,0.4*micrometer , ThetaTHSurface); 
+//    ELab = LightAl.EvaluateInitialEnergy( Energy ,0.4*micrometer , ThetaTHSurface); 
+    ELab = Energy;
     // Target Correction
     ELab = LightTarget.EvaluateInitialEnergy( ELab ,TargetThickness/2., ThetaNormalTarget); 
 
diff --git a/Projects/T40/run.mac b/Projects/T40/run.mac
new file mode 100644
index 0000000000000000000000000000000000000000..506556ed2882f22f4f6013ba592297aba63dec25
--- /dev/null
+++ b/Projects/T40/run.mac
@@ -0,0 +1 @@
+/run/beamOn 100000