constrainedpolymodel.h

#ifndef CONSTRAINEDPOLYMODEL__H
#define CONSTRAINEDPOLYMODEL__H

#include "eigenstuff.h"
#include "simplepolymapping.h"

class CcdImage;

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

class CcdImageList;

 
typedef unsigned ShootIdType;


class ConstrainedPolyModel : 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<TwoTransfoMapping> > mappingMapType;
  mappingMapType _mappings;
  typedef std::map<unsigned, CountedRef<SimpleGtransfoMapping> > chipMapType;
  chipMapType _chipMap;
  typedef std::map<ShootIdType, CountedRef<SimpleGtransfoMapping> > shootMapType;
  shootMapType _shootMap;
  const ProjectionHandler* _sky2TP;
  bool _fittingChips, _fittingShoots;


public :
  ConstrainedPolyModel(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);

  void DumpT2Transfos() const;
  
  const Gtransfo* Sky2TP(const Mapping* M, const CcdImage &C) const
  { return _sky2TP->Sky2TP(M,C);}

};



#endif /* CONSTRAINEDPOLYMODEL__H */