twotransfomapping.h

#ifndef TWOTRANSFOMAPPING__H
#define TWOTRANSFOMAPPING__H

#include "memory"

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

class TwoTransfoMapping: public Mapping
{
  SimpleGtransfoMapping *_m1, *_m2;
  unsigned _nPar1, _nPar2;
  struct tmpVars // just there to get around constness issues
  {
    Eigen::MatrixX2d h1,h2;
    Eigen::Matrix2d dt2dx;
  };

  std::auto_ptr<tmpVars> tmp;

 public : 
 TwoTransfoMapping(SimpleGtransfoMapping *ChipM, 
           SimpleGtransfoMapping *ShootM);
 unsigned Npar() const;

 void GetMappingIndices(std::vector<unsigned> &Indices) const;

 void ComputeTransformAndDerivatives(const FatPoint &Where,
                     FatPoint &OutPos,
                     Eigen::MatrixX2d &H)   const;
 void TransformPosAndErrors(const FatPoint &Where,
                FatPoint &OutPos) const;

 const Gtransfo& T1() const
 {  return _m1->Transfo();}

 const Gtransfo& T2() const
 {  return _m2->Transfo();}



 void  PosDerivative(Point &Where, Eigen::Matrix2d &Der, const double & Eps) const;

 private:

 friend class ConstrainedPolyModel;
 void SetWhatToFit(const bool FittingT1, const bool FittingT2);

};


#endif