- Timestamp:
- Sep 4, 2013, 8:26:45 AM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/modules/Calorimeter.cc
r1274 r1275 382 382 Double_t energy, pt, eta, phi; 383 383 Double_t ecalEnergy, hcalEnergy; 384 Double_t ecalSigma, hcalSigma; 384 385 385 386 if(!fTower) return; 386 387 387 // ecalEnergy = gRandom->Gaus(fTowerECalEnergy, fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy)); 388 ecalSigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy); 389 390 // ecalEnergy = gRandom->Gaus(fTowerECalEnergy, ecalSigma); 388 391 // if(ecalEnergy < 0.0) ecalEnergy = 0.0; 389 392 390 ecalEnergy = LogNormal(fTowerECalEnergy, fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy)); 391 392 // hcalEnergy = gRandom->Gaus(fTowerHCalEnergy, fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy)); 393 ecalEnergy = LogNormal(fTowerECalEnergy, ecalSigma); 394 395 hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy); 396 397 // hcalEnergy = gRandom->Gaus(fTowerHCalEnergy, hcalSigma); 393 398 // if(hcalEnergy < 0.0) hcalEnergy = 0.0; 394 399 395 hcalEnergy = LogNormal(fTowerHCalEnergy, fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy));400 hcalEnergy = LogNormal(fTowerHCalEnergy, hcalSigma); 396 401 397 402 energy = ecalEnergy + hcalEnergy; … … 436 441 437 442 ecalEnergy -= fTrackECalEnergy; 438 if(ecalEnergy < 0.0 ) ecalEnergy = 0.0;443 if(ecalEnergy < 0.0 || ecalEnergy < ecalSigma) ecalEnergy = 0.0; 439 444 440 445 hcalEnergy -= fTrackHCalEnergy; 441 if(hcalEnergy < 0.0 ) hcalEnergy = 0.0;446 if(hcalEnergy < 0.0 || hcalEnergy < hcalSigma) hcalEnergy = 0.0; 442 447 443 448 energy = ecalEnergy + hcalEnergy;
Note:
See TracChangeset
for help on using the changeset viewer.