ProfilePCA.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2008 by Paul Demorest
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
#ifndef __Pulsar_ProfilePCA_h
#define __Pulsar_ProfilePCA_h
 
#include "Pulsar/Algorithm.h"
 
namespace Pulsar
{
  class Profile;
  class Integration;
 
  class ProfilePCA : public Algorithm
  {   
  public:
 
    ProfilePCA ();
 
    ~ProfilePCA();
 
    void set_nharm_pca(unsigned nharm);
 
    void set_nharm_cov(unsigned nharm);
 
    unsigned get_nharm_pca();
 
    unsigned get_nharm_cov();
 
    void reset();
 
    void add_Profile (const Profile* p);
 
    void add_Integration (const Integration* i);
 
    void compute();
 
    double get_cov_value(unsigned i, unsigned j);
 
    double get_pc_value(unsigned i);
 
    Profile* get_pc_vector(int i, unsigned nbin);
 
    std::vector<double> decompose(const Profile *p, unsigned n_pc);
 
    void unload(const std::string& filename);
 
    static ProfilePCA *load(const std::string& filename);
 
  protected:
 
    unsigned nharm_cov;
 
    unsigned nharm_pca;
 
    double *cov;
 
    double *mean;
 
    unsigned nprof;
 
    double wt_sum;
 
    double wt2_sum;
 
    std::vector<double> pc_values;
 
    double *pc_vectors;
 
  };
 
}
 
 
#endif