#include <DuneDPhaseRoiBuildingService.h>

Inheritance diagram for DuneDPhaseRoiBuildingService:

Public Member Functions
	DuneDPhaseRoiBuildingService (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 Member Functions
AdcSignalVector	getLowFreqFiltered (const AdcSignalVector &adc) const

Private Attributes
bool	m_UseFilter

std::vector< float >	m_FltCoeffs

int	m_LogLevel

AdcSignal	m_NSigmaStart

AdcSignal	m_NSigmaEnd

AdcIndex	m_PadLow

AdcIndex	m_PadHigh

Detailed Description

Definition at line 47 of file DuneDPhaseRoiBuildingService.h.

Constructor & Destructor Documentation

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

Definition at line 23 of file DuneDPhaseRoiBuildingService_service.cc.

                                                                                   :
     m_UseFilter( pset.get<bool>("UseFilter") ),
     m_FltCoeffs( pset.get<std::vector<float>>("FltCoeffs") ),
     m_LogLevel(1)
 {
   const string myname = "DuneDPhaseRoiBuildingService::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 DuneDPhaseRoiBuildingService::build ( AdcChannelData & data ) const

virtual

Implements AdcRoiBuildingService.

Definition at line 40 of file DuneDPhaseRoiBuildingService_service.cc.

                                                                   {
   const string myname = "DuneDPhaseRoiBuildingService:build: ";
   if ( m_LogLevel >= 2 ) cout << myname << "Building ROIs for channel "
                               << data.channel() << "." << endl;
   data.rois.clear();
   // Get signal shaping service.
   art::ServiceHandle<util::SignalShapingServiceDUNEDPhase> hsss;
   AdcSignal sigma = hsss->GetDeconNoise(data.channel());
   
   // Use filtered or raw ADC.
   AdcSignalVector sigs;
   if (m_UseFilter) { sigs = getLowFreqFiltered(data.samples); }
   else { 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;
 }

AdcSignalVector DuneDPhaseRoiBuildingService::getLowFreqFiltered ( const AdcSignalVector & adc ) const

private

Definition at line 123 of file DuneDPhaseRoiBuildingService_service.cc.

 {
   art::ServiceHandle<util::LArFFT> fft;
   std::vector< TComplex > ch_spectrum(fft->FFTSize() / 2 + 1);
   std::vector< float > ch_waveform(fft->FFTSize(), 0);
 
   size_t n_samples = adc.size();
 
   std::copy(adc.begin(), adc.end(), ch_waveform.begin());
   for (size_t s = n_samples; s < ch_waveform.size(); ++s)
   {
       ch_waveform[s] = ch_waveform[s-1];
   }
   fft->DoFFT(ch_waveform, ch_spectrum);
   for (size_t c = 0; c < m_FltCoeffs.size(); ++c)
   {
       ch_spectrum[c] *= m_FltCoeffs[c];
   }
   fft->DoInvFFT(ch_spectrum, ch_waveform);
 
   AdcSignalVector flt_out(n_samples);
   std::copy(ch_waveform.begin(), ch_waveform.begin()+n_samples, flt_out.begin());
   return flt_out;
 }

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

virtual

Implements AdcRoiBuildingService.

Definition at line 150 of file DuneDPhaseRoiBuildingService_service.cc.

                                          {
   out << prefix << "DuneDPhaseRoiBuildingService:" << 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;
   out << prefix << "    UseFilter: " << m_UseFilter << endl;
   return out;
 }