HcalQIEShape.cc

Go to the documentation of this file.

#include "CondFormats/HcalObjects/interface/HcalQIEShape.h"

HcalQIEShape::HcalQIEShape() : nbins_(0) {}

HcalQIEShape::~HcalQIEShape() {}

void HcalQIEShape::expand() {
  int scale = 1;
  for (unsigned range = 1; range < 4; range++) {
    int factor = nbins_ == 32 ? 5 : 8;  // QIE8/QIE10 -> 5/8
    scale *= factor;
    unsigned index = range * nbins_;
    unsigned overlap = (nbins_ == 32) ? 2 : 3;  // QIE10 -> 3 bin overlap
    mValues[index] = mValues[index - overlap];  // link to previous range
    for (unsigned i = 1; i < nbins_; i++) {
      mValues[index + i] = mValues[index + i - 1] + scale * (mValues[i] - mValues[i - 1]);
    }
  }
  mValues[nbins_ * 4] = 2 * mValues[nbins_ * 4 - 1] - mValues[nbins_ * 4 - 2];  // extrapolate
}

float HcalQIEShape::lowEdge(unsigned fAdc) const {
  if (fAdc < 4 * nbins_)
    return mValues[fAdc];
  return 0.;
}

float HcalQIEShape::center(unsigned fAdc) const {
  if (fAdc < 4 * nbins_) {
    if (fAdc % nbins_ == nbins_ - 1)
      return 0.5 * (3 * mValues[fAdc] - mValues[fAdc - 1]);  // extrapolate
    else
      return 0.5 * (mValues[fAdc] + mValues[fAdc + 1]);  // interpolate
  }
  return 0.;
}

float HcalQIEShape::highEdge(unsigned fAdc) const {
  if (fAdc < 4 * nbins_)
    return mValues[fAdc + 1];
  return 0.;
}

bool HcalQIEShape::setLowEdge(float fValue, unsigned fAdc) {
  if (fAdc >= nbins_)
    return false;
  mValues[fAdc] = fValue;
  return true;
}

bool HcalQIEShape::setLowEdges(unsigned int nbins, const float *fValue) {
  nbins_ = nbins;
  mValues.clear();
  mValues.resize(4 * nbins_ + 1);
  bool result = true;
  for (unsigned int adc = 0; adc < nbins_; adc++)
    result = result && setLowEdge(fValue[adc], adc);
  expand();
  return result;
}