EcalCoder.h

Go to the documentation of this file.

#ifndef EcalSimAlgos_EcalCoder_h
#define EcalSimAlgos_EcalCoder_h 1

#include "CalibFormats/CaloObjects/interface/CaloTSamples.h"
#include "CondFormats/EcalObjects/interface/EcalPedestals.h"
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstantsMC.h"
#include "CondFormats/EcalObjects/interface/EcalGainRatios.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EcalCorrelatedNoiseMatrix.h"

template <typename M>
class CorrelatedNoisifier;
class EcalMGPASample;
class EcalDataFrame;
class DetId;

#include <vector>

namespace CLHEP {
  class HepRandomEngine;
}

/* \class EEDigitizerTraits
 * \brief Converts CaloDataFrame in CaloTimeSample and vice versa.
 *
 */
class EcalCoder {
public:
  typedef CaloTSamples<float, 10> EcalSamples;

  typedef CorrelatedNoisifier<EcalCorrMatrix> Noisifier;

  enum {
    NBITS = 12,            // number of available bits
    MAXADC = 4095,         // 2^12 -1,  adc max range
    ADCGAINSWITCH = 4079,  // adc gain switch
    NGAINS = 3             // number of electronic gains
  };

  EcalCoder(bool addNoise,
            bool PreMix1,
            Noisifier* ebCorrNoise0,
            Noisifier* eeCorrNoise0 = nullptr,
            Noisifier* ebCorrNoise1 = nullptr,
            Noisifier* eeCorrNoise1 = nullptr,
            Noisifier* ebCorrNoise2 = nullptr,
            Noisifier* eeCorrNoise2 = nullptr);  // make EE version optional for tb compatibility
  virtual ~EcalCoder();

  void setPedestals(const EcalPedestals* pedestals);

  void setGainRatios(const EcalGainRatios* gainRatios);

  void setFullScaleEnergy(double EBscale, double EEscale);

  void setIntercalibConstants(const EcalIntercalibConstantsMC* ical);

  virtual void analogToDigital(CLHEP::HepRandomEngine*, const EcalSamples& clf, EcalDataFrame& df) const;

private:
  double fullScaleEnergy(const DetId& did) const;

  void encode(const EcalSamples& ecalSamples, EcalDataFrame& df, CLHEP::HepRandomEngine*) const;

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

  //      void noisify( const EcalIntercalibConstantsMC* values ,
  //            int                              size     ) const ;

  void findPedestal(const DetId& detId, int gainId, double& pedestal, double& width) const;

  void findGains(const DetId& detId, double theGains[]) const;

  void findIntercalibConstant(const DetId& detId, double& icalconst) const;

  const EcalPedestals* m_peds;

  const EcalGainRatios* m_gainRatios;  // the electronics gains

  const EcalIntercalibConstantsMC* m_intercals;  //record specific for simulation of gain variation in MC

  double m_maxEneEB;  // max attainable energy in the ecal barrel
  double m_maxEneEE;  // max attainable energy in the ecal endcap

  bool m_addNoise;  // whether add noise to the pedestals and the gains
  bool m_PreMix1;   // Follow necessary steps for PreMixing input

  const Noisifier* m_ebCorrNoise[3];
  const Noisifier* m_eeCorrNoise[3];
};

#endif