photomfit.h

#ifndef PHOTOMFIT__H
#define PHOTOMFIT__H

#include <iostream>

#include "ccdimage.h"
#include "fittedstar.h"
#include <map>

#include "refzp.h"

class Associations;


class PhotomParam : public RefCount
{
    public :
    double fluxScale;
    int index;

    PhotomParam(const double &FS) : fluxScale(FS), index(-1) {};
    void SetIndex(const int I) { index=I;}
    int Index() const { return index;}
    
};

#include <list>
typedef std::list<CountedRef<PhotomParam> > PhotomParamList;

class PhotomFit 
{
private :
  RefZP refZP ;
 std::string photomBand; 
  Associations &assoc;
  double zpRef;
  CcdImageList refCcds;
  CcdImageList otherCcds;
  CcdImageList ccdImageList;
  FittedStarList &fittedStarList;
  PhotomParamList photomParamList;
  int nParTot;
  double chi2_measured, chi2_ref;
  int n_measured, n_ref;

  // pseudo pointer from CcdImage's to PhotomParam's
  std::map<const CcdImage *, PhotomParam*> imageParamMap;
  void SetPhotomParam(const CcdImage &C, PhotomParam &P);

  void init(); // routine to be called by different constructors


public :

  PhotomFit(Associations &A, const RefZP &R);

  PhotomFit(Associations &A);

  bool Minimize(const int MaxIter);

  PhotomParam& GetPhotomParam(const CcdImage &C) const;
  
  CcdImageList const& GetCcdImageList() const { return ccdImageList; }
  
  double ComputeChi2();
  void PrintChi2(const std::string &Message) const;

  void RemoveOutliers(const double &NSigCut);
  void DumpFluxScales(const std::string &FileNamePattern) const;
  void WriteCatalog(const std::string &FileName) const;
  
  void SaveResults(const std::string &Pattern);

  private :

  template <class Op> void loop_on_measurements(CcdImageList &List, Op &);

  void update_fluxes();
  void SetMags();

};


#endif /* PHOTOMFIT__H */