#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 &toa1, const mtd::MTDSimHitData &toa2, 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 ETLPulseShape	etlPulseShape_

const reco::FormulaEvaluator	formulaLandauNoise_

const MTDGeometry *	geom_

const edm::ESGetToken< MTDGeometry, MTDDigiGeometryRecord >	geomToken_

const float	integratedLum_

const float	iThreshold_MIP_

const float	noiseLevel_

const float	referenceChargeColl_

const float	sigmaDistorsion_

const float	sigmaTDC_

const uint32_t	tdcBitSaturation_

const uint32_t	tdcNbits_

const float	toaLSB_ns_

Detailed Description

Definition at line 27 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")),
       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")),
       iThreshold_MIP_(pset.getParameter<double>("iThreshold_MIP")),
       toaLSB_ns_(pset.getParameter<double>("toaLSB_ns")),
       tdcBitSaturation_(std::pow(2, tdcNbits_) - 1),
       referenceChargeColl_(pset.getParameter<double>("referenceChargeColl")),
       noiseLevel_(pset.getParameter<double>("noiseLevel")),
       sigmaDistorsion_(pset.getParameter<double>("sigmaDistorsion")),
       sigmaTDC_(pset.getParameter<double>("sigmaTDC")),
       formulaLandauNoise_(pset.getParameter<std::string>("formulaLandauNoise")) {}

Member Function Documentation

◆ getEvent()

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

inline

Definition at line 31 of file ETLElectronicsSim.h.

31 {}

◆ getEventSetup()

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

Definition at line 30 of file ETLElectronicsSim.cc.

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

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

ETLElectronicsSim::geomToken_

const edm::ESGetToken< MTDGeometry, MTDDigiGeometryRecord > geomToken_

Definition: ETLElectronicsSim.h:49

ETLElectronicsSim::geom_

const MTDGeometry * geom_

Definition: ETLElectronicsSim.h:50

◆ run()

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

Definition at line 32 of file ETLElectronicsSim.cc.

References adcThreshold_MIP_, bxTime_, TauDecayModes::dec, hcalRecHitTable_cff::detId, etlPulseShape_, reco::FormulaEvaluator::evaluate(), Exception, f, MTDShapeBase::fallTime(), formulaLandauNoise_, geom_, mps_fire::i, MTDGeometry::idToDet(), input, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, iThreshold_MIP_, PixelTopology::localX(), PixelTopology::localY(), MTDShapeBase::maximum(), noiseLevel_, CosmicsPD_Skims::radius, DetId::rawId(), referenceChargeColl_, runTrivialShaper(), mtdDigitizer_cfi::sigmaDistorsion, sigmaDistorsion_, mtdDigitizer_cfi::sigmaTDC, sigmaTDC_, mathSSE::sqrt(), MTDShapeBase::timeAtThr(), MTDShapeBase::timeOfMax(), GeomDet::topology(), compareTotals::tot, and updateOutput().

                                                              {
   MTDSimHitData chargeColl, toa1, toa2, tot;
 
   std::vector<double> emptyV;
   std::vector<double> radius(1);
   std::vector<double> fluence(1);
   std::vector<double> chOverMPV(1);
 
   for (MTDSimHitDataAccumulator::const_iterator it = input.begin(); it != input.end(); it++) {
     chargeColl.fill(0.f);
     toa1.fill(0.f);
     toa2.fill(0.f);
     tot.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.);
 
     //Here we are sampling over all the buckets. However for the ETL there should be only one bucket at i = 0
     for (size_t i = 0; i < it->second.hit_info[0].size(); i++) {
       if ((it->second).hit_info[0][i] == 0) {
         continue;
       }
       if ((it->second).hit_info[0][i] < adcThreshold_MIP_)
         continue;
       // time of arrival
       double finalToA = (it->second).hit_info[1][i];
       double finalToC = (it->second).hit_info[1][i];
 
       // fill the time and charge arrays
       const unsigned int ibucket = std::floor(finalToA / bxTime_);
       if ((i + ibucket) >= chargeColl.size())
         continue;
 
       //Calculate the jitter
       double SignalToNoise =
           etlPulseShape_.maximum() * ((it->second).hit_info[0][i] / referenceChargeColl_) / noiseLevel_;
       double sigmaJitter1 = etlPulseShape_.timeOfMax() / SignalToNoise;
       double sigmaJitter2 = (etlPulseShape_.fallTime() - etlPulseShape_.timeOfMax()) / SignalToNoise;
       //Calculate the distorsion
       double sigmaDistorsion = sigmaDistorsion_;
       //Calculate the TDC
       double sigmaTDC = sigmaTDC_;
       //Calculate the Landau Noise
       chOverMPV[0] = (it->second).hit_info[0][i] / (it->second).hit_info[2][i];
       double sigmaLN = formulaLandauNoise_.evaluate(chOverMPV, emptyV);
       double sigmaToA = sqrt(sigmaJitter1 * sigmaJitter1 + sigmaDistorsion * sigmaDistorsion + sigmaTDC * sigmaTDC +
                              sigmaLN * sigmaLN);
       double sigmaToC = sqrt(sigmaJitter2 * sigmaJitter2 + sigmaDistorsion * sigmaDistorsion + sigmaTDC * sigmaTDC +
                              sigmaLN * sigmaLN);
       double smearing1 = 0.0;
       double smearing2 = 0.0;
 
       if (sigmaToA > 0. && sigmaToC > 0.) {
         smearing1 = CLHEP::RandGaussQ::shoot(hre, 0., sigmaToA);
         smearing2 = CLHEP::RandGaussQ::shoot(hre, 0., sigmaToC);
       }
 
       finalToA += smearing1;
       finalToC += smearing1 + smearing2;
 
       std::array<float, 3> times = etlPulseShape_.timeAtThr(
           (it->second).hit_info[0][i] / referenceChargeColl_, iThreshold_MIP_, iThreshold_MIP_);
 
       //The signal is below the threshold
       if (times[0] == 0 && times[1] == 0 && times[2] == 0) {
         continue;
       }
       //The signal is considered to be below the threshold
       finalToA += times[0];
       finalToC += times[2];
       if (finalToA >= finalToC)
         continue;
       chargeColl[i + ibucket] += (it->second).hit_info[0][i];
 
       if (toa1[i + ibucket] == 0. || (finalToA - ibucket * bxTime_) < toa1[i + ibucket]) {
         toa1[i + ibucket] = finalToA - ibucket * bxTime_;
         toa2[i + ibucket] = finalToC - ibucket * bxTime_;
       }
 
       tot[i + ibucket] = finalToC - finalToA;
     }
     // Run the shaper to create a new data frame
     ETLDataFrame rawDataFrame(it->first.detid_);
     runTrivialShaper(rawDataFrame, chargeColl, toa1, tot, it->first.row_, it->first.column_);
     updateOutput(output, rawDataFrame);
   }
 }

◆ runTrivialShaper()

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

Definition at line 129 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_, toaLSB_ns_, and compareTotals::tot.

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. Notice that we are only interested in the first element of the array for the ETL
   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_time1 = std::min((uint32_t)std::floor(toa[it] / toaLSB_ns_), tdcBitSaturation_);
     const uint32_t tdc_time2 = std::min((uint32_t)std::floor(tot[it] / toaLSB_ns_), tdcBitSaturation_);
     //If time over threshold is 0 the event is assumed to not pass the threshold
     bool thres = true;
     if (tdc_time2 == 0 || chargeColl[it] < adcThreshold_MIP_)
       thres = false;
 
     ETLSample newSample;
     newSample.set(thres, false, tdc_time1, tdc_time2, adc, row, col);
 
     //ETLSample newSample;
     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 172 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 64 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

◆ adcLSB_MIP_

const float ETLElectronicsSim::adcLSB_MIP_

private

Definition at line 63 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

◆ adcNbits_

const uint32_t ETLElectronicsSim::adcNbits_

private

Definition at line 59 of file ETLElectronicsSim.h.

◆ adcSaturation_MIP_

const float ETLElectronicsSim::adcSaturation_MIP_

private

Definition at line 62 of file ETLElectronicsSim.h.

◆ adcThreshold_MIP_

const float ETLElectronicsSim::adcThreshold_MIP_

private

Definition at line 65 of file ETLElectronicsSim.h.

Referenced by run(), and runTrivialShaper().

◆ bxTime_

const float ETLElectronicsSim::bxTime_

private

Definition at line 53 of file ETLElectronicsSim.h.

Referenced by run().

◆ debug_

const bool ETLElectronicsSim::debug_

private

Definition at line 52 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

◆ dfSIZE

constexpr int ETLElectronicsSim::dfSIZE = 5

static

Definition at line 46 of file ETLElectronicsSim.h.

Referenced by updateOutput().

◆ etlPulseShape_

const ETLPulseShape ETLElectronicsSim::etlPulseShape_

private

Definition at line 56 of file ETLElectronicsSim.h.

Referenced by run().

◆ formulaLandauNoise_

const reco::FormulaEvaluator ETLElectronicsSim::formulaLandauNoise_

private

Definition at line 73 of file ETLElectronicsSim.h.

Referenced by run().

◆ geom_

const MTDGeometry* ETLElectronicsSim::geom_

private

Definition at line 50 of file ETLElectronicsSim.h.

Referenced by getEventSetup(), and run().

◆ geomToken_

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

private

Definition at line 49 of file ETLElectronicsSim.h.

Referenced by getEventSetup().

◆ integratedLum_

const float ETLElectronicsSim::integratedLum_

private

Definition at line 54 of file ETLElectronicsSim.h.

◆ iThreshold_MIP_

const float ETLElectronicsSim::iThreshold_MIP_

private

Definition at line 66 of file ETLElectronicsSim.h.

Referenced by run().

◆ noiseLevel_

const float ETLElectronicsSim::noiseLevel_

private

Definition at line 70 of file ETLElectronicsSim.h.

Referenced by run().

◆ referenceChargeColl_

const float ETLElectronicsSim::referenceChargeColl_

private

Definition at line 69 of file ETLElectronicsSim.h.

Referenced by run().

◆ sigmaDistorsion_

const float ETLElectronicsSim::sigmaDistorsion_

private

Definition at line 71 of file ETLElectronicsSim.h.

Referenced by run().

◆ sigmaTDC_

const float ETLElectronicsSim::sigmaTDC_

private

Definition at line 72 of file ETLElectronicsSim.h.

Referenced by run().

◆ tdcBitSaturation_

const uint32_t ETLElectronicsSim::tdcBitSaturation_

private

Definition at line 68 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

◆ tdcNbits_

const uint32_t ETLElectronicsSim::tdcNbits_

private

Definition at line 59 of file ETLElectronicsSim.h.

◆ toaLSB_ns_

const float ETLElectronicsSim::toaLSB_ns_

private

Definition at line 67 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

Public Member Functions

Static Public Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ETLElectronicsSim()

Member Function Documentation

◆ getEvent()

◆ getEventSetup()

◆ run()

◆ runTrivialShaper()

◆ updateOutput()

Member Data Documentation

◆ adcBitSaturation_

◆ adcLSB_MIP_

◆ adcNbits_

◆ adcSaturation_MIP_

◆ adcThreshold_MIP_

◆ bxTime_

◆ debug_

◆ dfSIZE

◆ etlPulseShape_

◆ formulaLandauNoise_

◆ geom_

◆ geomToken_

◆ integratedLum_

◆ iThreshold_MIP_

◆ noiseLevel_

◆ referenceChargeColl_

◆ sigmaDistorsion_

◆ sigmaTDC_

◆ tdcBitSaturation_

◆ tdcNbits_

◆ toaLSB_ns_