HGCalVFELinearizationImpl.cc

Go to the documentation of this file.

#include "L1Trigger/L1THGCal/interface/veryfrontend/HGCalVFELinearizationImpl.h"

HGCalVFELinearizationImpl::
HGCalVFELinearizationImpl(const edm::ParameterSet& conf):
    linLSB_(conf.getParameter<double>("linLSB")),
    adcsaturation_(conf.getParameter<double>("adcsaturation")),
    tdcnBits_(conf.getParameter<uint32_t>("tdcnBits")),
    tdcOnsetfC_(conf.getParameter<double>("tdcOnsetfC")),
    adcnBits_(conf.getParameter<uint32_t>("adcnBits")),
    tdcsaturation_(conf.getParameter<double>("tdcsaturation")),
    linnBits_(conf.getParameter<uint32_t>("linnBits"))
{
  adcLSB_ = adcsaturation_/pow(2.,adcnBits_);
  tdcLSB_ =  tdcsaturation_/pow(2.,tdcnBits_);
  linMax_ = (0x1<<linnBits_)-1;
}

void
HGCalVFELinearizationImpl::
linearize(const std::vector<HGCDataFrame<DetId,HGCSample>>& dataframes,
        std::vector<std::pair<DetId, uint32_t > >& linearized_dataframes)
{
    double amplitude = 0.; 
    uint32_t amplitude_int = 0;
    const int kIntimeSample = 2;

    for(const auto& frame : dataframes) {//loop on DIGI
        unsigned det = frame.id().det();
        if(det==DetId::Forward || det==DetId::HGCalEE || det==DetId::HGCalHSi) {
            if (frame[kIntimeSample].mode()) {//TOT mode
                amplitude =( floor(tdcOnsetfC_/adcLSB_) + 1.0 )* adcLSB_ + double(frame[kIntimeSample].data()) * tdcLSB_;
            }
            else {//ADC mode
                amplitude = double(frame[kIntimeSample].data()) * adcLSB_;
            }

            amplitude_int = uint32_t (floor(amplitude/linLSB_+0.5)); 
        }
        else if(det==DetId::Hcal || det==DetId::HGCalHSc) {
            // no linearization here. Take the raw ADC data
            amplitude_int = frame[kIntimeSample].data();
        }
        if (amplitude_int>linMax_) amplitude_int = linMax_;

        linearized_dataframes.push_back(std::make_pair (frame.id(), amplitude_int));
    }
}