#include <HcalFeatureHFEMBit.h>

Inheritance diagram for HcalFeatureHFEMBit:

Public Member Functions
bool	fineGrainbit (const QIE10DataFrame &short1, const QIE10DataFrame &short2, const QIE10DataFrame &long1, const QIE10DataFrame &long2, bool validShort1, bool validShort2, bool validLong1, bool validLong2, int idx) const override

bool	fineGrainbit (const HFDataFrame &shortDigi, const HFDataFrame &longDigi, int idx) const override

	HcalFeatureHFEMBit (double ShortMinE, double LongMinE, double ShortLongCutSlope, double ShortLongCutOffset, const HcalDbService &conditions)

	~HcalFeatureHFEMBit () override

Public Member Functions inherited from HcalFeatureBit
	HcalFeatureBit ()

virtual	~HcalFeatureBit ()

Private Member Functions
template<typename T >
float	getE (const T &f, int idx) const

Private Attributes
const HcalDbService &	conditions_

double	LongMinE_

double	ShortLongCutOffset_

double	ShortLongCutSlope_

double	ShortMinE_

Detailed Description

Definition at line 8 of file HcalFeatureHFEMBit.h.

Constructor & Destructor Documentation

HcalFeatureHFEMBit::HcalFeatureHFEMBit	(	double	ShortMinE,
		double	LongMinE,
		double	ShortLongCutSlope,
		double	ShortLongCutOffset,
		const HcalDbService &	conditions
	)

Definition at line 8 of file HcalFeatureHFEMBit.cc.

References LongMinE_, ShortLongCutOffset_, ShortLongCutSlope_, and ShortMinE_.

     : conditions_(conditions) {
   ShortMinE_ = ShortMinE;  //minimum energy deposited
   LongMinE_ = LongMinE;
   ShortLongCutSlope_ =
       ShortLongCutSlope;  // this is a the slope of the cut line related to energy deposited in short fibers vrs long fibers
   ShortLongCutOffset_ = ShortLongCutOffset;  // this is the offset of said line.
 }

HcalFeatureHFEMBit::~HcalFeatureHFEMBit ( )

override

Definition at line 21 of file HcalFeatureHFEMBit.cc.

21 {}

Member Function Documentation

bool HcalFeatureHFEMBit::fineGrainbit	(	const QIE10DataFrame &	short1,
		const QIE10DataFrame &	short2,
		const QIE10DataFrame &	long1,
		const QIE10DataFrame &	long2,
		bool	validShort1,
		bool	validShort2,
		bool	validLong1,
		bool	validLong2,
		int	idx
	)		const

overridevirtual

Implements HcalFeatureBit.

Definition at line 35 of file HcalFeatureHFEMBit.cc.

References getE(), LongMinE_, ShortLongCutOffset_, ShortLongCutSlope_, and ShortMinE_.

                                                      {
   float shortE = 0;
   if (validShort1)
     shortE += getE(short1, idx);
   if (validShort2)
     shortE += getE(short2, idx);
   if (validShort1 and validShort2)
     shortE *= .5;
 
   float longE = 0;
   if (validLong1)
     longE += getE(long1, idx);
   if (validLong2)
     longE += getE(long2, idx);
   if (validLong1 and validLong2)
     longE *= .5;
 
   if (shortE < ShortMinE_)
     return false;
   if (longE < LongMinE_)
     return false;
 
   return (shortE < (longE - ShortLongCutOffset_) * ShortLongCutSlope_);
 }

bool HcalFeatureHFEMBit::fineGrainbit	(	const HFDataFrame &	shortDigi,
		const HFDataFrame &	longDigi,
		int	idx
	)		const

overridevirtual

Implements HcalFeatureBit.

Definition at line 23 of file HcalFeatureHFEMBit.cc.

References getE(), LongMinE_, ShortLongCutOffset_, ShortLongCutSlope_, and ShortMinE_.

                                                                                                               {
   float shortE = getE(shortDigi, idx);
   float longE = getE(longDigi, idx);
 
   if (shortE < ShortMinE_)
     return false;
   if (longE < LongMinE_)
     return false;
 
   return (shortE < (longE - ShortLongCutOffset_) * ShortLongCutSlope_);
 }

template<typename T >

float HcalFeatureHFEMBit::getE	(	const T &	f,
		int	idx
	)		const

private

Definition at line 38 of file HcalFeatureHFEMBit.h.

References HcalCoderDb::adc2fC(), conditions_, alignCSCRings::corr, EcalCondDB::db, HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), and gpuClustering::id.

Referenced by fineGrainbit().

                                                         {
   const HcalDetId id(f.id());
   const HcalCalibrations& calibrations = conditions_.getHcalCalibrations(id);
   const auto* coder = conditions_.getHcalCoder(id);
   const auto* shape = conditions_.getHcalShape(coder);
 
   HcalCoderDb db(*coder, *shape);
   CaloSamples samples;
   db.adc2fC(f, samples);
 
   auto ped = calibrations.effpedestal(f[idx].capid());
   auto corr = calibrations.respcorrgain(f[idx].capid());
 
   return (samples[idx] - ped) * corr;
 }