simplepolymodel.h

#ifndef SIMPLEPOLYMODEL__H
#define SIMPLEPOLYMODEL__H

#include "eigenstuff.h"

class CcdImage;

#include "distortionmodel.h"
#include "gtransfo.h"
#include "simplepolymapping.h"
#include "projectionhandler.h"
#include <map>

class CcdImageList;

/* We deal here with coordinate transforms which are fitted
   and/or necessary to AstromFit. The classes SimplePolyModel 
and SimplePolyMapping implement a model where  there is one 
separate transfrom per CcdImage. One could chose other setups.

*/

/* This modeling of distortions can even accommodate images set mixing instruments */
class SimplePolyModel : public DistortionModel
{
  /* using ref counts here allows us to not write a destructor nor a copy
     constructor. I could *not* get it to work using std::auto_ptr. */
  typedef std::map<const CcdImage*, CountedRef<SimpleGtransfoMapping> > mapType;
  mapType _myMap;
  const ProjectionHandler* _sky2TP;

public :

  SimplePolyModel(const CcdImageList &L, 
          const ProjectionHandler* ProjH, 
          bool InitFromWCS,
          unsigned NNotFit=0);

  // The following routines are the interface to AstromFit
  const Mapping* GetMapping(const CcdImage &) const;

  unsigned AssignIndices(unsigned FirstIndex, std::string &WhatToFit);

  // dispaches the offsets after a fit step into the actual locations of parameters 
  void OffsetParams(const Eigen::VectorXd &Delta);
  
  const Gtransfo* Sky2TP(const Mapping* M, const CcdImage &C) const
  { return _sky2TP->Sky2TP(M,C);}

};



#endif /* SIMPLEPOLYMODEL__H */