WireCell::SigProc::Microboone::CoherentNoiseSub Class Reference

#include <Microboone.h>

Inheritance diagram for WireCell::SigProc::Microboone::CoherentNoiseSub:

Public Member Functions
	CoherentNoiseSub (const std::string &anode="AnodePlane", const std::string &noisedb="OmniChannelNoiseDB")
virtual	~CoherentNoiseSub ()
virtual WireCell::Waveform::ChannelMaskMap	apply (int channel, signal_t &sig) const
virtual WireCell::Waveform::ChannelMaskMap	apply (channel_signals_t &chansig) 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::SigProc::Microboone::ConfigFilterBase
	ConfigFilterBase (const std::string &anode="AnodePlane", const std::string &noisedb="OmniChannelNoiseDB")
virtual	~ConfigFilterBase ()
virtual void	configure (const WireCell::Configuration &config)
	IConfigurable configuration interface. More...
virtual WireCell::Configuration	default_configuration () const
	Optional, override to return a hard-coded default configuration. More...
void	set_channel_noisedb (WireCell::IChannelNoiseDatabase::pointer ndb)
Public Member Functions inherited from WireCell::IConfigurable
virtual	~IConfigurable ()
Public Member Functions inherited from WireCell::IComponent< IConfigurable >
virtual	~IComponent ()

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...
Protected Attributes inherited from WireCell::SigProc::Microboone::ConfigFilterBase
std::string	m_anode_tn
std::string	m_noisedb_tn
IAnodePlane::pointer	m_anode
IChannelNoiseDatabase::pointer	m_noisedb

Detailed Description

Microboone style coherent noise subtraction.

Fixme: in principle, this class could be general purpose for other detectors. However, it uses the functions above which hard code microboone-isms. If those microboone-specific parameters can be pulled out to a higher layer then this class can become generic and move outside of this file.

Definition at line 67 of file Microboone.h.

Constructor & Destructor Documentation

Microboone::CoherentNoiseSub::CoherentNoiseSub	(	const std::string &	anode = "AnodePlane",
		const std::string &	noisedb = "OmniChannelNoiseDB"
	)

Definition at line 871 of file Microboone.cxx.

    : ConfigFilterBase(anode, noisedb) {}

Microboone::CoherentNoiseSub::~CoherentNoiseSub ( )

virtual

Definition at line 873 of file Microboone.cxx.

{}

Member Function Documentation

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

virtual

Filter in place the signal sig from given channel.

Implements WireCell::IChannelFilter.

Definition at line 947 of file Microboone.cxx.

{
    return WireCell::Waveform::ChannelMaskMap();                // not implemented
}

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

virtual

Filter in place a group of signals together.

Implements WireCell::IChannelFilter.

Definition at line 876 of file Microboone.cxx.

{
    //std::cout << "Xin2: " << " " << chansig.size() << std::endl;
    // find the median among all 
    WireCell::Waveform::realseq_t medians = Derivations::CalcMedian(chansig);

    // std::cout << medians.size() << " " << medians.at(100) << " " << medians.at(101) << std::endl;
  

    // For Xin: here is how you can get the response spectrum for this group.
    const int achannel = chansig.begin()->first;

    //std::cerr << "CoherentNoiseSub: ch=" << achannel << " response offset:" << m_noisedb->response_offset(achannel) << std::endl;

    const Waveform::compseq_t& respec = m_noisedb->response(achannel);
    const int res_offset = m_noisedb->response_offset(achannel);
    const int pad_f = m_noisedb->pad_window_front(achannel);
    const int pad_b = m_noisedb->pad_window_back(achannel);

    // need to move these to data base, consult with Brett ...
    // also need to be time dependent ... 
    const float decon_limit = m_noisedb->coherent_nf_decon_limit(achannel);// 0.02;
    const float decon_lf_cutoff = m_noisedb->coherent_nf_decon_lf_cutoff(achannel);
    const float adc_limit = m_noisedb->coherent_nf_adc_limit(achannel);//15;
    const float decon_limit1 = m_noisedb->coherent_nf_decon_limit1(achannel);// 0.08; // loose filter
    const float roi_min_max_ratio = m_noisedb->coherent_nf_roi_min_max_ratio(achannel);// 0.8 default

    const float protection_factor = m_noisedb->coherent_nf_protection_factor(achannel);
    const float min_adc_limit = m_noisedb->coherent_nf_min_adc_limit(achannel);
    
    // std::cout << decon_limit << " " << adc_limit << " " << protection_factor << " " << min_adc_limit << std::endl;
    
    // if (respec.size()) {
    // now, apply the response spectrum to deconvolve the median
    // and apply the special protection or pass respec into
    // SignalProtection().
    //}

    //std::cout << achannel << std::endl;

    // do the signal protection and adaptive baseline
    std::vector< std::vector<int> > rois = Microboone::SignalProtection(medians,respec,res_offset,pad_f,pad_b,decon_limit, decon_lf_cutoff, adc_limit, protection_factor, min_adc_limit);

    // if (achannel == 3840){
    //  std::cout << "Xin1: " << rois.size() << std::endl;
    //  for (size_t i=0;i!=rois.size();i++){
    //      std::cout << "Xin1: " << rois.at(i).front() << " " << rois.at(i).back() << std::endl;
    //  }
    // }
    
    //std::cerr <<"\tSigprotection done: " << chansig.size() << " " << medians.size() << " " << medians.at(100) << " " << medians.at(101) << std::endl;

    // calculate the scaling coefficient and subtract
    Microboone::Subtract_WScaling(chansig, medians, respec, res_offset, rois, decon_limit1, roi_min_max_ratio);

    
    // WireCell::IChannelFilter::signal_t& signal = chansig.begin()->second;
    // for (size_t i=0;i!=signal.size();i++){
    //  signal.at(i) = medians.at(i);
    // }
    
    //std::cerr <<"\tSubtrace_WScaling done" << std::endl;

    // for (auto it: chansig){
    //  std::cout << "Xin3 " << it.first << std::endl;
    //  break;
    // }
  
    return WireCell::Waveform::ChannelMaskMap();                // not implemented
}

---

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