old_astromfit.h

#ifndef ASTROMFIT__H
#define ASTROMFIT__H

#include <string>
#include <iostream>
#include <map>

#include "ccdimage.h"
#include "chip.h"

class Associations;
class Vect;


enum FitKind { Total=0, FlipFlop=1};


class AstromFit {
  private : 

      Associations &assoc; 
  const CcdImageList &ccdImageList;
  FittedStarList &fittedStarList;

  unsigned int nParChips;
  unsigned int nParStars;
  unsigned int nParTot;
  // chi2 stuff
  double   chi2_measured, chi2_ref;
  int n_measured, n_ref;
  int maxNGroups;

  std::map<const CcdImage*, GeomParam*> imageParamMap; 

  void SetGeomParam(const CcdImage &C, GeomParam &G);

  void update_geom_params(Vect &B);

      void minimize_independant_geom_params();

  bool DoOneStep();

  bool FlipFlopStep();

  bool UpdateParams();

  public :
    AstromFit (Associations &A);
  
  bool Minimize(const FitKind, const int MaxIter =0);

  void SetStarIndexes();
  double ComputeChi2();
  void PrintChi2(const std::string &Message) const;

  GeomParam & GetGeomParam(const CcdImage & ccdImage) const;



  void SaveResults(const std::string &Tag);

};


#endif /* ASTROMFIT__H */