#include <ESElectronicsSim.h>

Public Types
enum	{ MINADC = 0 }

enum	{ MAXADC = 4095 }

Public Member Functions
virtual void	analogToDigital (CLHEP::HepRandomEngine *, const CaloSamples &cs, ESDataFrame &df) const

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

	ESElectronicsSim (bool addNoise)

void	newEvent (CLHEP::HepRandomEngine *)
	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, CLHEP::HepRandomEngine *) const

Private Attributes
bool	addNoise_

int	gain_

const ESIntercalibConstants *	mips_

double	MIPToGeV_

const ESPedestals *	peds_

Detailed Description

Definition at line 16 of file ESElectronicsSim.h.

Member Enumeration Documentation

anonymous enum

Enumerator
MINADC

Definition at line 19 of file ESElectronicsSim.h.

19 { MINADC = 0 };

ESElectronicsSim::MINADC

Definition: ESElectronicsSim.h:19

anonymous enum

Enumerator
MAXADC

Definition at line 18 of file ESElectronicsSim.h.

18 { MAXADC = 4095 };

ESElectronicsSim::MAXADC

Definition: ESElectronicsSim.h:18

Constructor & Destructor Documentation

ESElectronicsSim::ESElectronicsSim ( bool addNoise )

Definition at line 11 of file ESElectronicsSim.cc.

                                                 : addNoise_(addNoise), peds_(nullptr), mips_(nullptr) {
   // 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 18 of file ESElectronicsSim.cc.

18 {}

Member Function Documentation

void ESElectronicsSim::analogToDigital	(	CLHEP::HepRandomEngine *	engine,
		const CaloSamples &	cs,
		ESDataFrame &	df
	)		const

virtual

Definition at line 20 of file ESElectronicsSim.cc.

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

                                                                                                                  {
   std::vector<ESSample> essamples = encode(cs, engine);
 
   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 84 of file ESElectronicsSim.cc.

Referenced by digitalToAnalog().

84 { return 0.; }

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

virtual

Definition at line 29 of file ESElectronicsSim.cc.

References decode(), mps_fire::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,
		CLHEP::HepRandomEngine *	engine
	)		const

private

Definition at line 35 of file ESElectronicsSim.cc.

References ecalMGPA::adc(), addNoise_, Exception, ESCondObjectContainer< T >::find(), ESCondObjectContainer< T >::getMap(), mps_fire::i, CaloSamples::id(), createfilelist::int, edm::Service< T >::isAvailable(), MAXADC, MINADC, mips_, MIPToGeV_, peds_, 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_) {
       noi = CLHEP::RandGaussQ::shoot(engine, 0., sigma_);
     }
 
     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.emplace_back(adc);
   }
 
   return results;
 }