#include <CastorChannelCoder.h>

Public Member Functions
int	adc (const reco::castor::QieShape &fShape, double fCharge, int fCapId) const
	fC + capid[0..3] -> ADC conversion More...

	CastorChannelCoder (const float fOffset[16], const float fSlope[16])

double	charge (const reco::castor::QieShape &fShape, int fAdc, int fCapId) const
	ADC[0..127]+capid[0..3]->fC conversion. More...

int	index (int fCapId, int Range)

Private Attributes
double	mOffset [4][4]

double	mSlope [4][4]

Detailed Description

Container for ADC<->fQ conversion constants for HCAL/Castor QIE

Container for ADC<->fQ conversion constants for HCAL/Castor QIE $Original author: ratnikov

Definition at line 16 of file CastorChannelCoder.h.

Constructor & Destructor Documentation

◆ CastorChannelCoder()

CastorChannelCoder::CastorChannelCoder	(	const float	fOffset[16],
		const float	fSlope[16]
	)

Definition at line 12 of file CastorChannelCoder.cc.

                                                                                       {  // [CapId][Range]
   for (int range = 0; range < 4; range++) {
     for (int capId = 0; capId < 4; capId++) {
       mOffset[capId][range] = fOffset[index(capId, range)];
       mSlope[capId][range] = fSlope[index(capId, range)];
     }
   }
 }

References index(), mOffset, mSlope, and FastTimerService_cff::range.

Member Function Documentation

◆ adc()

int CastorChannelCoder::adc	(	const reco::castor::QieShape &	fShape,
		double	fCharge,
		int	fCapId
	)		const

fC + capid[0..3] -> ADC conversion

Definition at line 29 of file CastorChannelCoder.cc.

                                                                                               {
   int adc = -1;  //nothing found yet
   // search for the range
   for (int range = 0; range < 4; range++) {
     double qieCharge = (fCharge - mOffset[fCapId][range]) * mSlope[fCapId][range];
     double qieChargeMax = fShape.linearization(32 * range + 31) + 0.5 * fShape.binSize(32 * range + 31);
     if (range == 3 && qieCharge > qieChargeMax)
       adc = 127;  // overflow
     if (qieCharge > qieChargeMax)
       continue;  // next range
     for (int bin = 32 * range; bin < 32 * (range + 1); bin++) {
       if (qieCharge < fShape.linearization(bin) + 0.5 * fShape.binSize(bin)) {
         adc = bin;
         break;
       }
     }
     if (adc >= 0)
       break;  // found
   }
   if (adc < 0)
     adc = 0;  // underflow
  
   //   std::cout << "CastorChannelCoder::adc-> " << fCharge << '/' << fCapId
   //        << " result: " << adc << std::endl;
   return adc;
 }

References newFWLiteAna::bin, reco::castor::QieShape::binSize(), reco::castor::QieShape::linearization(), mOffset, mSlope, and FastTimerService_cff::range.

◆ charge()

double CastorChannelCoder::charge	(	const reco::castor::QieShape &	fShape,
		int	fAdc,
		int	fCapId
	)		const

ADC[0..127]+capid[0..3]->fC conversion.

Definition at line 21 of file CastorChannelCoder.cc.

                                                                                               {
   int range = (fAdc >> 6) & 0x3;
   double charge = fShape.linearization(fAdc) / mSlope[fCapId][range] + mOffset[fCapId][range];
   //   std::cout << "CastorChannelCoder::charge-> " << fAdc << '/' << fCapId
   //        << " result: " << charge << std::endl;
   return charge;
 }