AdcTimingMitigator Class Reference

#include <AdcTimingMitigator.h>

Inheritance diagram for AdcTimingMitigator:

Public Types
using	Index = unsigned int
using	IndexVector = std::vector< Index >
using	IndexSet = std::set< Index >
Public Types inherited from AdcChannelTool
using	Index = unsigned int

Public Member Functions
	AdcTimingMitigator (fhicl::ParameterSet const &ps)
DataMap	update (AdcChannelData &acd) const override
Public Member Functions inherited from TpcDataTool
virtual DataMap	updateTpcData (TpcData &) const
virtual DataMap	viewTpcData (const TpcData &) const
virtual int	forwardTpcData () const
Public Member Functions inherited from AdcChannelTool
virtual	~AdcChannelTool ()=default
virtual DataMap	view (const AdcChannelData &acd) const
virtual DataMap	updateMap (AdcChannelDataMap &acds) const
virtual DataMap	viewMap (const AdcChannelDataMap &acds) const
virtual bool	updateWithView () const
virtual bool	viewWithUpdate () const
virtual DataMap	beginEvent (const DuneEventInfo &) const
virtual DataMap	endEvent (const DuneEventInfo &) const
virtual DataMap	close (const DataMap *dmin=nullptr)

Private Attributes
int	m_LogLevel
float	m_SamplingRatio
float	m_SamplingOffset
IndexVector	m_FEMBs
IndexSet	m_fembSet

Additional Inherited Members
Static Public Member Functions inherited from AdcChannelTool
static int	interfaceNotImplemented ()

Detailed Description

Definition at line 41 of file AdcTimingMitigator.h.

Member Typedef Documentation

using AdcTimingMitigator::Index = unsigned int

Definition at line 46 of file AdcTimingMitigator.h.

using AdcTimingMitigator::IndexSet = std::set<Index>

Definition at line 48 of file AdcTimingMitigator.h.

using AdcTimingMitigator::IndexVector = std::vector<Index>

Definition at line 47 of file AdcTimingMitigator.h.

Constructor & Destructor Documentation

AdcTimingMitigator::AdcTimingMitigator

(

fhicl::ParameterSet const &

ps

)

Definition at line 16 of file AdcTimingMitigator_tool.cc.

: m_LogLevel(ps.get<int>("LogLevel")),
  m_SamplingRatio(ps.get<float>("SamplingRatio")),
  m_SamplingOffset(ps.get<float>("SamplingOffset")),
  m_FEMBs(ps.get<IndexVector>("FEMBs")) {
  const string myname = "AdcTimingMitigator::ctor: ";
  for ( Index ifmb : m_FEMBs ) m_fembSet.insert(ifmb);
  if ( m_LogLevel >= 1 ) {
    cout << myname << "Configuration parameters:" << endl;
    cout << myname << "        LogLevel: " << m_LogLevel << endl;
    cout << myname << "   SamplingRatio: " << m_SamplingRatio << endl;
    cout << myname << "  SamplingOffset: " << m_SamplingOffset << endl;
    cout << myname << "           FEMBs: [";
    bool first = true;
    for ( Index ifmb : m_fembSet ) {
      if ( first ) first = false;
      else cout << ", ";
      cout << ifmb;
    }
    cout << "]" << endl;
  }
}

Member Function Documentation

DataMap AdcTimingMitigator::update ( AdcChannelData &  acd ) const

overridevirtual

Reimplemented from AdcChannelTool.

Definition at line 41 of file AdcTimingMitigator_tool.cc.

                                                            {
  const string myname = "AdcTimingMitigator::update: ";
  DataMap ret;
  if ( m_fembSet.count(acd.fembID()) == 0 ) {
    if ( m_LogLevel >= 2 ) {
      cout << myname << "Skipping channel " << acd.channel() << " in FEMB " << acd.fembID() << endl;
    }
    return ret;
  }
  Index nsam = acd.samples.size();
  if ( nsam == 0 ) {
    if ( m_LogLevel >= 2 ) {
      cout << myname << "Skipping channel " << acd.channel() << " with no samples." << endl;
    }
    return ret;
  }
  if ( m_SamplingRatio < 0.01 ) return ret;
  AdcSignalVector newsams;
  newsams.reserve(nsam/m_SamplingRatio+1);
  // Loop over indices for new sample.
  if ( m_LogLevel >= 3 ) {
    cout << myname << "Updating samples for channel " << acd.channel() << "." << endl;
  }
  for ( Index j=0; ; ++j ) {
    double xj = j*m_SamplingRatio - m_SamplingOffset;
    double yj = acd.samples[0];
    if ( xj > 0.0 ) {
      Index i1 = xj;
      Index i2 = i1 + 1;
      if ( i2 >= acd.samples.size() ) break;
      double y1 = acd.samples[i1];
      double y2 = acd.samples[i2];
      yj = (y2 - y1)*(xj - i1) + y1;
      if ( m_LogLevel >= 3 ) {
        cout << myname << "  " << j << ": " << y1 << " , " << y2 << " ==> " << yj << endl;
      }
    } else {
      if ( m_LogLevel >= 3 ) {
        cout << myname << "  " << j << ": " << " Keeping first sample: " << yj << endl;
      }
    }
    newsams.push_back(yj);
  }
  if ( m_LogLevel >= 2 ) {
    cout << myname << "Updating samples for channel " << acd.channel() << ": size "
         << nsam << " --> " << newsams.size() << endl;
  }
  acd.samples = newsams;
  acd.signal.clear();
  acd.rois.clear();
  return ret;
}

Member Data Documentation

IndexVector AdcTimingMitigator::m_FEMBs

private

Definition at line 60 of file AdcTimingMitigator.h.

IndexSet AdcTimingMitigator::m_fembSet

private

Definition at line 62 of file AdcTimingMitigator.h.

int AdcTimingMitigator::m_LogLevel

private

Definition at line 57 of file AdcTimingMitigator.h.

float AdcTimingMitigator::m_SamplingOffset

private

Definition at line 59 of file AdcTimingMitigator.h.

float AdcTimingMitigator::m_SamplingRatio

private

Definition at line 58 of file AdcTimingMitigator.h.

---

The documentation for this class was generated from the following files:

• dunedataprep/dunedataprep/DataPrep/Tool/AdcTimingMitigator.h
• dunedataprep/dunedataprep/DataPrep/Tool/AdcTimingMitigator_tool.cc