AdcThresholdSignalFinder Class Reference

#include <AdcThresholdSignalFinder.h>

Inheritance diagram for AdcThresholdSignalFinder:

Public Member Functions
	AdcThresholdSignalFinder (fhicl::ParameterSet const &ps)
DataMap	view (const AdcChannelData &acd) const override
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	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 Types
using	Name = std::string

Private Attributes
int	m_LogLevel
ParFormula *	m_Threshold
unsigned int	m_BinsBefore
unsigned int	m_BinsAfter
bool	m_FlagPositive
bool	m_FlagNegative
RunDataTool *	m_prdtool

Additional Inherited Members
Public Types inherited from AdcChannelTool
using	Index = unsigned int
Static Public Member Functions inherited from AdcChannelTool
static int	interfaceNotImplemented ()

Detailed Description

Definition at line 37 of file AdcThresholdSignalFinder.h.

Member Typedef Documentation

using AdcThresholdSignalFinder::Name = std::string

private

Definition at line 49 of file AdcThresholdSignalFinder.h.

Constructor & Destructor Documentation

AdcThresholdSignalFinder::AdcThresholdSignalFinder

(

fhicl::ParameterSet const &

ps

)

Definition at line 29 of file AdcThresholdSignalFinder.cc.

: m_LogLevel(ps.get<int>("LogLevel")),
  m_Threshold(new RootParFormula("Threshold", ps.get<Name>("Threshold"))),
  m_BinsBefore(ps.get<unsigned int>("BinsBefore")),
  m_BinsAfter(ps.get<unsigned int>("BinsAfter")),
  m_FlagPositive(ps.get<bool>("FlagPositive")),
  m_FlagNegative(ps.get<bool>("FlagNegative")),
  m_prdtool(nullptr) {
  const string myname = "AdcThresholdSignalFinder::ctor: ";
  DuneToolManager* pdtm = DuneToolManager::instance();
  if ( pdtm == nullptr ) {
    cout << myname << "ERROR: Unable to retrieve tool manager." << endl;
  } else {
    string stnam = "runDataTool";
    if ( m_Threshold->npar() ) {
      m_prdtool = pdtm->getShared<RunDataTool>(stnam);
      if ( m_prdtool == nullptr ) {
        cout << myname << "ERROR: RunDataTool " << stnam
             << " not found. Scale factor formula will not be evaluated." << endl;
      } else {
        cout << myname << "RunDataTool retrieved." << endl;
      }
    }
  }
  if ( m_LogLevel >= 1 ) {
    cout << myname << "Configuration parameters:" << endl;
    cout << myname << "      LogLevel: " << m_LogLevel << endl;
    cout << myname << "     Threshold: " << m_Threshold->formulaString() << endl;
    cout << myname << "    BinsBefore: " << m_BinsBefore << endl;
    cout << myname << "     BinsAfter: " << m_BinsAfter << endl;
    cout << myname << "  FlagPositive: " << boolToString(m_FlagPositive) << endl;
    cout << myname << "  FlagNegative: " << boolToString(m_FlagNegative) << endl;
  }
}

Member Function Documentation

DataMap AdcThresholdSignalFinder::update ( AdcChannelData &  acd ) const

overridevirtual

Reimplemented from AdcChannelTool.

Definition at line 66 of file AdcThresholdSignalFinder.cc.

                                                                  {
  const string myname = "AdcThresholdSignalFinder::update: ";
  DataMap ret;
  AdcIndex nsam = acd.samples.size();
  if ( nsam == 0 ) {
    cout << myname << "ERROR: No samples found in channel " << acd.channel() << endl;
    acd.signal.clear();
    acd.rois.clear();
    return ret.setStatus(1);
  }
  if ( m_prdtool != nullptr ) {
    RunData rdat = m_prdtool->runData(acd.run());
    if ( ! rdat.isValid() ) cout << myname << "WARNING: RunData not found." << endl;
    else rdat.setFormulaPars(*m_Threshold);
  }
  if ( ! m_Threshold->ready() ) {
    cout << myname << "WARNING: Using default scale factor " << m_Threshold->defaultEval() << endl;
  } 
  float thr = m_Threshold->eval();
  if ( m_LogLevel >= 2 ) cout << myname << "Finding ROIs for channel " << acd.channel() << endl;
  AdcIndex nsamlo = m_BinsBefore;
  AdcIndex nsamhi = m_BinsAfter;
  acd.signal.clear();
  acd.signal.resize(nsam, false);
  AdcIndex nbinAbove = 0;
  AdcIndex isamUnknown = 0;  // First sample which is not known to be inside or ouside a ROI.
  for ( AdcIndex isam=0; isam<nsam; ++isam ) {
    float val = acd.samples[isam];
    bool keep = ( m_FlagPositive && val >  thr ) ||
                ( m_FlagNegative && val < -thr );
    if ( keep ) {
      ++nbinAbove;
      AdcIndex jsam1 = isam > nsamlo ? isam - nsamlo : 0;
      if ( jsam1 < isamUnknown ) jsam1 = isamUnknown;
      AdcIndex jsam2 = isam + nsamhi + 1;
      if ( jsam2 > nsam ) jsam2 = nsam;
      if ( m_LogLevel >= 4 ) cout << myname << "Trigger: sample[" << isam
                                  << "] = " << val << " ==> range: ["
                                  << jsam1 << ", " << jsam2 << ")" << endl;
      for ( AdcIndex jsam=jsam1; jsam<jsam2; ++jsam ) acd.signal[jsam] = true;
      isamUnknown = jsam2;
    }
  }
  acd.roisFromSignal();
  if ( m_LogLevel >= 3 ) {
    cout << myname << "  # ticks above threshold: " << nbinAbove << endl;
    cout << myname << "  # ROI: " << acd.rois.size() << endl;
  }
  ret.setFloat("threshold", thr);
  ret.setInt("nThresholdBins", nbinAbove);
  ret.setInt("nroi", acd.rois.size());
  return ret;
}

DataMap AdcThresholdSignalFinder::view ( const AdcChannelData &  acd ) const

overridevirtual

Reimplemented from AdcChannelTool.

Definition at line 122 of file AdcThresholdSignalFinder.cc.

                                                                      {
  AdcChannelData acdtmp;
  acdtmp.samples = acd.samples;
  return update(acdtmp);
}

Member Data Documentation

unsigned int AdcThresholdSignalFinder::m_BinsAfter

private

Definition at line 55 of file AdcThresholdSignalFinder.h.

unsigned int AdcThresholdSignalFinder::m_BinsBefore

private

Definition at line 54 of file AdcThresholdSignalFinder.h.

bool AdcThresholdSignalFinder::m_FlagNegative

private

Definition at line 57 of file AdcThresholdSignalFinder.h.

bool AdcThresholdSignalFinder::m_FlagPositive

private

Definition at line 56 of file AdcThresholdSignalFinder.h.

int AdcThresholdSignalFinder::m_LogLevel

private

Definition at line 52 of file AdcThresholdSignalFinder.h.

RunDataTool* AdcThresholdSignalFinder::m_prdtool

private

Definition at line 60 of file AdcThresholdSignalFinder.h.

ParFormula* AdcThresholdSignalFinder::m_Threshold

private

Definition at line 53 of file AdcThresholdSignalFinder.h.

---

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

• dunedataprep/dunedataprep/DataPrep/Tool/AdcThresholdSignalFinder.h
• dunedataprep/dunedataprep/DataPrep/Tool/AdcThresholdSignalFinder.cc