#include <ESElectronicsSimFast.h>

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

typedef CaloTSamples< float, 3 >	ESSamples

Public Member Functions
void	analogToDigital (CLHEP::HepRandomEngine *, 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

Detailed Description

Definition at line 16 of file ESElectronicsSimFast.h.

Member Typedef Documentation

typedef CaloTSamples<float, 3> ESElectronicsSimFast::ESSamples

Definition at line 18 of file ESElectronicsSimFast.h.

Member Enumeration Documentation

anonymous enum

Enumerator
MAXADC
MINADC

Definition at line 20 of file ESElectronicsSimFast.h.

20 { MAXADC = 4095, MINADC = 0 };

ESElectronicsSimFast::MAXADC

Definition: ESElectronicsSimFast.h:20

ESElectronicsSimFast::MINADC

Definition: ESElectronicsSimFast.h:20

Constructor & Destructor Documentation

ESElectronicsSimFast::ESElectronicsSimFast	(	bool	addNoise,
		bool	PreMix1
	)

Definition at line 8 of file ESElectronicsSimFast.cc.

     : m_addNoise(addNoise), m_PreMix1(PreMix1), m_MIPToGeV(0), m_peds(nullptr), m_mips(nullptr) {
   // 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
 }

ESElectronicsSimFast::~ESElectronicsSimFast ( )

Definition at line 16 of file ESElectronicsSimFast.cc.

16 {}

Member Function Documentation

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

Definition at line 24 of file ESElectronicsSimFast.cc.

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

                                                                {
   assert(nullptr != m_peds && nullptr != m_mips && 0 < m_MIPToGeV);  // 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 * CLHEP::RandGaussQ::shoot(engine, 0, 1));
     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 ;
 }