HcalZSAlgoEnergy.cc

Go to the documentation of this file.

#include "CondFormats/HcalObjects/interface/HcalPedestal.h"
#include "CondFormats/HcalObjects/interface/HcalQIECoder.h"
#include "CondFormats/HcalObjects/interface/HcalQIEShape.h"
#include "HcalZSAlgoEnergy.h"
#include <iostream>

HcalZSAlgoEnergy::HcalZSAlgoEnergy(bool mp, int level, int start, int samples, bool twosided)
    : HcalZeroSuppressionAlgo(mp), threshold_(level), firstsample_(start), samplecount_(samples), twosided_(twosided) {}

namespace ZSEnergy_impl {

  template <class DIGI>
  bool keepMe(const HcalDbService &db, const DIGI &inp, int threshold, int firstSample, int samples, bool twosided) {
    bool keepIt = false;
    const HcalQIECoder *channelCoder = db.getHcalCoder(inp.id());
    const HcalQIEShape *shape = db.getHcalShape(channelCoder);

    // determine average pedestal for channel
    float pedsum = 0, pedave = 0, offset = 0, slope = 0;
    const HcalPedestal *ped = db.getPedestal(inp.id());
    for (int j = 0; j < 4; j++) {
      pedave += ped->getValue(j) / 4.0;
      offset += channelCoder->charge(*shape, 0, j) / 4.0;
      slope += channelCoder->charge(*shape, 1, j) / 4.0;
    }
    slope -= offset;
    pedave -= offset;
    pedave /= slope;

    int sum = 0;

    for (int j = 0; j < samples && j + firstSample < (int)inp.size(); j++) {
      sum += inp[j + firstSample].adc();
      pedsum += pedave;
    }
    //    int presum=sum;
    sum -= (int)(pedsum);

    if (sum >= threshold)
      keepIt = true;
    else if (sum <= (-threshold) && twosided)
      keepIt = true;
    /*
    else
   std::cout << inp.id() << " " << sum << ":" << presum << " " << threshold
    << " " << pedsum << " " << pedave
    << " " << offset << " " << slope
    << std::endl;
  */
    return keepIt;
  }
}  // namespace ZSEnergy_impl

bool HcalZSAlgoEnergy::shouldKeep(const HBHEDataFrame &digi) const {
  return ZSEnergy_impl::keepMe<HBHEDataFrame>(*db_, digi, threshold_, firstsample_, samplecount_, twosided_);
}
bool HcalZSAlgoEnergy::shouldKeep(const HODataFrame &digi) const {
  return ZSEnergy_impl::keepMe<HODataFrame>(*db_, digi, threshold_, firstsample_, samplecount_, twosided_);
}
bool HcalZSAlgoEnergy::shouldKeep(const HFDataFrame &digi) const {
  return ZSEnergy_impl::keepMe<HFDataFrame>(*db_, digi, threshold_, firstsample_, samplecount_, twosided_);
}
bool HcalZSAlgoEnergy::shouldKeep(const QIE10DataFrame &digi) const {
  return ZSEnergy_impl::keepMe<QIE10DataFrame>(*db_, digi, threshold_, firstsample_, samplecount_, twosided_);
}
bool HcalZSAlgoEnergy::shouldKeep(const QIE11DataFrame &digi) const {
  return ZSEnergy_impl::keepMe<QIE11DataFrame>(*db_, digi, threshold_, firstsample_, samplecount_, twosided_);
}

void HcalZSAlgoEnergy::prepare(const HcalDbService *db) { db_ = db; }
void HcalZSAlgoEnergy::done() { db_ = nullptr; }