associations.h

// -*- C++ -*-
// 
#ifndef ASSOCIATIONS__H
#define ASSOCIATIONS__H

#include <string>
#include <iostream>

#include "stringlist.h"
#include "refstar.h"
#include "ccdimage.h"
#include "measuredstar.h" // for CatalogLoader
#include "point.h"


class Associations {
  public: 

  StringList imageNames;//  the names of individual ccd images
  CcdImageList ccdImageList; // the catalog handlers
  RefStarList refStarList;// the (e.g.) USNO stars
  RefStarList photRefStarList; // the (e.g) Landolt stars
  FittedStarList fittedStarList; //  the std::list of stars that are going to be fitted

    // fit cuts and stuff:
  Point commonTangentPoint;

  const CatalogLoader * load_it;

  void AssignMags();

public:


  Associations(const CatalogLoader  & LoadIt = *CatalogLoader::DefaultCatalogLoader);

  void SetCommonTangentPoint(const Point &CommonTangentPoint)
    { commonTangentPoint = CommonTangentPoint;};

  Point CommonTangentPoint() const { return commonTangentPoint;}

  bool AddImage(const std::string &ReducedImageName);

  bool AddImage(const ReducedImage &Ri);

  void AssociateCatalogs(const double MatchCutInArcSec = 0, 
               const bool UseFittedList = false,
               const bool EnlargeFittedList = true);

  void RefitWCSs(const unsigned Order);


  void CollectRefStars(const bool ProjectOnTP=true);



    

  void CollectMCStars(int realization = -1);
  //    void CheckMCStars(); // DELETE THIS METHOD
  
  void CollectPhotometricRefStars(std::string const& catalogname);
  void AssociatePhotometricRefStars(double MatchCutInArcSec);
  unsigned int GetNbAssociatedPhotometricRefStars() const { return nb_photref_associations; }

  void DeprojectFittedStars();

  /* If Color is "g-i", then the color is assigned from columns "g" and "i" of the colored catalog. */
  void SetFittedStarColors(std::string DicStarListName, 
                       std::string Color,
               const double &MatchCutArcSec);
  
  //    void SetRefPhotFactor(int chip, double photfact);
  
  void SelectFittedStars();
  
  void RefineFittedStars();

  void check_refcounts();
  
  void PrintStats() const;
  
  const CcdImageList& TheCcdImageList() const {return ccdImageList;}
  
  void SetFittedStarList(const FittedStarList &L);
  
  void AssignIndexes();
  
  unsigned int NShoots() const { return nshoots_; }

  unsigned NBands() const {return 1;}

  
private:
  void AssociateRefStars(const double &MatchCutInArcSec, const Gtransfo *T);
  unsigned int nshoots_;
  unsigned int nb_photref_associations;
};

#endif /* ASSOCIATIONS__H */