balloon.hh

#ifndef BALLOON_H
#define BALLOON_H


//class Balloon:

#include <iostream>

class Position;
class Vector;
class TChain;
class TTree;
class TFile;
class TTreeIndex;
class IceModel;
class Anita;
class Settings;
class Interaction;
class Ray;
class TH1F;
class TGraph;

using std::string;
using std::cout;


class Balloon {

private:
  //  string anitaliteflight; // the gps path of the anita-lite flight
  //string anitaflight;// gps path of anita flight
 
public:
  Balloon();

  
  // GPS positions of Anita-lite balloon flight
  int igps;                                                                  
  int ibnposition;
  double dtryingposition;                                                    
  static const int MAX_POSITIONS=50;                                         
  Vector slacpositions[MAX_POSITIONS];
  string sslacpositions[MAX_POSITIONS];
  int islacposition;
  TChain *flightdatachain;
  TTreeIndex *tindex;
  unsigned int realTime_flightdata_temp;                                     
  unsigned int realTime_flightdata;                                          
  float flatitude,flongitude,faltitude,fheading,froll, fpitch;
  double latitude,longitude,altitude,heading,roll,pitch;
  double MINALTITUDE;                                                        
  int igps_previous;                                                         
  int REDUCEBALLOONPOSITIONS;                                                
  int WHICHPATH;                                                             
  int RANDOMIZE_BN_ORIENTATION;                                              
  double BN_ALTITUDE;                                                        
  unsigned short surfTrigBandMask[9][2];                                     
  float powerthresh[9][32];                                                  
  float meanp[9][32];                                                        
  double altitude_bn;
  double theta_bn;
  double phi_bn;                                                             
  Position r_bn;                                                             
  double horizcoord_bn;                                                      
  double vertcoord_bn;                                                       
  Position r_boresights[Anita::NLAYERS_MAX][Anita::NPHI_MAX];                
  Vector x_axis_rot;
  Vector y_axis_rot;
  Vector z_axis_rot;
  Vector n_bn;                                                               
  Vector n_east;                                                             
  Vector n_north;                                                            
  double surface_under_balloon;                                              
  Position r_bn_shadow;                                                      
  double MAXHORIZON;                                                         
  double phi_spin;                                                           
  int NPOINTS;                                                               
  int NPOINTS_MIN;                                                           
  int NPOINTS_MAX;
  double latitude_bn_anitalite[100000];                                      
  double longitude_bn_anitalite[100000];                                     
  double altitude_bn_anitalite[100000];                                      
  double heading_bn_anitalite[100000];                                       
  double realtime_bn_anitalite[100000];                                      
  double BN_LONGITUDE;                                                       
  double BN_LATITUDE;                                                        
  double min_time; 
  double max_time; 



  void setObservationLocation(Interaction *interaction1,int inu,IceModel *antarctic,Settings *settings1);
    

  void GetBoresights(Settings *settings1,Anita *anita1,Position r_boresights[Anita::NLAYERS_MAX][Anita::NPHI_MAX]);
    

  void PickDownwardInteractionPoint(Interaction *interaction1,Anita *anita1,Settings *settings1,IceModel *antarctica1,
                    Ray *ray1, int &beyondhorizon, double len_int_km2, Vector * force_dir = 0); 
  

  void InitializeBalloon();
    

  void ReadAnitaliteFlight();
    

  void CenterPayload(double& hitangle_e);

    

  void PickBalloonPosition(Vector straightup,IceModel *antarctica,Settings *settings1, Anita *anita1);
    

  void PickBalloonPosition(IceModel *antarctica1,Settings *settings1,int inu,Anita *anita1, double randomNumber);


  int Getibnposition();


  double GetBalloonSpin(double heading);

    

  void GetAntennaOrientation(Settings *settings1,
                 Anita *anita1, 
                 int ilayer, int ifold, 
                 Vector& n_eplane,
                 Vector& n_hplane, 
                 Vector& n_normal); 


  void GetEcompHcompEvector(Settings *settings1,
                Vector n_eplane, 
                Vector n_hplane, 
                const Vector n_pol, 
                double& e_component, 
                double& h_component, 
                double& n_component);
        

  void GetEcompHcompkvector(Vector n_eplane,
                Vector n_hplane, 
                Vector n_normal, 
                const Vector n_exit2bn, 
                double& e_component_kvector, 
                double& h_component_kvector, 
                double& n_component_kvector);
    

  void GetHitAngles(
            double e_component_kvector,
            double h_component_kvector,
            double n_component_kvector, 
            double& hitangle_e,
            double& hitangle_h); 
 

  void SetDefaultBalloonPosition(IceModel *antarctica1);


  void setr_bn(double latitude,double longitude);


  void AdjustSlacBalloonPosition(int inu);
    

  void GetSlacPositions(Anita *anita1);


  void GetBoresights(Settings *settings1,Anita *anita1);
    

  void calculate_antenna_positions(Settings *settings1,Anita *anita1);
    

  Vector RotatePayload(Vector ant_pos);
    

  Vector unRotatePayload(Vector ant_pos);


   
}; //class Balloon


#endif