Skip to content
Snippets Groups Projects
Commit 4447e337 authored by Greg Christian's avatar Greg Christian
Browse files

Added commented code setting energy resolution of Ge detectors

parent bb6b6928
No related branches found
No related tags found
No related merge requests found
......@@ -574,7 +574,7 @@ void GeTAMU::ReadSensitive(const G4Event* event){
G4double* Info = *(HPGE_itr->second);
G4double Energy = Info[0];
G4double Energy = Info[0]; // RandGauss::shoot(Info[0], ResoEnergy/2.334);
G4double Time = Info[1];
G4int CloverNbr = (int)Info[7];
G4int CrystalNbr = (int)Info[8];
......@@ -586,12 +586,6 @@ void GeTAMU::ReadSensitive(const G4Event* event){
else if(CrystalNbr == 1 || CrystalNbr == 4) { SegmentNbr = 1; } // RIGHT
else { SegmentNbr = 3; } // LEFT
// These are just a guess at resolution, based on p. 430 in G.F. Knoll
// G4double FWHM = 0.00083*Energy + 0.0009;
// Energy = RandGauss::shoot(Info[0], FWHM/2.334);
m_GeTAMUData->SetCoreE(CloverNbr, CrystalNbr, Energy/keV);
m_GeTAMUData->SetCoreT(CloverNbr, CrystalNbr, Time/ns);
m_GeTAMUData->SetSegmentE(CloverNbr, SegmentNbr, Energy/keV);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment