#include <EcalUncalibRecHitMaxSampleAlgo.h>

Inheritance diagram for EcalUncalibRecHitMaxSampleAlgo< C >:

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

virtual	~EcalUncalibRecHitMaxSampleAlgo ()

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>

virtual EcalUncalibRecHitMaxSampleAlgo< C >::~EcalUncalibRecHitMaxSampleAlgo ( )

inlinevirtual

Definition at line 20 of file EcalUncalibRecHitMaxSampleAlgo.h.

20 { };

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
	)

virtual

compute rechits

Implements EcalUncalibRecHitRecAbsAlgo< C >.

Definition at line 34 of file EcalUncalibRecHitMaxSampleAlgo.h.

References flags, 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 );
 }