HiggsObservable.h

Go to the documentation of this file.

/*
 * Copyright (C) 2014 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef HIGGSOBSERVABLE_H
#define HIGGSOBSERVABLE_H
 
#include "Observable.h"
#include <TMatrixD.h>
#include "gslpp.h"
#include <boost/tokenizer.hpp>
 
class ThObsFactory;
class HiggsObservable : public Observable {
public:
    HiggsObservable(const Observable& Obs);
               
    HiggsObservable(const HiggsObservable& orig);
 
//    virtual ~HiggsObservable();
    
    // Read the necessary information from the config file. First row should 
    // contain the list of categories used in the analysis
 
    virtual void setParametricLikelihood(std::string filename, std::vector<ThObservable*> thObsV);
 
    virtual double computeWeight();
    
    void setIsnew(bool isnew)
    {
        this->isnew = isnew;
    }
    
    bool isNew()
    {
        return isnew;
    }
    
    void getTheoryValues(std::vector<double>& theoryValues_i) 
    {
        theoryValues_i = theoryValues;
    }
    
    int getNTheoryValues() 
    {
        return thobsvsize;
    }
    
    int getNChannels() 
    {
        return channelsize;
    }
    
    bool getIsCorrelated()
    {
        return isCorrelated;
    }
    
    void setIsCorrelated(bool correlated)
    {
        isCorrelated = correlated;
    }
    
    
    boost::tokenizer<boost::char_separator<char> >::iterator & ParseHiggsObservable(boost::tokenizer<boost::char_separator<char> >::iterator & beg, 
                                                                                    ThObsFactory& myObsFactory,
                                                                                    StandardModel * myModel,
                                                                                    int rank);
 
    private:
        TMatrixD channels;
        bool isnew;
        std::vector<ThObservable*> thObsV;
        std::vector<double> theoryValues;
        double thobsvsize;
        double channelsize;
        bool isCorrelated;
        bool covarianceFromConfig;
        gslpp::matrix<double>* InvCov;
        int rank;
};
 
#endif  /* HIGGSOBSERVABLE_H */