#include <SimpleChannelNoiseDB.h>

Inheritance diagram for WireCell::SigProc::SimpleChannelNoiseDB:

Public Types
typedef std::tuple< double, int, int >	mask_t

typedef std::vector< mask_t >	multimask_t

Public Types inherited from WireCell::IChannelNoiseDatabase
typedef WireCell::Waveform::compseq_t	filter_t
	The data type for all frequency-space, multiplicative filters. More...

typedef std::vector< int >	channel_group_t

Public Types inherited from WireCell::IComponent< IChannelNoiseDatabase >
typedef std::shared_ptr< IChannelNoiseDatabase >	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 Member Functions
	SimpleChannelNoiseDB (double tick=0.5 *units::us, int nsamples=9600)

virtual	~SimpleChannelNoiseDB ()

virtual double	sample_time () const
	FIXME: how to handle state changes? More...

virtual double	nominal_baseline (int channel) const
	Return a nominal baseline correction (additive offset) More...

virtual double	gain_correction (int channel) const

virtual double	response_offset (int channel) const
	Return a time offset associated with the response(). More...

virtual int	pad_window_front (int channel) const

virtual int	pad_window_back (int channel) const

virtual float	coherent_nf_decon_limit (int channel) const

virtual float	coherent_nf_decon_lf_cutoff (int channel) const

virtual float	coherent_nf_decon_limit1 (int channel) const

virtual float	coherent_nf_adc_limit (int channel) const

virtual float	coherent_nf_protection_factor (int channel) const

virtual float	coherent_nf_min_adc_limit (int channel) const

virtual float	coherent_nf_roi_min_max_ratio (int channel) const

virtual double	min_rms_cut (int channel) const

virtual double	max_rms_cut (int channel) const

virtual const filter_t &	rcrc (int channel) const
	Return the filter for the RC+RC coupling response function. More...

virtual const filter_t &	config (int channel) const
	Return the filter to correct any wrongly configured channels. More...

virtual const filter_t &	noise (int channel) const
	Return the filter to attenuate noise. More...

virtual const filter_t &	response (int channel) const
	A nominal detector response spectrum for a given channel. More...

virtual std::vector< channel_group_t >	coherent_channels () const
	Return channel grouping for coherent noise subtraction. More...

virtual channel_group_t	bad_channels () const
	Return channels which are considered a'priori "bad". More...

void	set_sampling (double tick=0.5 *units::us, int nsamples=9600)

void	set_nominal_baseline (const std::vector< int > &channels, double baseline)
	Set nominal baseline in units of ADC (eg uB is -2048 for U/V, -400 for W) More...

void	set_gains_shapings (const std::vector< int > &channels, double from_gain_mVfC=7.8, double to_gain_mVfC=14.0, double from_shaping=1.0 units::us, double to_shaping=2.0 units::us)

void	set_response_offset (const std::vector< int > &channels, double offset)
	Set a response offset for the given set of channels. More...

void	set_pad_window_front (const std::vector< int > &channels, int pad_f)

void	set_pad_window_back (const std::vector< int > &channels, int pad_b)

void	set_coherent_nf_decon_limit (const std::vector< int > &channels, float decon_limit)

void	set_coherent_nf_decon_lf_cutoff (const std::vector< int > &channels, float decon_lf_cutoff)

void	set_coherent_nf_decon_limit1 (const std::vector< int > &channels, float decon_limit)

void	set_coherent_nf_adc_limit (const std::vector< int > &channels, float adc_limit)

void	set_coherent_nf_protection_factor (const std::vector< int > &channels, float protection_factor)

void	set_coherent_nf_min_adc_limit (const std::vector< int > &channels, float min_adc_limit)

void	set_coherent_nf_roi_min_max_ratio (const std::vector< int > &channels, float roi_min_max_ratio)

void	set_min_rms_cut (const std::vector< int > &channels, double min_rms)

void	set_min_rms_cut_one (int channel, double min_rms)

void	set_max_rms_cut (const std::vector< int > &channels, double max_rms)

void	set_max_rms_cut_one (int channel, double max_rms)

void	set_rcrc_constant (const std::vector< int > &channels, double rcrc=2000.0)

void	set_response (const std::vector< int > &channels, const filter_t &spectrum)
	Set a detector response spectrum for the set of channels. More...

void	set_filter (const std::vector< int > &channels, const multimask_t &mask)

void	set_channel_groups (const std::vector< channel_group_t > &channel_groups)
	Set the channel groups. More...

void	set_bad_channels (const channel_group_t &bc)
	Set "bad" channels. More...

Public Member Functions inherited from WireCell::IChannelNoiseDatabase
virtual	~IChannelNoiseDatabase ()

virtual channel_group_t	miscfg_channels () const
	Return channels which are considered a'priori "misconfigured". More...

Public Member Functions inherited from WireCell::IComponent< IChannelNoiseDatabase >
virtual	~IComponent ()

Public Member Functions inherited from WireCell::Interface
virtual	~Interface ()

Private Types
typedef std::shared_ptr< filter_t >	shared_filter_t

typedef std::vector< shared_filter_t >	filter_vector_t

Private Member Functions
int	chind (int ch) const

const IChannelNoiseDatabase::filter_t &	get_filter (int channel, const filter_vector_t &fv) const

Private Attributes
double	m_tick

int	m_nsamples

double	m_default_baseline

double	m_default_gain

double	m_default_offset

double	m_default_min_rms

double	m_default_max_rms

int	m_default_pad_f

int	m_default_pad_b

float	m_default_decon_limit

float	m_default_decon_lf_cutoff

float	m_default_adc_limit

float	m_default_decon_limit1

float	m_default_protection_factor

float	m_default_min_adc_limit

float	m_default_roi_min_max_ratio

std::vector< double >	m_baseline

std::vector< double >	m_gain

std::vector< double >	m_offset

std::vector< double >	m_min_rms

std::vector< double >	m_max_rms

std::vector< int >	m_pad_f

std::vector< int >	m_pad_b

std::vector< float >	m_decon_limit

std::vector< float >	m_decon_lf_cutoff

std::vector< float >	m_adc_limit

std::vector< float >	m_decon_limit1

std::vector< float >	m_protection_factor

std::vector< float >	m_min_adc_limit

std::vector< float >	m_roi_min_max_ratio

filter_vector_t	m_rcrc

filter_vector_t	m_config

filter_vector_t	m_masks

filter_vector_t	m_response

shared_filter_t	m_default_filter

shared_filter_t	m_default_response

std::unordered_map< int, int >	m_ch2ind

std::vector< channel_group_t >	m_channel_groups

channel_group_t	m_bad_channels

Detailed Description

Definition at line 17 of file SimpleChannelNoiseDB.h.

Member Typedef Documentation

typedef std::vector<shared_filter_t> WireCell::SigProc::SimpleChannelNoiseDB::filter_vector_t

private

Definition at line 145 of file SimpleChannelNoiseDB.h.

typedef std::tuple<double, int, int> WireCell::SigProc::SimpleChannelNoiseDB::mask_t

Set a constant scaling to a band covering the given frequency bins (inclusively) for the given channels. Frequency bin "i" is from i*f to (i+1)*f where f is 1.0/(nsamples*tick). The largest meaningful frequency bin is nsamples/2. The frequency band is inclusive of both min and max frequency bins. Note, it's up to caller to appropriately segment multiple masks across multiple channels. For any given channel, last call to this method wins.

Definition at line 119 of file SimpleChannelNoiseDB.h.

typedef std::vector<mask_t> WireCell::SigProc::SimpleChannelNoiseDB::multimask_t

Definition at line 120 of file SimpleChannelNoiseDB.h.

typedef std::shared_ptr<filter_t> WireCell::SigProc::SimpleChannelNoiseDB::shared_filter_t

private

Definition at line 144 of file SimpleChannelNoiseDB.h.

Constructor & Destructor Documentation

SimpleChannelNoiseDB::SimpleChannelNoiseDB	(	double	tick = `0.5*units::us`,
		int	nsamples = `9600`
	)

Create a simple channel noise DB for digitized waveforms with the given size and number of samples. Default is for microboone.

Definition at line 15 of file SimpleChannelNoiseDB.cxx.

     : m_tick(-1)
     , m_nsamples(-1)
     , m_default_baseline(0.0)
     , m_default_gain(1.0)
     , m_default_offset(0.0)
     , m_default_min_rms(0.5)
     , m_default_max_rms(10)
     , m_default_pad_f(0)
     , m_default_pad_b(0)
     , m_default_decon_limit(0.02)
     , m_default_decon_lf_cutoff(0.08)
     , m_default_adc_limit(15.0)
     , m_default_decon_limit1(0.08)
     , m_default_protection_factor(5.0)
     , m_default_min_adc_limit(50)
     , m_default_roi_min_max_ratio(0.8)
 {
     set_sampling(tick, nsamples);
 }

SimpleChannelNoiseDB::~SimpleChannelNoiseDB ( )

virtual

Definition at line 35 of file SimpleChannelNoiseDB.cxx.

36 {

37 }

Member Function Documentation

virtual channel_group_t WireCell::SigProc::SimpleChannelNoiseDB::bad_channels ( ) const

inlinevirtual

Return channels which are considered a'priori "bad".

Reimplemented from WireCell::IChannelNoiseDatabase.

Definition at line 55 of file SimpleChannelNoiseDB.h.

                                                          {
                 return m_bad_channels;
             }

int SimpleChannelNoiseDB::chind ( int ch ) const

private

Definition at line 506 of file SimpleChannelNoiseDB.cxx.

 {
     auto it = m_ch2ind.find(ch);
     if (it == m_ch2ind.end()) {
         int ind = m_ch2ind.size();
         m_ch2ind[ch] = ind;
         return ind;
     }
     return it->second;
 }

virtual std::vector<channel_group_t> WireCell::SigProc::SimpleChannelNoiseDB::coherent_channels ( ) const

inlinevirtual

Return channel grouping for coherent noise subtraction.

Implements WireCell::IChannelNoiseDatabase.

Definition at line 52 of file SimpleChannelNoiseDB.h.

                                                                          {
                 return m_channel_groups;
             }

float SimpleChannelNoiseDB::coherent_nf_adc_limit ( int channel ) const

virtual

Implements WireCell::IChannelNoiseDatabase.

Definition at line 136 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_adc_limit.size()) {
         return m_adc_limit[ind];
     }
     return m_default_adc_limit;
 }

float SimpleChannelNoiseDB::coherent_nf_decon_lf_cutoff ( int channel ) const

virtual

Implements WireCell::IChannelNoiseDatabase.

Definition at line 118 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_decon_lf_cutoff.size()) {
         return m_decon_lf_cutoff[ind];
     }
     return m_default_decon_lf_cutoff;
 }

float SimpleChannelNoiseDB::coherent_nf_decon_limit ( int channel ) const

virtual

Implements WireCell::IChannelNoiseDatabase.

Definition at line 109 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_decon_limit.size()) {
         return m_decon_limit[ind];
     }
     return m_default_decon_limit;
 }

float SimpleChannelNoiseDB::coherent_nf_decon_limit1 ( int channel ) const

virtual

Implements WireCell::IChannelNoiseDatabase.

Definition at line 127 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_decon_limit1.size()) {
         return m_decon_limit1[ind];
     }
     return m_default_decon_limit1;
 }

float SimpleChannelNoiseDB::coherent_nf_min_adc_limit ( int channel ) const

virtual

Implements WireCell::IChannelNoiseDatabase.

Definition at line 154 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_min_adc_limit.size()) {
     return m_min_adc_limit[ind];
     }
     return m_default_min_adc_limit;
 }

float SimpleChannelNoiseDB::coherent_nf_protection_factor ( int channel ) const

virtual

Implements WireCell::IChannelNoiseDatabase.

Definition at line 145 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_protection_factor.size()) {
         return m_protection_factor[ind];
     }
     return m_default_protection_factor;
 }

float SimpleChannelNoiseDB::coherent_nf_roi_min_max_ratio ( int channel ) const

virtual

Implements WireCell::IChannelNoiseDatabase.

Definition at line 164 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_roi_min_max_ratio.size()) {
     return m_roi_min_max_ratio[ind];
     }
     return m_default_roi_min_max_ratio;
 }

const IChannelNoiseDatabase::filter_t & SimpleChannelNoiseDB::config ( int channel ) const

virtual

Return the filter to correct any wrongly configured channels.

Implements WireCell::IChannelNoiseDatabase.

Definition at line 197 of file SimpleChannelNoiseDB.cxx.

 {
     return get_filter(channel, m_config);
 }

double SimpleChannelNoiseDB::gain_correction ( int channel ) const

virtual

Return simple gain correction (a multiplicative, unitless scaling) to apply to a given channel.

Implements WireCell::IChannelNoiseDatabase.

Definition at line 50 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_gain.size()) {
         return m_gain[ind];
     }
     return m_default_gain;
 }

const IChannelNoiseDatabase::filter_t & SimpleChannelNoiseDB::get_filter	(	int	channel,
		const filter_vector_t &	fv
	)		const

private

Definition at line 174 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     //std::cerr << "ch=" << channel << " ind=" << ind << " " << fv.size() << std::endl;
     if (0 <= ind && ind < (int)fv.size()) {
         const shared_filter_t sf = fv[ind];
         if (sf == nullptr) {
             return *(m_default_filter.get());
         }
 
         const filter_t* filtp = sf.get();
         return *filtp;
     }
     const filter_t* filtp = m_default_filter.get();
     //std::cerr << "Filter: "<< (void*)filtp << std::endl;
     return *filtp;
 }

double SimpleChannelNoiseDB::max_rms_cut ( int channel ) const

virtual

Implements WireCell::IChannelNoiseDatabase.

Definition at line 78 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_max_rms.size()) {
         return m_max_rms[ind];
     }
     return m_default_max_rms;
 }

double SimpleChannelNoiseDB::min_rms_cut ( int channel ) const

virtual

Implements WireCell::IChannelNoiseDatabase.

Definition at line 69 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_min_rms.size()) {
         return m_min_rms[ind];
     }
     return m_default_min_rms;
 }

const IChannelNoiseDatabase::filter_t & SimpleChannelNoiseDB::noise ( int channel ) const

virtual

Return the filter to attenuate noise.

Implements WireCell::IChannelNoiseDatabase.

Definition at line 202 of file SimpleChannelNoiseDB.cxx.

 {
     return get_filter(channel, m_masks);
 }

double SimpleChannelNoiseDB::nominal_baseline ( int channel ) const

virtual

Return a nominal baseline correction (additive offset)

Implements WireCell::IChannelNoiseDatabase.

Definition at line 41 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_baseline.size()) {
         return m_baseline[ind];
     }
     return m_default_baseline;
 }

int SimpleChannelNoiseDB::pad_window_back ( int channel ) const

virtual

Implements WireCell::IChannelNoiseDatabase.

Definition at line 99 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_pad_b.size()) {
         return m_pad_b[ind];
     }
     return m_default_pad_b;
 
 }

int SimpleChannelNoiseDB::pad_window_front ( int channel ) const

virtual

Implements WireCell::IChannelNoiseDatabase.

Definition at line 89 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_pad_f.size()) {
         return m_pad_f[ind];
     }
     return m_default_pad_f;
 
 }

const IChannelNoiseDatabase::filter_t & SimpleChannelNoiseDB::rcrc ( int channel ) const

virtual

Return the filter for the RC+RC coupling response function.

Implements WireCell::IChannelNoiseDatabase.

Definition at line 192 of file SimpleChannelNoiseDB.cxx.

 {
     return get_filter(channel, m_rcrc);
 }

const IChannelNoiseDatabase::filter_t & SimpleChannelNoiseDB::response ( int channel ) const

virtual

A nominal detector response spectrum for a given channel.

Implements WireCell::IChannelNoiseDatabase.

Definition at line 207 of file SimpleChannelNoiseDB.cxx.

 {
     return get_filter(channel, m_response);
 }

double SimpleChannelNoiseDB::response_offset ( int channel ) const

virtual

Return a time offset associated with the response().

Implements WireCell::IChannelNoiseDatabase.

Definition at line 59 of file SimpleChannelNoiseDB.cxx.

 {
     const int ind = chind(channel);
     if (0 <= ind && ind < (int)m_offset.size()) {
         return m_offset[ind];
     }
     return m_default_offset;
 
 }

virtual double WireCell::SigProc::SimpleChannelNoiseDB::sample_time ( ) const

inlinevirtual

FIXME: how to handle state changes?

Return the time-domain sample period (time in system of units) which is was used in producing the filter response spectral functions (filter_t).

Warning: take care that the number of frequency samples (filter_t::size()) is fixed and may not match the number of ticks in your waveform.

Implements WireCell::IChannelNoiseDatabase.

Definition at line 27 of file SimpleChannelNoiseDB.h.

27 { return m_tick; }

WireCell::SigProc::SimpleChannelNoiseDB::m_tick

double m_tick

Definition: SimpleChannelNoiseDB.h:131

void SimpleChannelNoiseDB::set_bad_channels ( const channel_group_t & bc )

Set "bad" channels.

Definition at line 500 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB:: bad channels: " << bc.size() << std::endl;
     m_bad_channels = bc;
 }

void SimpleChannelNoiseDB::set_channel_groups ( const std::vector< channel_group_t > & channel_groups )

Set the channel groups.

Definition at line 493 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: channel groups: " << channel_groups.size() << " X " << channel_groups[0].size() << std::endl;
     m_channel_groups = channel_groups;
 }

void SimpleChannelNoiseDB::set_coherent_nf_adc_limit	(	const std::vector< int > &	channels,
		float	adc_limit
	)

Definition at line 433 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set pad window back on " << channels.size() << " channels: " << pad_b << std::endl;
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, adc_limit, m_adc_limit, m_default_adc_limit);
     }
 }

void SimpleChannelNoiseDB::set_coherent_nf_decon_lf_cutoff	(	const std::vector< int > &	channels,
		float	decon_lf_cutoff
	)

Definition at line 414 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set pad window back on " << channels.size() << " channels: " << pad_b << std::endl;
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, decon_lf_cutoff, m_decon_lf_cutoff, m_default_decon_lf_cutoff);
     }
 }

void SimpleChannelNoiseDB::set_coherent_nf_decon_limit	(	const std::vector< int > &	channels,
		float	decon_limit
	)

Definition at line 405 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set pad window back on " << channels.size() << " channels: " << pad_b << std::endl;
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, decon_limit, m_decon_limit, m_default_decon_limit);
     }
 }

void SimpleChannelNoiseDB::set_coherent_nf_decon_limit1	(	const std::vector< int > &	channels,
		float	decon_limit
	)

Definition at line 424 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set pad window back on " << channels.size() << " channels: " << pad_b << std::endl;
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, decon_limit1, m_decon_limit1, m_default_decon_limit1);
     }
 }

void SimpleChannelNoiseDB::set_coherent_nf_min_adc_limit	(	const std::vector< int > &	channels,
		float	min_adc_limit
	)

Definition at line 451 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set pad window back on " << channels.size() << " channels: " << pad_b << std::endl;
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, min_adc_limit, m_min_adc_limit, m_default_min_adc_limit);
     }
 }

void SimpleChannelNoiseDB::set_coherent_nf_protection_factor	(	const std::vector< int > &	channels,
		float	protection_factor
	)

Definition at line 442 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set pad window back on " << channels.size() << " channels: " << pad_b << std::endl;
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, protection_factor, m_protection_factor, m_default_protection_factor);
     }
 }

void SimpleChannelNoiseDB::set_coherent_nf_roi_min_max_ratio	(	const std::vector< int > &	channels,
		float	roi_min_max_ratio
	)

Definition at line 460 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set pad window back on " << channels.size() << " channels: " << pad_b << std::endl;
     for (auto ch : channels) {
     int ind = chind(ch);
     set_one(ind, roi_min_max_ratio, m_roi_min_max_ratio, m_default_roi_min_max_ratio);
     }
 }

void SimpleChannelNoiseDB::set_filter	(	const std::vector< int > &	channels,
		const multimask_t &	mask
	)

Definition at line 469 of file SimpleChannelNoiseDB.cxx.

 {
     Waveform::compseq_t spectrum;
     spectrum.resize(m_nsamples,std::complex<float>(1,0));
 
     for (auto m : masks) {
         for (int ind=get<1>(m); ind <= get<2>(m); ++ind) {
             //filt->assign(ind, get<0>(m));
             spectrum.at(ind) = get<0>(m);
         }
         //std::cerr << "set freqmasks: " << get<1>(m) << " " << get<2>(m) << " " << get<0>(m) << " " << spectrum.at(0) << " " << spectrum.at(169) << " " << spectrum.at(170)<< std::endl;
     }
     auto filt = std::make_shared<filter_t>(spectrum);
   
 
     for (auto ch : channels) {
         set_one(chind(ch), filt, m_masks, m_default_filter);
     }
   
 }

void SimpleChannelNoiseDB::set_gains_shapings	(	const std::vector< int > &	channels,
		double	from_gain_mVfC = `7.8`,
		double	to_gain_mVfC = `14.0`,
		double	from_shaping = `1.0*units::us`,
		double	to_shaping = `2.0*units::us`
	)

Set gain/shaping corrections for cnofig_correction. Gains are assumed to be in mV/fC. Shaping times should be in the system of units. Defaults are microboone (but you need to give channels).

Definition at line 295 of file SimpleChannelNoiseDB.cxx.

 {
     const double gain_ratio = to_gain/from_gain;
     // std::cerr << "SimpleChannelNoiseDB: set gain/shaping on " << channels.size() << " chans to: "
     //           << "g=" << from_gain << "->" << to_gain << ", "
     //           << "s=" << from_shaping << "->" << to_shaping << ", "
     //           << "rat=" << gain_ratio
     //           << "m_tick=" << m_tick/units::us << " us"
     //           << "\n";
 
 
     Response::ColdElec from_ce(from_gain, from_shaping);
     Response::ColdElec to_ce(to_gain, to_shaping);
     // auto to_sig   =   to_ce.generate(WireCell::Binning(m_nsamples, 0, m_nsamples*m_tick));
     // auto from_sig = from_ce.generate(WireCell::Binning(m_nsamples, 0, m_nsamples*m_tick));
     auto to_sig   =   to_ce.generate(WireCell::Waveform::Domain(0, m_nsamples*m_tick), m_nsamples);
     auto from_sig = from_ce.generate(WireCell::Waveform::Domain(0, m_nsamples*m_tick), m_nsamples);
     
     auto to_filt   = Waveform::dft(to_sig);
     auto from_filt = Waveform::dft(from_sig);
 
     // auto from_filt_sum = Waveform::sum(from_filt);
     // auto to_filt_sum   = Waveform::sum(to_filt);
     
     Waveform::shrink(to_filt, from_filt); // divide
     auto filt = std::make_shared<filter_t>(to_filt);
 
     //std::cout << filt->at(0) << " " << filt->at(1) << std::endl;
     // std::cerr << "SimpleChannelNoiseDB: "
     //           << " from_sig sum=" << Waveform::sum(from_sig)
     //           << " to_sig sum=" << Waveform::sum(to_sig)
     //           << " from_filt sum=" << from_filt_sum
     //           << " to_filt sum=" << to_filt_sum
     //           << " rat_filt sum=" << Waveform::sum(to_filt)
     //           << std::endl;
 
 
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, filt, m_config, m_default_filter);
         set_one(ind, gain_ratio, m_gain, m_default_gain);
     }
 }

void SimpleChannelNoiseDB::set_max_rms_cut	(	const std::vector< int > &	channels,
		double	max_rms
	)

Definition at line 370 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set max rms cut on "<<channels.size()<<":[" << channels.front() << "," << channels.back() << "] to: " << max_rms << std::endl;
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, max_rms, m_max_rms, m_default_max_rms);
     }
 }

void SimpleChannelNoiseDB::set_max_rms_cut_one	(	int	channel,
		double	max_rms
	)

Definition at line 379 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set max rms cut on " << ch << " to: " << max_rms << std::endl;
     int ind = chind(ch);
     set_one(ind, max_rms, m_max_rms, m_default_max_rms);
 }

void SimpleChannelNoiseDB::set_min_rms_cut	(	const std::vector< int > &	channels,
		double	min_rms
	)

Definition at line 352 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set min rms cut on "<<channels.size()<<":[" << channels.front() << "," << channels.back() << "] to: " << min_rms << std::endl;
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, min_rms, m_min_rms, m_default_min_rms);
     }
 }

void SimpleChannelNoiseDB::set_min_rms_cut_one	(	int	channel,
		double	min_rms
	)

Definition at line 361 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set min rms cut on " << ch << " to: " << min_rms << std::endl;
     int ind = chind(ch);
     set_one(ind, min_rms, m_min_rms, m_default_min_rms);
 }

void SimpleChannelNoiseDB::set_nominal_baseline	(	const std::vector< int > &	channels,
		double	baseline
	)

Set nominal baseline in units of ADC (eg uB is -2048 for U/V, -400 for W)

Definition at line 245 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set baseline to " << channels.size() << " chans: " << baseline << std::endl;
     for (auto ch : channels) {
         set_one(chind(ch), baseline, m_baseline, m_default_baseline);
     }
 }

void SimpleChannelNoiseDB::set_pad_window_back	(	const std::vector< int > &	channels,
		int	pad_b
	)

Definition at line 396 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set pad window back on " << channels.size() << " channels: " << pad_b << std::endl;
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, pad_b, m_pad_b, m_default_pad_b);
     }
 }

void SimpleChannelNoiseDB::set_pad_window_front	(	const std::vector< int > &	channels,
		int	pad_f
	)

Definition at line 387 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set pad window front on " << channels.size() << " channels: " << pad_f << std::endl;
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, pad_f, m_pad_f, m_default_pad_f);
     }
 }

void SimpleChannelNoiseDB::set_rcrc_constant	(	const std::vector< int > &	channels,
		double	rcrc = `2000.0`
	)

Set the RC+RC time constant in the system of units for the digitization sample time ("tick"). Default is for microboone.

Definition at line 252 of file SimpleChannelNoiseDB.cxx.

 {
     Response::SimpleRC rcres(rcrc,m_tick);
     //auto signal = rcres.generate(WireCell::Waveform::Domain(0, m_nsamples*m_tick), m_nsamples);
     //    auto signal = rcres.generate(WireCell::Binning(m_nsamples, 0, m_nsamples*m_tick));
     auto signal = rcres.generate(WireCell::Waveform::Domain(0, m_nsamples*m_tick), m_nsamples);
     
     Waveform::compseq_t spectrum = Waveform::dft(signal);
     
     //std::cout << rcrc << " " << m_tick << " " << m_nsamples << " " << signal.front() << " " << signal.at(1) << " " << signal.at(2) << std::endl;
 
     // get the square of it because there are two RC filters
     Waveform::compseq_t spectrum2 = spectrum;
     Waveform::scale(spectrum2,spectrum);
     // for (auto it : spectrum){
     //   float real_part = it.real();
     //   float imag_part = it.imag();
     //   std::complex<float> temp(real_part*real_part-imag_part*imag_part,2*real_part*imag_part);
     //   spectrum2.push_back(temp);
     // }
 
     // std::cerr << "SimpleChannelNoiseDB:: get rcrc as: " << rcrc 
     //           << " sum=" << Waveform::sum(spectrum2)
     //           << std::endl;
 
     auto filt = std::make_shared<filter_t>(spectrum2);
     
     for (auto ch : channels) {
         set_one(chind(ch), filt, m_rcrc, m_default_filter);
     }
 }

void SimpleChannelNoiseDB::set_response	(	const std::vector< int > &	channels,
		const filter_t &	spectrum
	)

Set a detector response spectrum for the set of channels.

Definition at line 284 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set respnose on " << channels.size() << " chans with: " << spectrum.size() << " samples\n";
 
     auto filt = std::make_shared<filter_t>(spectrum);
     for (auto ch : channels) {
         set_one(chind(ch), filt, m_response, m_default_response);
     }
 
 }

void SimpleChannelNoiseDB::set_response_offset	(	const std::vector< int > &	channels,
		double	offset
	)

Set a response offset for the given set of channels.

Definition at line 341 of file SimpleChannelNoiseDB.cxx.

 {
     //std::cerr << "SimpleChannelNoiseDB: set response offset " << channels.size() << " chans to: " << offset << std::endl;
     for (auto ch : channels) {
         int ind = chind(ch);
         set_one(ind, offset, m_offset, m_default_offset);
     }
 }

void SimpleChannelNoiseDB::set_sampling	(	double	tick = `0.5*units::us`,
		int	nsamples = `9600`
	)

Set the size and number of samples of a channel's waveform, default is for microboone.

Warning: calling this will reset any settings for gains+shaping and rcrc as they depend on knowing the sampling.

Definition at line 213 of file SimpleChannelNoiseDB.cxx.

 {
 
     if (m_nsamples == nsamples && tick == m_tick) {
         //std::cerr << "Sampling unchanged: " << nsamples << " @ " << m_tick/units::ms << " ms" << std::endl;
         return;
     }
     m_nsamples = nsamples;
     m_tick = tick;
 
     m_rcrc.clear();
     m_config.clear();
     m_response.clear();
 
     Waveform::compseq_t spectrum;
     spectrum.resize(nsamples,std::complex<float>(1,0));
     m_default_filter = std::make_shared<filter_t>(spectrum);
     Waveform::compseq_t empty;
     m_default_response = std::make_shared<filter_t>(empty);
 }

Member Data Documentation

std::vector<float> WireCell::SigProc::SimpleChannelNoiseDB::m_adc_limit

private

Definition at line 141 of file SimpleChannelNoiseDB.h.

channel_group_t WireCell::SigProc::SimpleChannelNoiseDB::m_bad_channels

private

Definition at line 156 of file SimpleChannelNoiseDB.h.

std::vector<double> WireCell::SigProc::SimpleChannelNoiseDB::m_baseline

private

Definition at line 139 of file SimpleChannelNoiseDB.h.

std::unordered_map<int,int> WireCell::SigProc::SimpleChannelNoiseDB::m_ch2ind

mutableprivate

Definition at line 150 of file SimpleChannelNoiseDB.h.

std::vector< channel_group_t > WireCell::SigProc::SimpleChannelNoiseDB::m_channel_groups

private

Definition at line 155 of file SimpleChannelNoiseDB.h.

filter_vector_t WireCell::SigProc::SimpleChannelNoiseDB::m_config

private

Definition at line 146 of file SimpleChannelNoiseDB.h.

std::vector<float> WireCell::SigProc::SimpleChannelNoiseDB::m_decon_lf_cutoff

private

Definition at line 141 of file SimpleChannelNoiseDB.h.

std::vector<float> WireCell::SigProc::SimpleChannelNoiseDB::m_decon_limit

private

Definition at line 141 of file SimpleChannelNoiseDB.h.

std::vector<float> WireCell::SigProc::SimpleChannelNoiseDB::m_decon_limit1

private

Definition at line 141 of file SimpleChannelNoiseDB.h.

float WireCell::SigProc::SimpleChannelNoiseDB::m_default_adc_limit

private

Definition at line 137 of file SimpleChannelNoiseDB.h.

double WireCell::SigProc::SimpleChannelNoiseDB::m_default_baseline

private

Definition at line 134 of file SimpleChannelNoiseDB.h.

float WireCell::SigProc::SimpleChannelNoiseDB::m_default_decon_lf_cutoff

private

Definition at line 137 of file SimpleChannelNoiseDB.h.

float WireCell::SigProc::SimpleChannelNoiseDB::m_default_decon_limit

private

Definition at line 137 of file SimpleChannelNoiseDB.h.

float WireCell::SigProc::SimpleChannelNoiseDB::m_default_decon_limit1

private

Definition at line 137 of file SimpleChannelNoiseDB.h.

shared_filter_t WireCell::SigProc::SimpleChannelNoiseDB::m_default_filter

private

Definition at line 147 of file SimpleChannelNoiseDB.h.

double WireCell::SigProc::SimpleChannelNoiseDB::m_default_gain

private

Definition at line 134 of file SimpleChannelNoiseDB.h.

double WireCell::SigProc::SimpleChannelNoiseDB::m_default_max_rms

private

Definition at line 135 of file SimpleChannelNoiseDB.h.

float WireCell::SigProc::SimpleChannelNoiseDB::m_default_min_adc_limit

private

Definition at line 137 of file SimpleChannelNoiseDB.h.

double WireCell::SigProc::SimpleChannelNoiseDB::m_default_min_rms

private

Definition at line 135 of file SimpleChannelNoiseDB.h.

double WireCell::SigProc::SimpleChannelNoiseDB::m_default_offset

private

Definition at line 134 of file SimpleChannelNoiseDB.h.

int WireCell::SigProc::SimpleChannelNoiseDB::m_default_pad_b

private

Definition at line 136 of file SimpleChannelNoiseDB.h.

int WireCell::SigProc::SimpleChannelNoiseDB::m_default_pad_f

private

Definition at line 136 of file SimpleChannelNoiseDB.h.

float WireCell::SigProc::SimpleChannelNoiseDB::m_default_protection_factor

private

Definition at line 137 of file SimpleChannelNoiseDB.h.

shared_filter_t WireCell::SigProc::SimpleChannelNoiseDB::m_default_response

private

Definition at line 148 of file SimpleChannelNoiseDB.h.

float WireCell::SigProc::SimpleChannelNoiseDB::m_default_roi_min_max_ratio

private

Definition at line 137 of file SimpleChannelNoiseDB.h.

std::vector<double> WireCell::SigProc::SimpleChannelNoiseDB::m_gain

private

Definition at line 139 of file SimpleChannelNoiseDB.h.

filter_vector_t WireCell::SigProc::SimpleChannelNoiseDB::m_masks

private

Definition at line 146 of file SimpleChannelNoiseDB.h.

std::vector<double> WireCell::SigProc::SimpleChannelNoiseDB::m_max_rms

private

Definition at line 139 of file SimpleChannelNoiseDB.h.

std::vector<float> WireCell::SigProc::SimpleChannelNoiseDB::m_min_adc_limit

private

Definition at line 141 of file SimpleChannelNoiseDB.h.

std::vector<double> WireCell::SigProc::SimpleChannelNoiseDB::m_min_rms

private

Definition at line 139 of file SimpleChannelNoiseDB.h.

int WireCell::SigProc::SimpleChannelNoiseDB::m_nsamples

private

Definition at line 132 of file SimpleChannelNoiseDB.h.

std::vector<double> WireCell::SigProc::SimpleChannelNoiseDB::m_offset

private

Definition at line 139 of file SimpleChannelNoiseDB.h.

std::vector<int> WireCell::SigProc::SimpleChannelNoiseDB::m_pad_b

private

Definition at line 140 of file SimpleChannelNoiseDB.h.

std::vector<int> WireCell::SigProc::SimpleChannelNoiseDB::m_pad_f

private

Definition at line 140 of file SimpleChannelNoiseDB.h.

std::vector<float> WireCell::SigProc::SimpleChannelNoiseDB::m_protection_factor

private

Definition at line 141 of file SimpleChannelNoiseDB.h.

filter_vector_t WireCell::SigProc::SimpleChannelNoiseDB::m_rcrc

private

Definition at line 146 of file SimpleChannelNoiseDB.h.

filter_vector_t WireCell::SigProc::SimpleChannelNoiseDB::m_response

private

Definition at line 146 of file SimpleChannelNoiseDB.h.

std::vector<float> WireCell::SigProc::SimpleChannelNoiseDB::m_roi_min_max_ratio

private

Definition at line 141 of file SimpleChannelNoiseDB.h.

double WireCell::SigProc::SimpleChannelNoiseDB::m_tick

private

Definition at line 131 of file SimpleChannelNoiseDB.h.

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

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

Public Types

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation