---

EcalCoder.h

Go to the documentation of this file.

#ifndef EcalSimAlgos_EcalCoder_h
#define EcalSimAlgos_EcalCoder_h 1

#include "FWCore/Framework/interface/EDProducer.h"
#include "DataFormats/Common/interface/EDProduct.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CalibFormats/CaloObjects/interface/CaloSamples.h"
#include "SimGeneral/NoiseGenerators/interface/CorrelatedNoisifier.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EcalCorrelatedNoiseMatrix.h"

#include "CalibCalorimetry/EcalTrivialCondModules/interface/EcalTrivialConditionRetriever.h"


class EcalMGPASample;
class EBDataFrame;
class EEDataFrame;
class DetId;
#include<vector>

/* \class EEDigitizerTraits
 * \brief Converts CaloDataFrame in CaloTimeSample and vice versa.
 *
 */
class EcalCoder
{
 public:
  enum {NBITS = 12};
  // 2^12 -1
  enum {MAXADC = 4095}; 
  enum {ADCGAINSWITCH = 4079};
  enum {NGAINS = 3};

  EcalCoder(bool addNoise, CorrelatedNoisifier<EcalCorrMatrix> * theCorrNoise) ;
  virtual ~EcalCoder() {}

  void setPedestals(const EcalPedestals * pedestals);

  void setGainRatios(const EcalGainRatios * gainRatios);

  void setFullScaleEnergy(const double EBscale , const double EEscale) {m_maxEneEB = EBscale; m_maxEneEE = EEscale; }

 
  virtual void digitalToAnalog(const EBDataFrame& df, CaloSamples& lf) const;
  virtual void digitalToAnalog(const EEDataFrame& df, CaloSamples& lf) const;
  virtual void analogToDigital(const CaloSamples& clf, EcalDataFrame& df) const;
 
  void newEvent() {}

 private:

  double fullScaleEnergy (const DetId & ) const ;

  void encode(const CaloSamples& caloSamples, EcalDataFrame& df) const;

  double decode(const EcalMGPASample & sample, const DetId & detId) const;

  void noisify(float * values, int size) const;

  void findPedestal(const DetId & detId, int gainId, 
                    double & pedestal, double & width) const;
    
  double theGains[NGAINS+1];
   
  void findGains(const DetId & detId, double theGains[] ) const;
   
  const EcalPedestals * thePedestals;
  const EcalGainRatios * theGainRatios;
  double m_maxEneEB ;
  double m_maxEneEE ;
  bool addNoise_;
  CorrelatedNoisifier<EcalCorrMatrix> * theCorrNoise_;

};


#endif

---