chip.h

#ifndef CHIP__H
#define CHIP__H



#include <gtransfo.h>
#include <frame.h>
#include <countedref.h>

#include <geomparam.h>

class CcdImage;



class Chip : public GeomParam{
  private :
    // what we gonna fit
    Gtransfo *fittedTransfo;

  const Gtransfo *CTP2TP;
  Gtransfo *ctp2Pix;
  Gtransfo *tP2Pix;

  /* in principle useless. but reduces numerical problems 
     when fitting polynomials  */
  GtransfoLinShift shiftInTangentPlane; 


  void update_transfos();

  int index; // provided and retrieved by AstromFit.


  public :

    Chip(const CcdImage &ccdImage, const int DistortionDegree);

  int NGroups() const {return 1;}

  int Npar(const int IGroup) const
    { if (IGroup) {}; return (fittedTransfo)? fittedTransfo->Npar() : 0;}
      

  int Index(const int IGroup) const
    { if (IGroup) {}; return index;}

  void SetIndex(const int &I, const int IGroup) 
      { if (IGroup) {};  index = I;};

  void ComputeH(const Point &Where, const int IGroup, double *H) const;

  void OffsetParams(const double *ParOffsets, const int IGroup);


  Point CTP2Pix(const Point &PointInCTP) const;

  virtual const Gtransfo* CTP2PixTransfo() const { return ctp2Pix;};



  void CTP2PixDerivative(const Point &Where, double *HxHy) const;

  ~Chip()
      {    
      delete fittedTransfo;
      delete ctp2Pix;
      delete tP2Pix;
      }


};


#include <list>

typedef CountedRef<Chip> ChipRef;
typedef std::list<CountedRef<Chip> > ChipList;
typedef ChipList::iterator       ChipIterator;
typedef ChipList::const_iterator ChipCIterator;




#endif /* CHIP__H */