EcalTBTDCRecInfoAlgo.cc

Go to the documentation of this file.

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

#include <list>
#include <iostream>

EcalTBTDCRecInfoAlgo::EcalTBTDCRecInfoAlgo(const std::vector<EcalTBTDCRanges>& tdcRanges) :
  tdcRanges_(tdcRanges), tdcRangeErrorMessageAlreadyDisplayed_(false), actualRun_(-1) , actualRange_(-1)
{
  
}

EcalTBTDCRecInfo EcalTBTDCRecInfoAlgo::reconstruct(const EcalTBTDCRawInfo& TDCRawInfo,const EcalTBEventHeader& eventHeader, bool use2004OffsetConvention) const 
{
  if (actualRun_ != eventHeader.runNumber() )
    {
      actualRun_ = eventHeader.runNumber();
      actualRange_ = -1;
      for (unsigned int i=0; i<tdcRanges_.size();i++)
    if ( eventHeader.runNumber() <= tdcRanges_[i].runRanges.second && eventHeader.runNumber() >= tdcRanges_[i].runRanges.first )
      actualRange_ = i;
      
      if (actualRange_ == -1)
    {
      edm::LogError("TDCRange not found") << "TDC range not found";
      return EcalTBTDCRecInfo(-1);
    }
    }
  
  int eventType;
  eventType=( (eventHeader.dbEventType() == 0) ? 0 : (eventHeader.dbEventType()-1));

  int tdcd = TDCRawInfo[0].tdcValue();

  if( !tdcRangeErrorMessageAlreadyDisplayed_ 
      && (tdcd < tdcRanges_[actualRange_].tdcMin[eventType] || tdcd > tdcRanges_[actualRange_].tdcMax[eventType]) ){
    std::cout << " ============================" <<std::endl;
    std::cout << " Event type: " << eventType << std::endl;
    std::cout << " TDC values not correctly set Raw TDC = "<< tdcd 
          << " (Min/Max = "<< tdcRanges_[actualRange_].tdcMin[eventType] << "/" << tdcRanges_[actualRange_].tdcMax[eventType]
          << std::endl;
    std::cout << " ============================" <<std::endl;
    tdcRangeErrorMessageAlreadyDisplayed_ = true;
  }

  double offset = ( (double)tdcd - (double)tdcRanges_[actualRange_].tdcMin[eventType] )
    / ((double)tdcRanges_[actualRange_].tdcMax[eventType]-(double)tdcRanges_[actualRange_].tdcMin[eventType]);
  if (use2004OffsetConvention)
    offset = (1. - offset) ;
  return EcalTBTDCRecInfo(offset); 
}