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