diff --git a/Inputs/EventGenerator/gamma.source b/Inputs/EventGenerator/gamma.source
index 877c06c326b9b0b11b0bc3c7e0805457e8b116d9..2e6472d36109149de8c84a031d2bfdac9cf64bb6 100644
--- a/Inputs/EventGenerator/gamma.source
+++ b/Inputs/EventGenerator/gamma.source
@@ -4,8 +4,8 @@
 %			   Energy are given in MeV , Position in mm				  %	
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Isotropic
-	EnergyLow=  .1
-	EnergyHigh= 5.
+	EnergyLow=  1
+	EnergyHigh= 1.
 	HalfOpenAngleMin= 0
 	HalfOpenAngleMax= 180
 	x0= 0	
diff --git a/NPSimulation/Paris/Paris.cc b/NPSimulation/Paris/Paris.cc
index 6b61ffcd8a6ef1dd7679e1258e77fbb9ab4e0e27..00705f42162c69dd2700ef62dfa59d17e77eeaa6 100644
--- a/NPSimulation/Paris/Paris.cc
+++ b/NPSimulation/Paris/Paris.cc
@@ -528,7 +528,7 @@ void Paris::ReadSensitive(const G4Event* event){
 
     G4double* Info = *(LaBr3_itr->second);
 
-    double Energy = Info[0];
+    double Energy = RandGauss::shoot(Info[0],ResoFirstStage);
 
     if(Energy>EnergyThreshold){
       double Time = Info[1];
@@ -555,7 +555,7 @@ void Paris::ReadSensitive(const G4Event* event){
 
     G4double* Info = *(NaI_itr->second);
 
-    double Energy = Info[0];
+    double Energy = RandGauss::shoot(Info[0],ResoSecondStage);
 
     if(Energy>EnergyThreshold){
       double Time = Info[1];