AcdWireReader_tool.cc

Go to the documentation of this file.

// AcdWireReader_tool.cc

#include "AcdWireReader.h"
#include <iostream>
#include "lardataobj/RecoBase/Wire.h"

using std::string;
using std::cout;
using std::endl;
using raw::RawDigit;
using recob::Wire;

//**********************************************************************

AcdWireReader::AcdWireReader(fhicl::ParameterSet const& ps)
: m_LogLevel(ps.get<int>("LogLevel")) {
  const string myname = "AcdWireReader::ctor: ";
  if ( m_LogLevel > 0 ) {
    cout << myname << "    LogLevel: " << m_LogLevel << endl;
  }
}

//**********************************************************************

DataMap AcdWireReader::update(AdcChannelData& acd) const {
  const string myname = "AcdWireReader::update: ";
  // Take the wire from the channel data.
  const Wire* pwir = acd.wire;
  if ( pwir == nullptr ) {
    cout << myname << "ERROR: Wire is null." << endl;
    return DataMap(1);
  }
  const Wire& wir = *pwir;
  // Check the input data is empty.
  if ( acd.samples.size() ) {
    cout << myname << "ERROR: ADC channel has prepared data." << endl;
    return DataMap(2);
  }
  // Check the signal flags record is empty.
  if ( acd.signal.size() ) {
    cout << myname << "ERROR: ADC channel has signal flags." << endl;
    return DataMap(3);
  }
  // Check the ROI data is empty.
  if ( acd.rois.size() ) {
    cout << myname << "ERROR: ADC channel has ROIs." << endl;
    return DataMap(4);
  }
  // Set or check the wire index.
  if ( acd.wireIndex != AdcChannelData::badIndex() ) {
    cout << myname << "ERROR: ADC channel has a wire index." << endl;
    return DataMap(5);
  }
  // Set or check the channel.
  if ( acd.channel() == AdcChannelData::badChannel() ) {
    acd.setChannelInfo(wir.Channel());
  } else {
    if ( acd.channel() != wir.Channel() ) {
      cout << myname << "ERROR: Wire has inconsistent channel number." << endl;
      return DataMap(6);
    }
  }
  // Resize the signal array.
  unsigned int nsig = pwir->NSignal();
  // Loop over ROIs and fetch the samples and signals.
  acd.samples.resize(nsig, 0.0);
  acd.signal.resize(nsig, false);
  const lar::sparse_vector<float>& inRois = pwir->SignalROI();
  for ( const lar::sparse_vector<double>::datarange_t& range : inRois.get_ranges()) {
    unsigned int isig1 = range.begin_index();
    unsigned int isig2 = isig1 + range.size();
    for ( unsigned int isig=isig1; isig<isig2; ++isig ) {
      acd.samples[isig] = inRois[isig];
      acd.signal[isig] = true;
    }
  }
  acd.roisFromSignal();   // Fill in the channel data ROIs

  return DataMap(0);
}

//**********************************************************************

DEFINE_ART_CLASS_TOOL(AcdWireReader)