source: Sophya/trunk/SigPredictor/levsptsrc.h@ 3730

#ifndef LevSPtSrc_H
#define LevSPtSrc_H

#include <ostream>
#include "unitvector.h"

struct LevSPtSrc 
{
        /*      // Coordinates in Healpix convention, radians
        double theta;   // 0 ---->M_PI
        double phi;             // 0 ----> 2*M_PI
        */
        UnitVector VecPoint;
        float flux30G;
        float flux44G;
        float flux70G;
        float flux100G;
        float flux143G;
        float flux217G;
        float flux352G;
        float flux545G;
        float flux857G;
        LevSPtSrc(      double Phi, double Theta, 
                                        float F30G, float F44G, float F70G, 
                                        float F100G, float F143G, float F217G, float F352G, float F545G, float F857G
                                        );                                      
        LevSPtSrc();
        LevSPtSrc(double theta, double phi);
        ~LevSPtSrc() { }
        double Phi() const      {       return VecPoint.Phi();  }
        double Theta() const {  return VecPoint.Theta();}
        LevSPtSrc& operator= (LevSPtSrc const & a); 
        
        void print(ostream& flux) const
        {       flux<<VecPoint.Theta()<<'\t'<<VecPoint.Phi()<<'\t'<<flux30G;
                flux<<'\t'<<flux44G<<'\t'<<flux70G<<'\t'<<flux100G<<'\t'<<flux143G;
                flux<<'\t'<<flux217G<<'\t'<<flux352G<<'\t'<<flux545G<<'\t'<<flux857G<<endl;
        }
                
        void printHeader(ostream& flux) const
        {       flux<<"theta\tphi\tflux30G\tflux44G\tflux70G\tflux100G\tflux143G\tflux217G\tflux352G\tflux545G\tflux857G"<<endl;
                return;
        }       
};

// Relation d'ordre sur la sphere a utiliser dans les contenaires.
inline bool operator < (const LevSPtSrc& Vec1, const LevSPtSrc& Vec2)
{       if ( Vec1.Theta() < Vec2.Theta() )  return true;
/*      if ( Vec1.Theta() == Vec2.Theta() ) 
        {       if ( Vec1.Phi() < Vec2.Phi() ) return true;
        }
*/      
        return false;
}

/*
// Pour utiliser dans le conteneur

struct lessThLevSPtSrc
{
        inline bool operator () (const LevSPtSrc& Vec1, const LevSPtSrc& Vec2) const
        {       if ( Vec1.Theta() < Vec2.Theta() )  return true;
                if ( Vec1.Theta() == Vec2.Theta() ) 
                {       if ( Vec1.Phi() < Vec2.Phi() ) return true;
                }
                return false;
        }
};

*/

#endif