source: Sophya/trunk/SigPredictor/lobecartographmoyen.cc@ 975

                        // Dominique YVON, CEA/DAPNIA/SPP 02/2000

#include <math.h>
#ifdef __MWERKS__
   #include "mwerksmath.h"
   #include "unixmac.h"
   #include "macenvvariables.h"
#endif
#include "unitvector.h"
#include "lobecartographmoyen.h"

LobeCartoMoyen::LobeCartoMoyen(CarteLobe* pcart,double freqMin, double freqMax)
  :MeanFreqLobe(freqMin,freqMax), pCarte(pcart) {
  
        angleMax=pCarte->halfAngAperture();
        cosanglemax=cos(angleMax);
        sprintf(Name,"Lobe CartographiŽ TICRA");
}

LobeCartoMoyen::~LobeCartoMoyen(){      }
double LobeCartoMoyen::lobeResol() const{
  // Du fait des extrapolations paraboliques, on peut calculer le lobe en
  // tout point sans reconter de probleme de precision
  return pCarte->resol();
}

UnitVector LobeCartoMoyen::VecShift(const UnitVector& Vpointe, const UnitVector& VY) const {
  double alpha,beta;
  pCarte->UVToAng(pCarte->UCentre(),pCarte->VCentre(),alpha,beta);
  UnitVector Vec;
  Vec=Vpointe.Rotate(VY,alpha);
  Vec=Vec.Rotate(Vpointe,beta);
  return Vec;
}


double LobeCartoMoyen::weigthAmpl(const UnitVector& Vcourant,const UnitVector& VP, 
                const UnitVector& VY) const{
  UnitVector VCentreLobe;
  VCentreLobe=VecShift(VP, VY);
  double cosbidon=(Vcourant*VCentreLobe);
  if (cosanglemax>cosbidon) return 0.;
  else 
  //return 1.;  //Debug
  {   double cosalpha=Vcourant*VP;
        //  cout<<1-cosalpha*cosalpha<<endl;
          double sinalpha=sqrt(1-cosalpha*cosalpha);    // Approx. des petits angles
  
          Vector3d VProjVcourant;
          VProjVcourant=Vcourant-cosalpha*VP;
          VProjVcourant.Normalize();
          
          UnitVector VX=VY^VP;          
          // Produit vectoriel
          // On peut optimiser en sortant le calcul de VX de la fonction.
          double cosbeta=VProjVcourant*VX;
          double sinbeta=VProjVcourant*VY;

          return pCarte->Value(sinalpha*cosbeta,sinalpha*sinbeta);
          // Il faut comprendre les homogeneites de la carte BUGGGGGG  XXXXXXX
         
  }
  
}