mc_antenna.h

Go to the documentation of this file.

 
/*===========================================================================
  === External declarations
  ===========================================================================*/
#ifndef mc_antenna_h
#define mc_antenna_h
 
#include "arts.h"
#include "matpack_data.h"
#include "rng.h"
#include <cmath>
#include <stdexcept>
 
enum AntennaType {
  ANTENNA_TYPE_PENCIL_BEAM = 1,
  ANTENNA_TYPE_GAUSSIAN = 2,
  ANTENNA_TYPE_LOOKUP = 3
};
 
struct MCAntenna {
  AntennaType atype{};
  Numeric sigma_aa{}, sigma_za{};
  Vector aa_grid, za_grid;
  Matrix G_lookup;
 
  MCAntenna()
      : aa_grid(),
        za_grid(),
        G_lookup() { /* Nothing to do here */
  }
 
  void set_pencil_beam();
  
  void set_gaussian(const Numeric& za_sigma, const Numeric& aa_sigma);
    
  void set_gaussian_fwhm(const Numeric& za_fwhm, const Numeric& aa_fwhm);
 
  void set_lookup(ConstVectorView za_grid,
                  ConstVectorView aa_grid,
                  ConstMatrixView G_lookup);
    
  void return_los(Numeric& wgt,
                  ConstMatrixView R_return,
                  ConstMatrixView R_enu2ant) const;
    
  void draw_los(VectorView sampled_rte_los,
                MatrixView R_los,
                RandomNumberGenerator<>& rng,
                ConstMatrixView R_ant2enu,
                ConstVectorView bore_sight_los) const;
 
  friend ostream& operator<<(ostream& os, const MCAntenna& mca);
};
 
void rotmat_enu(MatrixView R_ant2enu, ConstVectorView prop_los);
 
void rotmat_stokes(MatrixView R_pra,
                   const Index& stokes_dim,
                   const Numeric& f1_dir,
                   const Numeric& f2_dir,
                   ConstMatrixView R_f1,
                   ConstMatrixView R_f2);
 
#endif  // mc_antenna_h