ETLElectronicsSim Class Reference

#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 uint8_t row, const uint8_t column) 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 float	bxTime_
const bool	debug_
const reco::FormulaEvaluator	sigmaEta_
const uint32_t	tdcNbits_
const float	toaLSB_ns_

Detailed Description

Definition at line 21 of file ETLElectronicsSim.h.

Constructor & Destructor Documentation

ETLElectronicsSim::ETLElectronicsSim

(

const edm::ParameterSet &

pset

)

Definition at line 8 of file ETLElectronicsSim.cc.

                                                                :
  debug_( pset.getUntrackedParameter<bool>("debug",false) ),
  bxTime_(pset.getParameter<double>("bxTime") ),
  sigmaEta_( pset.getParameter<std::string>("etaResolution") ),
  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 25 of file ETLElectronicsSim.h.

{ }

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

inline

Definition at line 27 of file ETLElectronicsSim.h.

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

{ }

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

Definition at line 21 of file ETLElectronicsSim.cc.

References adcThreshold_MIP_, bxTime_, PVValHelper::eta, reco::FormulaEvaluator::evaluate(), f, mps_fire::i, runTrivialShaper(), sigmaEta_, and updateOutput().

                                                 {
  
  MTDSimHitData chargeColl, toa;
  

  std::vector<double> emptyV;
  std::vector<double> eta(1);

  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++) {

      if ( (it->second).hit_info[0][i] < adcThreshold_MIP_ ) continue;

      // time of arrival
      float finalToA = (it->second).hit_info[1][i];

      // Gaussian smearing of the time of arrival
      eta[0] = 2.;  // This is just temporary. Once the RECO geometry is 
                    // available, the actual module eta will be used.
      double sigmaToA = sigmaEta_.evaluate(eta, emptyV);

      if ( sigmaToA > 0. )
    finalToA += CLHEP::RandGaussQ::shoot(hre, 0., sigmaToA);

      // fill the time and charge arrays
      const unsigned int ibucket = std::floor( finalToA/bxTime_ );
      if ( (i+ibucket) >= chargeColl.size() ) continue;
      
      chargeColl[i+ibucket] += (it->second).hit_info[0][i];      

      if ( toa[i+ibucket] == 0. || (finalToA-ibucket*bxTime_) < toa[i+ibucket] )
    toa[i+ibucket] = finalToA - ibucket*bxTime_;
      
    }

    // run the shaper to create a new data frame
    ETLDataFrame rawDataFrame( it->first.detid_ );    
    runTrivialShaper(rawDataFrame,chargeColl,toa, it->first.row_, it->first.column_);
    updateOutput(output,rawDataFrame);
    
  }
    
}

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

Definition at line 73 of file ETLElectronicsSim.cc.

References ecalMGPA::adc(), adcLSB_MIP_, adcSaturation_MIP_, adcThreshold_MIP_, debug, debug_, createfilelist::int, min(), mps_check::msg, FTLDataFrameT< D, S, DECODE >::print(), ETLSample::set(), FTLDataFrameT< D, S, DECODE >::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,row,col);
      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 104 of file ETLElectronicsSim.cc.

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

Member Data Documentation

const float ETLElectronicsSim::adcLSB_MIP_

private

Definition at line 55 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

const uint32_t ETLElectronicsSim::adcNbits_

private

Definition at line 51 of file ETLElectronicsSim.h.

const float ETLElectronicsSim::adcSaturation_MIP_

private

Definition at line 54 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

const float ETLElectronicsSim::adcThreshold_MIP_

private

Definition at line 56 of file ETLElectronicsSim.h.

Referenced by run(), and runTrivialShaper().

const float ETLElectronicsSim::bxTime_

private

Definition at line 47 of file ETLElectronicsSim.h.

Referenced by run().

const bool ETLElectronicsSim::debug_

private

Definition at line 46 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().

int ETLElectronicsSim::dfSIZE = 5

static

Definition at line 42 of file ETLElectronicsSim.h.

Referenced by updateOutput().

const reco::FormulaEvaluator ETLElectronicsSim::sigmaEta_

private

Definition at line 48 of file ETLElectronicsSim.h.

Referenced by run().

const uint32_t ETLElectronicsSim::tdcNbits_

private

Definition at line 51 of file ETLElectronicsSim.h.

const float ETLElectronicsSim::toaLSB_ns_

private

Definition at line 57 of file ETLElectronicsSim.h.

Referenced by runTrivialShaper().