Fork me on GitHub

Changeset 49234af in git for modules/SimpleCalorimeter.cc


Ignore:
Timestamp:
Dec 9, 2014, 1:27:13 PM (10 years ago)
Author:
Michele <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
37deb3b, 9e991f8
Parents:
f6b6ee7 (diff), e7e90df (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.
Message:

Merge branch 'TestFastJet310b1'

File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/SimpleCalorimeter.cc

    rf6b6ee7 r49234af  
    147147  }
    148148*/
     149
     150  // read min E value for towers to be saved
     151  fEnergyMin = GetDouble("TowerMinEnergy", 0.0);
     152  fSigmaMin  = GetDouble("TowerMinSignificance", 0.0);
     153 
    149154  // read resolution formulas
    150155  fResolutionFormula->Compile(GetString("ResolutionFormula", "0"));
     
    402407  time = (fTowerWeightTime < 1.0E-09 ) ? 0 : fTowerTime/fTowerWeightTime;
    403408
     409  sigma = fResolutionFormula->Eval(0.0, fTowerEta, 0.0, energy);
     410 
     411  energy = (energy < fEnergyMin || energy < fSigmaMin*sigma) ? 0 : energy;
     412 
    404413  eta = gRandom->Uniform(fTowerEdges[0], fTowerEdges[1]);
    405414  phi = gRandom->Uniform(fTowerEdges[2], fTowerEdges[3]);
     
    423432  // fill energy flow candidates
    424433  energy -= fTrackEnergy;
    425   if(energy < 0.0) energy = 0.0;
     434  if(energy < fEnergyMin || energy < fSigmaMin*fResolutionFormula->Eval(0.0, fTowerEta, 0.0, energy)) energy = 0.0;
    426435   
    427436  // save energy excess as an energy flow tower
Note: See TracChangeset for help on using the changeset viewer.