Changes in / [bba5aca:b25184c] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/Calorimeter.cc
rbba5aca rb25184c 212 212 Double_t ecalFraction, hcalFraction; 213 213 Double_t ecalEnergy, hcalEnergy; 214 Double_t ecalSigma, hcalSigma;215 Double_t energyGuess;216 214 Int_t pdgCode; 217 215 … … 400 398 { 401 399 fECalTrackEnergy += ecalEnergy; 402 ecalSigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, momentum.E()); 403 if(ecalSigma/momentum.E() < track->TrackResolution) energyGuess = ecalEnergy; 404 else energyGuess = momentum.E(); 405 406 fECalTrackSigma += (track->TrackResolution)*energyGuess*(track->TrackResolution)*energyGuess; 400 fECalTrackSigma += (track->TrackResolution)*momentum.E()*(track->TrackResolution)*momentum.E(); 407 401 fECalTowerTrackArray->Add(track); 408 402 } … … 411 405 { 412 406 fHCalTrackEnergy += hcalEnergy; 413 hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, momentum.E()); 414 if(hcalSigma/momentum.E() < track->TrackResolution) energyGuess = hcalEnergy; 415 else energyGuess = momentum.E(); 416 417 fHCalTrackSigma += (track->TrackResolution)*energyGuess*(track->TrackResolution)*energyGuess; 407 fHCalTrackSigma += (track->TrackResolution)*momentum.E()*(track->TrackResolution)*momentum.E(); 418 408 fHCalTowerTrackArray->Add(track); 419 409 }
Note:
See TracChangeset
for help on using the changeset viewer.