Fork me on GitHub

Changeset f85a40e in git


Ignore:
Timestamp:
Sep 4, 2013, 3:07:05 PM (11 years ago)
Author:
pavel <pavel@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
a1ac20e
Parents:
852bc18
Message:

remove sigma cut and switch to gaussian distribution for calorimeter smearing

File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/Calorimeter.cc

    r852bc18 rf85a40e  
    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;
     
    441441
    442442  ecalEnergy -= fTrackECalEnergy;
    443   if(ecalEnergy < 0.0 || ecalEnergy < ecalSigma) ecalEnergy = 0.0;
     443  if(ecalEnergy < 0.0) ecalEnergy = 0.0;
    444444
    445445  hcalEnergy -= fTrackHCalEnergy;
    446   if(hcalEnergy < 0.0 || hcalEnergy < hcalSigma) hcalEnergy = 0.0;
     446  if(hcalEnergy < 0.0) hcalEnergy = 0.0;
    447447
    448448  energy = ecalEnergy + hcalEnergy;
Note: See TracChangeset for help on using the changeset viewer.