#include <BTLElectronicsSim.h>

Public Member Functions
	BTLElectronicsSim (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, BTLDigiCollection &output, CLHEP::HepRandomEngine *hre) const

void	runTrivialShaper (BTLDataFrame &dataFrame, const mtd::MTDSimHitData &chargeColl, const mtd::MTDSimHitData &toa1, const mtd::MTDSimHitData &toa2, const uint8_t row, const uint8_t col) const

void	updateOutput (BTLDigiCollection &coll, const BTLDataFrame &rawDataFrame) const

Static Public Attributes
static constexpr int	dfSIZE = 2

Private Member Functions
float	sigma2_pe (const float &Q, const float &R) const

Private Attributes
const uint32_t	adcBitSaturation_

const float	adcLSB_MIP_

const uint32_t	adcNbits_

const float	adcSaturation_MIP_

const float	adcThreshold_MIP_

const BTLPulseShape	btlPulseShape_

const float	bxTime_

const float	ChannelTimeOffset_

const float	CorrCoeff_

const float	cosPhi_

const float	DarkCountRate_

const float	DCRxRiseTime_

const bool	debug_

const float	EnergyThreshold_

const float	Npe_to_pC_

const float	Npe_to_V_

const float	ReferencePulseNpe_

const float	ScintillatorDecayTime2_

const float	ScintillatorDecayTime_

const float	ScintillatorDecayTimeInv_

const float	ScintillatorRiseTime_

const float	SigmaClock2_

const float	SigmaClock_

const float	SigmaElectronicNoise2_

const float	SigmaElectronicNoise_

const float	SinglePhotonTimeResolution_

const float	sinPhi_

const float	smearChannelTimeOffset_

const bool	smearTimeForOOTtails_

const float	SPTR2_

const uint32_t	tdcBitSaturation_

const uint32_t	tdcNbits_

const float	testBeamMIPTimeRes_

const float	TimeThreshold1_

const float	TimeThreshold2_

const float	toaLSB_ns_

Detailed Description

Definition at line 21 of file BTLElectronicsSim.h.

Constructor & Destructor Documentation

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

Definition at line 11 of file BTLElectronicsSim.cc.

     : debug_(pset.getUntrackedParameter<bool>("debug", false)),
       bxTime_(pset.getParameter<double>("bxTime")),
       testBeamMIPTimeRes_(pset.getParameter<double>("TestBeamMIPTimeRes")),
       ScintillatorRiseTime_(pset.getParameter<double>("ScintillatorRiseTime")),
       ScintillatorDecayTime_(pset.getParameter<double>("ScintillatorDecayTime")),
       ChannelTimeOffset_(pset.getParameter<double>("ChannelTimeOffset")),
       smearChannelTimeOffset_(pset.getParameter<double>("smearChannelTimeOffset")),
       EnergyThreshold_(pset.getParameter<double>("EnergyThreshold")),
       TimeThreshold1_(pset.getParameter<double>("TimeThreshold1")),
       TimeThreshold2_(pset.getParameter<double>("TimeThreshold2")),
       ReferencePulseNpe_(pset.getParameter<double>("ReferencePulseNpe")),
       SinglePhotonTimeResolution_(pset.getParameter<double>("SinglePhotonTimeResolution")),
       DarkCountRate_(pset.getParameter<double>("DarkCountRate")),
       SigmaElectronicNoise_(pset.getParameter<double>("SigmaElectronicNoise")),
       SigmaClock_(pset.getParameter<double>("SigmaClock")),
       smearTimeForOOTtails_(pset.getParameter<bool>("SmearTimeForOOTtails")),
       Npe_to_pC_(pset.getParameter<double>("Npe_to_pC")),
       Npe_to_V_(pset.getParameter<double>("Npe_to_V")),
       adcNbits_(pset.getParameter<uint32_t>("adcNbits")),
       tdcNbits_(pset.getParameter<uint32_t>("tdcNbits")),
       adcSaturation_MIP_(pset.getParameter<double>("adcSaturation_MIP")),
       adcBitSaturation_(std::pow(2, adcNbits_) - 1),
       adcLSB_MIP_(adcSaturation_MIP_ / adcBitSaturation_),
       adcThreshold_MIP_(pset.getParameter<double>("adcThreshold_MIP")),
       toaLSB_ns_(pset.getParameter<double>("toaLSB_ns")),
       tdcBitSaturation_(std::pow(2, tdcNbits_) - 1),
       CorrCoeff_(pset.getParameter<double>("CorrelationCoefficient")),
       cosPhi_(0.5 * (sqrt(1. + CorrCoeff_) + sqrt(1. - CorrCoeff_))),
       sinPhi_(0.5 * CorrCoeff_ / cosPhi_),
       ScintillatorDecayTime2_(ScintillatorDecayTime_ * ScintillatorDecayTime_),
       ScintillatorDecayTimeInv_(1. / ScintillatorDecayTime_),
       SPTR2_(SinglePhotonTimeResolution_ * SinglePhotonTimeResolution_),
       DCRxRiseTime_(DarkCountRate_ * ScintillatorRiseTime_),
       SigmaElectronicNoise2_(SigmaElectronicNoise_ * SigmaElectronicNoise_),
       SigmaClock2_(SigmaClock_ * SigmaClock_) {}

Member Function Documentation

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

inline

Definition at line 25 of file BTLElectronicsSim.h.

25 {}

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

inline

Definition at line 27 of file BTLElectronicsSim.h.

27 {}

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

Definition at line 48 of file BTLElectronicsSim.cc.

References btlPulseShape_, bxTime_, ChannelTimeOffset_, cosPhi_, DCRxRiseTime_, EnergyThreshold_, funct::exp(), validate-o2o-wbm::f, mtd_digitizer::kInTimeBX, Npe_to_pC_, Npe_to_V_, ReferencePulseNpe_, runTrivialShaper(), ScintillatorDecayTime2_, ScintillatorDecayTime_, ScintillatorDecayTimeInv_, ScintillatorRiseTime_, sigma2_pe(), SigmaClock2_, SigmaElectronicNoise2_, sinPhi_, smearChannelTimeOffset_, smearTimeForOOTtails_, SPTR2_, mathSSE::sqrt(), testBeamMIPTimeRes_, MTDShapeBase::timeAtThr(), TimeThreshold1_, TimeThreshold2_, and updateOutput().

                                                              {
   MTDSimHitData chargeColl, toa1, toa2;
 
   for (MTDSimHitDataAccumulator::const_iterator it = input.begin(); it != input.end(); it++) {
     // --- Digitize only the in-time bucket:
     const unsigned int iBX = mtd_digitizer::kInTimeBX;
 
     chargeColl.fill(0.f);
     toa1.fill(0.f);
     toa2.fill(0.f);
     for (size_t iside = 0; iside < 2; iside++) {
       // --- Fluctuate the total number of photo-electrons
       float Npe = CLHEP::RandPoissonQ::shoot(hre, (it->second).hit_info[2 * iside][iBX]);
       if (Npe < EnergyThreshold_)
         continue;
 
       // --- Get the time of arrival and add a channel time offset
       float finalToA1 = (it->second).hit_info[1 + 2 * iside][iBX] + ChannelTimeOffset_;
 
       if (smearChannelTimeOffset_ > 0.) {
         float timeSmearing = CLHEP::RandGaussQ::shoot(hre, 0., smearChannelTimeOffset_);
         finalToA1 += timeSmearing;
       }
 
       // --- Calculate and add the time walk: the time of arrival is read in correspondence
       //                                      with two thresholds on the signal pulse
       std::array<float, 3> times =
           btlPulseShape_.timeAtThr(Npe / ReferencePulseNpe_, TimeThreshold1_ * Npe_to_V_, TimeThreshold2_ * Npe_to_V_);
 
       // --- If the pulse amplitude is smaller than TimeThreshold2, the trigger does not fire
       if (times[1] == 0.)
         continue;
 
       float finalToA2 = finalToA1 + times[1];
       finalToA1 += times[0];
 
       // --- Estimate the time uncertainty due to photons from earlier OOT hits in the current BTL cell
       if (smearTimeForOOTtails_) {
         float rate_oot = 0.;
         // Loop on earlier OOT hits
         for (int ibx = 0; ibx < mtd_digitizer::kInTimeBX; ++ibx) {
           if ((it->second).hit_info[2 * iside][ibx] > 0.) {
             float hit_time = (it->second).hit_info[1 + 2 * iside][ibx] + bxTime_ * (ibx - mtd_digitizer::kInTimeBX);
             float npe_oot = CLHEP::RandPoissonQ::shoot(hre, (it->second).hit_info[2 * iside][ibx]);
             rate_oot += npe_oot * exp(hit_time * ScintillatorDecayTimeInv_) * ScintillatorDecayTimeInv_;
           }
         }  // ibx loop
 
         if (rate_oot > 0.) {
           float sigma_oot = sqrt(rate_oot * ScintillatorRiseTime_) * ScintillatorDecayTime_ / Npe;
           float smearing_oot = CLHEP::RandGaussQ::shoot(hre, 0., sigma_oot);
           finalToA1 += smearing_oot;
           finalToA2 += smearing_oot;
         }
       }  // if smearTimeForOOTtails_
 
       // --- Uncertainty due to the fluctuations of the n-th photon arrival time:
       if (testBeamMIPTimeRes_ > 0.) {
         // In this case the time resolution is parametrized from the testbeam.
         // The same parameterization is used for both thresholds.
         float sigma = testBeamMIPTimeRes_ / sqrt(Npe);
         float smearing_stat_thr1 = CLHEP::RandGaussQ::shoot(hre, 0., sigma);
         float smearing_stat_thr2 = CLHEP::RandGaussQ::shoot(hre, 0., sigma);
 
         finalToA1 += smearing_stat_thr1;
         finalToA2 += smearing_stat_thr2;
 
       } else {
         // In this case the time resolution is taken from the literature.
         // The fluctuations due to the first TimeThreshold1_ p.e. are common to both times
         float smearing_stat_thr1 =
             CLHEP::RandGaussQ::shoot(hre, 0., ScintillatorDecayTime_ * sqrt(sigma2_pe(TimeThreshold1_, Npe)));
         float smearing_stat_thr2 = CLHEP::RandGaussQ::shoot(
             hre, 0., ScintillatorDecayTime_ * sqrt(sigma2_pe(TimeThreshold2_ - TimeThreshold1_, Npe)));
         finalToA1 += smearing_stat_thr1;
         finalToA2 += smearing_stat_thr1 + smearing_stat_thr2;
       }
 
       // --- Add in quadrature the uncertainties due to the SiPM timing resolution, the SiPM DCR,
       //     the electronic noise and the clock distribution:
       float slew2 = ScintillatorDecayTime2_ / Npe / Npe;
 
       float sigma2_tot_thr1 =
           SPTR2_ / TimeThreshold1_ + (DCRxRiseTime_ + SigmaElectronicNoise2_) * slew2 + SigmaClock2_;
       float sigma2_tot_thr2 =
           SPTR2_ / TimeThreshold2_ + (DCRxRiseTime_ + SigmaElectronicNoise2_) * slew2 + SigmaClock2_;
 
       // --- Smear the arrival times using the correlated uncertainties:
       float smearing_thr1_uncorr = CLHEP::RandGaussQ::shoot(hre, 0., sqrt(sigma2_tot_thr1));
       float smearing_thr2_uncorr = CLHEP::RandGaussQ::shoot(hre, 0., sqrt(sigma2_tot_thr2));
 
       finalToA1 += cosPhi_ * smearing_thr1_uncorr + sinPhi_ * smearing_thr2_uncorr;
       finalToA2 += sinPhi_ * smearing_thr1_uncorr + cosPhi_ * smearing_thr2_uncorr;
 
       chargeColl[iside] = Npe * Npe_to_pC_;  // the p.e. number is here converted to pC
 
       toa1[iside] = finalToA1;
       toa2[iside] = finalToA2;
 
     }  // iside loop
 
     //run the shaper to create a new data frame
     BTLDataFrame rawDataFrame(it->first.detid_);
     runTrivialShaper(rawDataFrame, chargeColl, toa1, toa2, it->first.row_, it->first.column_);
     updateOutput(output, rawDataFrame);
 
   }  // MTDSimHitDataAccumulator loop
 }

void BTLElectronicsSim::runTrivialShaper	(	BTLDataFrame &	dataFrame,
		const mtd::MTDSimHitData &	chargeColl,
		const mtd::MTDSimHitData &	toa1,
		const mtd::MTDSimHitData &	toa2,
		const uint8_t	row,
		const uint8_t	col
	)		const

Definition at line 159 of file BTLElectronicsSim.cc.

References gpuClustering::adc, adcBitSaturation_, adcLSB_MIP_, adcThreshold_MIP_, debug, debug_, SiStripPI::min, mps_check::msg, FTLDataFrameT< D, S, DECODE >::print(), BTLSample::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("BTLElectronicsSim") << "[runTrivialShaper]" << std::endl;
 
   //set new ADCs
   for (int it = 0; it < (int)(chargeColl.size()); it++) {
     BTLSample newSample;
     newSample.set(false, false, 0, 0, 0, row, col);
 
     //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(toa1[it] / toaLSB_ns_), tdcBitSaturation_);
     const uint32_t tdc_time2 = std::min((uint32_t)std::floor(toa2[it] / toaLSB_ns_), tdcBitSaturation_);
 
     newSample.set(
         chargeColl[it] > adcThreshold_MIP_, tdc_time1 == tdcBitSaturation_, tdc_time2, tdc_time1, adc, row, col);
     dataFrame.setSample(it, newSample);
 
     if (debug)
       edm::LogVerbatim("BTLElectronicsSim") << adc << " (" << chargeColl[it] << "/" << adcLSB_MIP_ << ") ";
   }
 
   if (debug) {
     std::ostringstream msg;
     dataFrame.print(msg);
     edm::LogVerbatim("BTLElectronicsSim") << msg.str() << std::endl;
   }
 }

float BTLElectronicsSim::sigma2_pe	(	const float &	Q,
		const float &	R
	)		const

private

Definition at line 214 of file BTLElectronicsSim.cc.

References dttmaxenums::R.

Referenced by run().

                                                                        {
   float OneOverR = 1. / R;
   float OneOverR2 = OneOverR * OneOverR;
 
   // --- This is Eq. (17) from Nucl. Instr. Meth. A 564 (2006) 185
   float sigma2 = Q * OneOverR2 *
                  (1. + 2. * (Q + 1.) * OneOverR + (Q + 1.) * (6. * Q + 11) * OneOverR2 +
                   (Q + 1.) * (Q + 2.) * (2. * Q + 5.) * OneOverR2 * OneOverR);
 
   return sigma2;
 }

void BTLElectronicsSim::updateOutput	(	BTLDigiCollection &	coll,
		const BTLDataFrame &	rawDataFrame
	)		const

Definition at line 199 of file BTLElectronicsSim.cc.

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

Referenced by run().

                                                                                                     {
   BTLDataFrame dataFrame(rawDataFrame.id());
   dataFrame.resize(dfSIZE);
   bool putInEvent(false);
   for (int it = 0; it < dfSIZE; ++it) {
     dataFrame.setSample(it, rawDataFrame[it]);
     if (it == 0)
       putInEvent = rawDataFrame[it].threshold();
   }
 
   if (putInEvent) {
     coll.push_back(dataFrame);
   }
 }

Member Data Documentation

const uint32_t BTLElectronicsSim::adcBitSaturation_

private

Definition at line 72 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const float BTLElectronicsSim::adcLSB_MIP_

private

Definition at line 73 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const uint32_t BTLElectronicsSim::adcNbits_

private

Definition at line 68 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::adcSaturation_MIP_

private

Definition at line 71 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::adcThreshold_MIP_

private

Definition at line 74 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const BTLPulseShape BTLElectronicsSim::btlPulseShape_

private

Definition at line 89 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::bxTime_

private

Definition at line 47 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ChannelTimeOffset_

private

Definition at line 51 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::CorrCoeff_

private

Definition at line 78 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::cosPhi_

private

Definition at line 79 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::DarkCountRate_

private

Definition at line 60 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::DCRxRiseTime_

private

Definition at line 85 of file BTLElectronicsSim.h.

Referenced by run().

const bool BTLElectronicsSim::debug_

private

Definition at line 45 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

constexpr int BTLElectronicsSim::dfSIZE = 2

static

Definition at line 40 of file BTLElectronicsSim.h.

Referenced by updateOutput().

const float BTLElectronicsSim::EnergyThreshold_

private

Definition at line 54 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::Npe_to_pC_

private

Definition at line 64 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::Npe_to_V_

private

Definition at line 65 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ReferencePulseNpe_

private

Definition at line 57 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorDecayTime2_

private

Definition at line 82 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorDecayTime_

private

Definition at line 50 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorDecayTimeInv_

private

Definition at line 83 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorRiseTime_

private

Definition at line 49 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::SigmaClock2_

private

Definition at line 87 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::SigmaClock_

private

Definition at line 62 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::SigmaElectronicNoise2_

private

Definition at line 86 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::SigmaElectronicNoise_

private

Definition at line 61 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::SinglePhotonTimeResolution_

private

Definition at line 59 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::sinPhi_

private

Definition at line 80 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::smearChannelTimeOffset_

private

Definition at line 52 of file BTLElectronicsSim.h.

Referenced by run().

const bool BTLElectronicsSim::smearTimeForOOTtails_

private

Definition at line 63 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::SPTR2_

private

Definition at line 84 of file BTLElectronicsSim.h.

Referenced by run().

const uint32_t BTLElectronicsSim::tdcBitSaturation_

private

Definition at line 76 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const uint32_t BTLElectronicsSim::tdcNbits_

private

Definition at line 68 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::testBeamMIPTimeRes_

private

Definition at line 48 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::TimeThreshold1_

private

Definition at line 55 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::TimeThreshold2_

private

Definition at line 56 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::toaLSB_ns_

private

Definition at line 75 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

Public Member Functions

Static Public Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation