A class for Higgs experimental analyses. More...

#include <HiggsObservable.h>

Inheritance diagram for HiggsObservable:

Detailed Description

A class for Higgs experimental analyses.

Author: HEPfit Collaboration

Copyright: GNU General Public License

The class for building observables encoding Higgs experimental analyses, storing the parameters read from a file specified in the SomeModel.conf file or by the user. The names (thnames) of the observables have to correspond to the allowed names of observables listed in the ThFactory class.

Definition at line 28 of file HiggsObservable.h.

Public Member Functions
virtual double	computeWeight ()
	A method to compute the weight associated with the observable. More...

bool	getIsCorrelated ()

int	getNChannels ()

int	getNTheoryValues ()

void	getTheoryValues (std::vector< double > &theoryValues_i)

	HiggsObservable (const HiggsObservable &orig)

	HiggsObservable (const Observable &Obs)

bool	isNew ()
	A method to get the observable to the new parametric form. More...

boost::tokenizer< boost::char_separator< char > >::iterator &	ParseHiggsObservable (boost::tokenizer< boost::char_separator< char > >::iterator &beg, ThObsFactory &myObsFactory, StandardModel *myModel, int rank)
	the parser for HiggsObservables More...

void	setIsCorrelated (bool correlated)

void	setIsnew (bool isnew)
	A method to set the observable to the new parametric form. More...

virtual void	setParametricLikelihood (std::string filename, std::vector< ThObservable * > thObsV)
	Set the parametric likelihood describing one Higgs decay channel from a config file. More...

Public Member Functions inherited from Observable
double	computeTheoryValue ()
	A method to access the computed theory value of the observable. More...

virtual double	computeWeight (double th)
	A method to compute the weight associated with the observable. More...

virtual double	computeWeight (double th, double ave_i, double errg_i, double errf_i)
	A method to compute the weight associated with the observable. More...

virtual double	computeWeight (double th1, double th2)
	A method to compute the weight associated with the observable. More...

double	getAve () const
	A get method to access the average value of the observable. More...

std::string	getDistr () const
	A get method to access the name of the distribution of the observable. More...

double	getErrf () const
	A get method to access the flat error of the observable. More...

double	getErrg () const
	A get method to access the Gaussian error of the observble. More...

std::string	getFilename () const
	A get method to access the filename of the observables experimental likelihood file. More...

std::string	getHistoname () const
	A get method to access the name for the histogram of the observable. More...

std::string	getLabel () const
	A get method to access the label for the observable. More...

double	getMax () const
	A get method to access the maximum value of the observable. More...

double	getMin () const
	A get method to access the minimum value of the observable. More...

std::string	getName () const
	A get method to access the name of the observable. More...

std::string	getObsType () const
	A get method to get the Observable type. More...

std::string	getThname () const
	A get method to access the thname of the observable as defined in ThFactory class. More...

ThObservable *	getTho () const
	A get method to access the pointer to the object of the ThObservable class. More...

bool	isTMCMC () const
	A method to check if the observable is listed for MCMC. More...

bool	isWriteChain () const
	A method to check if the observable is listed for writing chains. More...

double	LogGaussian (double x, double ave, double sigma)

double	LogSplitGaussian (double x, double ave, double errl, double errr)

	Observable ()
	The default constructor. More...

	Observable (const Observable &orig)
	The copy constructor. More...

	Observable (const std::string name_i, const std::string thname_i, const std::string label_i, const bool tMCMC_i, const double min_i, const double max_i, ThObservable *tho_i)
	Constructor. More...

boost::tokenizer< boost::char_separator< char > >::iterator &	ParseObservable (std::string &type, boost::tokenizer< boost::char_separator< char > > *tok, boost::tokenizer< boost::char_separator< char > >::iterator &beg, std::string &filepath, std::string &infilename, int rank)
	The parser for Observables. More...

void	setAve (double ave)
	A set method to fix the average value of the observable. More...

void	setDistr (std::string distr)
	A set method to fix the name of the distribution of the observable. More...

void	setErrf (double errf)
	A set method to fix the flat error of the observable. More...

void	setErrg (double errg)
	A set method to fix the gaussian error of the observable. More...

void	setErrgl (double errgl)
	A set method to fix the left Gaussian error of the observable. More...

void	setErrgr (double errgr)
	A set method to fix the right Gaussian error of the observable. More...

void	setFilename (std::string filename_i)

void	setHasInverseCovariance (bool hasInverseCovariance)
	A set method to state that the Observable is a part of ObservablesWithInverseCovariance. More...

void	setHistoname (std::string histoname_i)
	A set method to set the name of the histogram containing the likelihood. More...

void	setLabel (std::string label)
	A set method to fix the label for the observable. More...

virtual void	setLikelihoodFromHisto (std::string filename, std::string histoname)
	A set method to set the likelihood from which the experimental likelihood of the observable will be read. More...

void	setMax (double max)
	A set method to fix the maximum value for the observable. More...

void	setMin (double min)
	A set method to fix the minimum value for the observable. More...

void	setName (std::string name)
	A set method to fix the name for the observable. More...

void	setObsType (std::string &obsType_s)
	A set method to set the Observable type. More...

virtual void	setParametricLikelihood (std::string filename)
	A set method to set a parametric likelihood reading parameters from a file. More...

void	setThname (std::string thname)
	A set method to fix the name of the observable as listed in ThFactory class. More...

void	setTho (ThObservable *tho_i)
	A set method to fix the pointer to object of type ThObservable. More...

void	setTho (ThObservable *tho_i, double bmin, double bmax)
	A set method to fix the pointer to object of type ThObservable. More...

void	setTMCMC (bool tMCMC)
	A set method to fix the observable's inclusion in the MCMC listing. More...

virtual	~Observable ()
	The default destructor. More...

Private Attributes
TMatrixD	channels
	A matrix holding the information of all the channels. More...

double	channelsize
	The number of channels in the Higgs Observable. More...

bool	covarianceFromConfig
	A flag for reading inverse covariance from config file. More...

gslpp::matrix< double > *	InvCov
	The inverse covariance matrix. More...

bool	isCorrelated
	A flag for correlated Higgs Observable. More...

bool	isnew
	A boolean which is true if the parametrization is new. More...

int	rank
	The rank of the process initializing this observable. More...

std::vector< double >	theoryValues
	A vector to contain the theoryvalues. More...

std::vector< ThObservable * >	thObsV
	A vector of ThObservables. More...

double	thobsvsize
	The size of the theory observables vector. More...

Additional Inherited Members
Protected Attributes inherited from Observable
double	ave
	The average value of the observable. More...

double	bin_max
	The maximum valus of the observable bin. More...

double	bin_min
	The minimum value of the observable bin. More...

std::string	distr
	The name of the distribution of the the observable. More...

double	errf
	The flat error of the observable. More...

double	errg
	The gaussian error of the observable. More...

double	errgl
	The lower gaussian error of the observable. More...

double	errgr
	The upper gaussian error of the observable. More...

std::string	filename
	The name of the file containing the experimental likelihood for the observable. More...

bool	hasInverseCovariance

std::string	histoname
	The name of the histogram for the observable. More...

TH1D *	inhisto
	1D Histogram containing the experimental likelihood for the observable More...

int	iterationNo
	A counter for the interation that helps with the observable caching. More...

std::string	label
	A label for the observable. More...

double	max
	The maximum valus of the observable. More...

double	min
	The minimum value of the observable. More...

std::string	name
	A name for the observable. More...

std::string	obsType
	Type of the Observable. 0: Observable, 1: HiggsObservable, 2: BinnedObservable, 3: FunctionObservable. More...

std::string	thname
	The name for the observable as fixed in the ThObservable class. More...

ThObservable *	tho
	A pointer of to the object of the ThObservables class. More...

double	thValue
	The theory value of the first observable. More...

bool	tMCMC
	The flag to include or exclude the observable from the MCMC run. More...

bool	writeChain
	The flag to write the chain for the observable from the MCMC run. More...

Constructor & Destructor Documentation

◆ HiggsObservable() [1/2]

HiggsObservable::HiggsObservable ( const Observable & Obs )

Definition at line 20 of file HiggsObservable.cpp.

                                                      :
Observable(Obs)
{
    isnew = true;
    isCorrelated = false;
    covarianceFromConfig = false;
};

◆ HiggsObservable() [2/2]

HiggsObservable::HiggsObservable ( const HiggsObservable & orig )

Definition at line 28 of file HiggsObservable.cpp.

{
    channels = orig.channels;
    thObsV = orig.thObsV;
    thobsvsize = orig.thobsvsize;
    channelsize = orig.channelsize;
    isnew = orig.isnew;
    isCorrelated = orig.isCorrelated;
    covarianceFromConfig = orig.covarianceFromConfig;
    InvCov = orig.InvCov;
    rank = orig.rank;
}

Member Function Documentation

◆ computeWeight()

double HiggsObservable::computeWeight ( )

virtual

A method to compute the weight associated with the observable.

Reimplemented from Observable.

Definition at line 154 of file HiggsObservable.cpp.

{
    double logprob = 0;
 
    if (isnew) {
        for (int i = 0; i < channels.GetNrows(); i++) {
            double mu = 0;
            for (unsigned int j = 0; j < thobsvsize; j++) {
                theoryValues.at(j) = thObsV.at(j)->computeThValue();
                mu += channels(i, j) * theoryValues.at(j);
            }
            theoryValues.at(thobsvsize) = tho->computeThValue();
            if (thObsV.at(0)->getModel().isModelLinearized())
                theoryValues.at(thobsvsize + i + 1) = -1. + mu + theoryValues.at(thobsvsize);
            else
                theoryValues.at(thobsvsize + i + 1) = mu * theoryValues.at(thobsvsize);
            if (!isCorrelated) logprob += LogSplitGaussian(theoryValues.at(thobsvsize + i + 1), channels(i, thobsvsize), channels(i, thobsvsize + 1), channels(i, thobsvsize + 2));
        }
        if (isCorrelated) {
            gslpp::vector<double> theObsVal(channels.GetNrows(), 0.);
            for (int i = 0; i < channels.GetNrows(); i++) {
                theObsVal(i) = theoryValues.at(thobsvsize + i + 1);
            }
            logprob =  (-0.5 * theObsVal * ((*InvCov) * theObsVal));
        }
    } else {
        for (int i = 0; i < channels.GetNrows(); i++) {
            double mu = 0, sum = 0.;
            for (unsigned int j = 0; j < thobsvsize - 1; j++) {
                theoryValues.at(j) = thObsV.at(j)->computeThValue();
                mu += channels(i, j) * theoryValues.at(j);
                sum += channels(i, j);
            }
            mu += (1. - sum) * thObsV.at(thobsvsize - 1)->computeThValue();
            theoryValues.at(thobsvsize) = tho->computeThValue();
            theoryValues.at(thobsvsize + i + 1) = mu * theoryValues.at(thobsvsize);
            logprob += LogSplitGaussian(theoryValues.at(thobsvsize + i + 1), channels(i, 3), channels(i, 4), channels(i, 5));
        }
    }
 
    return (logprob);
}

◆ getIsCorrelated()

bool HiggsObservable::getIsCorrelated ( )

inline

Definition at line 84 of file HiggsObservable.h.

    {
        return isCorrelated;
    }

◆ getNChannels()

int HiggsObservable::getNChannels ( )

inlinevirtual

Reimplemented from Observable.

Definition at line 79 of file HiggsObservable.h.

    {
        return channelsize;
    }

◆ getNTheoryValues()

int HiggsObservable::getNTheoryValues ( )

inlinevirtual

Reimplemented from Observable.

Definition at line 74 of file HiggsObservable.h.

    {
        return thobsvsize;
    }

◆ getTheoryValues()

void HiggsObservable::getTheoryValues ( std::vector< double > & theoryValues_i )

inlinevirtual

Reimplemented from Observable.

Definition at line 69 of file HiggsObservable.h.

    {
        theoryValues_i = theoryValues;
    }

◆ isNew()

bool HiggsObservable::isNew ( )

inlinevirtual

A method to get the observable to the new parametric form.

Returns: a boolean which is true if the observable is in the new parametric form

Reimplemented from Observable.

Definition at line 64 of file HiggsObservable.h.

    {
        return isnew;
    }

◆ ParseHiggsObservable()

boost::tokenizer< boost::char_separator< char > >::iterator & HiggsObservable::ParseHiggsObservable	(	boost::tokenizer< boost::char_separator< char > >::iterator &	beg,
		ThObsFactory &	myObsFactory,
		StandardModel *	myModel,
		int	rank
	)

the parser for HiggsObservables

Parameters

[in]	beg	the iterator that parses a line in the config file
[in]	myObsFactory	a reference to the Observables Factory
[in]	myModel	a pointer to the model
[in]	rank	the rank of the process that is using the parser

Returns: the iterator for the line being parsed

Definition at line 197 of file HiggsObservable.cpp.

{
    this->rank = rank;
    std::string type = "HiggsObservable";
    setObsType(type);
    std::vector<ThObservable*> hthobs;
    ++beg;
    distr = *beg;
    if (distr.compare("parametric") == 0) {
        setIsnew(false);
        ++beg;
        distr = *beg;
        if (distr.compare("LHC7") == 0) {
            hthobs.push_back(myObsFactory.CreateThMethod("ggH7", *myModel));
            hthobs.push_back(myObsFactory.CreateThMethod("VBF7", *myModel));
            hthobs.push_back(myObsFactory.CreateThMethod("VH7", *myModel));
            hthobs.push_back(myObsFactory.CreateThMethod("ttH7", *myModel));
        } else if (distr.compare("LHC8") == 0) {
            hthobs.push_back(myObsFactory.CreateThMethod("ggH8", *myModel));
            hthobs.push_back(myObsFactory.CreateThMethod("VBF8", *myModel));
            hthobs.push_back(myObsFactory.CreateThMethod("VH8", *myModel));
            hthobs.push_back(myObsFactory.CreateThMethod("ttH8", *myModel));
        } else if (distr.compare("TeV196") == 0) {
            hthobs.push_back(myObsFactory.CreateThMethod("ggH196", *myModel));
            hthobs.push_back(myObsFactory.CreateThMethod("VBF196", *myModel));
            hthobs.push_back(myObsFactory.CreateThMethod("VH196", *myModel));
            hthobs.push_back(myObsFactory.CreateThMethod("ttH196", *myModel));
        } else if (rank == 0)
            throw std::runtime_error("ERROR: wrong keyword " + distr + " in " + getName());
        setParametricLikelihood(*(++beg), hthobs);
    } else if (distr.compare("new_parametric") == 0) {
        std::string filename_h = *(++beg);
        std::ifstream ifile(filename_h.c_str());
        if (!ifile.is_open()) {
            if (rank == 0) throw std::runtime_error("\nERROR: " + filename_h + " does not exist. Make sure to specify a valid Higgs parameters configuration file.\n");
            else sleep(2);
        }
        std::string line_h;
        bool IsEOF_h;
        do {
            IsEOF_h = getline(ifile, line_h).eof();
            if (line_h.compare(0, 10, "CATEGORIES") == 0) {
                boost::char_separator<char> sep(" \t");
                boost::tokenizer<boost::char_separator<char> > mytok(line_h, sep);
                boost::tokenizer<boost::char_separator<char> >::iterator beg2 = mytok.begin();
                beg2++;
                while (beg2 != mytok.end()) {
                    std::string cat = *beg2;
                    hthobs.push_back(myObsFactory.CreateThMethod(cat, *myModel));
                    beg2++;
                }
            }
        } while (!IsEOF_h);
        if (hthobs.size() > 0)
            setParametricLikelihood(filename_h, hthobs);
        else {
            if (rank == 0) throw std::runtime_error("\nERROR: " + getName() + " does not provide at least one category\n");
            else sleep(2);
        }
    } else {
        if (rank == 0) throw std::runtime_error("ERROR: wrong distribution flag " + distr + " in " + getName());
        else sleep(2);
    }
 
    return beg;
}

◆ setIsCorrelated()

void HiggsObservable::setIsCorrelated ( bool correlated )

inline

Definition at line 89 of file HiggsObservable.h.

    {
        isCorrelated = correlated;
    }

◆ setIsnew()

void HiggsObservable::setIsnew ( bool isnew )

inline

A method to set the observable to the new parametric form.

Definition at line 55 of file HiggsObservable.h.

    {
        this->isnew = isnew;
    }

◆ setParametricLikelihood()

void HiggsObservable::setParametricLikelihood	(	std::string	filename,
		std::vector< ThObservable * >	thObsV
	)

virtual

Set the parametric likelihood describing one Higgs decay channel from a config file.

Parameters

filename	the name of the config file
thObsV	a vector of ThObservables containing the ratio of the production cross section for the specified categories in the model analyzed over the SM prediction

Reimplemented from Observable.

Definition at line 43 of file HiggsObservable.cpp.

{
    this->filename = filename;
    this->thObsV = thObsV;
    std::ifstream ifile(filename.c_str());
    if (!ifile.is_open())
        throw std::runtime_error("\nERROR: " + filename + " does not exist. Make sure to specify a valid Higgs parameters configuration file.\n");
    std::string line;
    bool IsEOF;
    bool readCorrelations = false;
    unsigned int i = 0; 
    unsigned int nrows = 0;
    do {
        IsEOF = getline(ifile, line).eof();
        if (line.compare(0, 11, "MEASUREMENT") == 0) nrows++;
    } while (!IsEOF);
    if (nrows == 0) {
        std::stringstream ss;
        ss << "\nERROR: " << filename << " should contain at least 1 measurement.\n";
        throw std::runtime_error(ss.str());
    }
    int implicit_tth = (isnew ? 0 : -1);
    channels.ResizeTo(nrows, thObsV.size() + 3 + implicit_tth);
    ifile.clear();
    ifile.seekg(0, ifile.beg);
    int lineNo = 0;
    do {
        IsEOF = getline(ifile, line).eof();
        lineNo++;
        if (*line.rbegin() == '\r') line.erase(line.length() - 1); // for CR+LF
        if (line.compare(0, 11, "MEASUREMENT") != 0 && !readCorrelations) {
            continue;
        } else {
            boost::char_separator<char> sep(" |\t");
            boost::tokenizer<boost::char_separator<char> > tok(line, sep);
            boost::tokenizer<boost::char_separator<char> >::iterator beg = tok.begin();
            if (isCorrelated) readCorrelations = true;
            // Read the necessary information from the config file. Each row contains:
            // ggH fraction
            // VBF fraction
            // VH fraction (ttH is computed as 1-ggH-VBF-VH)
            // average value of mu
            // left-side error
            // right-side error
            beg++; // skip label
            for (int j = 0; j < channels.GetNcols(); j++)
                channels(i, j) = atof((*(++beg)).c_str());
            if (isCorrelated && (std::fabs(channels(i, thObsV.size() + 1)) != std::fabs(channels(i, thObsV.size() + 2)))) {
                if (rank == 0) throw std::runtime_error("ERROR: Gaussian errors must be symmetric for CorrelatedHiggsObservables in " + filename + " at line number:" + boost::lexical_cast<std::string>(lineNo) + ".\n");
                else sleep(2);
            }
            i++;
            if (readCorrelations && (i == nrows)) {
                gslpp::matrix<double> myCorr(gslpp::matrix<double>::Id(nrows));
                int ni = 0;
                for (unsigned int irow = 0; irow < nrows; irow++) {
                    IsEOF = getline(ifile, line).eof();
                    if (line.empty() || line.at(0) == '#') {
                        if (rank == 0) throw std::runtime_error("ERROR: no comments or empty lines in CorrelatedHiggsObservables please! In file " + filename + " at line number:" + boost::lexical_cast<std::string>(lineNo) + ".\n");
                        else sleep(2);
                    }
                    lineNo++;
                    boost::tokenizer<boost::char_separator<char> > mytok(line, sep);
                    beg = mytok.begin();
                    int nj = 0;
                    for (unsigned int jcol = 0; jcol < nrows; jcol++) {
                        if ((*beg).compare(0, 1, "0") == 0
                                || (*beg).compare(0, 1, "1") == 0
                                || (*beg).compare(0, 1, "-") == 0 || (covarianceFromConfig == true)) {
                            if (std::distance(mytok.begin(), mytok.end()) < nrows) {
                                if (rank == 0) throw std::runtime_error(("ERROR: Correlation matrix is of wrong size in CorrelatedHiggsObservables: " + name + +" at line number:" + boost::lexical_cast<std::string>(lineNo) + ".\n").c_str());
                                else sleep(2);
                            }
                            myCorr(ni, nj) = atof((*beg).c_str());
                            nj++;
                            beg++;
                        } else {
                            if (rank == 0) throw std::runtime_error("ERROR: invalid correlation matrix for " + name + ". Check element (" + boost::lexical_cast<std::string>(ni + 1) + "," + boost::lexical_cast<std::string>(nj + 1) + ") in line number " + boost::lexical_cast<std::string>(lineNo) + " in file " + filename + ".\n");
                            else sleep(2);
                        }
                    }
                    ni++;
                }
                if (myCorr.size_i() != nrows || myCorr.size_j() != nrows)
                    throw std::runtime_error("The size of the correlated observables in " + name + " does not match the size of the correlation matrix!");
                InvCov = new gslpp::matrix<double>(nrows, nrows, 0.);
                if (covarianceFromConfig) {
                    for (unsigned int i = 0; i < nrows; i++)
                        for (unsigned int j = 0; j < nrows; j++)
                            (*InvCov)(i, j) = myCorr(i, j);
                } else {
                    for (unsigned int i = 0; i < nrows; i++)
                        for (unsigned int j = 0; j < nrows; j++)
                            (*InvCov)(i, j) = channels(i, thObsV.size() + 2) * myCorr(i, j) * channels(i, thObsV.size() + 2); // left = right imposed above
                    *InvCov = InvCov->inverse();
                }
                readCorrelations = false;
            }
        }
    } while (!IsEOF);
    if (i != nrows) {
        std::stringstream ss;
        ss << "\nERROR: " << filename << " should contain " << nrows << " measurements, but I have read " << i << " ones instead.\n";
        throw std::runtime_error(ss.str());
    }
 
    thobsvsize = thObsV.size();
    channelsize = channels.GetNrows();
    theoryValues.resize(thobsvsize + 1 + channelsize);
}

Member Data Documentation

◆ channels

TMatrixD HiggsObservable::channels

private

A matrix holding the information of all the channels.

Definition at line 109 of file HiggsObservable.h.

◆ channelsize

double HiggsObservable::channelsize

private

The number of channels in the Higgs Observable.

Definition at line 114 of file HiggsObservable.h.

◆ covarianceFromConfig

bool HiggsObservable::covarianceFromConfig

private

A flag for reading inverse covariance from config file.

Definition at line 116 of file HiggsObservable.h.

◆ InvCov

gslpp::matrix<double>* HiggsObservable::InvCov

private

The inverse covariance matrix.

Definition at line 117 of file HiggsObservable.h.

◆ isCorrelated

bool HiggsObservable::isCorrelated

private

A flag for correlated Higgs Observable.

Definition at line 115 of file HiggsObservable.h.

◆ isnew

bool HiggsObservable::isnew

private

A boolean which is true if the parametrization is new.

Definition at line 110 of file HiggsObservable.h.

◆ rank

int HiggsObservable::rank

private

The rank of the process initializing this observable.

Definition at line 118 of file HiggsObservable.h.

◆ theoryValues

std::vector<double> HiggsObservable::theoryValues

private

A vector to contain the theoryvalues.

Definition at line 112 of file HiggsObservable.h.

◆ thObsV

std::vector<ThObservable*> HiggsObservable::thObsV

private

A vector of ThObservables.

Definition at line 111 of file HiggsObservable.h.

◆ thobsvsize

double HiggsObservable::thobsvsize

private

The size of the theory observables vector.

Definition at line 113 of file HiggsObservable.h.

The documentation for this class was generated from the following files:

Detailed Description

Public Member Functions

Private Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ HiggsObservable() [1/2]

◆ HiggsObservable() [2/2]

Member Function Documentation

◆ computeWeight()

◆ getIsCorrelated()

◆ getNChannels()

◆ getNTheoryValues()

◆ getTheoryValues()

◆ isNew()

◆ ParseHiggsObservable()

◆ setIsCorrelated()

◆ setIsnew()

◆ setParametricLikelihood()

Member Data Documentation

◆ channels

◆ channelsize

◆ covarianceFromConfig

◆ InvCov

◆ isCorrelated

◆ isnew

◆ rank

◆ theoryValues

◆ thObsV

◆ thobsvsize