EcalElectronicsSim.h

Go to the documentation of this file.

#ifndef EcalSimAlgos_EcalElectronicsSim_h
#define EcalSimAlgos_EcalElectronicsSim_h

#include "SimCalorimetry/CaloSimAlgos/interface/CaloVNoiseSignalGenerator.h"
#include "CLHEP/Random/RandGaussQ.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EcalSimParameterMap.h"

class EcalSimParameterMap;

namespace CLHEP {
  class HepRandomEngine;
}

/* \class EcalElectronicsSim
 * \brief Converts CaloDataFrame in CaloTimeSample and vice versa.
 * 
 */

template <typename CoderType, typename SamplesType, typename DataFrameType>
class EcalElectronicsSim {
public:
  EcalElectronicsSim(const EcalSimParameterMap* parameterMap,
                     CoderType* coder,
                     bool applyConstantTerm,
                     double rmsConstantTerm)
      : m_simMap(parameterMap), m_theCoder(coder), m_thisCT(rmsConstantTerm), m_applyConstantTerm(applyConstantTerm) {}


  void analogToDigital(CLHEP::HepRandomEngine* engine, SamplesType& clf, DataFrameType& df) const {
    // input signal is in pe.  Converted in GeV
    amplify(clf, engine);
    m_theCoder->analogToDigital(engine, clf, df);
  }

  void newEvent() {}

  void setNoiseSignalGenerator(const CaloVNoiseSignalGenerator* noiseSignalGenerator) {
    theNoiseSignalGenerator = noiseSignalGenerator;
  }

private:
  void amplify(SamplesType& clf, CLHEP::HepRandomEngine* engine) const {
    const double fac(m_simMap->simParameters(clf.id()).photoelectronsToAnalog());
    if (m_applyConstantTerm) {
      clf *= fac * CLHEP::RandGaussQ::shoot(engine, 1.0, m_thisCT);
    } else {
      clf *= fac;
    }
  }


  const EcalSimParameterMap* m_simMap;

  const CaloVNoiseSignalGenerator* theNoiseSignalGenerator;

  CoderType* m_theCoder;

  const double m_thisCT;
  const bool m_applyConstantTerm;
};

#endif