EcalElectronicsSim.cc

Go to the documentation of this file.

#include "SimCalorimetry/EcalSimAlgos/interface/EcalElectronicsSim.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EcalCoder.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EcalSimParameterMap.h"

#include "CLHEP/Random/RandGaussQ.h"

#include <cstring>
#include <sstream>
#include <iostream>
#include <unistd.h>
#include <fstream>

EcalElectronicsSim::EcalElectronicsSim( const EcalSimParameterMap* parameterMap      , 
                    EcalCoder*                 coder             , 
                    bool                       applyConstantTerm , 
                    double                     rmsConstantTerm     ) :
   m_simMap             ( parameterMap ) ,
   m_theCoder           ( coder        ) ,
   m_thisCT             ( rmsConstantTerm ),
   m_applyConstantTerm  ( applyConstantTerm )
{
}

EcalElectronicsSim::~EcalElectronicsSim()
{  
}

void 
EcalElectronicsSim::analogToDigital( CLHEP::HepRandomEngine* engine,
                                     EcalElectronicsSim::EcalSamples& clf ,
                     EcalDataFrame&                   df    ) const 
{

   //PG input signal is in pe.  Converted in GeV
  amplify( clf, engine ) ;

  m_theCoder->analogToDigital( engine, clf, df ) ;
}

void 
EcalElectronicsSim::amplify( EcalElectronicsSim::EcalSamples& 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 ;
   }
}