/data/refman/pasoursint/CMSSW_6_1_2_SLHC4_patch1/src/RecoLocalCalo/HcalRecAlgos/src/HcalTDCReco.cc

Go to the documentation of this file.

#include "RecoLocalCalo/HcalRecAlgos/src/HcalTDCReco.h"
#include "DataFormats/HcalDigi/interface/HcalUpgradeDataFrame.h"
#include "DataFormats/HcalRecHit/interface/HBHERecHit.h"

HcalTDCReco::HcalTDCReco() 
{}

void HcalTDCReco::reconstruct(const HcalUpgradeDataFrame & digi, 
                              HBHERecHit & recHit) const
{
  int n = digi.size();
  double risingTime = -999.;
  double fallingTime = -999.;
  int signalBX = 4;              // NB: HARDWIRED !!!
  // start at bunch crossing 3
  for(int i=3; i < n; ++i)
  {
    unsigned tdc = digi.tdc(i);
    unsigned rising = tdc & 0x7F;
    // altered packing to tdcBins*2=7
    unsigned falling = (tdc >> 7) & 0x7F;
    // only set the first time
    if(risingTime < -998. && rising != 64 && rising != 65) {
      risingTime = rising*25./64. + (i-signalBX)*25.;
    }
    if(fallingTime < -998. && falling != 64 && falling != 65) {
      fallingTime = falling*25./64. + (i-signalBX)*25.;
    }    
  }
  recHit = HBHERecHit(recHit.id(), recHit.energy(), risingTime, fallingTime);
}