#include <HGCalVFELinearizationImpl.h>

Public Member Functions
	HGCalVFELinearizationImpl (const edm::ParameterSet &conf)

void	linearize (const std::vector< HGCDataFrame< DetId, HGCSample >> &, std::vector< std::pair< DetId, uint32_t > > &)

uint32_t	linnBits () const

Private Attributes
double	adcLSB_

uint32_t	adcnBits_

double	adcsaturation_

double	linLSB_

uint32_t	linMax_

uint32_t	linnBits_

double	tdcLSB_

uint32_t	tdcnBits_

double	tdcOnsetfC_

double	tdcsaturation_

Detailed Description

Definition at line 12 of file HGCalVFELinearizationImpl.h.

Constructor & Destructor Documentation

HGCalVFELinearizationImpl::HGCalVFELinearizationImpl ( const edm::ParameterSet & conf )

Definition at line 4 of file HGCalVFELinearizationImpl.cc.

References adcLSB_, adcnBits_, adcsaturation_, linearize(), linMax_, linnBits_, funct::pow(), tdcLSB_, tdcnBits_, and tdcsaturation_.

                                                       :
     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;
 }

Member Function Documentation

void HGCalVFELinearizationImpl::linearize	(	const std::vector< HGCDataFrame< DetId, HGCSample >> &	dataframes,
		std::vector< std::pair< DetId, uint32_t > > &	linearized_dataframes
	)

Definition at line 20 of file HGCalVFELinearizationImpl.cc.

References adcLSB_, CustomPhysics_cfi::amplitude, data, DetId::Forward, DetId::Hcal, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, linLSB_, linMax_, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, tdcLSB_, and tdcOnsetfC_.

Referenced by HGCalVFELinearizationImpl(), and HGCalVFEProcessorSums::run().

 {
     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));
     }
 }