- Timestamp:
- Sep 4, 2013, 3:07:05 PM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/modules/Calorimeter.cc
r1276 r1277 388 388 ecalSigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy); 389 389 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); 394 394 395 395 hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy); 396 396 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); 401 401 402 402 energy = ecalEnergy + hcalEnergy; … … 441 441 442 442 ecalEnergy -= fTrackECalEnergy; 443 if(ecalEnergy < 0.0 || ecalEnergy < ecalSigma) ecalEnergy = 0.0;443 if(ecalEnergy < 0.0) ecalEnergy = 0.0; 444 444 445 445 hcalEnergy -= fTrackHCalEnergy; 446 if(hcalEnergy < 0.0 || hcalEnergy < hcalSigma) hcalEnergy = 0.0;446 if(hcalEnergy < 0.0) hcalEnergy = 0.0; 447 447 448 448 energy = ecalEnergy + hcalEnergy;
Note:
See TracChangeset
for help on using the changeset viewer.