#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 16 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 42 of file SimpleElectronicsSimInMIPs.cc.

References ecalMGPA::adc(), adcLSB_MIP_, adcSaturation_MIP_, adcThreshold_MIP_, debug, debug_, createfilelist::int, min(), mps_check::msg, FTLDataFrameT< D, S, DECODE >::print(), FTLSample::set(), FTLDataFrameT< D, S, DECODE >::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 74 of file SimpleElectronicsSimInMIPs.cc.

References FTLDataFrameT< D, S, DECODE >::id(), edm::SortedCollection< T, SORT >::push_back(), FTLDataFrameT< D, S, DECODE >::resize(), and FTLDataFrameT< D, S, DECODE >::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);
   }
 }