A class for Higgs experimental analyses. More...

#include <HiggsObservable.h>

Inheritance diagram for HiggsObservable:

[legend]

Collaboration diagram for HiggsObservable:

[legend]

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 26 of file HiggsObservable.h.

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

int	getNChannels ()

int	getNTheoryValues ()

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

	HiggsObservable (const Observable &Obs)

	HiggsObservable (const HiggsObservable &orig)

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	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...

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

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

	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...

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

	Observable ()
	The default 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	setFilename (std::string filename_i)

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

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

std::vector< double >	theoryValues

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

double	thobsvsize

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...

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

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 oservable 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...

Constructor & Destructor Documentation

HiggsObservable::HiggsObservable ( const Observable & Obs )

inline

Definition at line 28 of file HiggsObservable.h.

                                           : 
                Observable(Obs)
                {
                    isnew = true;
                };

HiggsObservable::HiggsObservable ( const HiggsObservable & orig )

Member Function Documentation

double HiggsObservable::computeWeight ( )

virtual

A method to compute the weight associated with the observable.

Reimplemented from Observable.

Definition at line 77 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();
                 theoryValues.at(thobsvsize + i + 1) = mu * theoryValues.at(thobsvsize);
             logprob += LogSplitGaussian(theoryValues.at(thobsvsize + i), channels(i, thobsvsize), channels(i, thobsvsize + 1), channels(i, thobsvsize + 2));
         }
       
     } 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), channels(i, 3), channels(i, 4), channels(i, 5));
         }
     }
     
     return (logprob);
 }

int HiggsObservable::getNChannels ( )

inlinevirtual

Reimplemented from Observable.

Definition at line 81 of file HiggsObservable.h.

     {
         return channelsize;
     }

int HiggsObservable::getNTheoryValues ( )

inlinevirtual

Reimplemented from Observable.

Definition at line 76 of file HiggsObservable.h.

     {
         return thobsvsize;
     }

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

inlinevirtual

Reimplemented from Observable.

Definition at line 71 of file HiggsObservable.h.

     {
         theoryValues_i = theoryValues;
     }

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 66 of file HiggsObservable.h.

     {
         return isnew;
     }

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 112 of file HiggsObservable.cpp.

 {
     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;
 }

void HiggsObservable::setIsnew ( bool isnew )

inline

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

Definition at line 57 of file HiggsObservable.h.

     {
         this->isnew = isnew;
     }

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 22 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;
     int i = 0, 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);
     do {
         IsEOF = getline(ifile, line).eof();
         if (*line.rbegin() == '\r') line.erase(line.length() - 1); // for CR+LF
         if (line.compare(0, 11, "MEASUREMENT") != 0)
             continue;
         boost::char_separator<char> sep(" |");
         boost::tokenizer<boost::char_separator<char> > tok(line, sep);
         boost::tokenizer<boost::char_separator<char> >::iterator beg = tok.begin();
         // 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());
         i++;
 
     } 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

TMatrixD HiggsObservable::channels

private

A matrix holding the information of all the channels.

Definition at line 101 of file HiggsObservable.h.

double HiggsObservable::channelsize

private

Definition at line 106 of file HiggsObservable.h.

bool HiggsObservable::isnew

private

A boolean which is true if the parametrization is new.

Definition at line 102 of file HiggsObservable.h.

std::vector<double> HiggsObservable::theoryValues

private

Definition at line 104 of file HiggsObservable.h.

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

private

A vector of ThObservables.

Definition at line 103 of file HiggsObservable.h.

double HiggsObservable::thobsvsize

private

Definition at line 105 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

Member Function Documentation

Member Data Documentation