Skip to content
Snippets Groups Projects
Commit 166b97c6 authored by adrien-matta's avatar adrien-matta
Browse files

* adding resolution to Paris detector

parent 2e4e7115
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Inputs/EventGenerator/gamma.source
+ 2
− 2
View file @ 166b97c6
...@@ -4,8 +4,8 @@ ...@@ -4,8 +4,8 @@
% Energy are given in MeV , Position in mm % % Energy are given in MeV , Position in mm %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Isotropic Isotropic
EnergyLow= .1 EnergyLow= 1
EnergyHigh= 5. EnergyHigh= 1.
HalfOpenAngleMin= 0 HalfOpenAngleMin= 0
HalfOpenAngleMax= 180 HalfOpenAngleMax= 180
x0= 0 x0= 0
......
NPSimulation/Paris/Paris.cc
+ 2
− 2
View file @ 166b97c6
...@@ -528,7 +528,7 @@ void Paris::ReadSensitive(const G4Event* event){ ...@@ -528,7 +528,7 @@ void Paris::ReadSensitive(const G4Event* event){
G4double* Info = *(LaBr3_itr->second); G4double* Info = *(LaBr3_itr->second);
double Energy = Info[0]; double Energy = RandGauss::shoot(Info[0],ResoFirstStage);
if(Energy>EnergyThreshold){ if(Energy>EnergyThreshold){
double Time = Info[1]; double Time = Info[1];
...@@ -555,7 +555,7 @@ void Paris::ReadSensitive(const G4Event* event){ ...@@ -555,7 +555,7 @@ void Paris::ReadSensitive(const G4Event* event){
G4double* Info = *(NaI_itr->second); G4double* Info = *(NaI_itr->second);
double Energy = Info[0]; double Energy = RandGauss::shoot(Info[0],ResoSecondStage);
if(Energy>EnergyThreshold){ if(Energy>EnergyThreshold){
double Time = Info[1]; double Time = Info[1];
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment