#include <Microboone.h>

Inheritance diagram for WireCell::SigProc::Microboone::ADCBitShift:

Public Member Functions
	ADCBitShift (int nbits=12, int exam_nticks=500, double threshold_sigma=7.5, double threshold_fix=0.8)

virtual	~ADCBitShift ()

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...

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
int	m_nbits

int	m_exam_nticks

double	m_threshold_sigma

double	m_threshold_fix

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 148 of file Microboone.h.

Constructor & Destructor Documentation

Microboone::ADCBitShift::ADCBitShift	(	int	nbits = `12`,
		int	exam_nticks = `500`,
		double	threshold_sigma = `7.5`,
		double	threshold_fix = `0.8`
	)

Definition at line 1144 of file Microboone.cxx.

     : m_nbits(nbits)
     , m_exam_nticks(exam_nticks)
     , m_threshold_sigma(threshold_sigma)
     , m_threshold_fix(threshold_fix)
 {
 }

Microboone::ADCBitShift::~ADCBitShift ( )

virtual

Definition at line 1151 of file Microboone.cxx.

1152 {

1153 }

Member Function Documentation

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

virtual

Filter in place the signal sig from given channel.

Implements WireCell::IChannelFilter.

Definition at line 1183 of file Microboone.cxx.

 {
     std::vector<int> counter(m_nbits,0);
     std::vector<int> s,s1(m_nbits,0);
     int nbin = m_exam_nticks;
 
     for (int i=0;i!=nbin;i++){
         int x = signal.at(i);
         s.clear();
         do
             {
                 s.push_back( (x & 1));
             } while (x >>= 1);
         s.resize(m_nbits);
         
         for (int j=0;j!=m_nbits;j++){
             counter.at(j) += abs(s.at(j) - s1.at(j));
         }
         s1=s;
     }
     
     int nshift = 0;
     for (int i=0;i!=m_nbits;i++){
         if (counter.at(i) < nbin/2. - nbin/2. *sqrt(1./nbin) * m_threshold_sigma){
             nshift ++;
         }else{
             break;
         }
     }
     
     WireCell::Waveform::ChannelMaskMap ret;
     if (nshift!=0 && nshift < 11){
         WireCell::Waveform::BinRange ADC_bit_shifts;
         ADC_bit_shifts.first = nshift;
         ADC_bit_shifts.second = nshift;
 
         ret["ADCBitShift"][ch].push_back(ADC_bit_shifts);
 
         // do the correction ...
         const int nl = signal.size();
         std::vector<int> x(nl,0), x_orig(nl,0);
         for (int i=0;i!=nl;i++){
             x_orig[i] = signal.at(i);
             x[i] = signal.at(i);
         }
         
         std::set<int> fillings;
         int filling;
         int mean = 0;
         for (int i=0;i!=nl;i++){
             filling = WireCell::Bits::lowest_bits(x_orig[i], nshift);
             int y = WireCell::Bits::shift_right(x_orig[i], nshift, filling, 12);
             fillings.insert(filling);
             // cout << y << " ";
             // filling = lowest_bits(x[i], nshift);
             // filling = x[i] & int(pow(2, nshift)-1);
             x[i] = y;
             mean += x[i];
         }
         mean = mean/nl;
 
         int exp_diff = pow(2,m_nbits-nshift)*m_threshold_fix;
 
         // examine the results ... 
         int prev1_bin_content = mean;
         int prev_bin_content = mean;
         int next_bin_content = mean;
         int next1_bin_content = mean;
         int curr_bin_content;
         
         for (int i=0;i<nl;i=i+1){
             curr_bin_content = x[i];
             // when to judge if one is likely bad ... 
             if (abs(curr_bin_content-mean)>exp_diff && 
                 (abs(curr_bin_content - prev_bin_content) > exp_diff ||
                  abs(curr_bin_content - next_bin_content) > exp_diff)
                 ){
                 int exp_value = ( (2*prev_bin_content - prev1_bin_content) +
                                     (prev_bin_content + next_bin_content)/2. + 
                                     (prev_bin_content * 2./3. + next1_bin_content/3.))/3.;
                 for (auto it = fillings.begin(); it!=fillings.end(); it++){
                     int y = WireCell::Bits::shift_right(x_orig[i], nshift, (*it), 12);
                     // when to switch ... 
                     if (fabs(y-exp_value) < fabs(x[i] - exp_value)){
                         x[i] = y;//hs->SetBinContent(i+1,y);
                     }
                 }
             }
             
             prev1_bin_content = prev_bin_content;
             prev_bin_content = x[i];
             if (i+2 < nl){
                 next_bin_content = x[i+2];
             }else{
                 next_bin_content = mean;
             }
             if (i+3 < nl){
                 next1_bin_content = x[i+3];
             }else{
                 next1_bin_content = mean;
             }
         }
         
         // change the histogram ...
         for (int i=0;i!=nl;i++){
             signal.at(i) = x[i];
         }
     }
     return ret;
 }

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

virtual

Filter in place a group of signals together.

Implements WireCell::IChannelFilter.

Definition at line 1176 of file Microboone.cxx.

 {
     // No CMM's are produced.
     return WireCell::Waveform::ChannelMaskMap();
 }

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

virtual

Accept a configuration.

Implements WireCell::IConfigurable.

Definition at line 1155 of file Microboone.cxx.

 {
     m_nbits = get<int>(cfg, "Number_of_ADC_bits", m_nbits);
     m_exam_nticks = get<int>(cfg, "Exam_number_of_ticks_test", m_exam_nticks);
     m_threshold_sigma = get<double>(cfg, "Threshold_sigma_test", m_threshold_sigma);
     m_threshold_fix = get<double>(cfg, "Threshold_fix", m_threshold_fix);
     //std::cerr << "ADCBitShift: \n" << cfg << "\n";
 }

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

virtual

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

Reimplemented from WireCell::IConfigurable.

Definition at line 1163 of file Microboone.cxx.

 {
     Configuration cfg;
     cfg["Number_of_ADC_bits"] = m_nbits;
     cfg["Exam_number_of_ticks_test"] = m_exam_nticks;
     cfg["Threshold_sigma_test"] = m_threshold_sigma;
     cfg["Threshold_fix"] = m_threshold_fix;
         
     return cfg;
 }

Member Data Documentation

int WireCell::SigProc::Microboone::ADCBitShift::m_exam_nticks

private

Definition at line 167 of file Microboone.h.

int WireCell::SigProc::Microboone::ADCBitShift::m_nbits

private

Definition at line 165 of file Microboone.h.

double WireCell::SigProc::Microboone::ADCBitShift::m_threshold_fix

private

Definition at line 171 of file Microboone.h.

double WireCell::SigProc::Microboone::ADCBitShift::m_threshold_sigma

private

Definition at line 169 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

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation