AnalysisConfig.h

#ifndef UCORRELATOR_ANALYSIS_CONFIG_H
#define UCORRELATOR_ANALYSIS_CONFIG_H

#include "AnitaConventions.h" 
class TH2; 


namespace AnitaResponse
{
class DeconvolutionMethod; 
}

namespace UCorrelator
{


  class AnalysisConfig
  {
    public: 


      AnalysisConfig(const char * config_file = 0); 
      ~AnalysisConfig();

      void loadFromFile(const char * config_file); 

      unsigned correlator_nphi; 
      unsigned correlator_ntheta;  
      double correlator_theta_lowest;  
      double correlator_theta_highest;  
      bool enable_group_delay;  

      int zoomed_nphi; 
      int zoomed_ntheta; 
      double zoomed_dphi;  
      double zoomed_dtheta;  
      int zoomed_nant;  

      unsigned combine_nantennas;  
      unsigned combine_npad;  

      double saturation_threshold; 

      AnitaPol::AnitaPol_t start_pol;  
      AnitaPol::AnitaPol_t end_pol;  

      double peak_isolation_requirement; 
      double max_peak_trigger_angle; 
      int min_peak_distance_from_unmasked; 

      enum FinePeakFindingOption_t
      {
        FinePeakFindingAbby,  
        FinePeakFindingBicubic, 
        FinePeakFindingGaussianFit, 
        FinePeakFindingQuadraticFit9, 
        FinePeakFindingQuadraticFit16, 
        FinePeakFindingQuadraticFit25, 
        FinePeakFindingHistogram, 
      } fine_peak_finding_option; 

      static const char * getPeakFindingString(FinePeakFindingOption_t opt); 

      enum ResponseOption_t
      {
        ResponseNone ,  
        ResponseSingleBRotter, 
        ResponseIndividualBRotter, 
        ResponseHarmSignalOnly, 
        ResponseTUFF, 
        ResponseA4NoNotches, 
        ResponseA4, 
        ResponseA4Old, // the old A4 impulse response pre UH 2019 measurements.
        ResponseCustomString // make sure you set response_string; 
      } response_option;  

      const char * response_string; 

      static const char * getResponseString(ResponseOption_t opt); 
      int response_npad; //number of times to pad (in freq domain) the response 


      int nmaxima; 
      bool use_bin_center; 

      double bw_ndb; 
      double spectral_fit_start; 
      double spectral_fit_stop; 

      double noise_estimate_t0;  
      double noise_estimate_t1;  

      bool scale_by_cos_theta;  
      bool use_offline_mask; 
      double baseline_weight; 
      double correlation_gain_correction ; // off-boresight gain correction sigma used in making maps. Default is 0, which disables it


      //payload blast cuts (if 0, then ignored) 
      double max_mean_power_filtered;  
      double max_median_power_filtered; 
      double max_bottom_to_top_ratio; 

      //maximum amount theta can be adjusted downward (default, 3 degrees) 
      double max_theta_adjustment; 

      class Pulser
      {
        public: 
          Pulser(double offset, double distance, double dt) : GPS_offset(offset), max_distance(distance), max_dt(dt)  { ; } 
          double GPS_offset; 
          double max_distance;
          double max_dt; 
      }; 

      Pulser wais_hpol, wais_vpol, siple_hpol, siple_vpol; 

      TH2* ldb_hist() const; 
      int ldb_max_run; 

      bool fill_blast_fraction; 

      bool compute_shape_parameters;
      
      bool trace_to_continent;

      bool set_bottom_first;
      
      bool delay_to_center;

      bool r_time_shift_correction; //if delaying to center, apply a correction for the effective different R's 
      
      bool simulation_time_shift_correction; //make the correction for different R's correct for simulated events (to fix polarization calculations).  Need r_time_shift_correction ON for this to work 
      
      double cross_correlate_hv; 

      AnitaResponse::DeconvolutionMethod * deconvolution_method; 
      
      bool only_use_usable;
      
      int use_antenna_level_snr;

      double stokes_fracI; 
      

      bool use_forced_trigger_rms; 
      
      bool use_coherent_spectra; 

      bool use_hilbert_for_antenna_average; 

      double combine_t0;
      double combine_t1;

      bool fill_channel_info; 
    private: 
      mutable TH2 * the_ldb_hist; 
  };
}


#endif