causticdata.h

//
//  causticdata.h
//  SLsimLib
//
//  Created by bmetcalf on 5/6/13.
//
//
 
#ifndef __SLsimLib__causticdata__
#define __SLsimLib__causticdata__
 
#include <iostream>
#include "standard.h"
#include "utilities_slsim.h"
#include "simpleTree.h"
#include "simpleTreeVec.h"
#include "grid_maintenance.h"
 
/*
struct CausticSummary{
  CausticSummary(){}
  CausticSummary(const CausticSummary &tmp){
    redshift = tmp.redshift;
    crit_center = tmp.crit_center;
    crit_radius[0] = tmp.crit_radius[0];
    crit_radius[1] = tmp.crit_radius[1];
    crit_radius[2] = tmp.crit_radius[2];
    crit_area = tmp.crit_area;
    crit_type = tmp.crit_type;
    
    caustic_center = tmp.caustic_center;
    caustic_radius[0] = tmp.caustic_radius[0];
    caustic_radius[1] = tmp.caustic_radius[1];
    caustic_radius[2] = tmp.caustic_radius[2];
    caustic_area = tmp.caustic_area;
    
  };
  
  CausticSummary &operator=(const CausticSummary &tmp){
    if(&tmp==this) return *this;
    
    redshift = tmp.redshift;
    crit_center = tmp.crit_center;
    crit_radius[0] = tmp.crit_radius[0];
    crit_radius[1] = tmp.crit_radius[1];
    crit_radius[2] = tmp.crit_radius[2];
    crit_area = tmp.crit_area;
    crit_type = tmp.crit_type;
    
    caustic_center = tmp.caustic_center;
    caustic_radius[0] = tmp.caustic_radius[0];
    caustic_radius[1] = tmp.caustic_radius[1];
    caustic_radius[2] = tmp.caustic_radius[2];
    caustic_area = tmp.caustic_area;
    
    return *this;
  }
  double redshift;
  Point_2d crit_center;
  double crit_radius[3];
  double crit_area;
  CritType crit_type;
  Point_2d caustic_center;
  double caustic_radius[3];
  double caustic_area;
};
 

  
public:
  CausticDataStore(std::string filename,bool verbose = false);
  CausticDataStore(std::vector<ImageFinding::CriticalCurve> &crticurve_vec);
  CausticDataStore(const CausticDataStore &input);
  ~CausticDataStore();
 
 
  void addcrits(std::vector<ImageFinding::CriticalCurve> &crticurve_vec);
  
  double getTotalCritArea(){return totalcritarea;}
  double getTotalCausticArea(){return totalcausticarea;}
  
  void printfile(std::string filename,std::string paramfile,double fieldofview,double minscale);

  void resize(size_t size){ data.resize(size);}
 
  size_t numberOfCaustics(){return data.size();}
  
  CausticSummary & operator[](size_t index){return data[index];}
  
  CausticSummary & CritAreaOrder(size_t n){return data[crit_area_index[n]];}
  CausticSummary & CausticAreaOrder(size_t n){return data[caus_area_index[n]];}
  CausticSummary & CritRadiusOrder(size_t n){return data[caus_area_index[n]];}
  
  size_t getNthIndexCritArea(size_t n){return crit_area_index[n];}
  size_t getNthIndexCaustArea(size_t n){return caus_area_index[n];}
  size_t getNthIndexCritRadius(size_t n){return crit_radius_index[n];}
 
  
  //void SortByCritSize();
  //void SortByCritArea();
  //void SortByCausticArea();
  
  size_t init_for_random(
        short type              /// select according to: (1) critical curve area, (2) caustic curve area
        ,double limit = 0.0     /// minimum accepted area
  );
  */
 
//std::ostream &operator<<(std::ostream &os, CausticSummary const &caust);
 
#endif /* defined(__SLsimLib__causticdata__) */