#include <ESElectronicsSim.h>

Public Types
enum	{ MAXADC = 4095 }

enum	{ MINADC = 0 }

Public Member Functions
virtual void	analogToDigital (const CaloSamples &cs, ESDataFrame &df) const

virtual void	digitalToAnalog (const ESDataFrame &df, CaloSamples &cs) const

	ESElectronicsSim (bool addNoise)

void	newEvent ()
	anything that needs to be done once per event More...

void	setGain (const int gain)

void	setMIPs (const ESIntercalibConstants *mips)

void	setMIPToGeV (const double MIPToGeV)

void	setPedestals (const ESPedestals *peds)

virtual	~ESElectronicsSim ()

Private Member Functions
double	decode (const ESSample &sample, const DetId &detId) const

std::vector< ESSample >	encode (const CaloSamples &timeframe) const

Private Attributes
bool	addNoise_

int	gain_

const ESIntercalibConstants *	mips_

double	MIPToGeV_

const ESPedestals *	peds_

Detailed Description

Definition at line 12 of file ESElectronicsSim.h.

Member Enumeration Documentation

anonymous enum

Enumerator
MAXADC

Definition at line 16 of file ESElectronicsSim.h.

16 {MAXADC = 4095};

ESElectronicsSim::MAXADC

Definition: ESElectronicsSim.h:16

anonymous enum

Enumerator
MINADC

Definition at line 17 of file ESElectronicsSim.h.

17 {MINADC = 0};

ESElectronicsSim::MINADC

Definition: ESElectronicsSim.h:17

Constructor & Destructor Documentation

ESElectronicsSim::ESElectronicsSim ( bool addNoise )

Definition at line 11 of file ESElectronicsSim.cc.

                                                 :
   addNoise_(addNoise), peds_(0), mips_(0)
 {
   // Preshower 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
 }

ESElectronicsSim::~ESElectronicsSim ( )

virtual

Definition at line 20 of file ESElectronicsSim.cc.

21 {}

Member Function Documentation

void ESElectronicsSim::analogToDigital	(	const CaloSamples &	cs,
		ESDataFrame &	df
	)		const

virtual

Definition at line 23 of file ESElectronicsSim.cc.

References encode(), i, ESDataFrame::setSample(), ESDataFrame::setSize(), ESDataFrame::size(), and CaloSamples::size().

 {
 
   std::vector<ESSample> essamples = encode(cs);
 
   df.setSize(cs.size());
   for(int i=0; i<df.size(); i++) {
     df.setSample(i, essamples[i]);
   }
 
 }

double ESElectronicsSim::decode	(	const ESSample &	sample,
		const DetId &	detId
	)		const

private

Definition at line 96 of file ESElectronicsSim.cc.

Referenced by digitalToAnalog().

 {
   return 0. ;
 }

void ESElectronicsSim::digitalToAnalog	(	const ESDataFrame &	df,
		CaloSamples &	cs
	)		const

virtual

Definition at line 35 of file ESElectronicsSim.cc.

References decode(), i, ESDataFrame::id(), and ESDataFrame::size().

 {
 
   for(int i = 0; i < df.size(); i++) {
     cs[i] = decode(df[i], df.id());
   }
 
 }

std::vector< ESSample > ESElectronicsSim::encode ( const CaloSamples & timeframe ) const

private

Definition at line 45 of file ESElectronicsSim.cc.

References ecalMGPA::adc(), addNoise_, edm::hlt::Exception, ESCondObjectContainer< T >::find(), edm::RandomNumberGenerator::getEngine(), ESCondObjectContainer< T >::getMap(), i, CaloSamples::id(), edm::Service< T >::isAvailable(), MAXADC, MINADC, mips_, MIPToGeV_, peds_, python.entryComment::results, and CaloSamples::size().

Referenced by analogToDigital().

 {
   edm::Service<edm::RandomNumberGenerator> rng;
   if ( ! rng.isAvailable()) {
     throw cms::Exception("Configuration")
       << "ESElectroncSim 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.";
   }
 
 
   std::vector<ESSample> results;
   results.reserve(timeframe.size());
 
   ESPedestals::const_iterator it_ped = peds_->find(timeframe.id());
   ESIntercalibConstantMap::const_iterator it_mip = mips_->getMap().find(timeframe.id());
   int baseline_  = (int) it_ped->getMean();
   double sigma_  = (double) it_ped->getRms();
   double MIPADC_ = (double) (*it_mip);
 
   int adc = 0; 
   double ADCGeV = MIPADC_/MIPToGeV_;
 
   for (int i=0; i<timeframe.size(); i++) {
 
     double noi = 0;
     double signal = 0;    
 
     if (addNoise_) {
       CLHEP::RandGaussQ gaussQDistribution(rng->getEngine(), 0., sigma_);
       noi = gaussQDistribution.fire();
     }
 
     signal = timeframe[i]*ADCGeV + noi + baseline_;
     
     if (signal>0) 
       signal += 0.5;
     else if (signal<0)
       signal -= 0.5;
     
     adc = int(signal);
 
     if (adc>MAXADC) adc = MAXADC;
     if (adc<MINADC) adc = MINADC;
 
     results.push_back(ESSample(adc));
   }
 
   return results;
 }