PedAlgoRmsSlider.h

Go to the documentation of this file.

#ifndef larana_OPTICALDETECTOR_PEDALGORMSSLIDER_H
#define larana_OPTICALDETECTOR_PEDALGORMSSLIDER_H

#include "larana/OpticalDetector/OpHitFinder/OpticalRecoTypes.h"

#include "PMTPedestalBase.h"
namespace fhicl {
  class ParameterSet;
}

#include <fstream>
#include <string>
#include <vector>

namespace pmtana {

  class PedAlgoRmsSlider : public PMTPedestalBase {

  public:
    PedAlgoRmsSlider(const std::string name = "PedRmsSlider");

    PedAlgoRmsSlider(const fhicl::ParameterSet& pset, const std::string name = "PedRmsSlider");

    void PrintInfo();

  protected:
    bool ComputePedestal(const pmtana::Waveform_t& wf,
                         pmtana::PedestalMean_t& mean_v,
                         pmtana::PedestalSigma_t& sigma_v);

  private:
    size_t _sample_size; 
    double _threshold;   

    float _max_sigma;     
    float _ped_range_max; 
    float _ped_range_min; 

    bool _verbose;        
    int _n_wf_to_csvfile; 
    int _wf_saved = 0;
    int _num_presample;  
    int _num_postsample; 
    std::ofstream _csvfile;

    double CalcMean(const std::vector<double>& wf, size_t start, size_t nsample);

    double CalcStd(const std::vector<double>& wf,
                   const double ped_mean,
                   size_t start,
                   size_t nsample);

    bool CheckSanity(pmtana::PedestalMean_t& mean_v, pmtana::PedestalSigma_t& sigma_v);
  };
}
#endif
 // end of doxygen group