Changeset 7da1826 in git for modules/SimpleCalorimeter.cc

Timestamp:

Aug 26, 2016, 2:04:18 PM (8 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

4122b0e, ebf057e

Parents:

46fb86c

Message:

apply better energy estimate also for simple calo

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);
-      
       }