HcalTDC.cc

Go to the documentation of this file.

#include "SimCalorimetry/HcalSimAlgos/interface/HcalTDC.h"
#include "CalibFormats/HcalObjects/interface/HcalCalibrations.h"
#include "CalibFormats/HcalObjects/interface/HcalDbService.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CLHEP/Random/RandGaussQ.h"

HcalTDC::HcalTDC(double threshold_currentTDC)
    : theTDCParameters(), theDbService(nullptr), threshold_currentTDC_(threshold_currentTDC), lsb(3.74) {}

HcalTDC::~HcalTDC() {}

//template <class Digi>
void HcalTDC::timing(const CaloSamples& lf, QIE11DataFrame& digi) const {
  std::vector<int> packedTDCs = leadingEdgeTDC(lf);

  for (int ibin = 0; ibin < lf.size(); ++ibin) {
    digi.setSample(ibin, digi[ibin].adc(), packedTDCs[ibin], digi[ibin].soi());

  }  // loop over bunch crossing bins
}

void HcalTDC::timing(const CaloSamples& lf, QIE10DataFrame& digi) const {
  std::vector<int> packedTDCs = leadingEdgeTDC(lf);

  for (int ibin = 0; ibin < lf.size(); ++ibin) {
    QIE10DataFrame::Sample sam = digi[ibin];
    digi.setSample(ibin, sam.adc(), packedTDCs[ibin] /*LE TDC*/, 0. /*TE TDC*/, sam.capid(), sam.soi(), sam.ok());

  }  // loop over bunch crossing bins
}

//leading edge TDC
std::vector<int> HcalTDC::leadingEdgeTDC(const CaloSamples& lf) const {
  std::vector<int> result;

  float const TDC_Threshold(getThreshold());
  bool risingReady(true);
  int tdcBins = theTDCParameters.nbins();
  bool hasTDCValues = true;
  if (lf.preciseSize() == 0)
    hasTDCValues = false;

  for (int ibin = 0; ibin < lf.size(); ++ibin) {
    /*
    If in a given 25ns bunch/time sample, the pulse is above
    TDC_Thresh, then TDC_RisingEdge set to time when threshold
    was crossed.
    TDC_RisingEdge=0 if it was low in the last precision bin 
    on the previous bunch crossing, but above by first precision
    bin in current bunch crossing.
    TDC_RisingEdge=62 if pulse starts above threshold by end of
    previous bunch crossing and stays above threshold in current 
    bunch crossing. 
    TDC_RisingEdge=63 if the pulse never crosses the threshold.
    */
    // special codes
    int TDC_RisingEdge = (risingReady) ? theTDCParameters.noTransitionCode() : theTDCParameters.alreadyTransitionCode();
    int preciseBegin = ibin * tdcBins;
    int preciseEnd = preciseBegin + tdcBins;

    if (hasTDCValues) {
      for (int i = preciseBegin; i < preciseEnd; ++i) {
        if ((!risingReady) && (i == preciseBegin) && (i != 0)) {
          if (lf.preciseAt(i) / theTDCParameters.deltaT() > TDC_Threshold) {
            TDC_RisingEdge = theTDCParameters.alreadyTransitionCode();
          } else {
            risingReady = true;
          }
        }

        if (risingReady) {
          if (lf.preciseAt(i) / theTDCParameters.deltaT() > TDC_Threshold) {
            TDC_RisingEdge = i - preciseBegin;
            risingReady = false;
          } else if (i == (preciseEnd - 1)) {
            TDC_RisingEdge = theTDCParameters.noTransitionCode();
          }
        }

        if ((!risingReady) && (i == (preciseEnd - 1))) {
          if (lf.preciseAt(i) / theTDCParameters.deltaT() < TDC_Threshold) {
            risingReady = true;
          }
        }

      }  //end of looping precise bins
    }

    // change packing to allow for special codes
    int packedTDC = TDC_RisingEdge;

    result.push_back(packedTDC);

  }  // loop over bunch crossing bins

  return result;
}

void HcalTDC::setDbService(const HcalDbService* service) { theDbService = service; }