A class for analyzing observables pairwise. More...

#include <Observable2D.h>

Inheritance diagram for Observable2D:

Detailed Description

A class for analyzing observables pairwise.

Author: HEPfit Collaboration

Copyright: GNU General Public License

The class for building a pair of observables and storing their different parameters read from the SomeModel.conf file or specified by the user. The names (thname) of the observables have to correspond to the allowed names of observables listed in the ThFactory class.

Definition at line 24 of file Observable2D.h.

Public Member Functions
double	computeTheoryValue2 ()
	A method to access the computed theory value of the second observable.

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

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

double	getAve2 () const
	A get method to access the average value of the second observable.

double	getErrf2 () const
	A get method to access the flat error of the second observable.

double	getErrg2 () const
	A get method to access the Gaussian error of the second observble.

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

double	getMax2 () const
	A get method to access the maximum value of the second observable.

double	getMin2 () const
	A get method to access the minimum value of the second observable.

std::string	getObsType2 () const

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

ThObservable *	getTho2 () const
	A get method to access the pointer to the object of the ThObservable class for the second observable.

bool	isEOF ()
	A method to check if the end of file has been reached.

	Observable2D ()
	The default constructor.

	Observable2D (const Observable &o1d)
	A conversion constructor. Constructs Observable2D with just one observable.

	Observable2D (const Observable2D &orig)
	The copy constructor.

	Observable2D (const std::string name_i, const std::string thname_i, const std::string thname2_i, const std::string label_i, const std::string label2_i, const bool tMCMC_i, const double min_i, const double max_i, const double min2_i, const double max2_i, ThObservable tho_i, ThObservable tho2_i)
	Constructor.

int	ParseObservable2D (std::string &type, boost::tokenizer< boost::char_separator< char > > *tok, boost::tokenizer< boost::char_separator< char > >::iterator &beg, std::string &infilename, std::ifstream &ifile, int lineNo, int rank)

void	setAve2 (double ave2)
	A set method to fix the average value of the second observable.

void	setErrf2 (double errf2)
	A set method to fix the flat error of the second observable.

void	setErrg2 (double errg2)
	A set method to fix the gaussian error of the observable.

void	setLabel2 (std::string label2)
	A set method to fix the label for the second observable.

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.

void	setMax2 (double max2)
	A set method to fix the maximum value for the second observable.

void	setMin2 (double min2)
	A set method to fix the minimum value for the second observable.

void	setObsType2 (std::string &obsType2_i)

void	setThname2 (std::string thname2)
	A set method to fix the name of the second observable as listed in ThFactory class.

void	setTho1Tho2 (ThObservable tho1_i, ThObservable tho2_i)
	A set method to fix the pointer to object of type ThObservable class for the second observable.

void	setTho2 (ThObservable *tho2_i)
	A set method to fix the pointer to object of type ThObservable class for the second observable.

virtual	~Observable2D ()
	The default destructor.

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

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

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

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

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

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

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

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

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

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

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

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

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

virtual int	getNChannels ()

virtual int	getNTheoryValues ()

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

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

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

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

virtual bool	isNew ()

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

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

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

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

	Observable ()
	The default constructor.

	Observable (const Observable &orig)
	The copy constructor.

	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.

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.

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

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

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

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

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

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

void	setFilename (std::string filename_i)

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

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

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

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

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

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

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

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

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

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

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

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

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

virtual	~Observable ()
	The default destructor.

Private Attributes
double	ave2
	The average value of the second observable.

std::vector< double >	bin_max
	The maximum value of the bin.

std::vector< double >	bin_min
	The minimum value of the bin.

double	errf2
	the flat error of the second observable.

double	errg2
	The Gaussian error of the second observable.

double	errgl2
	The left Gaussian error of the second observable.

double	errgr2
	The right Gaussian error of the second observable.

std::string	filepath
	The path to the file being parsed.

TH2D *	inhisto2d
	2D Histogram containing the experimental likelihood for the observable.

bool	IsEOF
	A bolean that is true if the end of file is reached.

int	iterationNo2
	Counts the iteration to help with caching.

std::string	label2
	A label for the second observable.

double	max2
	The maximum valus of the second observable.

double	min2
	The minimum value of the second observable.

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

std::string	thname2
	The name for the second observable as fixed in the ThObservable() class.

ThObservable *	tho2
	A pointer to an object of the ThObservable class.

double	thValue2
	The theory value of the second observable.

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

double	bin_max
	The maximum valus of the observable bin.

double	bin_min
	The minimum value of the observable bin.

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

double	errf
	The flat error of the observable.

double	errg
	The gaussian error of the observable.

double	errgl
	The lower gaussian error of the observable.

double	errgr
	The upper gaussian error of the observable.

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

bool	hasInverseCovariance

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

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

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

std::string	label
	A label for the observable.

double	max
	The maximum valus of the observable.

double	min
	The minimum value of the observable.

std::string	name
	A name for the observable.

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

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

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

double	thValue
	The theory value of the first observable.

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

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

Constructor & Destructor Documentation

◆ Observable2D() [1/4]

Observable2D::Observable2D	(	const std::string	name_i,
		const std::string	thname_i,
		const std::string	thname2_i,
		const std::string	label_i,
		const std::string	label2_i,
		const bool	tMCMC_i,
		const double	min_i,
		const double	max_i,
		const double	min2_i,
		const double	max2_i,
		ThObservable *	tho_i,
		ThObservable *	tho2_i
	)

Constructor.

Parameters

[in]	name_i	a given name for the observable pair
[in]	thname_i	the thname for the first observable fixed in ThFactory
[in]	thname2_i	the thname for the second observable fixed in ThFactory
[in]	label_i	the label assigned to the first observable
[in]	label2_i	the label assigned to the second observable
[in]	tMCMC_i	boolean flag to indicate inclusion in MCMC
[in]	min_i	minimum value for the first observable
[in]	max_i	maximum value for the first observable
[in]	min2_i	minimum value for the second observable
[in]	max2_i	maximum value for the second observable
[in]	tho_i	a pointer to an object of type ThObservable for the first observable
[in]	tho2_i	a pointer to an object of type ThObservable for the second observable

Definition at line 15 of file Observable2D.cpp.

: Observable(name_i, thname_i, label_i, tMCMC_i, min_i, max_i, tho_i), 
  bin_min(2, 0.),
  bin_max(2, 0.)
{
    thname2 = thname2_i;
    label2 = label2_i;
    min2 = min2_i;
    max2 = max2_i;
    tho2 = tho2_i;
    obsType2 = "";
    filepath = "";
    iterationNo2 = std::numeric_limits<int>::max();
    IsEOF = false;
}

◆ Observable2D() [2/4]

Observable2D::Observable2D ( )

The default constructor.

Definition at line 58 of file Observable2D.cpp.

: Observable(),
  bin_min(2, 0.),
  bin_max(2, 0.)
{
    thname2 = "";
    label2 = "";
    min2 = 0.;
    max2 = 0.;
    tho2 = NULL;
    obsType2 = "";
    filepath = "";
    iterationNo2 = std::numeric_limits<int>::max();
    IsEOF = false;
}

◆ Observable2D() [3/4]

Observable2D::Observable2D ( const Observable & o1d )

A conversion constructor. Constructs Observable2D with just one observable.

Definition at line 42 of file Observable2D.cpp.

: Observable(o1d),
  bin_min(2, 0.),
  bin_max(2, 0.)
{
    thname2 = "";
    label2 = "";
    min2 = 0.;
    max2 = 0.;
    tho2 = NULL;
    obsType2 = "";
    filepath = "";
    iterationNo2 = std::numeric_limits<int>::max();
    IsEOF = false;
}

◆ Observable2D() [4/4]

Observable2D::Observable2D ( const Observable2D & orig )

The copy constructor.

Definition at line 74 of file Observable2D.cpp.

: Observable(orig.name, orig.thname, orig.label, orig.tMCMC, orig.min, orig.max, orig.tho),
  bin_min(orig.bin_min),
  bin_max(orig.bin_max)
{
 
    distr = orig.distr;
    filename = orig.filename;
    histoname = orig.histoname;
    ave = orig.ave;
    errg = orig.errg;
    errg2 = orig.errg2;
    errf = orig.errf;
    errf2 = orig.errf2;
    errgl = orig.errgl;
    errgl2 = orig.errgl2;
    errgr = orig.errgr;
    errgr2 = orig.errgr2;
 
    thname2 = orig.thname2;
    label2 = orig.label2;
    min2 = orig.min2;
    max2 = orig.max2;
    tho2 = orig.tho2;
    filepath = orig.filepath;
    iterationNo2 = orig.iterationNo2;
    IsEOF = orig.IsEOF;
    inhisto2d = orig.inhisto2d;
}

◆ ~Observable2D()

Observable2D::~Observable2D ( )

virtual

The default destructor.

Definition at line 104 of file Observable2D.cpp.

105{}

Member Function Documentation

◆ computeTheoryValue2()

double Observable2D::computeTheoryValue2 ( )

A method to access the computed theory value of the second observable.

Definition at line 107 of file Observable2D.cpp.

{
   if (tho2->getModel().getIterationNo() == iterationNo2) {
        return thValue2;
    } else {
       iterationNo2 = tho2->getModel().getIterationNo();
        thValue2 = tho2->computeThValue();
        return thValue2;
    }
}

◆ computeWeight() [1/2]

virtual double Observable2D::computeWeight ( )

inlinevirtual

A method to compute the weight associated with the observable.

Reimplemented from Observable.

Definition at line 212 of file Observable2D.h.

    {
        return computeWeight(computeTheoryValue(),computeTheoryValue2());
    }

◆ computeWeight() [2/2]

double Observable2D::computeWeight	(	double	th1,
		double	th2
	)

virtual

A method to compute the weight associated with the observable.

Parameters

[in]	th1	the theoretical value of the first observable
[in]	th2	the theoretical value of the second observable

Reimplemented from Observable.

Definition at line 139 of file Observable2D.cpp.

{
    double logprob;
    if (distr.compare("file") == 0) {
        int i = inhisto2d->FindBin(th1, th2);
        if (inhisto2d->IsBinOverflow(i) || inhisto2d->IsBinUnderflow(i))
            logprob = log(0.);
        else
            logprob = log(inhisto2d->GetBinContent(i));
        //logprob = log(h->GetBinContent(h->GetXaxis()->FindBin(th1),h->GetYaxis()->FindBin(th2)));
    } else if (distr.compare("weight") == 0) {
        double wt2;
        if (obsType2.compare("AsyGausObservable") == 0) wt2 = Observable::LogSplitGaussian(th2, ave2, errgl, errgr);
        else wt2 = Observable::computeWeight(th2, ave2, errg2, errf2);
        logprob = Observable::computeWeight(th1) + wt2;
    } else
        throw std::runtime_error("ERROR: 2D MonteCarloEngine::Weight() called without file for "
            + name);
    return (logprob);
}

◆ getAve2()

double Observable2D::getAve2 ( ) const

inline

A get method to access the average value of the second observable.

Returns: the average value of the observable

Definition at line 231 of file Observable2D.h.

    {
        return ave2;
    }

◆ getErrf2()

double Observable2D::getErrf2 ( ) const

inline

A get method to access the flat error of the second observable.

Returns: the flat error of the observable

Definition at line 249 of file Observable2D.h.

    {
        return errf2;
    }

◆ getErrg2()

double Observable2D::getErrg2 ( ) const

inline

A get method to access the Gaussian error of the second observble.

Returns: the Gauissian error of the observable

Definition at line 267 of file Observable2D.h.

    {
        return errg2;
    }

◆ getLabel2()

std::string Observable2D::getLabel2 ( ) const

inline

A get method to access the label for the second observable.

Returns: the label for the second observable

Definition at line 101 of file Observable2D.h.

    {
        return label2;
    }

◆ getMax2()

double Observable2D::getMax2 ( ) const

inline

A get method to access the maximum value of the second observable.

Returns: the maximum value of the second observable

Definition at line 119 of file Observable2D.h.

    {
        return max2;
    }

◆ getMin2()

double Observable2D::getMin2 ( ) const

inline

A get method to access the minimum value of the second observable.

Returns: the minimum value of the second observable

Definition at line 137 of file Observable2D.h.

    {
        return min2;
    }

◆ getObsType2()

std::string Observable2D::getObsType2 ( ) const

inline

Definition at line 222 of file Observable2D.h.

    {
        return obsType2;
    }

◆ getThname2()

std::string Observable2D::getThname2 ( ) const

inline

A get method to access the thname of the second observable as defined in ThFactory class.

Returns: thname the name of the second observable as listed in ThFactory() class

Definition at line 155 of file Observable2D.h.

    {
        return thname2;
    }

◆ getTho2()

ThObservable * Observable2D::getTho2 ( ) const

inline

A get method to access the pointer to the object of the ThObservable class for the second observable.

Returns: pointer to the object of type ThObservable() for the second observable

Definition at line 174 of file Observable2D.h.

    {
        return tho2;
    }

◆ isEOF()

bool Observable2D::isEOF ( )

inline

A method to check if the end of file has been reached.

Returns: a boolean which is true if the end of file has been reached

Definition at line 285 of file Observable2D.h.

    {
        return IsEOF;
    }

◆ ParseObservable2D()

int Observable2D::ParseObservable2D	(	std::string &	type,
		boost::tokenizer< boost::char_separator< char > > *	tok,
		boost::tokenizer< boost::char_separator< char > >::iterator &	beg,
		std::string &	infilename,
		std::ifstream &	ifile,
		int	lineNo,
		int	rank
	)

Definition at line 160 of file Observable2D.cpp.

{
    if (infilename.find("\\/") == std::string::npos) filepath = infilename.substr(0, infilename.find_last_of("\\/") + 1);
    if (std::distance(tok->begin(), tok->end()) < 12) {
        setName(*beg);
        ++beg;
        if (std::distance(tok->begin(), tok->end()) < 4) {
            if(rank == 0) throw std::runtime_error("ERROR: lack of information on " + name + " in " + infilename + " at line number" + boost::lexical_cast<std::string>(lineNo));
            else sleep (2);
        }
        std::string toMCMC = *beg;
        if (toMCMC.compare("MCMC") == 0)
            setTMCMC(true);
        else if (toMCMC.compare("noMCMC") == 0)
            setTMCMC(false);
        else {
            if (rank == 0) throw std::runtime_error("ERROR: wrong MCMC flag in Observable2D" + name + " at line number:" + boost::lexical_cast<std::string>(lineNo) + " of file " + infilename + ".\n");
            else sleep(2);
        }
        
        ++beg;
        setDistr(*beg);
        if (distr.compare("file") == 0) {
            if (std::distance(tok->begin(), tok->end()) < 6) {
                if(rank == 0) throw std::runtime_error("ERROR: lack of information on "+ *beg + " in " + infilename);
                else sleep (2);
            }
            setFilename(filepath + *(++beg));
            setHistoname(*(++beg));
        }
 
        std::vector<double> min(2, 0.);
        std::vector<double> max(2, 0.);
        std::vector<double> ave(2, 0.);
        std::vector<double> errg(2, 0.);
        std::vector<double> errf(2, 0.);
        std::vector<double> errgl(2, 0.);
        std::vector<double> errgr(2, 0.);
        std::vector<std::string> thname(2, "");
        std::vector<std::string> label(2, "");
        std::vector<std::string> type2D(2, "");
        std::string line;
        size_t pos = 0;
        boost::char_separator<char> sep(" \t");
        for (int i = 0; i < 2; i++) {
            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 Observable2D please! In file " + infilename + " 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();
            type2D[i] = *beg;
            if (type2D[i].compare("Observable") != 0 && type2D[i].compare("BinnedObservable") != 0 && type2D[i].compare("FunctionObservable") != 0 && type2D[i].compare("AsyGausObservable") != 0) {
                if (rank == 0) throw std::runtime_error("ERROR: in line no." + boost::lexical_cast<std::string>(lineNo) + " of file " + infilename + ", expecting an Observable or BinnedObservable or FunctionObservable type here...\n");
                else sleep(2);
            }
            ++beg;
            thname[i] = *beg;
            ++beg;
            label[i] = *beg;
            while ((pos = label[i].find("~", pos)) != std::string::npos)
                label[i].replace(pos++, 1, " ");
            ++beg;
            min[i] = atof((*beg).c_str());
            ++beg;
            max[i] = atof((*beg).c_str());
            if (distr.compare("weight") == 0) {
                ++beg;
                ave[i] = atof((*beg).c_str());
                ++beg;
                if (type.compare("AsyGausObservable") == 0) errgl[i] = atof((*beg).c_str());
                else errg[i] = atof((*beg).c_str());
                ++beg;
                if (type.compare("AsyGausObservable") == 0) errgr[i] = atof((*beg).c_str());
                else errf[i] = atof((*beg).c_str());
                if (type.compare("AsyGausObservable") != 0 && errg[i] == 0. && errf[i] == 0.) {
                    if (rank == 0) throw std::runtime_error("ERROR: The Gaussian and flat error in weight for " + name + " cannot both be 0. in the " + infilename + " file, line number:" + boost::lexical_cast<std::string>(lineNo) + ".\n");
                    else sleep(2);
                }
                if (type.compare("AsyGausObservable") == 0 && errgl[i] == 0. && errgr[i] == 0.) {
                    if (rank == 0) throw std::runtime_error("ERROR: The left Gaussian and right Gaussian error in weight for " + name + " cannot both be 0. in the " + infilename + " file, line number:" + boost::lexical_cast<std::string>(lineNo) + ".\n");
                    else sleep(2);
                }
            } else if (distr.compare("noweight") == 0 || distr.compare("file") == 0) {
                if (type2D[i].compare("BinnedObservable") == 0 || type2D[i].compare("FunctionObservable") == 0) {
                    ++beg;
                    ++beg;
                    ++beg;
                }
            } else {
                if (rank == 0) throw std::runtime_error("ERROR: wrong distribution flag in " + name + " in file " + infilename + ".\n");
                else sleep(2);
            }
            if (type2D[i].compare("BinnedObservable") == 0) {
                ++beg;
                bin_min[i] = atof((*beg).c_str());
                ++beg;
                bin_max[i] = atof((*beg).c_str());
            } else if (type2D[i].compare("FunctionObservable") == 0) {
                ++beg;
                bin_min[i] = atof((*beg).c_str());
                ++beg;
            }
        }
        setObsType(type2D[0]);
        obsType2 = type2D[1];
        setThname(thname[0]);
        thname2 = thname[1];
        setLabel(label[0]);
        label2 = label[1];
        setMin(min[0]);
        min2 = min[1];
        setMax(max[0]);
        max2= max[1];
        setAve(ave[0]);
        ave2 = ave[1];
        setErrg(errg[0]);
        errg2 = errg[1];
        setErrf(errf[0]);
        errf2 = errf[1];
        setErrgl(errgl[0]);
        errgl2 = errgl[1];
        setErrgr(errgr[0]);
        errgr2 = errgr[1];
        if (distr.compare("file") == 0) {
            setLikelihoodFromHisto(filename, histoname);
            if (rank == 0) std::cout << "added input histogram " << filename << "/" << histoname << std::endl;
        }
        return lineNo;
    } else {
        beg = ParseObservable(type, tok, beg, filepath, filename, rank);
        ++beg;
        std::string distr = *beg;
        if (distr.compare("file") == 0) {
            if (std::distance(tok->begin(), tok->end()) < 14) {
                if (rank == 0) throw std::runtime_error("ERROR: lack of information on " + *beg + " in " + infilename + ".\n");
                else sleep(2);
            }
            setFilename(filepath + *(++beg));
            setHistoname(*(++beg));
            setLikelihoodFromHisto(filename, histoname);
            if (rank == 0) std::cout << "added input histogram " << filename << "/" << histoname << std::endl;
        } else if (distr.compare("noweight") == 0) {
        } else { 
            if (rank == 0) throw std::runtime_error("ERROR: wrong distribution flag in " + name);
            else sleep(2);
        }
        setDistr(distr);
        ++beg;
        thname2 = *beg;
        ++beg;
        std::string label = *beg;
        size_t pos = 0;
        while ((pos = label.find("~", pos)) != std::string::npos)
            label.replace(pos, 1, " ");
        label2 = label;
        ++beg;
        min2 = atof((*beg).c_str());
        ++beg;
        max2 = atof((*beg).c_str());
        ++beg;
        if (beg != tok->end())
            if (rank == 0) std::cout << "WARNING: unread information in observable2D " << name << std::endl;
        return lineNo;
    }
}

◆ setAve2()

void Observable2D::setAve2 ( double ave2 )

inline

A set method to fix the average value of the second observable.

Parameters

[in] ave2 the average value of the observable

Definition at line 240 of file Observable2D.h.

    {
        this->ave2 = ave2;
    }

◆ setErrf2()

void Observable2D::setErrf2 ( double errf2 )

inline

A set method to fix the flat error of the second observable.

Parameters

[in] errf2 the flat error of the observable

Definition at line 258 of file Observable2D.h.

    {
        this->errf2 = errf2;
    }

◆ setErrg2()

void Observable2D::setErrg2 ( double errg2 )

inline

A set method to fix the gaussian error of the observable.

Parameters

[in] errg2 the Gaussian error of the observable

Definition at line 276 of file Observable2D.h.

    {
        this->errg2 = errg2;
    }

◆ setLabel2()

void Observable2D::setLabel2 ( std::string label2 )

inline

A set method to fix the label for the second observable.

Parameters

[in] label2 the label for the second observable

Definition at line 110 of file Observable2D.h.

    {
        this->label2 = label2;
    }

◆ setLikelihoodFromHisto()

void Observable2D::setLikelihoodFromHisto	(	std::string	filename,
		std::string	histoname
	)

virtual

A set method to set the likelihood from which the experimental likelihood of the observable will be read.

Parameters

filename	the name of the file
histoname	the name of the histogram

Reimplemented from Observable.

Definition at line 118 of file Observable2D.cpp.

{
    this->filename = filename;
    this->histoname = histoname;
    TFile *lik2 = new TFile((filename + ".root").c_str(), "read");
    TH2D *htmp2 = (TH2D*) (lik2->Get(histoname.c_str()));
    if (htmp2 == NULL)
        throw std::runtime_error("ERROR: nonexistent histogram called "
            + histoname
            + " in " + filename + ".root");
    inhisto2d = (TH2D *) htmp2->Clone((filename + "/" + histoname).c_str());
    inhisto2d->SetDirectory(gROOT);
    std::cout << "added 2D input histogram " << name << std::endl;
    setMin(inhisto2d->GetXaxis()->GetXmin());
    setMax(inhisto2d->GetXaxis()->GetXmax());
    setMin2(inhisto2d->GetYaxis()->GetXmin());
    setMax2(inhisto2d->GetYaxis()->GetXmax());
    lik2->Close();
    delete lik2;
}

◆ setMax2()

void Observable2D::setMax2 ( double max2 )

inline

A set method to fix the maximum value for the second observable.

Parameters

[in] max2 the maximum value for the second observable

Definition at line 128 of file Observable2D.h.

    {
        this->max2 = max2;
    }

◆ setMin2()

void Observable2D::setMin2 ( double min2 )

inline

A set method to fix the minimum value for the second observable.

Parameters

[in] min2 the minimum value for the second observable

Definition at line 146 of file Observable2D.h.

    {
        this->min2 = min2;
    }

◆ setObsType2()

void Observable2D::setObsType2 ( std::string & obsType2_i )

inline

Definition at line 217 of file Observable2D.h.

    {
        obsType2 = obsType2_i;
    }

◆ setThname2()

void Observable2D::setThname2 ( std::string thname2 )

inline

A set method to fix the name of the second observable as listed in ThFactory class.

Parameters

[in] thname2 the name of the second observable as listed in ThFactory class

Definition at line 164 of file Observable2D.h.

    {
        this->thname2 = thname2;
    }

◆ setTho1Tho2()

void Observable2D::setTho1Tho2	(	ThObservable *	tho1_i,
		ThObservable *	tho2_i
	)

inline

A set method to fix the pointer to object of type ThObservable class for the second observable.

Parameters

[in] tho2 pointer to the object of type ThObservable() for the second observable

Definition at line 196 of file Observable2D.h.

    {
        setTho(tho1_i, bin_min[0], bin_max[0]);
        setTho2(tho2_i);
    }

◆ setTho2()

void Observable2D::setTho2 ( ThObservable * tho2_i )

inline

A set method to fix the pointer to object of type ThObservable class for the second observable.

Parameters

[in] tho2 pointer to the object of type ThObservable() for the second observable

Definition at line 184 of file Observable2D.h.

    {
        tho2 = tho2_i;
        tho2->setBinMin(bin_min[1]);
        tho2->setBinMax(bin_max[1]);
    }

Member Data Documentation

◆ ave2

double Observable2D::ave2

private

The average value of the second observable.

Definition at line 295 of file Observable2D.h.

◆ bin_max

std::vector<double> Observable2D::bin_max

private

The maximum value of the bin.

Definition at line 304 of file Observable2D.h.

◆ bin_min

std::vector<double> Observable2D::bin_min

private

The minimum value of the bin.

Definition at line 303 of file Observable2D.h.

◆ errf2

double Observable2D::errf2

private

the flat error of the second observable.

Definition at line 297 of file Observable2D.h.

◆ errg2

double Observable2D::errg2

private

The Gaussian error of the second observable.

Definition at line 296 of file Observable2D.h.

◆ errgl2

double Observable2D::errgl2

private

The left Gaussian error of the second observable.

Definition at line 298 of file Observable2D.h.

◆ errgr2

double Observable2D::errgr2

private

The right Gaussian error of the second observable.

Definition at line 299 of file Observable2D.h.

◆ filepath

std::string Observable2D::filepath

private

The path to the file being parsed.

Definition at line 305 of file Observable2D.h.

◆ inhisto2d

TH2D* Observable2D::inhisto2d

private

2D Histogram containing the experimental likelihood for the observable.

Definition at line 302 of file Observable2D.h.

◆ IsEOF

bool Observable2D::IsEOF

private

A bolean that is true if the end of file is reached.

Definition at line 308 of file Observable2D.h.

◆ iterationNo2

int Observable2D::iterationNo2

private

Counts the iteration to help with caching.

Definition at line 306 of file Observable2D.h.

◆ label2

std::string Observable2D::label2

private

A label for the second observable.

Definition at line 292 of file Observable2D.h.

◆ max2

double Observable2D::max2

private

The maximum valus of the second observable.

Definition at line 294 of file Observable2D.h.

◆ min2

double Observable2D::min2

private

The minimum value of the second observable.

Definition at line 293 of file Observable2D.h.

◆ obsType2

std::string Observable2D::obsType2

private

Type of the second Observable. 0: Observable, 1: HiggsObservable, 2: BinnedObservable, 3: FunctionObservable.

Definition at line 300 of file Observable2D.h.

◆ thname2

std::string Observable2D::thname2

private

The name for the second observable as fixed in the ThObservable() class.

Definition at line 291 of file Observable2D.h.

◆ tho2

ThObservable* Observable2D::tho2

private

A pointer to an object of the ThObservable class.

Definition at line 301 of file Observable2D.h.

◆ thValue2

double Observable2D::thValue2

private

The theory value of the second observable.

Definition at line 307 of file Observable2D.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

◆ Observable2D() [1/4]

◆ Observable2D() [2/4]

◆ Observable2D() [3/4]

◆ Observable2D() [4/4]

◆ ~Observable2D()

Member Function Documentation

◆ computeTheoryValue2()

◆ computeWeight() [1/2]

◆ computeWeight() [2/2]

◆ getAve2()

◆ getErrf2()

◆ getErrg2()

◆ getLabel2()

◆ getMax2()

◆ getMin2()

◆ getObsType2()

◆ getThname2()

◆ getTho2()

◆ isEOF()

◆ ParseObservable2D()

◆ setAve2()

◆ setErrf2()

◆ setErrg2()

◆ setLabel2()

◆ setLikelihoodFromHisto()

◆ setMax2()

◆ setMin2()

◆ setObsType2()

◆ setThname2()

◆ setTho1Tho2()

◆ setTho2()

Member Data Documentation

◆ ave2

◆ bin_max

◆ bin_min

◆ errf2

◆ errg2

◆ errgl2

◆ errgr2

◆ filepath

◆ inhisto2d

◆ IsEOF

◆ iterationNo2

◆ label2

◆ max2

◆ min2

◆ obsType2

◆ thname2

◆ tho2

◆ thValue2