#include <HcalChannelCoder.h>

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

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

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

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 QIE $Author: ratnikov

Definition at line 12 of file HcalChannelCoder.h.

Constructor & Destructor Documentation

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

Definition at line 12 of file HcalChannelCoder.cc.

References index(), mOffset, and mSlope.

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

Member Function Documentation

int HcalChannelCoder::adc	(	const QieShape &	fShape,
		double	fCharge,
		int	fCapId
	)		const

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

Definition at line 29 of file HcalChannelCoder.cc.

References newFWLiteAna::bin, QieShape::binSize(), QieShape::linearization(), mOffset, and mSlope.

                                                                                    {
 
   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 << "HcalChannelCoder::adc-> " << fCharge << '/' << fCapId 
 //      << " result: " << adc << std::endl;
   return adc;
 }

double HcalChannelCoder::charge	(	const QieShape &	fShape,
		int	fAdc,
		int	fCapId
	)		const

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

Definition at line 21 of file HcalChannelCoder.cc.

References QieShape::linearization(), mOffset, and mSlope.

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