Fork me on GitHub

Changes in / [bba5aca:b25184c] in git


Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/Calorimeter.cc

    rbba5aca rb25184c  
    212212  Double_t ecalFraction, hcalFraction;
    213213  Double_t ecalEnergy, hcalEnergy;
    214   Double_t ecalSigma, hcalSigma;
    215   Double_t energyGuess;
    216214  Int_t pdgCode;
    217215
     
    400398      {
    401399        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();
    407401        fECalTowerTrackArray->Add(track);
    408402      }
     
    411405      {
    412406        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();
    418408        fHCalTowerTrackArray->Add(track);
    419409      }
Note: See TracChangeset for help on using the changeset viewer.