/data/refman/pasoursint/CMSSW_5_3_10_patch1/src/RecoLocalCalo/EcalRecAlgos/interface/EcalUncalibRecHitLeadingEdgeAlgo.h

Go to the documentation of this file.

#ifndef RecoLocalCalo_EcalRecAlgos_EcalUncalibRecHitLeadingEdgeAlgo_HH
#define RecoLocalCalo_EcalRecAlgos_EcalUncalibRecHitLeadingEdgeAlgo_HH

#include "Math/SVector.h"
#include "Math/SMatrix.h"
#include "RecoLocalCalo/EcalRecAlgos/interface/EcalUncalibRecHitRecAbsAlgo.h"

#include "DataFormats/EcalDetId/interface/EBDetId.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"


template < class C > class EcalUncalibRecHitLeadingEdgeAlgo : public EcalUncalibRecHitRecAbsAlgo < C > {
      public:
        // destructor
        EcalUncalibRecHitLeadingEdgeAlgo < C > () : leadingSample_(0), shape_(0) { };
        virtual ~ EcalUncalibRecHitLeadingEdgeAlgo < C > () { };

        void setLeadingEdgeSample( int isample ) { leadingSample_ = isample; }
        int getLeadingEdgeSample() { return leadingSample_; }

        void setPulseShape( std::vector<double> & shape ) { shape_ = shape; }
        std::vector<double> & getPulseShape() { return shape_; }

        virtual EcalUncalibratedRecHit makeRecHit(const C & dataFrame, 
                        const double *pedestals,
                        const double *gainRatios,
                        const EcalWeightSet::EcalWeightMatrix** weights, 
                        const EcalWeightSet::EcalChi2WeightMatrix** chi2Matrix
                        )
        {
                double amplitude_(-1.), pedestal_(-1.), jitter_(-1.), chi2_(-1.);

                // compute amplitude
                amplitude_ = double (((double) (dataFrame.sample( leadingSample_ ).adc() -
                                              pedestals[ dataFrame.sample( leadingSample_ ).gainId() - 1]) * 
                                              saturationCorrection( leadingSample_ ) *
                                              gainRatios[ dataFrame.sample( leadingSample_ ).gainId() - 1] )); 



                return EcalUncalibratedRecHit(dataFrame.id(), amplitude_, pedestal_, jitter_, chi2_);
        }

        // saturation correction  
        double saturationCorrection(int unsaturatedSample)
        {
                if ( unsaturatedSample > 0 && unsaturatedSample < (int)shape_.size() ) {
                        return 1./ shape_[ unsaturatedSample ];
                } else {
                        edm::LogError("EcalUncalibRecHitLeadingEdgeAlgo") << "Invalid sample " << unsaturatedSample 
                                << " for a shape vector of size " << shape_.size();
                        return 0;
                }
        }

        private:
                int leadingSample_;
                std::vector< double > shape_;
};
#endif