ETLElectronicsSim Class Reference

#include <ETLElectronicsSim.h>

Public Member Functions
	ETLElectronicsSim (const edm::ParameterSet &pset)
void	getEvent (const edm::Event &evt)
void	getEventSetup (const edm::EventSetup &evt)
void	run (const mtd::MTDSimHitDataAccumulator &input, ETLDigiCollection &output, CLHEP::HepRandomEngine *hre) const
void	runTrivialShaper (ETLDataFrame &dataFrame, const mtd::MTDSimHitData &chargeColl, const mtd::MTDSimHitData &toa, const uint8_t row, const uint8_t column) const
void	updateOutput (ETLDigiCollection &coll, const ETLDataFrame &rawDataFrame) const

Static Public Attributes
static constexpr int	dfSIZE = 5

Private Attributes
const float	adcLSB_MIP_
const uint32_t	adcNbits_
const float	adcSaturation_MIP_
const float	adcThreshold_MIP_
const float	bxTime_
const bool	debug_
const reco::FormulaEvaluator	sigmaEta_
const uint32_t	tdcNbits_
const float	toaLSB_ns_

Detailed Description

Definition at line 21 of file ETLElectronicsSim.h.

Constructor & Destructor Documentation

ETLElectronicsSim()

ETLElectronicsSim::ETLElectronicsSim

(

const edm::ParameterSet &

pset

)

Definition at line 8 of file ETLElectronicsSim.cc.

    : debug_(pset.getUntrackedParameter<bool>("debug", false)),
      bxTime_(pset.getParameter<double>("bxTime")),
      sigmaEta_(pset.getParameter<std::string>("etaResolution")),
      adcNbits_(pset.getParameter<uint32_t>("adcNbits")),
      tdcNbits_(pset.getParameter<uint32_t>("tdcNbits")),
      adcSaturation_MIP_(pset.getParameter<double>("adcSaturation_MIP")),
      adcLSB_MIP_(adcSaturation_MIP_ / std::pow(2., adcNbits_)),
      adcThreshold_MIP_(pset.getParameter<double>("adcThreshold_MIP")),
      toaLSB_ns_(pset.getParameter<double>("toaLSB_ns")) {}

Member Function Documentation

getEvent()

void ETLElectronicsSim::getEvent ( const edm::Event &  evt )

inline

Definition at line 25 of file ETLElectronicsSim.h.

{}

getEventSetup()

void ETLElectronicsSim::getEventSetup ( const edm::EventSetup &  evt )

inline

Definition at line 27 of file ETLElectronicsSim.h.

{}

run()

void ETLElectronicsSim::run	(	const mtd::MTDSimHitDataAccumulator &	input,
		ETLDigiCollection &	output,
		CLHEP::HepRandomEngine *	hre
	)		const

Definition at line 19 of file ETLElectronicsSim.cc.

                                                             {
  MTDSimHitData chargeColl, toa;
 
  std::vector<double> emptyV;
  std::vector<double> eta(1);
 
  for (MTDSimHitDataAccumulator::const_iterator it = input.begin(); it != input.end(); it++) {
    chargeColl.fill(0.f);
    toa.fill(0.f);
    for (size_t i = 0; i < it->second.hit_info[0].size(); i++) {
      if ((it->second).hit_info[0][i] < adcThreshold_MIP_)
        continue;
 
      // time of arrival
      float finalToA = (it->second).hit_info[1][i];
 
      // Gaussian smearing of the time of arrival
      eta[0] = 2.;  // This is just temporary. Once the RECO geometry is
                    // available, the actual module eta will be used.
      double sigmaToA = sigmaEta_.evaluate(eta, emptyV);
 
      if (sigmaToA > 0.)
        finalToA += CLHEP::RandGaussQ::shoot(hre, 0., sigmaToA);
 
      // fill the time and charge arrays
      const unsigned int ibucket = std::floor(finalToA / bxTime_);
      if ((i + ibucket) >= chargeColl.size())
        continue;
 
      chargeColl[i + ibucket] += (it->second).hit_info[0][i];
 
      if (toa[i + ibucket] == 0. || (finalToA - ibucket * bxTime_) < toa[i + ibucket])
        toa[i + ibucket] = finalToA - ibucket * bxTime_;
    }
 
    // run the shaper to create a new data frame
    ETLDataFrame rawDataFrame(it->first.detid_);
    runTrivialShaper(rawDataFrame, chargeColl, toa, it->first.row_, it->first.column_);
    updateOutput(output, rawDataFrame);
  }
}

References adcThreshold_MIP_, bxTime_, PVValHelper::eta, reco::FormulaEvaluator::evaluate(), f, mps_fire::i, input, convertSQLitetoXML_cfg::output, runTrivialShaper(), sigmaEta_, and updateOutput().

runTrivialShaper()

void ETLElectronicsSim::runTrivialShaper	(	ETLDataFrame &	dataFrame,
		const mtd::MTDSimHitData &	chargeColl,
		const mtd::MTDSimHitData &	toa,
		const uint8_t	row,
		const uint8_t	column
	)		const

Definition at line 63 of file ETLElectronicsSim.cc.

                                                                  {
  bool debug = debug_;
#ifdef EDM_ML_DEBUG
  for (int it = 0; it < (int)(chargeColl.size()); it++)
    debug |= (chargeColl[it] > adcThreshold_MIP_);
#endif
 
  if (debug)
    edm::LogVerbatim("ETLElectronicsSim") << "[runTrivialShaper]" << std::endl;
 
  //set new ADCs
  for (int it = 0; it < (int)(chargeColl.size()); it++) {
    //brute force saturation, maybe could to better with an exponential like saturation
    const uint32_t adc = std::floor(std::min(chargeColl[it], adcSaturation_MIP_) / adcLSB_MIP_);
    const uint32_t tdc_time = std::floor(toa[it] / toaLSB_ns_);
    ETLSample newSample;
    newSample.set(chargeColl[it] > adcThreshold_MIP_, false, tdc_time, adc, row, col);
    dataFrame.setSample(it, newSample);
 
    if (debug)
      edm::LogVerbatim("ETLElectronicsSim") << adc << " (" << chargeColl[it] << "/" << adcLSB_MIP_ << ") ";
  }
 
  if (debug) {
    std::ostringstream msg;
    dataFrame.print(msg);
    edm::LogVerbatim("ETLElectronicsSim") << msg.str() << std::endl;
  }
}

References ecalMGPA::adc(), adcLSB_MIP_, adcSaturation_MIP_, adcThreshold_MIP_, cuy::col, debug, debug_, createfilelist::int, min(), mps_check::msg, FTLDataFrameT< D, S, DECODE >::print(), ETLSample::set(), FTLDataFrameT< D, S, DECODE >::setSample(), and toaLSB_ns_.

Referenced by run().

updateOutput()

void ETLElectronicsSim::updateOutput	(	ETLDigiCollection &	coll,
		const ETLDataFrame &	rawDataFrame
	)		const

Definition at line 97 of file ETLElectronicsSim.cc.

                                                                                                    {
  int itIdx(9);
  if (rawDataFrame.size() <= itIdx + 2)
    return;
 
  ETLDataFrame dataFrame(rawDataFrame.id());
  dataFrame.resize(dfSIZE);
  bool putInEvent(false);
  for (int it = 0; it < dfSIZE; ++it) {
    dataFrame.setSample(it, rawDataFrame[itIdx - 2 + it]);
    if (it == 2)
      putInEvent = rawDataFrame[itIdx - 2 + it].threshold();
  }
 
  if (putInEvent) {
    coll.push_back(dataFrame);
  }
}

References dfSIZE, FTLDataFrameT< D, S, DECODE >::id(), edm::SortedCollection< T, SORT >::push_back(), FTLDataFrameT< D, S, DECODE >::resize(), and FTLDataFrameT< D, S, DECODE >::size().

Referenced by run().

Member Data Documentation

adcLSB_MIP_

const float ETLElectronicsSim::adcLSB_MIP_

private

Definition at line 51 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

adcNbits_

const uint32_t ETLElectronicsSim::adcNbits_

private

Definition at line 47 of file ETLElectronicsSim.h.

adcSaturation_MIP_

const float ETLElectronicsSim::adcSaturation_MIP_

private

Definition at line 50 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

adcThreshold_MIP_

const float ETLElectronicsSim::adcThreshold_MIP_

private

Definition at line 52 of file ETLElectronicsSim.h.

Referenced by run(), and runTrivialShaper().

bxTime_

const float ETLElectronicsSim::bxTime_

private

Definition at line 43 of file ETLElectronicsSim.h.

Referenced by run().

debug_

const bool ETLElectronicsSim::debug_

private

Definition at line 42 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

dfSIZE

constexpr int ETLElectronicsSim::dfSIZE = 5

staticconstexpr

Definition at line 39 of file ETLElectronicsSim.h.

Referenced by updateOutput().

sigmaEta_

const reco::FormulaEvaluator ETLElectronicsSim::sigmaEta_

private

Definition at line 44 of file ETLElectronicsSim.h.

Referenced by run().

tdcNbits_

const uint32_t ETLElectronicsSim::tdcNbits_

private

Definition at line 47 of file ETLElectronicsSim.h.

toaLSB_ns_

const float ETLElectronicsSim::toaLSB_ns_

private

Definition at line 53 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().