WireCell::SigProc::Microboone::OneChannelStatus Class Reference

#include <Microboone.h>

Inheritance diagram for WireCell::SigProc::Microboone::OneChannelStatus:

Public Member Functions
	OneChannelStatus (const std::string anode_tn="AnodePlane", double threshold=3.5, int window=5, int nbins=250, double cut=14)
virtual	~OneChannelStatus ()
virtual WireCell::Waveform::ChannelMaskMap	apply (int channel, signal_t &sig) const
virtual WireCell::Waveform::ChannelMaskMap	apply (channel_signals_t &chansig) const
virtual void	configure (const WireCell::Configuration &config)
	Accept a configuration. More...
virtual WireCell::Configuration	default_configuration () const
	Optional, override to return a hard-coded default configuration. More...
bool	ID_lf_noisy (signal_t &sig) const
Public Member Functions inherited from WireCell::IChannelFilter
virtual	~IChannelFilter ()
Public Member Functions inherited from WireCell::IComponent< IChannelFilter >
virtual	~IComponent ()
Public Member Functions inherited from WireCell::Interface
virtual	~Interface ()
Public Member Functions inherited from WireCell::IConfigurable
virtual	~IConfigurable ()
Public Member Functions inherited from WireCell::IComponent< IConfigurable >
virtual	~IComponent ()

Private Attributes
std::string	m_anode_tn
IAnodePlane::pointer	m_anode
double	m_threshold
int	m_window
int	m_nbins
double	m_cut

Additional Inherited Members
Public Types inherited from WireCell::IChannelFilter
typedef Waveform::realseq_t	signal_t
typedef std::map< int, signal_t >	channel_signals_t
Public Types inherited from WireCell::IComponent< IChannelFilter >
typedef std::shared_ptr< IChannelFilter >	pointer
	Access subclass facet by pointer. More...
typedef std::vector< pointer >	vector
	Vector of shared pointers. More...
Public Types inherited from WireCell::Interface
typedef std::shared_ptr< Interface >	pointer
Public Types inherited from WireCell::IComponent< IConfigurable >
typedef std::shared_ptr< IConfigurable >	pointer
	Access subclass facet by pointer. More...
typedef std::vector< pointer >	vector
	Vector of shared pointers. More...

Detailed Description

Definition at line 121 of file Microboone.h.

Constructor & Destructor Documentation

Microboone::OneChannelStatus::OneChannelStatus	(	const std::string	anode_tn = "AnodePlane",
		double	threshold = 3.5,
		int	window = 5,
		int	nbins = 250,
		double	cut = 14
	)

Definition at line 1297 of file Microboone.cxx.

    : m_anode_tn(anode_tn)
    , m_threshold(threshold)
    , m_window(window)
    , m_nbins(nbins)
    , m_cut(cut)
{
}

Microboone::OneChannelStatus::~OneChannelStatus ( )

virtual

Definition at line 1305 of file Microboone.cxx.

{
}

Member Function Documentation

WireCell::Waveform::ChannelMaskMap Microboone::OneChannelStatus::apply ( int  channel, signal_t &  sig  ) const

virtual

Filter in place the signal sig from given channel.

Implements WireCell::IChannelFilter.

Definition at line 1342 of file Microboone.cxx.

{
    WireCell::Waveform::ChannelMaskMap ret;
    auto wpid = m_anode->resolve(ch);
    const int iplane = wpid.index();
    if (iplane!=2){ // not collection
        //std::cout << ch << std::endl;
        if (ID_lf_noisy(signal)){
            WireCell::Waveform::BinRange temp_chirped_bins;
            temp_chirped_bins.first = 0;
            temp_chirped_bins.second = signal.size();
            ret["lf_noisy"][ch].push_back(temp_chirped_bins);
            // std::cout << "lf noisy " << ch << std::endl;
        }
    }
    return ret;
}

WireCell::Waveform::ChannelMaskMap Microboone::OneChannelStatus::apply ( channel_signals_t &  chansig ) const

virtual

Filter in place a group of signals together.

Implements WireCell::IChannelFilter.

Definition at line 1334 of file Microboone.cxx.

{

    
    return WireCell::Waveform::ChannelMaskMap();
}

void Microboone::OneChannelStatus::configure ( const WireCell::Configuration &  config )

virtual

Accept a configuration.

Implements WireCell::IConfigurable.

Definition at line 1309 of file Microboone.cxx.

{
    m_threshold = get<double>(cfg, "Threshold", m_threshold);
    m_window = get<int>(cfg, "Window", m_window);
    m_nbins = get<double>(cfg, "Nbins", m_nbins);
    m_cut = get<double>(cfg, "Cut", m_cut);

    m_anode_tn = get(cfg, "anode", m_anode_tn);
    m_anode = Factory::find_tn<IAnodePlane>(m_anode_tn);
    if (!m_anode) {
        THROW(KeyError() << errmsg{"failed to get IAnodePlane: " + m_anode_tn});
    }
    //std::cerr << "OneChannelStatus: \n" << cfg << "\n";
}

WireCell::Configuration Microboone::OneChannelStatus::default_configuration ( ) const

virtual

Optional, override to return a hard-coded default configuration.

Reimplemented from WireCell::IConfigurable.

Definition at line 1323 of file Microboone.cxx.

{
    Configuration cfg;
    cfg["Threshold"] = m_threshold;
    cfg["Window"] = m_window;
    cfg["Nbins"] = m_nbins;
    cfg["Cut"] = m_cut;
    cfg["anode"] = m_anode_tn; 
    return cfg;
}

bool Microboone::OneChannelStatus::ID_lf_noisy

(

signal_t &

sig

)

const

Definition at line 1361 of file Microboone.cxx.

                                                               {
    // do something ...
    //std::pair<double,double> results = Derivations::CalcRMS(sig);

    //Waveform::mean_rms(sig);
    //double mean = results.first;
    //double rms = results.second;

    double mean = Waveform::percentile_binned(sig,0.5);
    double val1 = Waveform::percentile_binned(sig,0.5-0.34);
    double val2 = Waveform::percentile_binned(sig,0.5+0.34);
    double rms = sqrt((pow(val1-mean,2)+pow(val2-mean,2))/2.);
    
    double valid = 0 ;
    for (int i=0;i<int(sig.size());i++){
        if (sig.at(i)!=0){
            valid ++;
        }
    }
    if (!valid) {
        return false;
    }

    
    
    signal_t temp_sig = sig;
    for (int i=0;i<int(sig.size());i++){
        if (fabs(sig.at(i)-mean) > m_threshold * rms){
            for (int k=-m_window;k!=m_window;k++){
                if (k+i>=0&& k+i < int(sig.size()))
                    temp_sig.at(k+i) = mean;
            }
        }
    }

    double content = 0;
    // for (int i=0;i!=temp_sig.size();i++){
    //     temp_sig.at(i)=i;
    // }
    // do FFT 
    Waveform::compseq_t sig_freq = Waveform::dft(temp_sig);     
    for (int i=0;i!=m_nbins;i++){
        content += abs(sig_freq.at(i+1));
    }

    //std::cout << mean << " " << rms << " " << content << " " << valid << std::endl;
    // std::cout << m_threshold << " " << m_window << " " << m_nbins << " " << m_cut << std::endl;
    
    if (content/valid>m_cut) {
        // std::cerr << "OneChannelStatus::ID_lf_noisy: content=" << content << " valid=" << valid << " m_cut="<<m_cut<< std::endl;
        return true;
    }
    
    return false;
    
    
}

Member Data Documentation

IAnodePlane::pointer WireCell::SigProc::Microboone::OneChannelStatus::m_anode

private

Definition at line 140 of file Microboone.h.

std::string WireCell::SigProc::Microboone::OneChannelStatus::m_anode_tn

private

Definition at line 139 of file Microboone.h.

double WireCell::SigProc::Microboone::OneChannelStatus::m_cut

private

Definition at line 144 of file Microboone.h.

int WireCell::SigProc::Microboone::OneChannelStatus::m_nbins

private

Definition at line 143 of file Microboone.h.

double WireCell::SigProc::Microboone::OneChannelStatus::m_threshold

private

Definition at line 141 of file Microboone.h.

int WireCell::SigProc::Microboone::OneChannelStatus::m_window

private

Definition at line 142 of file Microboone.h.

---

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

• wire-cell-build/sigproc/inc/WireCellSigProc/Microboone.h
• wire-cell-build/sigproc/src/Microboone.cxx