ETLElectronicsSim Class Reference

#include <ETLElectronicsSim.h>

Public Member Functions
	ETLElectronicsSim (const edm::ParameterSet &pset, edm::ConsumesCollector iC)
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 uint32_t	adcBitSaturation_
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	fluence_
const MTDGeometry *	geom_
const edm::ESGetToken< MTDGeometry, MTDDigiGeometryRecord >	geomToken_
const float	integratedLum_
const reco::FormulaEvaluator	lgadGain_
const uint32_t	tdcBitSaturation_
const uint32_t	tdcNbits_
const reco::FormulaEvaluator	timeRes2_
const float	toaLSB_ns_

Detailed Description

Definition at line 25 of file ETLElectronicsSim.h.

Constructor & Destructor Documentation

ETLElectronicsSim()

ETLElectronicsSim::ETLElectronicsSim	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector	iC
	)

Definition at line 9 of file ETLElectronicsSim.cc.

    : geomToken_(iC.esConsumes()),
      geom_(nullptr),
      debug_(pset.getUntrackedParameter<bool>("debug", false)),
      bxTime_(pset.getParameter<double>("bxTime")),
      integratedLum_(pset.getParameter<double>("IntegratedLuminosity")),
      fluence_(pset.getParameter<std::string>("FluenceVsRadius")),
      lgadGain_(pset.getParameter<std::string>("LGADGainVsFluence")),
      timeRes2_(pset.getParameter<std::string>("TimeResolution2")),
      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_)),
      adcBitSaturation_(std::pow(2, adcNbits_) - 1),
      adcThreshold_MIP_(pset.getParameter<double>("adcThreshold_MIP")),
      toaLSB_ns_(pset.getParameter<double>("toaLSB_ns")),
      tdcBitSaturation_(std::pow(2, tdcNbits_) - 1) {}

Member Function Documentation

getEvent()

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

inline

Definition at line 29 of file ETLElectronicsSim.h.

{}

getEventSetup()

void ETLElectronicsSim::getEventSetup

(

const edm::EventSetup &

evt

)

Definition at line 27 of file ETLElectronicsSim.cc.

References geom_, geomToken_, and edm::EventSetup::getData().

{ geom_ = &evs.getData(geomToken_); }

run()

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

Definition at line 29 of file ETLElectronicsSim.cc.

References adcThreshold_MIP_, bxTime_, TauDecayModes::dec, hcalRecHitTable_cff::detId, reco::FormulaEvaluator::evaluate(), Exception, f, fluence_, PedestalClient_cfi::gain, geom_, mps_fire::i, MTDGeometry::idToDet(), input, integratedLum_, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, lgadGain_, PixelTopology::localX(), PixelTopology::localY(), CosmicsPD_Skims::radius, DetId::rawId(), runTrivialShaper(), mathSSE::sqrt(), timeRes2_, GeomDet::toGlobal(), GeomDet::topology(), and updateOutput().

                                                             {
  MTDSimHitData chargeColl, toa;

  std::vector<double> emptyV;
  std::vector<double> radius(1);
  std::vector<double> fluence(1);
  std::vector<double> gain(1);

  for (MTDSimHitDataAccumulator::const_iterator it = input.begin(); it != input.end(); it++) {
    chargeColl.fill(0.f);
    toa.fill(0.f);

    ETLDetId detId = it->first.detid_;
    DetId geoId = detId.geographicalId();
    const MTDGeomDet* thedet = geom_->idToDet(geoId);
    if (thedet == nullptr)
      throw cms::Exception("EtlElectronicsSim") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                                << detId.rawId() << ") is invalid!" << std::dec << std::endl;
    const PixelTopology& topo = static_cast<const PixelTopology&>(thedet->topology());

    Local3DPoint local_point(topo.localX(it->first.row_), topo.localY(it->first.column_), 0.);
    const auto& global_point = thedet->toGlobal(local_point);

    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];

      // calculate the LGAD gain as a function of the fluence at R = radius
      radius[0] = global_point.perp();
      fluence[0] = integratedLum_ * fluence_.evaluate(radius, emptyV);
      gain[0] = lgadGain_.evaluate(fluence, emptyV);
      if (gain[0] <= 0.)
        throw cms::Exception("EtlElectronicsSim") << "Null or negative LGAD gain!" << std::endl;

      // Gaussian smearing of the time of arrival
      double sigmaToA = sqrt(timeRes2_.evaluate(gain, 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);
  }
}

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 92 of file ETLElectronicsSim.cc.

References gpuClustering::adc, adcBitSaturation_, adcLSB_MIP_, adcThreshold_MIP_, cuy::col, debug, debug_, createfilelist::int, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, SiStripPI::min, mps_check::msg, FTLDataFrameT< D, S, DECODE >::print(), ETLSample::set(), FTLDataFrameT< D, S, DECODE >::setSample(), tdcBitSaturation_, and toaLSB_ns_.

Referenced by run().

                                                                  {
  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::min((uint32_t)std::floor(chargeColl[it] / adcLSB_MIP_), adcBitSaturation_);
    const uint32_t tdc_time = std::min((uint32_t)std::floor(toa[it] / toaLSB_ns_), tdcBitSaturation_);
    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;
  }
}

updateOutput()

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

Definition at line 126 of file ETLElectronicsSim.cc.

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

Referenced by run().

                                                                                                    {
  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);
  }
}

Member Data Documentation

adcBitSaturation_

const uint32_t ETLElectronicsSim::adcBitSaturation_

private

Definition at line 62 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

adcLSB_MIP_

const float ETLElectronicsSim::adcLSB_MIP_

private

Definition at line 61 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

adcNbits_

const uint32_t ETLElectronicsSim::adcNbits_

private

Definition at line 57 of file ETLElectronicsSim.h.

adcSaturation_MIP_

const float ETLElectronicsSim::adcSaturation_MIP_

private

Definition at line 60 of file ETLElectronicsSim.h.

adcThreshold_MIP_

const float ETLElectronicsSim::adcThreshold_MIP_

private

Definition at line 63 of file ETLElectronicsSim.h.

Referenced by run(), and runTrivialShaper().

bxTime_

const float ETLElectronicsSim::bxTime_

private

Definition at line 50 of file ETLElectronicsSim.h.

Referenced by run().

debug_

const bool ETLElectronicsSim::debug_

private

Definition at line 49 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

dfSIZE

constexpr int ETLElectronicsSim::dfSIZE = 5

static

Definition at line 43 of file ETLElectronicsSim.h.

Referenced by updateOutput().

fluence_

const reco::FormulaEvaluator ETLElectronicsSim::fluence_

private

Definition at line 52 of file ETLElectronicsSim.h.

Referenced by run().

geom_

const MTDGeometry* ETLElectronicsSim::geom_

private

Definition at line 47 of file ETLElectronicsSim.h.

Referenced by getEventSetup(), and run().

geomToken_

const edm::ESGetToken<MTDGeometry, MTDDigiGeometryRecord> ETLElectronicsSim::geomToken_

private

Definition at line 46 of file ETLElectronicsSim.h.

Referenced by getEventSetup().

integratedLum_

const float ETLElectronicsSim::integratedLum_

private

Definition at line 51 of file ETLElectronicsSim.h.

Referenced by run().

lgadGain_

const reco::FormulaEvaluator ETLElectronicsSim::lgadGain_

private

Definition at line 53 of file ETLElectronicsSim.h.

Referenced by run().

tdcBitSaturation_

const uint32_t ETLElectronicsSim::tdcBitSaturation_

private

Definition at line 65 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

tdcNbits_

const uint32_t ETLElectronicsSim::tdcNbits_

private

Definition at line 57 of file ETLElectronicsSim.h.

timeRes2_

const reco::FormulaEvaluator ETLElectronicsSim::timeRes2_

private

Definition at line 54 of file ETLElectronicsSim.h.

Referenced by run().

toaLSB_ns_

const float ETLElectronicsSim::toaLSB_ns_

private

Definition at line 64 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().