ChanTrigger.h

#ifndef CHANTRIGGER_HH_
#define CHANTRIGGER_HH_

//class ChanTrigger:

class RX;
class Anita;
using std::vector;
class Balloon;
class IceModel;
class Settings;
class Screen;
class GlobalTrigger;

class ChanTrigger {
    
 private:

  double ADCCountstoPowerThreshold(Anita *anita1, int ipol, int iant);

  const static int NSURF=9;                                     
  const static int NSURFPLUSONE=10;                             
  const static int NSURFMINUSONE=8;                             
  const static int NCHANNELS=32;                                
  const static int NPOINTS=4073;                                
  static const unsigned NFOUR = 1024;                           
  static const unsigned HALFNFOUR = 512;                        
                                   
  double thisrate;                                              
  double thispowerthresh;                                       
  double e_component;                                           
  double h_component;                                           
  double n_component;                                           
  double e_component_kvector;                                   
  double h_component_kvector;                                   
  double n_component_kvector;                                   
  double hitangle_e;                                            
  double hitangle_h;                                            
      
 public:

  ChanTrigger(); 

  

  void InitializeEachBand(Anita *anita1);


  void ApplyAntennaGain(Settings *settings1, Anita *anita1, Balloon *bn1, Screen *panel1, int ant, Vector &n_eplane, Vector &n_hplane, Vector &n_normal);


  void TriggerPath(Settings *settings1, Anita *anita1, int ant, Balloon *bn1);


  void DigitizerPath(Settings *settings1, Anita *anita1, int ant, Balloon *bn1);


  void TimeShiftAndSignalFluct(Settings *settings1, Anita *anita1, int ilayer, int ifold, double volts_rx_rfcm_lab_e_all[48][512], double volts_rx_rfcm_lab_h_all[48][512]);
  

  static void ConvertEHtoLREfield(double,double,double&,double&);


  static void ConvertEHtoLREnergy(double,double,double&,double&);


  void L1Trigger(Anita *anita1,double timedomain_output_1[5][Anita::NFOUR],double timedomain_output_2[5][Anita::NFOUR],double powerthreshold[2][5],int *channels_passing_e_forglob,int *channels_passing_h_forglob,int &npass);    



  double getRate();


  double rateToThreshold(double rate, int band);


  static double GetNoise(Settings *settings1,double altitude_bn,double geoid,double theta,double bw,double temp);


  void WhichBandsPass(Settings *settings1, Anita *anita1, GlobalTrigger *globaltrig1, Balloon *bn1, int ilayer, int ifold, double dangle, double emfrac, double hadfrac, double thresholds[2][5]);


  void WhichBandsPassTrigger1(Settings *settings1, Anita *anita1, GlobalTrigger *globaltrig1, Balloon *bn1, int ilayer, int ifold, double thresholds[2][5]);


  void WhichBandsPassTrigger2(Settings *settings1, Anita *anita1, GlobalTrigger *globaltrig1, Balloon *bn1, int ilayer, int ifold, double dangle, double emfrac, double hadfrac, double thresholds[2][5]);


  static double FindPeak(double *waveform,int n);

  

  void GetThresholds(Settings *settings1,Anita *anita1,int ilayer,double thresholds[2][5]); // get thresholds for this layer


  void DiodeConvolution(Settings *settings1, Anita *anita1, GlobalTrigger *globaltrig1, int ilayer, int ifold, double mindiodeconvl[5], double onediodeconvl[5], double psignal[5][Anita::NFOUR],  double timedomain_output[5][Anita::NFOUR],  double timedomain_output_justNoise[5][Anita::NFOUR], int ibinshift, int ipol, double thresholds[2][5]);



  void addToChannelSums(Settings *settings1,Anita *anita1,int ibw,int k);


  static int IsItUnmasked(unsigned short surfTrigBandMask[9][2],int ibw,int ilayer, int ifold, int ipol);


  void applyImpulseResponseDigitizer(Settings *settings1, Anita *anita1, int nPoints, int ant, double *x, double y[512], bool pol);


  void applyImpulseResponseTrigger(Settings *settings1, Anita *anita1, int ant, double y[512], double *vhz, bool pol);


  void getNoiseFromFlight(Settings* settings1, Anita* anita1, int ant, bool also_digi = true);


  void injectImpulseAfterAntenna(Anita *anita1, int ant);


  TGraph *getPulserAtAMPA(Anita *anita1, int ant);


  void saveTriggerWaveforms(Anita *anita1, double sig0[48], double sig1[48], double noise0[48], double noise1[48]);


  void saveDigitizerWaveforms(Anita *anita1, double sig0[48], double sig1[48], double noise0[48], double noise1[48]);


  void injectImpulseAtSurf(Anita *anita1, double volts_triggerPath_e[Anita::HALFNFOUR], double volts_triggerPath_h[Anita::HALFNFOUR], int ant);


  void calculateCW(Anita *anita1, double frequency, double phase, double amplitude);

  double applyButterworthFilter(double ff, double ampl, int notchStatus[3]);

  
  double vhz_rx[2][5][Anita::NFREQ];                           
  double volts_rx_forfft[2][5][Anita::HALFNFOUR];              
  vector<int> flag_e[5];                                       
  vector<int> flag_h[5];                                       
  double bwslice_volts_pol0[5];                                
  double bwslice_volts_pol1[5];                                
  double bwslice_energy_pol0[5];                               
  double bwslice_energy_pol1[5];                               
  double bwslice_volts_pol0_em[5];                             
  double bwslice_volts_pol1_em[5];                             
  double bwslice_volts_pole[5];                                
  double bwslice_energy_pole[5];                               
  double bwslice_volts_polh[5];                                
  double bwslice_energy_polh[5];                               
  double volts_rx_rfcm_lab[2][Anita::HALFNFOUR];               
  double volts_rx_rfcm_lab_all[2][48][Anita::HALFNFOUR];       
  double volts_rx_rfcm[2][Anita::HALFNFOUR];                   
  double justNoise_digPath[2][Anita::HALFNFOUR];               
  double justNoise_trigPath[2][Anita::HALFNFOUR];              
  double cw_digPath[2][Anita::HALFNFOUR];                      
  double justSig_trigPath[2][Anita::HALFNFOUR];                
  double justSig_digPath[2][Anita::HALFNFOUR];                
  
  // these are filled for triggerscheme==0 and triggerscheme==1
  // frequency domain voltage and energy based
  double signal_eachband[2][Anita::NBANDS_MAX];               
  double threshold_eachband[2][Anita::NBANDS_MAX];            
  double noise_eachband[2][Anita::NBANDS_MAX];                
  int passes_eachband[2][Anita::NBANDS_MAX];                  
  vector<double> vsignal_eachband[2];                         
  vector<double> vthreshold_eachband[2];                      
  vector<double> vnoise_eachband[2];                          
  vector<int>    vpasses_eachband[2];                         
    
  double v_banding_rfcm[2][5][Anita::NFREQ];                  
  double v_banding_rfcm_forfft[2][5][HALFNFOUR];              
  double vm_banding_rfcm_forfft[2][5][HALFNFOUR];             
  double vm_banding_rfcm_forfft_justNoise[2][5][HALFNFOUR];   
  double v_banding_rfcm_forfft_temp[2][5][HALFNFOUR];         
  double integral_vmmhz;                                      
  int unwarned;                                               
}; //class ChanTrigger

#endif