#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	(	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

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

inline

Definition at line 29 of file ETLElectronicsSim.h.

29 {}

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

Definition at line 27 of file ETLElectronicsSim.cc.

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

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

ETLElectronicsSim::geomToken_

const edm::ESGetToken< MTDGeometry, MTDDigiGeometryRecord > geomToken_

Definition: ETLElectronicsSim.h:46

ETLElectronicsSim::geom_

const MTDGeometry * geom_

Definition: ETLElectronicsSim.h:47

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, reco::FormulaEvaluator::evaluate(), Exception, validate-o2o-wbm::f, fluence_, ETLDetId::geographicalId(), geom_, mps_fire::i, MTDGeometry::idToDet(), integratedLum_, lgadGain_, PixelTopology::localX(), PixelTopology::localY(), PV3DBase< T, PVType, FrameType >::perp(), 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);
   }
 }

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_, debug, debug_, 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;
   }
 }

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

Definition at line 126 of file ETLElectronicsSim.cc.

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().

                                                                                                     {
   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

const uint32_t ETLElectronicsSim::adcBitSaturation_

private

Definition at line 62 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

const float ETLElectronicsSim::adcLSB_MIP_

private

Definition at line 61 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

const uint32_t ETLElectronicsSim::adcNbits_

private

Definition at line 57 of file ETLElectronicsSim.h.

const float ETLElectronicsSim::adcSaturation_MIP_

private

Definition at line 60 of file ETLElectronicsSim.h.

const float ETLElectronicsSim::adcThreshold_MIP_

private

Definition at line 63 of file ETLElectronicsSim.h.

Referenced by run(), and runTrivialShaper().

const float ETLElectronicsSim::bxTime_

private

Definition at line 50 of file ETLElectronicsSim.h.

Referenced by run().

const bool ETLElectronicsSim::debug_

private

Definition at line 49 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

constexpr int ETLElectronicsSim::dfSIZE = 5

static

Definition at line 43 of file ETLElectronicsSim.h.

Referenced by updateOutput().

const reco::FormulaEvaluator ETLElectronicsSim::fluence_

private

Definition at line 52 of file ETLElectronicsSim.h.

Referenced by run().

const MTDGeometry* ETLElectronicsSim::geom_

private

Definition at line 47 of file ETLElectronicsSim.h.

Referenced by getEventSetup(), and run().

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

private

Definition at line 46 of file ETLElectronicsSim.h.

Referenced by getEventSetup().

const float ETLElectronicsSim::integratedLum_

private

Definition at line 51 of file ETLElectronicsSim.h.

Referenced by run().

const reco::FormulaEvaluator ETLElectronicsSim::lgadGain_

private

Definition at line 53 of file ETLElectronicsSim.h.

Referenced by run().

const uint32_t ETLElectronicsSim::tdcBitSaturation_

private

Definition at line 65 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

const uint32_t ETLElectronicsSim::tdcNbits_

private

Definition at line 57 of file ETLElectronicsSim.h.

const reco::FormulaEvaluator ETLElectronicsSim::timeRes2_

private

Definition at line 54 of file ETLElectronicsSim.h.

Referenced by run().

const float ETLElectronicsSim::toaLSB_ns_

private

Definition at line 64 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

Public Member Functions

Static Public Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation