#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

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

Static Public Attributes
static int	dfSIZE = 5

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

Private Attributes
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	tdcNbits_

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") ),
   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") ),
   adcLSB_MIP_( adcSaturation_MIP_/std::pow(2.,adcNbits_) ),
   adcThreshold_MIP_( pset.getParameter<double>("adcThreshold_MIP") ),
   toaLSB_ns_( pset.getParameter<double>("toaLSB_ns") ),
   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 coll, input, convertSQLitetoXML_cfg::output, and findQualityFiles::run.

27 { }

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

Definition at line 45 of file BTLElectronicsSim.cc.

References btlPulseShape_, bxTime_, 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(), 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:
       //     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;
 
 
       // --- Fill the time and charge arrays
       const unsigned int ibucket = std::floor( finalToA1/bxTime_ );
       if ( (i+ibucket) >= chargeColl.size() ) continue;
       
       chargeColl[i+ibucket] = Npe*Npe_to_pC_; // the p.e. number is here converted to pC
       
       if ( toa1[i+ibucket] == 0. || (finalToA1-ibucket*bxTime_) < toa1[i+ibucket] ){
     toa1[i+ibucket] = finalToA1 - ibucket*bxTime_;
     toa2[i+ibucket] = finalToA2 - ibucket*bxTime_;
       }
 
     }
 
     //run the shaper to create a new data frame
     BTLDataFrame rawDataFrame( it->first );    
     runTrivialShaper(rawDataFrame,chargeColl,toa1,toa2);
     updateOutput(output,rawDataFrame);
     
   }
     
 }

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

Definition at line 137 of file BTLElectronicsSim.cc.

References ecalMGPA::adc(), adcLSB_MIP_, adcSaturation_MIP_, adcThreshold_MIP_, debug, debug_, createfilelist::int, min(), mps_check::msg, FTLDataFrameT< D, S >::print(), BTLSample::set(), FTLDataFrameT< D, S >::setSample(), 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++) {
 
     if ( chargeColl[it] == 0. ) continue;
 
     //brute force saturation, maybe could to better with an exponential like saturation      
     const uint32_t adc=std::floor( std::min(chargeColl[it],adcSaturation_MIP_) / adcLSB_MIP_ );
     const uint32_t tdc_time1=std::floor( toa1[it] / toaLSB_ns_ );
     const uint32_t tdc_time2=std::floor( toa2[it] / toaLSB_ns_ );
     BTLSample newSample;
     newSample.set(chargeColl[it] > adcThreshold_MIP_,false,tdc_time2,tdc_time1,adc);
     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 191 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 173 of file BTLElectronicsSim.cc.

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

Referenced by run().

                                                              {
   int itIdx(9);
   if(rawDataFrame.size()<=itIdx+2) return;
   
   BTLDataFrame 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 float BTLElectronicsSim::adcLSB_MIP_

private

Definition at line 75 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const uint32_t BTLElectronicsSim::adcNbits_

private

Definition at line 71 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::adcSaturation_MIP_

private

Definition at line 74 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const float BTLElectronicsSim::adcThreshold_MIP_

private

Definition at line 76 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 51 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ChannelTimeOffset_

private

Definition at line 54 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::CorrCoeff_

private

Definition at line 79 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::cosPhi_

private

Definition at line 80 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::DarkCountRate_

private

Definition at line 63 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 49 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

int BTLElectronicsSim::dfSIZE = 5

static

Definition at line 42 of file BTLElectronicsSim.h.

Referenced by updateOutput().

const float BTLElectronicsSim::EnergyThreshold_

private

Definition at line 57 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::Npe_to_pC_

private

Definition at line 67 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::Npe_to_V_

private

Definition at line 68 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ReferencePulseNpe_

private

Definition at line 60 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorDecayTime2_

private

Definition at line 83 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorDecayTime_

private

Definition at line 53 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorRiseTime_

private

Definition at line 52 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::SigmaClock2_

private

Definition at line 87 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::SigmaClock_

private

Definition at line 65 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 64 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::SinglePhotonTimeResolution_

private

Definition at line 62 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::sinPhi_

private

Definition at line 81 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::smearChannelTimeOffset_

private

Definition at line 55 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::tdcNbits_

private

Definition at line 71 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::TimeThreshold1_

private

Definition at line 58 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::TimeThreshold2_

private

Definition at line 59 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::toaLSB_ns_

private

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