#include <HcalFeatureHFEMBit.h>

Inheritance diagram for HcalFeatureHFEMBit:

Public Member Functions
virtual 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

virtual 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 ()

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 9 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 ( )

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 40 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 24 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(), AlignmentProducer_cff::calibrations, conditions_, corr, HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), hcalTTPDigis_cfi::id, and hcalTTPDigis_cfi::samples.

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.pedestal(f[idx].capid());
    auto corr = calibrations.respcorrgain(f[idx].capid());
 
    return (samples[idx] - ped) * corr;
 }