/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/SimCalorimetry/EcalSimAlgos/src/ESElectronicsSimFast.cc

Go to the documentation of this file.

#include "SimCalorimetry/EcalSimAlgos/interface/ESElectronicsSimFast.h"
#include "DataFormats/DetId/interface/DetId.h"

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/RandomNumberGenerator.h"
#include "CLHEP/Random/RandGaussQ.h"
#include "FWCore/Utilities/interface/Exception.h"

#include <iostream>

ESElectronicsSimFast::ESElectronicsSimFast( bool addNoise ) :
   m_addNoise ( addNoise ) ,
   m_MIPToGeV (        0 ) ,
   m_peds     (        0 ) ,
   m_mips     (        0 ) ,
   m_ranGau   (        0 )
{
   // Preshower "Fast" Electronics Simulation
   // gain = 1 : low gain for data taking 
   // gain = 2 : high gain for calibration and low energy runs
   // For 300(310/320) um Si, the MIP is 78.47(81.08/83.7) keV

   if( m_addNoise )
   {
      edm::Service<edm::RandomNumberGenerator> rng;
      if( !rng.isAvailable() ) 
      {
         throw cms::Exception("Configuration")
            << "ESElectroncSimFast requires the RandomNumberGeneratorService\n"
            "which is not present in the configuration file.  You must add the service\n"
            "in the configuration file or remove the modules that require it.";
      }
      m_ranGau = new CLHEP::RandGaussQ( rng->getEngine(), 0, 1 ) ;
   }
}

ESElectronicsSimFast::~ESElectronicsSimFast()
{
   delete m_ranGau ;
}

void 
ESElectronicsSimFast::setPedestals( const ESPedestals* peds ) 
{
   m_peds = peds ;
} 

void 
ESElectronicsSimFast::setMIPs( const ESIntercalibConstants* mips ) 
{
   m_mips = mips ;
}

void 
ESElectronicsSimFast::setMIPToGeV( double MIPToGeV )
{
   m_MIPToGeV = MIPToGeV ;
}

void 
ESElectronicsSimFast::analogToDigital( ESSamples&   cs, 
                                       ESDataFrame& df, 
                                       bool         isNoise ) const
{
   assert( 0 != m_peds &&
           0 != m_mips &&
           0 < m_MIPToGeV ) ; // sanity check

   assert( ( !m_addNoise ) ||
           0 != m_ranGau ) ; // sanity check

   df.setSize( cs.size() ) ;

   const DetId id ( cs.id() ) ;
   ESPedestals::const_iterator it_ped ( m_peds->find( id ) ) ;
   ESIntercalibConstantMap::const_iterator it_mip (
      isNoise ? m_mips->getMap().end() : m_mips->getMap().find( id ) ) ;

   const double baseline ( (double) it_ped->getMean() ) ;
   const double sigma    ( isNoise ? 0. : (double) it_ped->getRms() ) ;
   const double MIPADC   ( isNoise ? 0. : (double) (*it_mip) ) ;
   const double ADCGeV   ( isNoise ? 1. : MIPADC/m_MIPToGeV ) ;

   int adc = 0 ;
//   std::cout<<"   **Id="<<ESDetId(df.id())<<", size="<<df.size();
   for( unsigned int i ( 0 ) ; i != cs.size(); ++i ) 
   {
      const double noi ( isNoise || (!m_addNoise) ? 0 :
                         sigma*m_ranGau->fire() ) ;
    
      double signal = cs[i]*ADCGeV + noi + baseline ;

      if( 0 <= signal )
      { 
         signal += 0.5 ;
      }
      else
      {
         signal -= 0.5 ;
      }
    
      adc = int( signal ) ;
      assert( 0 < adc ) ;

      if( 0.5 < signal - adc ) ++adc ;

      if( MAXADC < adc )
      {
         adc = MAXADC ;
      }
      else
      {
         if( MINADC > adc ) adc = MINADC ;
      }

      df.setSample( i, ESSample( adc ) ) ;
//      std::cout<<", "<<df[i];
   }
//   std::cout<<std::endl ;
}