PedestalSub.cc

Go to the documentation of this file.

#include <iostream>
#include <cmath>
#include <climits>
#include "RecoLocalCalo/HcalRecAlgos/interface/PedestalSub.h"

using namespace std;

PedestalSub::PedestalSub() : fThreshold(2.7),fQuantile(0.0),fCondition(0){
}

PedestalSub::~PedestalSub() { 
}

void PedestalSub::init(int runCond=0, float threshold=0.0, float quantile=0.0) {
  fThreshold=threshold;
  fQuantile=quantile;
  fCondition=runCond;
}

void PedestalSub::calculate(const std::vector<double> & inputCharge, const std::vector<double> & inputPedestal, const std::vector<double> & inputNoise, std::vector<double> & corrCharge, int soi, int nTS) const {

  double bseCorr=PedestalSub::getCorrection(inputCharge, inputPedestal, inputNoise, soi, nTS);
  for (auto i=0; i<nTS; i++) {
    corrCharge.push_back(inputCharge[i]-inputPedestal[i]-bseCorr);
  }
}

double PedestalSub::getCorrection(const std::vector<double> & inputCharge, const std::vector<double> & inputPedestal, const std::vector<double> & inputNoise, int soi, int nTS) const {

  double baseline=0;

    for (auto i=0; i<nTS; i++) {
      if (i==soi || i==(soi+1)) continue;
      if ( (inputCharge[i]-inputPedestal[i])<3*inputNoise[i]) {
    baseline+=(inputCharge[i]-inputPedestal[i]);
      }
      else {
    baseline+=3*inputNoise[i];
      }
    }
    baseline/=(nTS-2);
  return baseline;
  
}