/data/git/CMSSW_5_3_11_patch5/src/RecoLocalCalo/EcalRecAlgos/interface/EcalUncalibRecHitMaxSampleAlgo.h

Go to the documentation of this file.

#ifndef RecoLocalCalo_EcalRecAlgos_EcalUncalibRecHitMaxSampleAlgo_HH
#define RecoLocalCalo_EcalRecAlgos_EcalUncalibRecHitMaxSampleAlgo_HH

#include "RecoLocalCalo/EcalRecAlgos/interface/EcalUncalibRecHitRecAbsAlgo.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"


template<class C> class EcalUncalibRecHitMaxSampleAlgo : public EcalUncalibRecHitRecAbsAlgo<C>
{
  
 public:
  
  virtual ~EcalUncalibRecHitMaxSampleAlgo<C>() { };
  virtual EcalUncalibratedRecHit makeRecHit(const C& dataFrame, const double* pedestals,
                                            const double* gainRatios,
                                            const EcalWeightSet::EcalWeightMatrix** weights,
                                            const EcalWeightSet::EcalChi2WeightMatrix** chi2Matrix);

 private:
  int16_t  amplitude_,  pedestal_, jitter_, sampleAdc_, gainId_;
  double   chi2_;

};

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

  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 );
}

#endif