Changeset 4b9a2dc in git for modules/SimpleCalorimeter.cc

Timestamp:

Sep 24, 2014, 10:08:15 AM (10 years ago)

Author:

Michele <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

1c8d9db

Parents:

fcc959c

Message:

added EnergyMin and SigmaMin parameters

File:

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

modules/SimpleCalorimeter.cc

@@ -129 +129 @@
   }
 */
+
+  // read min E value for towers to be saved
+  fEnergyMin = GetDouble("TowerMinEnergy", 0.0); 
+  fSigmaMin  = GetDouble("TowerMinSignificance", 0.0); 
+ 
   // read resolution formulas
   fResolutionFormula->Compile(GetString("ResolutionFormula", "0"));
@@ -384 +389 @@
   time = (fTowerWeightTime < 1.0E-09 ) ? 0 : fTowerTime/fTowerWeightTime;
 
+  sigma = fResolutionFormula->Eval(0.0, fTowerEta, 0.0, energy);
+  
+  energy = (energy < fEnergyMin || energy < fSigmaMin*sigma) ? 0 : energy;
+  
   eta = gRandom->Uniform(fTowerEdges[0], fTowerEdges[1]);
   phi = gRandom->Uniform(fTowerEdges[2], fTowerEdges[3]);
@@ -405 +414 @@
   // fill energy flow candidates
   energy -= fTrackEnergy;
-  if(energy < 0.0) energy = 0.0;
+  if(energy < fEnergyMin || energy < fSigmaMin*fResolutionFormula->Eval(0.0, fTowerEta, 0.0, energy)) energy = 0.0;
    
   // save energy excess as an energy flow tower