SimpleElectronicsSimInMIPs.cc

Go to the documentation of this file.

#include "SimFastTiming/FastTimingCommon/interface/SimpleElectronicsSimInMIPs.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

using namespace ftl;

SimpleElectronicsSimInMIPs::SimpleElectronicsSimInMIPs(const edm::ParameterSet& pset) :
  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")) {    
}

void SimpleElectronicsSimInMIPs::run(const ftl::FTLSimHitDataAccumulator& input,
                     FTLDigiCollection& output) const {
  
  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 {
    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 {
  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);    
  }
}