Changeset 99437cb in git for modules/Calorimeter.cc

Timestamp:

Sep 4, 2013, 5:20:22 PM (11 years ago)

Author:

pavel <pavel@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

0390d62

Parents:

ce1cbf1

Message:

revert to log-normal distribution for calorimeter smearing

File:

• modules/Calorimeter.cc (modified) (1 diff)

modules/Calorimeter.cc

@@ -388 +388 @@
   ecalSigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy);
 
-  ecalEnergy = gRandom->Gaus(fTowerECalEnergy, ecalSigma);
-  if(ecalEnergy < 0.0) ecalEnergy = 0.0;
-
-//  ecalEnergy = LogNormal(fTowerECalEnergy, ecalSigma);
+//  ecalEnergy = gRandom->Gaus(fTowerECalEnergy, ecalSigma);
+//  if(ecalEnergy < 0.0) ecalEnergy = 0.0;
+
+  ecalEnergy = LogNormal(fTowerECalEnergy, ecalSigma);
 
   hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy);
 
-  hcalEnergy = gRandom->Gaus(fTowerHCalEnergy, hcalSigma);
-  if(hcalEnergy < 0.0) hcalEnergy = 0.0;
-
-//  hcalEnergy = LogNormal(fTowerHCalEnergy, hcalSigma);
+//  hcalEnergy = gRandom->Gaus(fTowerHCalEnergy, hcalSigma);
+//  if(hcalEnergy < 0.0) hcalEnergy = 0.0;
+
+  hcalEnergy = LogNormal(fTowerHCalEnergy, hcalSigma);
 
   energy = ecalEnergy + hcalEnergy;