Fork me on GitHub

Changeset 99437cb in git


Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/Calorimeter.cc

    rce1cbf1 r99437cb  
    388388  ecalSigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy);
    389389
    390   ecalEnergy = gRandom->Gaus(fTowerECalEnergy, ecalSigma);
    391   if(ecalEnergy < 0.0) ecalEnergy = 0.0;
    392 
    393 //  ecalEnergy = LogNormal(fTowerECalEnergy, ecalSigma);
     390//  ecalEnergy = gRandom->Gaus(fTowerECalEnergy, ecalSigma);
     391//  if(ecalEnergy < 0.0) ecalEnergy = 0.0;
     392
     393  ecalEnergy = LogNormal(fTowerECalEnergy, ecalSigma);
    394394
    395395  hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy);
    396396
    397   hcalEnergy = gRandom->Gaus(fTowerHCalEnergy, hcalSigma);
    398   if(hcalEnergy < 0.0) hcalEnergy = 0.0;
    399 
    400 //  hcalEnergy = LogNormal(fTowerHCalEnergy, hcalSigma);
     397//  hcalEnergy = gRandom->Gaus(fTowerHCalEnergy, hcalSigma);
     398//  if(hcalEnergy < 0.0) hcalEnergy = 0.0;
     399
     400  hcalEnergy = LogNormal(fTowerHCalEnergy, hcalSigma);
    401401
    402402  energy = ecalEnergy + hcalEnergy;
Note: See TracChangeset for help on using the changeset viewer.