#include <SimpleElectronicsSimInMIPs.h>

Public Member Functions
void	getEvent (const edm::Event &evt)

void	getEventSetup (const edm::EventSetup &evt)

void	run (const ftl::FTLSimHitDataAccumulator &input, FTLDigiCollection &output) const

void	runTrivialShaper (FTLDataFrame &dataFrame, const ftl::FTLSimHitData &chargeColl, const ftl::FTLSimHitData &toa) const

	SimpleElectronicsSimInMIPs (const edm::ParameterSet &pset)

void	updateOutput (FTLDigiCollection &coll, const FTLDataFrame &rawDataFrame) const

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 15 of file SimpleElectronicsSimInMIPs.h.

Constructor & Destructor Documentation

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

Definition at line 7 of file SimpleElectronicsSimInMIPs.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 SimpleElectronicsSimInMIPs::getEvent ( const edm::Event & evt )

inline

Definition at line 19 of file SimpleElectronicsSimInMIPs.h.

19 { }

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

inline

Definition at line 21 of file SimpleElectronicsSimInMIPs.h.

References coll, input, convertSQLitetoXML_cfg::output, run(), runTrivialShaper(), and updateOutput().

21 { }

void SimpleElectronicsSimInMIPs::run	(	const ftl::FTLSimHitDataAccumulator &	input,
		FTLDigiCollection &	output
	)		const

Definition at line 17 of file SimpleElectronicsSimInMIPs.cc.

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

Referenced by getEventSetup().

                                                       {
   
   FTLSimHitData chargeColl,toa;
   
   for(FTLSimHitDataAccumulator::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
     FTLDataFrame rawDataFrame( it->first );    
     runTrivialShaper(rawDataFrame,chargeColl,toa);
     updateOutput(output,rawDataFrame);
     
   }
     
 }

void SimpleElectronicsSimInMIPs::runTrivialShaper	(	FTLDataFrame &	dataFrame,
		const ftl::FTLSimHitData &	chargeColl,
		const ftl::FTLSimHitData &	toa
	)		const

Definition at line 49 of file SimpleElectronicsSimInMIPs.cc.

References ecalMGPA::adc(), adcLSB_MIP_, adcSaturation_MIP_, adcThreshold_MIP_, debug, debug_, createfilelist::int, min(), mps_alisetup::msg, FTLDataFrameT< D, S >::print(), FTLSample::set(), FTLDataFrameT< D, S >::setSample(), and toaLSB_ns_.

Referenced by getEventSetup(), and 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("FTLSimpleElectronicsSimInMIPs") << "[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_ );
       FTLSample newSample;
       newSample.set(chargeColl[it] > adcThreshold_MIP_,false,tdc_time,adc);
       dataFrame.setSample(it,newSample);
 
       if(debug) edm::LogVerbatim("FTLSimpleElectronicsSimInMIPs") << adc << " (" << chargeColl[it] << "/" << adcLSB_MIP_ << ") ";
     }
 
   if(debug) { 
     std::ostringstream msg;
     dataFrame.print(msg);
     edm::LogVerbatim("FTLSimpleElectronicsSimInMIPs") << msg.str() << std::endl;
   } 
 }

void SimpleElectronicsSimInMIPs::updateOutput	(	FTLDigiCollection &	coll,
		const FTLDataFrame &	rawDataFrame
	)		const

Definition at line 79 of file SimpleElectronicsSimInMIPs.cc.

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

Referenced by getEventSetup(), and run().

                                                                   {
   int itIdx(9);
   if(rawDataFrame.size()<=itIdx+2) return;
   
   FTLDataFrame dataFrame( rawDataFrame.id() );
   dataFrame.resize(5);
   bool putInEvent(false);
   for(int it=0;it<5; ++it) {    
     dataFrame.setSample(it, rawDataFrame[itIdx-2+it]);
     if(it==2) putInEvent = rawDataFrame[itIdx-2+it].threshold(); 
   }
 
   if(putInEvent) {
     coll.push_back(dataFrame);    
   }
 }