From 166b97c6c847185528a4d2ea8c63f1415c8804e7 Mon Sep 17 00:00:00 2001
From: adrien-matta <a.matta@surrey.ac.uk>
Date: Mon, 15 Dec 2014 16:03:30 +0000
Subject: [PATCH] * adding resolution to Paris detector
---
Inputs/EventGenerator/gamma.source | 4 ++--
NPSimulation/Paris/Paris.cc | 4 ++--
2 files changed, 4 insertions(+), 4 deletions(-)
diff --git a/Inputs/EventGenerator/gamma.source b/Inputs/EventGenerator/gamma.source
index 877c06c32..2e6472d36 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 6b61ffcd8..00705f421 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];
--
GitLab