#include <ESElectronicsSimFast.h>

Public Types
enum	{ MAXADC = 4095, MINADC = 0 }

typedef CaloTSamples< float, 3 >	ESSamples

Public Member Functions
void	analogToDigital (ESSamples &cs, ESDataFrame &df, bool isNoise=false) const

	ESElectronicsSimFast (bool addNoise, bool PreMix1)

void	newEvent ()

void	setMIPs (const ESIntercalibConstants *mips)

void	setMIPToGeV (double MIPToGeV)

void	setPedestals (const ESPedestals *peds)

	~ESElectronicsSimFast ()

Private Attributes
bool	m_addNoise

const ESIntercalibConstants *	m_mips

double	m_MIPToGeV

const ESPedestals *	m_peds

bool	m_PreMix1

CLHEP::RandGaussQ *	m_ranGau

Detailed Description

Definition at line 15 of file ESElectronicsSimFast.h.

Member Typedef Documentation

typedef CaloTSamples<float,3> ESElectronicsSimFast::ESSamples

Definition at line 19 of file ESElectronicsSimFast.h.

Member Enumeration Documentation

anonymous enum

Enumerator
MAXADC
MINADC

Definition at line 21 of file ESElectronicsSimFast.h.

21 { MAXADC = 4095,

22 MINADC = 0 } ;

ESElectronicsSimFast::MINADC

Definition: ESElectronicsSimFast.h:22

ESElectronicsSimFast::MAXADC

Definition: ESElectronicsSimFast.h:21

Constructor & Destructor Documentation

ESElectronicsSimFast::ESElectronicsSimFast	(	bool	addNoise,
		bool	PreMix1
	)

Definition at line 11 of file ESElectronicsSimFast.cc.

References edm::hlt::Exception, edm::RandomNumberGenerator::getEngine(), edm::Service< T >::isAvailable(), m_addNoise, and m_ranGau.

                                                                          :
    m_addNoise ( addNoise ) ,
    m_PreMix1  ( PreMix1  ) ,
    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 ( )

Definition at line 38 of file ESElectronicsSimFast.cc.

References m_ranGau.

 {
    delete m_ranGau ;
 }

Member Function Documentation

void ESElectronicsSimFast::analogToDigital	(	ESSamples &	cs,
		ESDataFrame &	df,
		bool	isNoise = `false`
	)		const

Definition at line 62 of file ESElectronicsSimFast.cc.

References ecalMGPA::adc(), ESCondObjectContainer< T >::end(), ESCondObjectContainer< T >::find(), ESCondObjectContainer< T >::getMap(), i, CaloTSamplesBase< Ttype >::id(), m_addNoise, m_mips, m_MIPToGeV, m_peds, m_PreMix1, m_ranGau, MAXADC, MINADC, ESDataFrame::setSample(), ESDataFrame::setSize(), and CaloTSamplesBase< Ttype >::size().

 {
    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;
 
       if(!m_PreMix1) signal = cs[i]*ADCGeV + noi + baseline ;
       else signal = cs[i]*ADCGeV ;
 
       if( 0 <= signal )
       { 
      signal += 0.5 ;
       }
       else
       {
      signal -= 0.5 ;
       }
     
       adc = int( signal ) ;
 
       if(!m_PreMix1) 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 ;
 }