ESElectronicsSimFast Class Reference

#include <ESElectronicsSimFast.h>

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

Public Member Functions
virtual void	analogToDigital (const CaloSamples &cs, ESDataFrame &df, bool wasEmpty, CLHEP::RandGeneral *histoDistribution, double hInf, double hSup, double hBin) const
void	digitalToAnalog (const ESDataFrame &df, CaloSamples &cs) const
	ESElectronicsSimFast (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)

Private Member Functions
double	decode (const ESSample &sample, const DetId &detId) const
std::vector< ESSample >	fastEncode (const CaloSamples &timeframe, CLHEP::RandGeneral *histoDistribution, double hInf, double hSup, double hBin) const
std::vector< ESSample >	standEncode (const CaloSamples &timeframe) const

Private Attributes
bool	addNoise_
int	gain_
const ESIntercalibConstants *	mips_
double	MIPToGeV_
const ESPedestals *	peds_

Detailed Description

Definition at line 14 of file ESElectronicsSimFast.h.

Member Enumeration Documentation

anonymous enum

Enumerator
MAXADC

Definition at line 18 of file ESElectronicsSimFast.h.

{MAXADC = 4095};

anonymous enum

Enumerator
MINADC

Definition at line 19 of file ESElectronicsSimFast.h.

{MINADC = 0};

Constructor & Destructor Documentation

ESElectronicsSimFast::ESElectronicsSimFast

(

bool

addNoise

)

Definition at line 13 of file ESElectronicsSimFast.cc.

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

Member Function Documentation

void ESElectronicsSimFast::analogToDigital ( const CaloSamples &  cs, ESDataFrame &  df, bool  wasEmpty, CLHEP::RandGeneral *  histoDistribution, double  hInf, double  hSup, double  hBin  ) const

virtual

Definition at line 22 of file ESElectronicsSimFast.cc.

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

Referenced by ESFastTDigitizer::run().

{
  std::vector<ESSample> essamples;
  if (!wasEmpty) essamples = standEncode(cs);
  if ( wasEmpty) essamples = fastEncode(cs, histoDistribution, hInf, hSup, hBin);
  
  df.setSize(cs.size());
  for(int i=0; i<df.size(); i++) {
    df.setSample(i, essamples[i]);
  }
}

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

private

Definition at line 131 of file ESElectronicsSimFast.cc.

Referenced by digitalToAnalog().

{
  return 0. ;
}

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

Definition at line 34 of file ESElectronicsSimFast.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 > ESElectronicsSimFast::fastEncode ( const CaloSamples &  timeframe, CLHEP::RandGeneral *  histoDistribution, double  hInf, double  hSup, double  hBin  ) const

private

Definition at line 93 of file ESElectronicsSimFast.cc.

References ecalMGPA::adc(), newFWLiteAna::bin, ESCondObjectContainer< T >::find(), CaloSamples::id(), MAXADC, MINADC, peds_, python.entryComment::results, CaloSamples::size(), and tablePrinter::width.

Referenced by analogToDigital().

{
  std::vector<ESSample> results;
  results.reserve(timeframe.size());

  int bin[3]; 
  double hBin2 = hBin*hBin;
  double hBin3 = hBin*hBin*hBin;
  double width = (hSup - hInf)/hBin;  

  double thisRnd  = histoDistribution->fire();  
  int thisRndCell = (int)((hBin3)*(thisRnd)/width);  
  bin[2] = (int)(thisRndCell/hBin2);                              // sample2 - bin [0,N-1]
  bin[1] = (int)((thisRndCell - hBin2*bin[2])/hBin);              // sample1
  bin[0] = (int)(thisRndCell - hBin*(bin[1] + hBin*bin[2]));      // sample0

  ESPedestals::const_iterator it_ped = peds_->find(timeframe.id());
  int baseline_  = (int) it_ped->getMean();

  int adc[3];
  double noi[3];
  for(int ii=0; ii<3; ii++){

    noi[ii] = hInf + bin[ii]*width; 
    if (noi[ii]>0)      noi[ii] += 0.5;
    else if (noi[ii]<0) noi[ii] -= 0.5;
    
    adc[ii] = int(noi[ii]) - 1000 + baseline_;      

    if (adc[ii]>MAXADC) adc[ii] = MAXADC;
    if (adc[ii]<MINADC) adc[ii] = MINADC;
    
    results.push_back(ESSample(adc[ii]));
  }
  
  return results;
}

void ESElectronicsSimFast::newEvent ( )

inline

anything that needs to be done once per event

Definition at line 33 of file ESElectronicsSimFast.h.

Referenced by ESFastTDigitizer::run().

{}

void ESElectronicsSimFast::setGain ( const int  gain )

inline

Definition at line 23 of file ESElectronicsSimFast.h.

References gain_.

Referenced by EcalDigiProducer::checkCalibrations().

{ gain_ = gain; }

void ESElectronicsSimFast::setMIPs ( const ESIntercalibConstants *  mips )

inline

Definition at line 25 of file ESElectronicsSimFast.h.

References mips_.

Referenced by EcalDigiProducer::checkCalibrations().

{ mips_ = mips; }

void ESElectronicsSimFast::setMIPToGeV ( const double  MIPToGeV )

inline

Definition at line 26 of file ESElectronicsSimFast.h.

References MIPToGeV_.

Referenced by EcalDigiProducer::checkCalibrations().

{ MIPToGeV_ = MIPToGeV; }

void ESElectronicsSimFast::setPedestals ( const ESPedestals *  peds )

inline

Definition at line 24 of file ESElectronicsSimFast.h.

References peds_.

Referenced by EcalDigiProducer::checkCalibrations().

{ peds_ = peds; } 

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

private

Definition at line 42 of file ESElectronicsSimFast.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")
      << "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.";
  }
  
  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;
}

Member Data Documentation

bool ESElectronicsSimFast::addNoise_

private

Definition at line 37 of file ESElectronicsSimFast.h.

Referenced by standEncode().

int ESElectronicsSimFast::gain_

private

Definition at line 38 of file ESElectronicsSimFast.h.

Referenced by setGain().

const ESIntercalibConstants* ESElectronicsSimFast::mips_

private

Definition at line 40 of file ESElectronicsSimFast.h.

Referenced by setMIPs(), and standEncode().

double ESElectronicsSimFast::MIPToGeV_

private

Definition at line 41 of file ESElectronicsSimFast.h.

Referenced by setMIPToGeV(), and standEncode().

const ESPedestals* ESElectronicsSimFast::peds_

private

Definition at line 39 of file ESElectronicsSimFast.h.

Referenced by fastEncode(), setPedestals(), and standEncode().