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

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	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) ),
   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") ),
   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 41 of file BTLElectronicsSim.cc.

References btlPulseShape_, ChannelTimeOffset_, DCRxRiseTime_, EnergyThreshold_, f, mps_fire::i, Npe_to_pC_, Npe_to_V_, ReferencePulseNpe_, runTrivialShaper(), ScintillatorDecayTime2_, sigma2_pe(), SigmaClock2_, SigmaElectronicNoise2_, 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:
       float sigma2_tot_thr1 = ScintillatorDecayTime2_*sigma2_pe(TimeThreshold1_,Npe);
       float sigma2_tot_thr2 = ScintillatorDecayTime2_*sigma2_pe(TimeThreshold2_,Npe);
 
 
       // --- Add in quadrature the uncertainty due to the SiPM timing resolution:
       sigma2_tot_thr1 += SPTR2_/TimeThreshold1_;
       sigma2_tot_thr2 += SPTR2_/TimeThreshold2_;
 
 
       // --- Add in quadrature the uncertainties due to the SiPM DCR, 
       //     the electronic noise and the clock distribution:
       float slew2 = ScintillatorDecayTime2_/Npe/Npe;
 
       sigma2_tot_thr1 += ( (DCRxRiseTime_ + SigmaElectronicNoise2_)*slew2 + SigmaClock2_ );
       sigma2_tot_thr2 += ( (DCRxRiseTime_ + SigmaElectronicNoise2_)*slew2 + SigmaClock2_ );
 
 
       // --- Smear the arrival times using the total uncertainty:
       finalToA1 += CLHEP::RandGaussQ::shoot(hre, 0., sqrt(sigma2_tot_thr1));
       finalToA2 += CLHEP::RandGaussQ::shoot(hre, 0., sqrt(sigma2_tot_thr2));
 
 
       while(finalToA1 < 0.f)  finalToA1 += 25.f;
       while(finalToA1 > 25.f) finalToA1 -= 25.f;
       toa1[i]=finalToA1;
 
       while(finalToA2 < 0.f)  finalToA2 += 25.f;
       while(finalToA2 > 25.f) finalToA2 -= 25.f;
       toa2[i]=finalToA2;
 
       chargeColl[i] = Npe*Npe_to_pC_;
       
     }
 
     //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 129 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 183 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 165 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 74 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const uint32_t BTLElectronicsSim::adcNbits_

private

Definition at line 70 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::adcSaturation_MIP_

private

Definition at line 73 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const float BTLElectronicsSim::adcThreshold_MIP_

private

Definition at line 75 of file BTLElectronicsSim.h.

Referenced by runTrivialShaper().

const BTLPulseShape BTLElectronicsSim::btlPulseShape_

private

Definition at line 84 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ChannelTimeOffset_

private

Definition at line 53 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::DarkCountRate_

private

Definition at line 62 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::DCRxRiseTime_

private

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

Referenced by run().

const float BTLElectronicsSim::Npe_to_pC_

private

Definition at line 66 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::Npe_to_V_

private

Definition at line 67 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ReferencePulseNpe_

private

Definition at line 59 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorDecayTime2_

private

Definition at line 78 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::ScintillatorDecayTime_

private

Definition at line 52 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::ScintillatorRiseTime_

private

Definition at line 51 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::SigmaClock2_

private

Definition at line 82 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::SigmaClock_

private

Definition at line 64 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::SigmaElectronicNoise2_

private

Definition at line 81 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::SigmaElectronicNoise_

private

Definition at line 63 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::SinglePhotonTimeResolution_

private

Definition at line 61 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::smearChannelTimeOffset_

private

Definition at line 54 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::SPTR2_

private

Definition at line 79 of file BTLElectronicsSim.h.

Referenced by run().

const uint32_t BTLElectronicsSim::tdcNbits_

private

Definition at line 70 of file BTLElectronicsSim.h.

const float BTLElectronicsSim::TimeThreshold1_

private

Definition at line 57 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::TimeThreshold2_

private

Definition at line 58 of file BTLElectronicsSim.h.

Referenced by run().

const float BTLElectronicsSim::toaLSB_ns_

private

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