HFPreRecAlgo.cc

Go to the documentation of this file.

#include <algorithm>

#include "RecoLocalCalo/HcalRecAlgos/interface/HFPreRecAlgo.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HcalChannelProperties.h"

#include "DataFormats/HcalDigi/interface/QIE10DataFrame.h"
#include "DataFormats/HcalRecHit/interface/HcalSpecialTimes.h"

#include "CalibFormats/HcalObjects/interface/HcalCoder.h"
#include "CalibFormats/CaloObjects/interface/CaloSamples.h"

HFQIE10Info HFPreRecAlgo::reconstruct(const QIE10DataFrame& digi,
                                      const int tsToUse,
                                      const HcalCoder& coder,
                                      const HcalChannelProperties& properties) const {
  // Scrap the trailing edge time for now -- until the front-end
  // FPGA firmware is finalized and the description of the FPGA
  // output becomes available
  static const float timeFalling = HcalSpecialTimes::UNKNOWN_T_NOTDC;

  HFQIE10Info result;

  CaloSamples cs;
  coder.adc2fC(digi, cs);
  const int nRead = cs.size();

  // Number of raw samples to store in HFQIE10Info
  const int nStore = std::min(nRead, static_cast<int>(HFQIE10Info::N_RAW_MAX));

  if (sumAllTS_) {
    // This branch is intended for use with cosmic runs
    double charge = 0.0, energy = 0.0;
    HFQIE10Info::raw_type raw[HFQIE10Info::N_RAW_MAX];

    for (int ts = 0; ts < nRead; ++ts) {
      const QIE10DataFrame::Sample s(digi[ts]);
      const int capid = s.capid();
      const HcalPipelinePedestalAndGain& pAndGain(properties.pedsAndGains.at(capid));
      const float q = cs[ts] - pAndGain.pedestal(false);
      charge += q;
      energy += q * pAndGain.gain();
      if (ts < nStore)
        raw[ts] = s.wideRaw();
    }

    // Timing measurement does not appear to be useful here
    const float timeRising = HcalSpecialTimes::UNKNOWN_T_NOTDC;

    // The following HFQIE10Info arguments correspond to SOI
    // not stored in the raw data. Essentially, only charge
    // and energy are meaningful.
    result = HFQIE10Info(digi.id(), charge, energy, timeRising, timeFalling, raw, nStore, nStore);
  } else if (0 <= tsToUse && tsToUse < nRead) {
    const QIE10DataFrame::Sample s(digi[tsToUse]);
    const int capid = s.capid();
    const HcalPipelinePedestalAndGain& pAndGain(properties.pedsAndGains.at(capid));
    const float charge = cs[tsToUse] - pAndGain.pedestal(false);
    const float energy = charge * pAndGain.gain();
    const float timeRising = HcalSpecialTimes::getTDCTime(s.le_tdc());

    // Figure out the window in the raw data
    // that we want to store. This window will
    // have the width given by "nStore" and
    // will start at "shift".
    int shift = 0;
    if (nRead > static_cast<int>(HFQIE10Info::N_RAW_MAX)) {
      // Try to center the window on "tsToUse"
      const int winCenter = nStore / 2;
      if (tsToUse > winCenter)
        shift = tsToUse - winCenter;
      if (shift + nStore > nRead)
        shift = nRead - nStore;
    }

    // Fill an array of raw values
    HFQIE10Info::raw_type raw[HFQIE10Info::N_RAW_MAX];
    for (int i = 0; i < nStore; ++i)
      raw[i] = digi[i + shift].wideRaw();

    result = HFQIE10Info(digi.id(), charge, energy, timeRising, timeFalling, raw, nStore, tsToUse - shift);
  }
  return result;
}