EcalLiteDTUCoder.h

Go to the documentation of this file.

#ifndef SimCalorimetry_EcalSimAlgos_EcalLiteDTUCoder_h
#define SimCalorimetry_EcalSimAlgos_EcalLiteDTUCoder_h
 
#include "CalibFormats/CaloObjects/interface/CaloTSamples.h"
#include "CondFormats/EcalObjects/interface/EcalLiteDTUPedestals.h"
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstantsMC.h"
#include "CondFormats/EcalObjects/interface/EcalCATIAGainRatios.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EcalCorrelatedNoiseMatrix.h"
#include "DataFormats/EcalDigi/interface/EcalConstants.h"
 
template <typename M>
class CorrelatedNoisifier;
class EcalLiteDTUSample;
class EcalDataFrame_Ph2;
class DetId;
class EcalLiteDTUPed;
 
#include <vector>
 
namespace CLHEP {
  class HepRandomEngine;
}
 
class EcalLiteDTUCoder {
public:
  typedef CaloTSamples<float, ecalPh2::sampleSize> EcalSamples;
 
  typedef CorrelatedNoisifier<EcalCorrMatrix_Ph2> Noisifier;
 
  EcalLiteDTUCoder(bool addNoise, bool PreMix1, Noisifier* ebCorrNoise0, Noisifier* ebCorrNoise1 = nullptr);
 
  virtual ~EcalLiteDTUCoder();
 
  void setPedestals(const EcalLiteDTUPedestalsMap* pedestals);
 
  void setGainRatios(float gainRatios);
 
  void setFullScaleEnergy(double EBscale);
 
  void setIntercalibConstants(const EcalIntercalibConstantsMC* ical);
 
  virtual void analogToDigital(CLHEP::HepRandomEngine*, const EcalSamples& clf, EcalDataFrame_Ph2& df) const;
 
private:
  double fullScaleEnergy(const DetId& did) const;
 
  void encode(const EcalSamples& ecalSamples, EcalDataFrame_Ph2& df, CLHEP::HepRandomEngine*) const;
 
  void findPedestal(const DetId& detId, int gainId, double& pedestal, double& width) const;
 
  void findGains(const DetId& detId, float theGains[]) const;
 
  void findIntercalibConstant(const DetId& detId, double& icalconst) const;
 
  const EcalLiteDTUPedestalsMap* m_peds;
 
  float 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
 
  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[ecalPh2::NGAINS];
};
 
#endif