BaselineDetrend Class Reference

#include <BaselineDetrend.h>

Inheritance diagram for BaselineDetrend:

Public Member Functions
	BaselineDetrend (fhicl::ParameterSet const &ps)
	~BaselineDetrend () override
DataMap	update (AdcChannelData &acd) const override
Public Member Functions inherited from TpcDataTool
virtual DataMap	updateTpcData (TpcData &) const
virtual DataMap	viewTpcData (const TpcData &) const
virtual int	forwardTpcData () const
Public Member Functions inherited from AdcChannelTool
virtual	~AdcChannelTool ()=default
virtual DataMap	view (const AdcChannelData &acd) const
virtual DataMap	updateMap (AdcChannelDataMap &acds) const
virtual DataMap	viewMap (const AdcChannelDataMap &acds) const
virtual bool	updateWithView () const
virtual bool	viewWithUpdate () const
virtual DataMap	beginEvent (const DuneEventInfo &) const
virtual DataMap	endEvent (const DuneEventInfo &) const
virtual DataMap	close (const DataMap *dmin=nullptr)

Private Member Functions
AdcSignalVector	Smoother (const AdcSignalVector &data, const std::vector< unsigned > &pedidx) const

Private Attributes
int	m_LogLevel
bool	m_UseBasicROI
float	m_Thresh
unsigned	m_Pad
float	m_MinFrac
float	m_Span
TGraphSmooth *	m_GS

Additional Inherited Members
Public Types inherited from AdcChannelTool
using	Index = unsigned int
Static Public Member Functions inherited from AdcChannelTool
static int	interfaceNotImplemented ()

Detailed Description

Definition at line 45 of file BaselineDetrend.h.

Constructor & Destructor Documentation

BaselineDetrend::BaselineDetrend

(

fhicl::ParameterSet const &

ps

)

Definition at line 18 of file BaselineDetrend_tool.cc.

                                                             :
  m_LogLevel( ps.get<int>("LogLevel") ),
  m_UseBasicROI( ps.get<bool>("UseBasicROI") ),  
  m_Thresh( ps.get<float>("Thresh") ), 
  m_Pad( ps.get<unsigned>("Pad") ), 
  m_MinFrac( ps.get<float>("MinFrac") ), 
  m_Span( ps.get<float>("Span") )
{
  const string myname = "BaselineDetrend::ctor: ";
  if( m_LogLevel >= 1 ){
    cout<<myname<<"Use basic ROI "<<m_UseBasicROI<<endl;
    cout<<myname<<"Threshold to drop samples from average "<<m_Thresh<<endl;
    cout<<myname<<"Number of ticks to pad around signals above threshold "<<m_Pad<<endl;
    cout<<myname<<"Min fraction of pedestal samples for detrending "<<m_MinFrac<<endl;
    cout<<myname<<"LOWESS span parameter "<<m_Span<<endl;
  }
  
  m_GS  = new TGraphSmooth( "smoother" );
}

BaselineDetrend::~BaselineDetrend ( )

override

Definition at line 39 of file BaselineDetrend_tool.cc.

                                 {
  m_GS->Delete();
}

Member Function Documentation

AdcSignalVector BaselineDetrend::Smoother ( const AdcSignalVector &  data, const std::vector< unsigned > &  pedidx  ) const

private

Definition at line 74 of file BaselineDetrend_tool.cc.

{
  const string myname = "BaselineDetrend::Smoother: ";
  AdcSignalVector trend( data.size(), 0 );
  
  vector<unsigned> idx;
  if( !pedidx.empty() ){
    idx = pedidx;
  }
  else { // simple ROI threshold finder
                          
    for( unsigned i = 0; i < data.size(); ++i ){
      if( data[i] > m_Thresh ) {
        unsigned count = 1;
        // remove last elements from the vector
        while( !idx.empty() ){
          if( count >= m_Pad ) break;
          //
          if( (int)idx.back() != ((int)i - (int)count) ) break;
          idx.pop_back();
          count++;
        }

        // move forward until we are below threshold + num of samples to skip after
        count = 0;
        for( unsigned j = i+1; j < data.size(); ++j ){
          if( data[j] > m_Thresh ) continue;
          if( count >= m_Pad ) { 
            i = j;
            break;
          }
          count++;
        }

        continue;
      }
   
      //
      idx.push_back( i );
    }
  } // simple roi finder

  // do nothing if we have only a limited number of ped samples
  if( ((float)idx.size())/data.size() < m_MinFrac ){
    if( m_LogLevel >= 2 ){
      cout<<myname<<"number of baseline samples is too low\n";
    }
    return trend;
  }
    
  //
  TGraph gin( idx.size() );
  for( unsigned i=0; i<idx.size(); ++i ){
    gin.SetPoint(i, idx[i], data[idx[i]] );
  }
  
  // do smoothing
  auto gout  = m_GS->SmoothLowess(&gin,"", m_Span);
  
  // our graph data are sorted in ticks
  gout->SetBit(TGraph::kIsSortedX);
  
  // compute detranding function interpolating inside ROIs
  unsigned nidx = idx.size();
  for( unsigned i = 0; i<nidx; ++i ){
    unsigned ii = idx[i];
    double x, y;
    gout->GetPoint(i, x, y);
    trend[ii]  = y;

    //
    if( i == 0 ){
      if( ii != 0 ){
        for(unsigned j=0;j<ii;++j)
          trend[j] = gout->Eval( j );
      }
    }
    else if( i == nidx - 1 ) {
      unsigned iii = data.size() - 1;
      if( ii != iii){
        for(unsigned j=ii;j<=iii;++j)
          trend[j] = gout->Eval( j );
      }
    }
    else {
      unsigned iii = idx[i-1];
      if( (ii - iii ) > 1 ){
        for(unsigned j=iii;j<ii;++j)
          trend[j] = gout->Eval( j );
      }
    }
    
  }
  
  return trend;
}

DataMap BaselineDetrend::update ( AdcChannelData &  acd ) const

overridevirtual

Reimplemented from AdcChannelTool.

Definition at line 44 of file BaselineDetrend_tool.cc.

                                                         {
  DataMap ret;
  AdcSignalVector& data = acd.samples;

  vector<unsigned> pedidx;
  if( !m_UseBasicROI ){
    // use info from ROI search already performed
    auto sig = acd.signal;
    for( unsigned i=0; i<data.size(); ++i ){
      if( sig[i] ) continue;
      pedidx.push_back( i );
    }
  }

  // computed baseline trend
  auto trend = Smoother( data, pedidx );

  // detrending
  for(size_t i=0;i<data.size();++i){
    data[i] -= trend[i];
  }

  // finished
  return ret;
}

Member Data Documentation

TGraphSmooth* BaselineDetrend::m_GS

private

Definition at line 64 of file BaselineDetrend.h.

int BaselineDetrend::m_LogLevel

private

Definition at line 55 of file BaselineDetrend.h.

float BaselineDetrend::m_MinFrac

private

Definition at line 59 of file BaselineDetrend.h.

unsigned BaselineDetrend::m_Pad

private

Definition at line 58 of file BaselineDetrend.h.

float BaselineDetrend::m_Span

private

Definition at line 60 of file BaselineDetrend.h.

float BaselineDetrend::m_Thresh

private

Definition at line 57 of file BaselineDetrend.h.

bool BaselineDetrend::m_UseBasicROI

private

Definition at line 56 of file BaselineDetrend.h.

---

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

• dunedataprep/dunedataprep/DataPrep/Tool/BaselineDetrend.h
• dunedataprep/dunedataprep/DataPrep/Tool/BaselineDetrend_tool.cc