ObjectLightCurve.h

// -*- mode: c++; mode: auto-fill; mode: flyspell-prog; -*-

#ifndef ObjectLightCurve_h
#define ObjectLightCurve_h 1

#include "luiza/Processor.h"

using namespace gloria;

namespace luiza {
  
  

  class ObjectLightCurve : public Processor {
    
  public:
    
    virtual Processor*  newProcessor() { return new ObjectLightCurve ; }
    
    
    ObjectLightCurve() ;
    
    virtual void init() ;


    virtual void startAnalysis( GloriaDataContainer* data ) ;
    
    virtual void processData( GloriaDataContainer* data ) ;
    
    
    virtual void check( GloriaDataContainer* data ) ; 
    
    

    virtual void concludeAnalysis( GloriaDataContainer* data ) ;

    virtual void end() ;
    
    
  protected:

    void printEndMessage() const ;

    void initializeLCTable( GloriaDataContainer* data );

    void initializeRefStars( GloriaDataContainer* data );

    string _inputCollectionName;

    string _outputCollectionName;

    std::vector<float> _objectPosition;

    double _matchingTolerance;

    std::vector<std::string> _positionColumns;


    std::string _brightnessColumn;


    bool _magnitudoGiven;

    std::string _polarAngleColumn;


    std::vector<float> _referenceStars;

    // Reference star file name(s) - no longer supported
    //    std::vector< std::string> _noninterference;

    string _referenceStarCatalog;


    string _referenceMagColumn;


    double _referenceMaxDist;

    int _maxReferenceCount;

    int _correctionFitOrder;

    bool _referenceUseMedian;

    bool _estimateError;

    std::vector<std::string> _copyColumns;

    std::vector<std::string> _headerRecords;


    int _nLoops ;
    int _nInputLists ;
    int _nObjectsMatched;
    int _nOutputPoints ;

    double xmin,xmax,ymin,ymax;

    GloriaFitsTable *_lightCurveTable;


    // Private code for matrix equation solving

    int GaussjSolve(double *alfa,double *beta,int n);

  } ;
  
} // end namespace luiza
#endif