HFPreRecAlgo Class Reference

#include <HFPreRecAlgo.h>

Public Member Functions
	HFPreRecAlgo (const bool sumAllTS)
HFQIE10Info	reconstruct (const QIE10DataFrame &digi, int tsToUse, const HcalCoder &coder, const HcalChannelProperties &prop) const
	~HFPreRecAlgo ()

Private Attributes
bool	sumAllTS_

Detailed Description

Definition at line 10 of file HFPreRecAlgo.h.

Constructor & Destructor Documentation

HFPreRecAlgo()

HFPreRecAlgo::HFPreRecAlgo ( const bool  sumAllTS )

inlineexplicit

Definition at line 12 of file HFPreRecAlgo.h.

: sumAllTS_(sumAllTS) {}

~HFPreRecAlgo()

HFPreRecAlgo::~HFPreRecAlgo ( )

inline

Definition at line 14 of file HFPreRecAlgo.h.

{}

Member Function Documentation

reconstruct()

HFQIE10Info HFPreRecAlgo::reconstruct	(	const QIE10DataFrame &	digi,
		int	tsToUse,
		const HcalCoder &	coder,
		const HcalChannelProperties &	prop
	)		const

Definition at line 12 of file HFPreRecAlgo.cc.

                                                                                     {
  // 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;
}

References HcalCoder::adc2fC(), ALCARECOTkAlJpsiMuMu_cff::charge, fwrapper::cs, HCALHighEnergyHPDFilter_cfi::energy, HcalSpecialTimes::getTDCTime(), mps_fire::i, QIE10DataFrame::id(), min(), HFQIE10Info::N_RAW_MAX, HcalChannelProperties::pedsAndGains, submitPVResolutionJobs::q, mps_fire::result, alignCSCRings::s, edm::shift, sumAllTS_, and HcalSpecialTimes::UNKNOWN_T_NOTDC.

Referenced by HFPreReconstructor::fillInfos().

Member Data Documentation

sumAllTS_

bool HFPreRecAlgo::sumAllTS_

private

Definition at line 22 of file HFPreRecAlgo.h.

Referenced by reconstruct().