#include <ETLElectronicsSim.h>

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

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 &toa) const

void	updateOutput (ETLDigiCollection &coll, const ETLDataFrame &rawDataFrame) const

Static Public Attributes
static int	dfSIZE = 5

Private Attributes
const float	adcLSB_MIP_

const uint32_t	adcNbits_

const float	adcSaturation_MIP_

const float	adcThreshold_MIP_

const bool	debug_

const uint32_t	tdcNbits_

const float	toaLSB_ns_

Detailed Description

Definition at line 19 of file ETLElectronicsSim.h.

Constructor & Destructor Documentation

ETLElectronicsSim::ETLElectronicsSim ( const edm::ParameterSet & pset )

Definition at line 7 of file ETLElectronicsSim.cc.

                                                                 :
   debug_( pset.getUntrackedParameter<bool>("debug",false) ),
   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")) {    
 }

Member Function Documentation

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

inline

Definition at line 23 of file ETLElectronicsSim.h.

23 { }

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

inline

Definition at line 25 of file ETLElectronicsSim.h.

References coll, input, convertSQLitetoXML_cfg::output, and findQualityFiles::run.

25 { }

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

Definition at line 18 of file ETLElectronicsSim.cc.

References f, mps_fire::i, runTrivialShaper(), and updateOutput().

                                                  {
   
   MTDSimHitData chargeColl,toa;
   
   for(MTDSimHitDataAccumulator::const_iterator it=input.begin();
       it!=input.end();
       it++) {
     
     chargeColl.fill(0.f); 
     toa.fill(0.f);
     for(size_t i=0; i<it->second.hit_info[0].size(); i++) {
       //time of arrival
       float finalToA = (it->second).hit_info[1][i];
       while(finalToA < 0.f)  finalToA+=25.f;
       while(finalToA > 25.f) finalToA-=25.f;
       toa[i]=finalToA;
       
       // collected charge (in this case in MIPs)
       chargeColl[i] = (it->second).hit_info[0][i];      
     }
     
     //run the shaper to create a new data frame
     ETLDataFrame rawDataFrame( it->first );    
     runTrivialShaper(rawDataFrame,chargeColl,toa);
     updateOutput(output,rawDataFrame);
     
   }
     
 }

void ETLElectronicsSim::runTrivialShaper	(	ETLDataFrame &	dataFrame,
		const mtd::MTDSimHitData &	chargeColl,
		const mtd::MTDSimHitData &	toa
	)		const

Definition at line 51 of file ETLElectronicsSim.cc.

References ecalMGPA::adc(), adcLSB_MIP_, adcSaturation_MIP_, adcThreshold_MIP_, debug, debug_, createfilelist::int, min(), mps_check::msg, FTLDataFrameT< D, S >::print(), ETLSample::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("ETLElectronicsSim") << "[runTrivialShaper]" << std::endl;
   
   //set new ADCs 
   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::floor( std::min(chargeColl[it],adcSaturation_MIP_) / adcLSB_MIP_ );
       const uint32_t tdc_time=std::floor( toa[it] / toaLSB_ns_ );
       ETLSample newSample;
       newSample.set(chargeColl[it] > adcThreshold_MIP_,false,tdc_time,adc);
       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;
   } 
 }

void ETLElectronicsSim::updateOutput	(	ETLDigiCollection &	coll,
		const ETLDataFrame &	rawDataFrame
	)		const

Definition at line 81 of file ETLElectronicsSim.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;
   
   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);    
   }
 }