DuneRoiBuildingService Class Reference

#include <DuneRoiBuildingService.h>

Inheritance diagram for DuneRoiBuildingService:

Public Member Functions
	DuneRoiBuildingService (fhicl::ParameterSet const &pset, art::ActivityRegistry &)
int	build (AdcChannelData &data) const
std::ostream &	print (std::ostream &out=std::cout, std::string prefix="") const
Public Member Functions inherited from AdcRoiBuildingService
virtual	~AdcRoiBuildingService ()=default

Private Attributes
int	m_LogLevel
AdcSignal	m_NSigmaStart
AdcSignal	m_NSigmaEnd
AdcIndex	m_PadLow
AdcIndex	m_PadHigh

Detailed Description

Definition at line 39 of file DuneRoiBuildingService.h.

Constructor & Destructor Documentation

DuneRoiBuildingService::DuneRoiBuildingService	(	fhicl::ParameterSet const &	pset,
		art::ActivityRegistry &
	)

Definition at line 22 of file DuneRoiBuildingService_service.cc.

: m_LogLevel(1) {
  const string myname = "DuneRoiBuildingService::ctor: ";
  pset.get_if_present<int>("LogLevel", m_LogLevel);
  m_NSigmaStart = pset.get<AdcSignal>("NSigmaStart");
  m_NSigmaEnd   = pset.get<AdcSignal>("NSigmaEnd");
  m_PadLow      = pset.get<AdcSignal>("PadLow");
  m_PadHigh     = pset.get<AdcSignal>("PadHigh");
  if ( m_LogLevel > 0 ) print(cout, myname);
}

Member Function Documentation

int DuneRoiBuildingService::build ( AdcChannelData &  data ) const

virtual

Implements AdcRoiBuildingService.

Definition at line 36 of file DuneRoiBuildingService_service.cc.

                                                            {
  const string myname = "DuneRoiBuildingService:build: ";
  if ( m_LogLevel >= 2 ) cout << myname << "Building ROIs for channel "
                              << data.channel() << "." << endl;
  data.rois.clear();
  // Get signal shaping service.
  art::ServiceHandle<util::SignalShapingServiceDUNE> hsss;
  AdcSignal sigma = hsss->GetDeconNoise(data.channel());
  if ( m_LogLevel >= 3 ) {
    cout << myname << "  Noise level: " << sigma << " " << data.sampleUnit << endl;
  }
  const AdcSignalVector& sigs = data.samples;
  // Build ROIS before padding and merging.
  AdcFilterVector& signal = data.signal;
  AdcRoiVector& rois = data.rois;
  signal.clear();
  signal.resize(sigs.size(), false);
  bool inroi = false;
  AdcSignal siglow = m_NSigmaEnd*sigma;
  AdcSignal sighigh = m_NSigmaStart*sigma;
  AdcIndex nsig = sigs.size();
  if ( nsig < 1 ) {
    if ( m_LogLevel >= 2 ) cout << myname << "Channel " << data.channel()
                                << " has no samples." << endl;
    return 0;
  }
  for ( AdcIndex isig=0; isig<sigs.size(); ++isig ) {
    AdcSignal sig = sigs[isig];
    if ( inroi ) {
      if ( sig > siglow ) {
        signal[isig] = true;
      } else  {
        inroi = false;
      }
    } else {
      if ( sig > sighigh ) {
        signal[isig] = true;
        inroi = true;
      }
    }
  }
  // Fill the unpadded ROIs.
  data.roisFromSignal();
  // Display ROIs before padding and merging.
  if ( m_LogLevel >= 3 ) {
    cout << myname << "  ROIs before merge (size = " << rois.size() << "):" << endl;
    for ( const AdcRoi& roi : rois ) {
      cout << myname << setw(8) << roi.first << " " << setw(8) << roi.second << endl;
    }
  } else if ( m_LogLevel >= 2 ) {
    cout << myname << "  ROI count before merge: " << data.rois.size() << endl;
  }
  if ( rois.size() == 0 ) return 0;
  // Pad ROIs.
  unsigned int isig1 = 0;
  unsigned int isig2 = 0;
  for ( AdcRoi roi : rois ) {
    isig2 = roi.first;
    isig1 = isig2 > m_PadLow ? isig2 - m_PadLow : 0;
    for ( unsigned int isig=isig1; isig<isig2; ++isig ) signal[isig] = true;
    isig1 = roi.second + 1;
    isig2 = isig1 + m_PadHigh;
    if ( isig2 > nsig ) isig2 = nsig;
    for ( unsigned int isig=isig1; isig<isig2; ++isig ) signal[isig] = true;
  }
  // Fill the final ROIs.
  data.roisFromSignal();
  // Display final ROIs.
  if ( m_LogLevel >= 3 ) {
    cout << myname << "  ROIs after merge (size = " << rois.size() << "):" << endl;
    for ( const AdcRoi& roi : rois ) {
      cout << myname << setw(8) << roi.first << " " << setw(8) << roi.second << endl;
    }
  } else if ( m_LogLevel >= 2 ) {
    cout << myname << "  ROI count after merge: " << data.rois.size() << endl;
  }
  return 0;
}

ostream & DuneRoiBuildingService::print ( std::ostream &  out = std::cout, std::string  prefix = ""  ) const

virtual

Implements AdcRoiBuildingService.

Definition at line 118 of file DuneRoiBuildingService_service.cc.

                                         {
  out << prefix << "DuneRoiBuildingService:" << endl;
  out << prefix << "    LogLevel: " << m_LogLevel << endl;
  out << prefix << " NSigmaStart: " << m_NSigmaStart << endl;
  out << prefix << "   NSigmaEnd: " << m_NSigmaEnd << endl;
  out << prefix << "      PadLow: " << m_PadLow << endl;
  out << prefix << "     PadHigh: " << m_PadHigh << endl;
  return out;
}

Member Data Documentation

int DuneRoiBuildingService::m_LogLevel

private

Definition at line 52 of file DuneRoiBuildingService.h.

AdcSignal DuneRoiBuildingService::m_NSigmaEnd

private

Definition at line 54 of file DuneRoiBuildingService.h.

AdcSignal DuneRoiBuildingService::m_NSigmaStart

private

Definition at line 53 of file DuneRoiBuildingService.h.

AdcIndex DuneRoiBuildingService::m_PadHigh

private

Definition at line 56 of file DuneRoiBuildingService.h.

AdcIndex DuneRoiBuildingService::m_PadLow

private

Definition at line 55 of file DuneRoiBuildingService.h.

---

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

• dunedataprep/dunedataprep/DataPrep/Service/DuneRoiBuildingService.h
• dunedataprep/dunedataprep/DataPrep/Service/DuneRoiBuildingService_service.cc