#include <EcalUncalibRecHitMaxSampleAlgo.h>

Inheritance diagram for EcalUncalibRecHitMaxSampleAlgo< C >:

Public Member Functions
EcalUncalibratedRecHit	makeRecHit (const C &dataFrame, const double pedestals, const double gainRatios, const EcalWeightSet::EcalWeightMatrix weights, const EcalWeightSet::EcalChi2WeightMatrix chi2Matrix) override
	compute rechits More...

	~EcalUncalibRecHitMaxSampleAlgo () override

Public Member Functions inherited from EcalUncalibRecHitRecAbsAlgo< C >
virtual	~EcalUncalibRecHitRecAbsAlgo ()=default
	Constructor. More...

Private Attributes
int16_t	amplitude_

double	chi2_

int16_t	gainId_

int16_t	jitter_

int16_t	pedestal_

int16_t	sampleAdc_

Additional Inherited Members
Public Types inherited from EcalUncalibRecHitRecAbsAlgo< C >
enum	{ nWeightsRows = 3, iAmplitude = 0, iPedestal = 1, iTime = 2 }

Detailed Description

template<class C>
class EcalUncalibRecHitMaxSampleAlgo< C >

Amplitude reconstucted by the difference MAX_adc - min_adc jitter is sample number of MAX_adc, pedestal is min_adc

Author: G. Franzoni, E. Di Marco

Definition at line 15 of file EcalUncalibRecHitMaxSampleAlgo.h.

Constructor & Destructor Documentation

template<class C>

EcalUncalibRecHitMaxSampleAlgo< C >::~EcalUncalibRecHitMaxSampleAlgo ( )

inlineoverride

Definition at line 17 of file EcalUncalibRecHitMaxSampleAlgo.h.

17 {};

Member Function Documentation

template<class C>

EcalUncalibratedRecHit EcalUncalibRecHitMaxSampleAlgo< C >::makeRecHit	(	const C &	dataFrame,
		const double *	pedestals,
		const double *	gainRatios,
		const EcalWeightSet::EcalWeightMatrix **	weights,
		const EcalWeightSet::EcalChi2WeightMatrix **	chi2Matrix
	)

overridevirtual

compute rechits

Implements EcalUncalibRecHitRecAbsAlgo< C >.

Definition at line 31 of file EcalUncalibRecHitMaxSampleAlgo.h.

References EcalUncalibratedRecHit::kSaturated, and min().

Referenced by EcalUncalibRecHitWorkerMaxSample::run().

                                                           {
   amplitude_ = std::numeric_limits<int16_t>::min();
   pedestal_ = 4095;
   jitter_ = -1;
   chi2_ = -1;
   //bool isSaturated = 0;
   uint32_t flags = 0;
   for (int16_t iSample = 0; iSample < C::MAXSAMPLES; iSample++) {
     gainId_ = dataFrame.sample(iSample).gainId();
 
     if (gainId_ == 0) {
       flags = EcalUncalibratedRecHit::kSaturated;
     }
 
     // ampli gain 12
     if (gainId_ == 1) {
       sampleAdc_ = dataFrame.sample(iSample).adc();
     }
 
     else {
       if (gainId_ == 2) {  // ampli gain 6
         sampleAdc_ = 200 + (dataFrame.sample(iSample).adc() - 200) * 2;
       } else {  // accounts for gainId_==3 or 0 - ampli gain 1 and gain0
         sampleAdc_ = 200 + (dataFrame.sample(iSample).adc() - 200) * 12;
       }
     }
 
     if (sampleAdc_ > amplitude_) {
       amplitude_ = sampleAdc_;
       jitter_ = iSample;
     }  // if statement
 
     if (sampleAdc_ < pedestal_)
       pedestal_ = sampleAdc_;
 
   }  // loop on samples
 
   return EcalUncalibratedRecHit(dataFrame.id(),
                                 static_cast<double>(amplitude_ - pedestal_),
                                 static_cast<double>(pedestal_),
                                 static_cast<double>(jitter_ - 5),
                                 chi2_,
                                 flags);
 }