DuneDeconvolutionService Class Reference

#include <DuneDeconvolutionService.h>

Inheritance diagram for DuneDeconvolutionService:

Public Member Functions
	DuneDeconvolutionService (fhicl::ParameterSet const &pset, art::ActivityRegistry &)
int	update (detinfo::DetectorClocksData const &clockData, AdcChannelData &data) const
std::ostream &	print (std::ostream &out=std::cout, std::string prefix="") const
Public Member Functions inherited from AdcDeconvolutionService
virtual	~AdcDeconvolutionService ()=default

Private Attributes
int	m_LogLevel

Detailed Description

Definition at line 26 of file DuneDeconvolutionService.h.

Constructor & Destructor Documentation

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

Definition at line 17 of file DuneDeconvolutionService_service.cc.

: m_LogLevel(1) {
  const string myname = "DuneDeconvolutionService::ctor: ";
  pset.get_if_present<int>("LogLevel", m_LogLevel);
  print(cout, myname);
}

Member Function Documentation

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

virtual

Implements AdcDeconvolutionService.

Definition at line 68 of file DuneDeconvolutionService_service.cc.

                                               {
  out << prefix << "DuneDeconvolutionService:"                   << endl;
  out << prefix << "               LogLevel: " << m_LogLevel              << endl;
  return out;
}

int DuneDeconvolutionService::update ( detinfo::DetectorClocksData const &  clockData, AdcChannelData &  data  ) const

virtual

Implements AdcDeconvolutionService.

Definition at line 27 of file DuneDeconvolutionService_service.cc.

                                   {
  const string myname = "DuneDeconvolutionService::update: ";
  AdcChannel chan = data.channel();
  AdcSignalVector& samples = data.samples;
  unsigned int nsam = samples.size();
  if ( m_LogLevel >= 2 ) cout << myname << "Deconvoluting channel " << chan
                             << " with " << nsam << " samples." << endl;
  // Fetch the FFT size.
  art::ServiceHandle<util::LArFFT> hFFT;
  unsigned int fftsize = hFFT->FFTSize();
  if ( nsam > fftsize ) {
    cout << myname << "ERROR: Data has too many ticks for FFT: "
         << nsam << " > " << fftsize << "." << endl;
    return 1;
  }
  // Pad the data to the FFT size.
  bool pad = fftsize > nsam;
  if ( pad ) {
    if ( m_LogLevel >= 3 ) cout << myname << "  Padding sample vector to " << fftsize << endl;
    samples.resize(fftsize);
    for ( unsigned int isam=nsam; isam<fftsize; ++isam ) {
      samples[isam] = samples[isam-nsam];
    }
  }
  // Deconvolute.
  art::ServiceHandle<util::SignalShapingServiceDUNE> hsss;
  if ( m_LogLevel >= 3 ) cout << myname << "  Deconvoluting." << endl;
  hsss->Deconvolute(clockData, chan, samples);
  if ( pad ) samples.resize(nsam);
  if ( m_LogLevel >= 3 ) cout << myname << "  Normalizing." << endl;
  float norm = 1.0/hsss->GetDeconNorm();
  if ( m_LogLevel >= 3 ) cout << myname << "  Scale factor: " << norm << endl;
  for ( float& sam : samples ) sam *= norm;
  // Done.
  return 0;
}

Member Data Documentation

int DuneDeconvolutionService::m_LogLevel

private

Definition at line 40 of file DuneDeconvolutionService.h.

---

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

• dunedataprep/dunedataprep/DataPrep/Service/DuneDeconvolutionService.h
• dunedataprep/dunedataprep/DataPrep/Service/DuneDeconvolutionService_service.cc