MedianPedestalService Class Reference

#include <MedianPedestalService.h>

Inheritance diagram for MedianPedestalService:

Public Member Functions
	MedianPedestalService (fhicl::ParameterSet const &pset, art::ActivityRegistry &)
int	evaluate (const AdcChannelData &data, AdcSignal *pped=nullptr, AdcSignal *prms=nullptr, AdcSignal *ppederr=nullptr, AdcSignal *prmserr=nullptr) const
std::ostream &	print (std::ostream &out=std::cout, std::string prefix="") const
Public Member Functions inherited from PedestalEvaluationService
virtual	~PedestalEvaluationService ()=default

Private Attributes
int	m_LogLevel
bool	m_UseMean
bool	m_SkipFlaggedSamples
bool	m_SkipSignals
std::vector< std::vector< AdcChannelVector > >	m_GroupChannels

Detailed Description

Definition at line 30 of file MedianPedestalService.h.

Constructor & Destructor Documentation

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

Definition at line 22 of file MedianPedestalService_service.cc.

: m_LogLevel(1),
  m_SkipFlaggedSamples(false),
  m_SkipSignals(false) {
  const string myname = "MedianPedestalService::ctor: ";
  pset.get_if_present<int>("LogLevel", m_LogLevel);
  m_UseMean = pset.get<bool>("UseMean");
  m_SkipFlaggedSamples = pset.get<bool>("SkipFlaggedSamples");
  m_SkipSignals = pset.get<bool>("SkipSignals");
  print(cout, myname);
}

Member Function Documentation

int MedianPedestalService::evaluate ( const AdcChannelData &  data, AdcSignal *  pped = nullptr, AdcSignal *  prms = nullptr, AdcSignal *  ppederr = nullptr, AdcSignal *  prmserr = nullptr  ) const

virtual

Implements PedestalEvaluationService.

Definition at line 37 of file MedianPedestalService_service.cc.

                                                       {
  const string myname = "MedianPedestalService::evaluate: ";
  AdcSignal ped = 0.0;
  AdcSignal rms = 0.0;
  AdcSignal pederr = 0.0;
  AdcSignal rmserr = 0.0;
  const AdcSignalVector sigsin = data.samples;
  const AdcFlagVector flags = data.flags;
  const AdcFilterVector signal = data.signal;
  if ( sigsin.size() == 0 ) return 0;
  AdcSignalVector sigs;
  unsigned int nsig = sigsin.size();
  double sigsum = 0.0;
  unsigned int nsigout = 0;
  for ( unsigned int isig=0; isig<nsig; ++isig ) {
    AdcSignal sig = sigsin[isig];
    if ( m_SkipFlaggedSamples && flags.size() ) {
      AdcFlag flag = flags[isig];
      bool skip = flag != AdcGood;
      if ( skip ) continue;
    }
    if ( m_SkipSignals && signal.size() ) {
      bool isSignal = signal[isig];
      if ( isSignal ) continue;
    }
    if ( m_UseMean ) sigsum += sig;
    else sigs.push_back(sig);
    ++nsigout;
  }
  if ( nsigout ) {
    if ( m_UseMean ) {
      ped = sigsum/nsigout;
    } else {
      sort(sigs.begin(), sigs.end());
      unsigned int isig = sigs.size()/2;
      bool isodd = sigs.size()%2;
      ped = isodd ? sigs[isig-1] : 0.5*(sigs[isig-1] + sigs[isig]);
    }
  }
  if ( pped != nullptr ) *pped = ped;
  if ( prms != nullptr ) *prms = rms;
  if ( ppederr != nullptr ) *ppederr = pederr;
  if ( prmserr != nullptr ) *prmserr = rmserr;
  return 0;
}

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

virtual

Implements PedestalEvaluationService.

Definition at line 87 of file MedianPedestalService_service.cc.

                                         {
  out << prefix << "MedianPedestalService:"                          << endl;
  out << prefix << "             LogLevel: " << m_LogLevel           << endl;
  out << prefix << "   SkipFlaggedSamples: " << m_SkipFlaggedSamples << endl;
  out << prefix << "          SkipSignals: " << m_SkipSignals        << endl;
  return out;
}

Member Data Documentation

std::vector<std::vector<AdcChannelVector> > MedianPedestalService::m_GroupChannels

private

Definition at line 55 of file MedianPedestalService.h.

int MedianPedestalService::m_LogLevel

private

Definition at line 45 of file MedianPedestalService.h.

bool MedianPedestalService::m_SkipFlaggedSamples

private

Definition at line 47 of file MedianPedestalService.h.

bool MedianPedestalService::m_SkipSignals

private

Definition at line 48 of file MedianPedestalService.h.

bool MedianPedestalService::m_UseMean

private

Definition at line 46 of file MedianPedestalService.h.

---

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

• dunedataprep/dunedataprep/DataPrep/Service/MedianPedestalService.h
• dunedataprep/dunedataprep/DataPrep/Service/MedianPedestalService_service.cc