Dune35tZeroSuppressService Class Reference

#include <Dune35tZeroSuppressService.h>

Inheritance diagram for Dune35tZeroSuppressService:

Public Types
typedef unsigned int	Index
Public Types inherited from AdcSuppressService
typedef unsigned int	Channel

Public Member Functions
	Dune35tZeroSuppressService (AdcCount ts, AdcCount tl, AdcCount td, Index ns, Index nl, Index nd, Index nt)
	Dune35tZeroSuppressService (fhicl::ParameterSet const &pset, art::ActivityRegistry &)
int	filter (const AdcCountVector &sigs, Channel chan, AdcPedestal ped, AdcFilterVector &keep) const
std::ostream &	print (std::ostream &out=std::cout, std::string prefix="") const
void	setDebug (int dbg)
Public Member Functions inherited from AdcSuppressService
virtual	~AdcSuppressService ()=default
int	filter (AdcCountSelection &acs)

Private Attributes
AdcCount	m_ts
AdcCount	m_tl
AdcCount	m_td
Index	m_ns
Index	m_nl
Index	m_nd
Index	m_nt
int	m_LogLevel

Detailed Description

Definition at line 30 of file Dune35tZeroSuppressService.h.

Member Typedef Documentation

typedef unsigned int Dune35tZeroSuppressService::Index

Definition at line 34 of file Dune35tZeroSuppressService.h.

Constructor & Destructor Documentation

Dune35tZeroSuppressService::Dune35tZeroSuppressService	(	AdcCount	ts,
		AdcCount	tl,
		AdcCount	td,
		Index	ns,
		Index	nl,
		Index	nd,
		Index	nt
	)

Definition at line 77 of file Dune35tZeroSuppressService_service.cc.

: m_ts(ts), m_tl(tl), m_td(td), m_ns(ns), m_nl(nl), m_nd(nd), m_nt(nt),
  m_LogLevel(1) { }

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

Member Function Documentation

int Dune35tZeroSuppressService::filter ( const AdcCountVector &  sigs, Channel  chan, AdcPedestal  ped, AdcFilterVector &  keep  ) const

virtual

Implements AdcSuppressService.

Definition at line 102 of file Dune35tZeroSuppressService_service.cc.

                                                                                               {
  const string myname = "ZeroSuppress35t::filter: ";
  if ( m_LogLevel >= 2 ) cout << myname << "Filtering channel " << chan << " with tick count " << sigs.size() << endl;
  bool m_skipStuck = false;
  AdcCount m_ts = 0; 
  unsigned int nsig = sigs.size();
  keep.clear();
  keep.resize(nsig, false);
  // Loop over signals.
  SigState state = OUT;
  unsigned int nlow = 0;
  for ( unsigned int isig=0; isig<nsig; ++isig ) {
    AdcCount sig = sigs[isig];
    // Evaluate a running signal sum of the preceding m_nl signals.
    RunningSum rs(sigs, ped, isig, m_ns, m_ts, m_skipStuck);
    AdcCount asigsum = std::abs(rs.sigsum);
    if ( m_LogLevel >= 3 ) cout << myname << setw(6) << isig << setw(6) << sig << setw(5) << sstate(state) << endl;
    // Last tick is outside a signal.
    if ( state == OUT || state == END ) {
      // If this tick is above TL, we are in the live region of a signal.
      // Keep the NL preceding signals.
      AdcCount sumthresh = m_tl*rs.count;
      bool keepit = asigsum > sumthresh;
      if ( m_LogLevel >= 3 ) cout << myname << " RS sum/thresh=" << setw(3) << rs.sigsum << "/" << setw(3) << sumthresh << endl;
      if ( keepit ) {
        if ( m_LogLevel == 2 ) cout << myname << setw(6) << isig << ": RS sum/thresh=" << setw(3) << rs.sigsum
                                    << "/" << setw(3) << sumthresh << endl;
        state = LIVE;
        unsigned int jsig1 = 0;
        if ( isig > m_nl ) jsig1 = isig - m_nl;
        unsigned int jsig2 = isig; 
        for ( unsigned int jsig=jsig1; jsig<jsig2; ++jsig) {
          keep[jsig] = true;
        }
      } else {
        state = OUT;
      }
    } else {
      // Last tick is is in the live region of a signal.
      AdcCount sumthresh = m_td*rs.count;
      if ( m_LogLevel >= 3 ) cout << myname << " RS sum/thresh=" << setw(3) << rs.sigsum << "/" << setw(3) << sumthresh << endl;
        if ( state == LIVE ) {
        // If this tick is below TD, we are in the dead region of a signal.
        if ( asigsum <= sumthresh ) {
          state = DEAD; 
          nlow = 1;
        }
      // Last tick is is in the dead region of a signal.
      } else if ( state == DEAD ) {
        // If signal is above TD, we are back in the live region.
        if ( asigsum > sumthresh ) {
          state = LIVE;
          nlow = 0;
        // If this is the ND'th consecutive signal in the dead region, we
        // have reached the end of the signal.
        // Keep this signal and a tail.
        } else if ( ++nlow >= m_nd ) {
          state = END;
          nlow = 0;
          // Protect the tail.
          unsigned int jsig1 = isig + 1;
          unsigned int jsig2 = jsig1 + m_nt;
          if ( jsig2 > nsig ) jsig2 = nsig;
          for ( unsigned int jsig=jsig1; jsig<jsig2; ++jsig) {
            keep[jsig] = true;
          }
        }
      } else {
        assert(false);
      }
    }
    if ( state != OUT ) keep[isig] = true;
  }
  return 0;
}

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

virtual

Implements AdcSuppressService.

Definition at line 180 of file Dune35tZeroSuppressService_service.cc.

                                                                            {
  out << prefix << "    TS = " << m_ts << endl;
  out << prefix << "    TL = " << m_tl << endl;
  out << prefix << "    TD = " << m_td << endl;
  out << prefix << "    NS = " << m_ns << endl;
  out << prefix << "    NL = " << m_nl << endl;
  out << prefix << "    ND = " << m_nd << endl;
  out << prefix << "    NT = " << m_nt << endl;
  out << prefix << "LogLevel: " << m_LogLevel << endl;
  return out;
}

void Dune35tZeroSuppressService::setDebug

(

int

dbg

)

Definition at line 194 of file Dune35tZeroSuppressService_service.cc.

                                                 {
  m_LogLevel = dbg;
}

Member Data Documentation

int Dune35tZeroSuppressService::m_LogLevel

private

Definition at line 61 of file Dune35tZeroSuppressService.h.

Index Dune35tZeroSuppressService::m_nd

private

Definition at line 59 of file Dune35tZeroSuppressService.h.

Index Dune35tZeroSuppressService::m_nl

private

Definition at line 58 of file Dune35tZeroSuppressService.h.

Index Dune35tZeroSuppressService::m_ns

private

Definition at line 57 of file Dune35tZeroSuppressService.h.

Index Dune35tZeroSuppressService::m_nt

private

Definition at line 60 of file Dune35tZeroSuppressService.h.

AdcCount Dune35tZeroSuppressService::m_td

private

Definition at line 56 of file Dune35tZeroSuppressService.h.

AdcCount Dune35tZeroSuppressService::m_tl

private

Definition at line 55 of file Dune35tZeroSuppressService.h.

AdcCount Dune35tZeroSuppressService::m_ts

private

Definition at line 54 of file Dune35tZeroSuppressService.h.

---

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

• dunesim/dunesim/DetSim/Service/Dune35tZeroSuppressService.h
• dunesim/dunesim/DetSim/Service/Dune35tZeroSuppressService_service.cc