Changes in modules/SimpleCalorimeter.cc [e0f8f99:7da1826] in git

File:

• modules/SimpleCalorimeter.cc (modified) (2 diffs)

modules/SimpleCalorimeter.cc

@@ -190 +190 @@
   Double_t fraction;
   Double_t energy;
+  Double_t sigma;
+  Double_t energyGuess;
+
   Int_t pdgCode;
 
@@ -364 +367 @@
        // compute total charged energy   
        fTrackEnergy += energy;
-       fTrackSigma += ((track->TrackResolution)*momentum.E())*((track->TrackResolution)*momentum.E());
-       
+       sigma = fResolutionFormula->Eval(0.0, fTowerEta, 0.0, momentum.E());
+       if(sigma/momentum.E() < track->TrackResolution) energyGuess = energy;
+       else energyGuess = momentum.E();
+
+       fTrackSigma += ((track->TrackResolution)*energyGuess)*((track->TrackResolution)*energyGuess);
        fTowerTrackArray->Add(track);
-      
       }