PulseFinder.hh

Go to the documentation of this file.

#ifndef PULSEFINDER_h
#define PULSEFINDER_h

#include "BaseModule.hh"

class PulseFinder : public BaseModule{
public:
  PulseFinder();
  ~PulseFinder();
   
  int Initialize();
  int Finalize();
  int Process(EventPtr evt);
  
  int EvaluatePulse(Pulse& pulse, ChannelData* chdata,
                    int start_index, int end_index) const;
  
  void DiscriminatorSearch(ChannelData* chdata,
                           std::vector<int>& start_index, 
                           std::vector<int>& end_index);
  
  void VarianceSearch(ChannelData* chdata,
                      std::vector<int>& start_index, 
                      std::vector<int>& end_index);
  
  void IntegralSearch(ChannelData* chdata,
                      std::vector<int>& start_index,
                      std::vector<int>& end_index);
  
  void CurvatureSearch(ChannelData* chdata,
                       std::vector<int>& start_index,
                       std::vector<int>& end_index);
  
  static const std::string GetDefaultName(){ return "PulseFinder"; }
  
  enum SEARCH_MODE { VARIANCE , DISCRIMINATOR , INTEGRAL , CURVATURE };

private:
  //parameters
  bool align_pulses;               
  SEARCH_MODE mode;                
  bool normalize;                  
  int start_window;                
  double min_start_variance;       
  int min_resolution;              
  int min_start_peak_sep;          
  bool discriminator_relative;     
  double discriminator_value;      
  int discriminator_start_add;     
  int discriminator_end_add;       
  double integral_start_time;      
  double integral_end_time;        
  double integral_start_threshold; 
  double integral_end_threshold;   
  double min_sep_time;             
  double multipulse_thresh_value;  
  double amplitude_start_threshold;
  double amplitude_end_threshold;  
  double min_pulse_time;           
  double lookback_time;            


  //parameters for curvature search
  int down_sample_factor;          
  double pulse_start_curvature;    
  int pile_up_curvature;           
  double pulse_end_slope;          
  double fixed_time1;      
  double fixed_time2;      
  
  double PredictedNPE(double t0, double t, double fprompt, 
                      double tau1, double tau2)
  {
    return fprompt*(1.-exp(-(t-t0)/tau1)) + (1-fprompt)*(1.-exp(-(t-t0)/tau2));
  }

};

//override stream ops for SearchMode
std::istream& operator>>(std::istream& in, PulseFinder::SEARCH_MODE& m);
std::ostream& operator<<(std::ostream& out, const PulseFinder::SEARCH_MODE& m);

#endif