HcalTDCReco Class Reference

#include <HcalTDCReco.h>

Public Member Functions
	HcalTDCReco ()
void	reconstruct (const HcalUpgradeDataFrame &digi, HBHERecHit &recHit) const

Detailed Description

Definition at line 7 of file HcalTDCReco.h.

Constructor & Destructor Documentation

HcalTDCReco::HcalTDCReco

(

)

Definition at line 5 of file HcalTDCReco.cc.

{}

Member Function Documentation

void HcalTDCReco::reconstruct	(	const HcalUpgradeDataFrame &	digi,
		HBHERecHit &	recHit
	)		const

Definition at line 8 of file HcalTDCReco.cc.

References CaloRecHit::energy(), i, HBHERecHit::id(), gen::n, pileupCalc::nbins, HcalUpgradeDataFrame::size(), and HcalUpgradeDataFrame::tdc().

{
  int n = digi.size();
  double risingTime = -999.;
  double fallingTime = -999.;
  int signalBX = 4;              // NB: HARDWIRED !!!
  int nbins = 50; // as in HcalTDCParameters.h (SimCalorimetry/HcalSimAlgos)
  int direction(-1), stepSize(1); // where to go after signalBX
  // start at bunch crossing 3 by default
  int i(signalBX);

  // for(int i = 3; i < n; ++i)
  while ((i > 2) && (i < 8) && (i < n) && 
     ((risingTime < -998.) || (fallingTime < 998.)))
  {
    unsigned tdc = digi.tdc(i);

    /* 
    unsigned rising = tdc & 0x7F;
    unsigned falling = (tdc >> 7) & 0x7F;
    */
    // temporary "unpacking" instead, which directly corresponds to 
    // SimCalorimetry/HcalSimAlgos/src/HcalTDC.cc  (nibs = 50...)
    // packedTDC = TDC_RisingEdge + (tdcBins*2) * TDC_FallingEdge;
    unsigned rising  =  tdc%100; 
    unsigned falling =  tdc/100;

    // only set the first time, avoiding "special" codes 
    if(risingTime < -998. && rising != 62 && rising != 63) {
      risingTime = rising*25./nbins + (i-signalBX)*25.;
    }
    if(((fallingTime < -998.) || (fallingTime < risingTime)) && 
       (falling != 62) && (falling != 63)) {
      fallingTime = falling*25./nbins + (i-signalBX)*25.;
    }  

    i += direction*stepSize;
    ++stepSize;
    direction *= -1;
    /*
    std::cout << " digi.tdc[" << i << "] = " << tdc 
          << "  rising = " << rising << "  falling = " << falling
          << "   Rt = " << risingTime
          << "   Ft = " << fallingTime
          << std::endl;
    */
  }
  recHit = HBHERecHit(recHit.id(), recHit.energy(), risingTime, fallingTime);
}