#include <BTLElectronicsSim.h>

Public Member Functions
	BTLElectronicsSim (const edm::ParameterSet &pset)

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 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	ScintillatorRiseTime_

const float	SigmaClock2_

const float	SigmaClock_

const float	SigmaElectronicNoise2_

const float	SigmaElectronicNoise_

const float	SinglePhotonTimeResolution_

const float	sinPhi_

const float	smearChannelTimeOffset_

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 )

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")),
   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_),
   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.

References cuy::col, coll, input, convertSQLitetoXML_cfg::output, and writedatasetfile::run.

27 { }

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

Definition at line 48 of file BTLElectronicsSim.cc.

References btlPulseShape_, ChannelTimeOffset_, cosPhi_, DCRxRiseTime_, EnergyThreshold_, f, mps_fire::i, Npe_to_pC_, Npe_to_V_, ReferencePulseNpe_, runTrivialShaper(), ScintillatorDecayTime2_, ScintillatorDecayTime_, sigma2_pe(), SigmaClock2_, SigmaElectronicNoise2_, sinPhi_, smearChannelTimeOffset_, SPTR2_, mathSSE::sqrt(), testBeamMIPTimeRes_, MTDShapeBase::timeAtThr(), create_public_lumi_plots::times, hcalSimParameters_cfi::timeSmearing, TimeThreshold1_, TimeThreshold2_, and updateOutput().

                                                  {
   
   MTDSimHitData chargeColl, toa1, toa2;
 
   for(MTDSimHitDataAccumulator::const_iterator it=input.begin();
       it!=input.end();
       it++) {
 
     chargeColl.fill(0.f);
     toa1.fill(0.f);
     toa2.fill(0.f);
     for(size_t i=0; i<it->second.hit_info[0].size(); i++) {
 
       // --- Fluctuate the total number of photo-electrons
       float Npe = CLHEP::RandPoissonQ::shoot(hre, (it->second).hit_info[0][i]);
       if ( Npe < EnergyThreshold_ ) continue;
 
 
       // --- Get the time of arrival and add a channel time offset
       float finalToA1 = (it->second).hit_info[1][i] + 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];
 
 
       // --- 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[i] = Npe*Npe_to_pC_; // the p.e. number is here converted to pC
       
       toa1[i] = finalToA1;
       toa2[i] = finalToA2;
 
     }
 
     //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);
     
   }
     
 }

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 147 of file BTLElectronicsSim.cc.

References ecalMGPA::adc(), adcBitSaturation_, adcLSB_MIP_, adcThreshold_MIP_, debug, debug_, createfilelist::int, 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_fC_);
 #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 202 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 186 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 78 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const float BTLElectronicsSim::adcLSB_MIP_

private

Definition at line 79 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const uint32_t BTLElectronicsSim::adcNbits_

private

Definition at line 74 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::adcSaturation_MIP_

private

Definition at line 77 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::adcThreshold_MIP_

private

Definition at line 80 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const BTLPulseShape BTLElectronicsSim::btlPulseShape_

private

Definition at line 94 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::bxTime_

private

Definition at line 53 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::ChannelTimeOffset_

private

Definition at line 57 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::CorrCoeff_

private

Definition at line 84 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::cosPhi_

private

Definition at line 85 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::DarkCountRate_

private

Definition at line 66 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::DCRxRiseTime_

private

Definition at line 90 of file BTLElectronicsSim.h.

Referenced by run().

const bool BTLElectronicsSim::debug_

private

Definition at line 51 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

int BTLElectronicsSim::dfSIZE = 2

static

Definition at line 44 of file BTLElectronicsSim.h.

Referenced by updateOutput().

const float BTLElectronicsSim::EnergyThreshold_

private

Definition at line 60 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::Npe_to_pC_

private

Definition at line 70 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::Npe_to_V_

private

Definition at line 71 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ReferencePulseNpe_

private

Definition at line 63 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorDecayTime2_

private

Definition at line 88 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorDecayTime_

private

Definition at line 56 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorRiseTime_

private

Definition at line 55 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::SigmaClock2_

private

Definition at line 92 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::SigmaClock_

private

Definition at line 68 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::SigmaElectronicNoise2_

private

Definition at line 91 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::SigmaElectronicNoise_

private

Definition at line 67 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::SinglePhotonTimeResolution_

private

Definition at line 65 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::sinPhi_

private

Definition at line 86 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::smearChannelTimeOffset_

private

Definition at line 58 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::SPTR2_

private

Definition at line 89 of file BTLElectronicsSim.h.

Referenced by run().

const uint32_t BTLElectronicsSim::tdcBitSaturation_

private

Definition at line 82 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const uint32_t BTLElectronicsSim::tdcNbits_

private

Definition at line 74 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::testBeamMIPTimeRes_

private

Definition at line 54 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::TimeThreshold1_

private

Definition at line 61 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::TimeThreshold2_

private

Definition at line 62 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::toaLSB_ns_

private

Definition at line 81 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