GloriaAverageMap.h

// -*- mode: c++; mode: auto-fill; mode: flyspell-prog; -*-


#ifndef GloriaAverageMap_h
#define GloriaAverageMap_h 1

#include "GloriaValueMap.h"
#include "gloria/GloriaFitsImage.h"

/*
 *   This source code is part of the LUIZA software package for GLORIA.
 *   You are free to use this source files for your own development as
 *   long as it stays in a public research context. You are not
 *   allowed to use it for commercial purpose. You must put this
 *   header with author names in all development based on this file.
 *
 */

namespace gloria {


  class GloriaAverageMap : public GloriaValueMap {

  public:

    GloriaAverageMap() : GloriaValueMap() { };


    GloriaAverageMap(int Nbx, int Nby, int Nvar = 1);



    GloriaAverageMap(const GloriaFitsImage* fitsImage, int Nbx, int Nby, int Nvar=1);


    GloriaAverageMap(GloriaFitsTable* mapTable) : GloriaValueMap(mapTable) { };


    void AddAverage();


    void AddAverage(string name);


    void SetAverageName(int ivar, string name);
 

    inline string GetAverageName(int ivar) {return _valueNames.at(2*ivar); } ;



    long CalculateMap(int iVar, int nVal, double* Xpos, double* Ypos, double* Val);


    long CalculateMap(int iVar, vector<double>* Xpos, vector<double>* Ypos, vector<double>* Val);



    long CalculateFitMap(int iVar, int nVal, double* Xpos, double* Ypos, double* Val, int Npol);


    long CalculateFitMap(int iVar, vector<double>* Xpos, vector<double>* Ypos, vector<double>* Val, int Npol);


    inline double GetAverage(int ivar, int ix, int iy, int ipol=1)
                         { return GetMapValue(2*ivar,ix,iy,ipol); } ;

    inline double GetAverage(int ivar, double x, double y, int ipol=1)
                         { return GetMapValue(2*ivar,x,y,ipol); } ;

    inline double GetAverage(int ivar, long ipix, int ipol=1)
                        { return GetMapValue(2*ivar,ipix,ipol); } ;

    inline double GetSpread(int ivar, int ix, int iy, int ipol=1)
                         { return GetMapValue(2*ivar+1,ix,iy,ipol); } ;

    inline double GetSpread(int ivar, double x, double y, int ipol=1)
                         { return GetMapValue(2*ivar+1,x,y,ipol); } ;

    inline double GetSpread(int ivar, long ipix, int ipol=1)
                        { return GetMapValue(2*ivar+1,ipix,ipol); } ;




    inline double* GetGloriaAverageMap(int ivar=0) {return _valueMaps.at(2*ivar); };

    inline double* GetSpreadMap(int ivar=0) {return _valueMaps.at(2*ivar+1); };


    inline  double GetAverageMean(int ivar=0){ return GetMapMean(2*ivar);} ;

    inline double GetSpreadMean(int ivar=0){ return GetMapMean(2*ivar+1);} ;


    inline GloriaFitsImage* GetAverageImage(int ivar=0, int ipol=1)
                    { return GetMapImage(2*ivar, ipol);}  ;

    inline GloriaFitsImage* GetSpreadImage(int ivar=0, int ipol=1)
                    { return GetMapImage(2*ivar+1, ipol);}  ;

    void InitializeMaps(int Nbx, int Nby, int Nvar);

  private:

    int  GaussjSolve(double *alfa,double *beta,int n);

  } ;


}

#endif