CorrelatedGaussianObservables.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2013 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef CORRELATEDGAUSSIANOBSERVABLES_H
#define CORRELATEDGAUSSIANOBSERVABLES_H
 
#include "Observable.h"
#include "gslpp.h"
#include <boost/ptr_container/ptr_vector.hpp>
#include <TMatrixDSym.h>
#include <TMatrixDSymEigen.h>
 
class ThObsFactory;
 
class CorrelatedGaussianObservables {
public:
    
    CorrelatedGaussianObservables(std::string name_i);
    
    CorrelatedGaussianObservables();
    
    CorrelatedGaussianObservables(const CorrelatedGaussianObservables& orig);
    
    virtual ~CorrelatedGaussianObservables();
 
    void ComputeCov(const TMatrixDSym& Corr);
 
    virtual double computeWeight();
 
    std::vector<double> computeLeaveOneOutWeights();
 
    void AddObs(Observable& Obs_i);
 
    std::vector<Observable> getObs() const
    {
        return Obs;
    }
    
    Observable getObs(int i) const
    {
        return (Obs.at(i));
    }
 
    std::string getName() const
    {
        return name;
    }
 
    const TMatrixDSym& getCov() const
    {
        return InvCov;
    }
    
    int ParseCGO(boost::ptr_vector<Observable>& Observables, 
                 std::ifstream& ifile, 
                 boost::tokenizer<boost::char_separator<char> >::iterator& beg, 
                 std::string& infilename, 
                 ThObsFactory& myObsFactory,
                 StandardModel * myModel,
                 int lineNo,
                 int rank);
    bool isEOF()
    {
        return IsEOF;
    }
    
    void setIsPrediction(bool IsPrediction_i)
    {
        IsPrediction = IsPrediction_i;
    }
    
    bool isPrediction()
    {
        return IsPrediction;
    }
    
    void setName(std::string name_i)
    {
        name = name_i;
    }
    
    void setCovarianceFromConfig (bool setInvCov)
    {
        covarianceFromConfig = setInvCov;
    }
    
private:
    std::vector<Observable> Obs;
    TMatrixDSym InvCov;
    TVectorD x;
    std::string name;
    std::string filepath;
    bool IsEOF;
    bool IsPrediction;
    bool covarianceFromConfig;
};
 
#endif  /* CORRELATEDGAUSSIANOBSERVABLES_H */