HcalSiPMShape.cc

Go to the documentation of this file.

#include "CalibCalorimetry/HcalAlgos/interface/HcalPulseShapes.h"
#include "SimCalorimetry/HcalSimAlgos/interface/HcalSiPMShape.h"
#include "SimCalorimetry/HcalSimAlgos/interface/HcalShapes.h"
#include <iostream>

HcalSiPMShape::HcalSiPMShape(unsigned int signalShape) : CaloVShape(), nBins_(HcalPulseShapes::nBinsSiPM_*HcalPulseShapes::invDeltaTSiPM_), nt_(nBins_, 0.) {
  computeShape(signalShape);
}

HcalSiPMShape::HcalSiPMShape(const HcalSiPMShape & other) : CaloVShape(other), nBins_(other.nBins_), nt_(other.nt_) {}

double HcalSiPMShape::operator () (double time) const {
  int jtime(time*HcalPulseShapes::invDeltaTSiPM_ + 0.5);
  if (jtime>=0 && jtime<nBins_) 
    return nt_[jtime];
  return 0.;
}

void HcalSiPMShape::computeShape(unsigned int signalShape) {
  //grab correct function pointer based on shape
  double (*analyticPulseShape)(double);
  if(signalShape==HcalShapes::ZECOTEK || signalShape==HcalShapes::HAMAMATSU) analyticPulseShape = &HcalPulseShapes::analyticPulseShapeSiPMHO;
  else if(signalShape==HcalShapes::HE2017 or signalShape==HcalShapes::HE2018) analyticPulseShape = &HcalPulseShapes::analyticPulseShapeSiPMHE;
  else return;

  double norm = 0.;
  for (int j = 0; j < nBins_; ++j) {
    nt_[j] = analyticPulseShape(j*HcalPulseShapes::deltaTSiPM_);
    norm += (nt_[j]>0) ? nt_[j] : 0.;
  }

  for (int j = 0; j < nBins_; ++j) {
    nt_[j] /= norm;
  }
}