EcalTBH2TDCRecInfoAlgo.cc

Go to the documentation of this file.

#include "RecoTBCalo/EcalTBTDCReconstructor/interface/EcalTBH2TDCRecInfoAlgo.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <list>
#include <vector>
#include <iostream>

EcalTBH2TDCRecInfoAlgo::EcalTBH2TDCRecInfoAlgo(const std::vector<EcalTBH2TDCRanges>& tdcZeros)
    : tdcZeros_(tdcZeros), actualRun_(-1), actualRange_(-1) {}

EcalTBTDCRecInfo EcalTBH2TDCRecInfoAlgo::reconstruct(const int& runNumber, const HcalTBTiming& TDCRawInfo) const {
  if (actualRun_ != runNumber) {
    actualRun_ = runNumber;
    actualRange_ = -1;
    for (unsigned int i = 0; i < tdcZeros_.size(); i++)
      if (runNumber <= tdcZeros_[i].runRanges.second && runNumber >= tdcZeros_[i].runRanges.first)
        actualRange_ = i;

    if (actualRange_ == -1) {
      edm::LogError("TDCRange not found") << "TDC range not found";
      return EcalTBTDCRecInfo(-1);
    }
  }

  int ntdc = TDCRawInfo.BeamCoincidenceCount();
  edm::LogInfo("") << "EcalTBH2TDCRecInfoAlgo::reconstruct # tdc hits: " << ntdc << std::endl;
  ;
  if (ntdc > 1) {
    for (int i = 0; i < ntdc; ++i) {
      edm::LogInfo("") << "hit i: " << i << " tdc: " << TDCRawInfo.BeamCoincidenceHits(i) << std::endl;
    }
  }

  if (ntdc == 0) {
    edm::LogError("NoTDCHits") << "no TDC hits. TDC info not reliable" << std::endl;
    return EcalTBTDCRecInfo(-999.);
  }

  //double tdcd = TDCRawInfo.ttcL1Atime() - TDCRawInfo.beamCoincidence();
  double tdcd = TDCRawInfo.ttcL1Atime() - TDCRawInfo.BeamCoincidenceHits(0);

  if (  
      (tdcd < tdcZeros_[actualRange_].tdcZero - 1 || tdcd > tdcZeros_[actualRange_].tdcZero + 26)) {
    edm::LogError("TDCOutOfRange") << " ============================\n"
                                   << " tdc value out of range = " << tdcd
                                   << " tdcZero = " << tdcZeros_[actualRange_].tdcZero << "\n"
                                   << " ============================\n"
                                   << std::endl;
    tdcRangeErrorMessageAlreadyDisplayed_ = true;
    return EcalTBTDCRecInfo(-999.);
  }

  double offset = ((double)tdcd - (double)tdcZeros_[actualRange_].tdcZero) / 25.;  //
  //   if (use2004OffsetConvention)
  //  offset = (1. - offset) ;
  return EcalTBTDCRecInfo(offset);
}