#include <Neg2LLMultiConstraint.h>

Inheritance diagram for MINT::Neg2LLMultiConstraint:

Public Member Functions
	Neg2LLMultiConstraint (MinuitParameterSet *mps=0, string label="")

virtual void	beginFit ()

virtual void	parametersChanged ()

virtual void	endFit ()

virtual double	getVal ()

virtual double	getNewVal ()

int	getNumberParams ()

void	smearInputValues ()

double	smearInputValuesChol (int index=0, int seed=0, int offset=1)

RooDataSet *	generateToys (int N=100)

TMatrixDSym *	getCovMatrix ()

vector< string >	split (string s, string delimiter)

vector< double >	splitValues (string s, string delimiter)

virtual	~Neg2LLMultiConstraint ()

Public Member Functions inherited from MINT::Minimisable
	Minimisable (MinuitParameterSet *mps=0)

	Minimisable (const Minimisable &other)

virtual	~Minimisable ()

void	setPset (MinuitParameterSet *mps)

MinuitParameterSet *	getParSet ()

Public Member Functions inherited from MINT::IMinimisable
virtual void	Gradient (std::vector< double > &grad)

virtual bool	useAnalyticGradient ()

virtual void	setUseAnalyticGradient (bool useAnalyticGradient)

virtual	~IMinimisable ()

Protected Attributes
RooArgList *	_x

RooArgList *	_mu

RooMultiVarGaussian *	_gauss_cov

TMatrixTSym< double > *	_cov

MinuitParameterSet *	_mps

TMatrixD *	_UT

Detailed Description

Definition at line 23 of file Neg2LLMultiConstraint.h.

Constructor & Destructor Documentation

◆ Neg2LLMultiConstraint()

MINT::Neg2LLMultiConstraint::Neg2LLMultiConstraint	(	MinuitParameterSet *	mps = `0`,
		string	label = `""`
	)

inline

Parameters to constrain

Cov. matrix

Sanity checks

Definition at line 32 of file Neg2LLMultiConstraint.h.

         : Minimisable(mps), _mps(mps) {
 
             NamedParameter<string> constrain(("ConstrainMulti"+label).c_str(), (string)"" , (char*) 0);
             string constrain_s = constrain;         
             vector<string> constrain_v = split(constrain_s, " ");
 
             NamedParameter<string> corr(("ConstrainMulti"+label+"_corr").c_str(),(string)"", (char*) 0);
             string corr_s = corr;           
             vector<double> corr_v = splitValues(corr_s, " ");
 
             TMatrixDSym cov(constrain_v.size());
             int index = 0;
             for(int i=0; i < cov.GetNcols(); i++){
                 for(int j=i; j < cov.GetNcols(); j++){    
                     cov(i,j) = corr_v[index];
                     cov(j,i) = corr_v[index];
                     index++;
                 }
             }
 
             //double* corr_a = &corr_v[0];
             //TMatrixDSym cov(constrain_v.size(),corr_a);
             cov.Print();
 
             if(constrain_v.size() * (constrain_v.size()+1)/2 != corr_v.size()) { 
                 cout << "ERROR in Neg2LLMultiConstraint::Inconsistent number of parameters and cov. matrix" << endl;
                 throw "ERROR";          
             }
 
             _x = new RooArgList();
             _mu = new RooArgList();
             vector<double> sigma_v;
 
             for(int j= 0; j < constrain_v.size();j++){
 
                 cout << "Adding gauss constraint for parameter: " << constrain_v[j] << endl;
 
                 double mean = -99999;
                 double sigma = -99999;
 
                 for(unsigned int i=0; i < _mps->size(); i++){
                         if(0 == _mps->getParPtr(i)) continue;
                     if(constrain_v[j] == _mps->getParPtr(i)->name()){   
                         mean = _mps->getParPtr(i)->mean();
                         sigma = ((FitParameter*)_mps->getParPtr(i))->stepInit();
                     }
                 }
                 if(mean == -99999) {
                     cout << "ERROR in Neg2LLMultiConstraint::Parameter not found" << endl;
                     throw "ERROR";          
                 }
                 else {
                     RooRealVar* x = new RooRealVar((constrain_v[j]).c_str(), (constrain_v[j]).c_str(),mean);
                     _x->add(*x);
                     _mu->add(RooRealConstant::value(mean));
                     sigma_v.push_back(sigma);
                 }
             }
                 
             for(int i=0; i < cov.GetNcols(); i++){
                 for(int j=0; j < cov.GetNcols(); j++){    
                     cov(i,j) = cov(i,j) * sigma_v[i] * sigma_v[j];
                 }
             }
             _cov = new TMatrixDSym(cov);
 
             _x->Print();
             _mu->Print();
             _cov->Print();
             _gauss_cov = new RooMultiVarGaussian("gauss_cov","gauss_cov",*_x, *_mu, *_cov);
 
             TDecompChol tdc(*_cov);
             tdc.Decompose();
             TMatrixD U = tdc.GetU();
             _UT = new TMatrixD(TMatrixD::kTransposed,U);
             _UT->Print();
 
     };

◆ ~Neg2LLMultiConstraint()

virtual MINT::Neg2LLMultiConstraint::~Neg2LLMultiConstraint ( )

inlinevirtual

Definition at line 207 of file Neg2LLMultiConstraint.h.

207 {}

Member Function Documentation

◆ beginFit()

virtual void MINT::Neg2LLMultiConstraint::beginFit ( )

inlinevirtual

Reimplemented from MINT::Minimisable.

Definition at line 115 of file Neg2LLMultiConstraint.h.

115 {};

◆ endFit()

virtual void MINT::Neg2LLMultiConstraint::endFit ( )

inlinevirtual

Reimplemented from MINT::Minimisable.

Definition at line 121 of file Neg2LLMultiConstraint.h.

121 {};

◆ generateToys()

RooDataSet* MINT::Neg2LLMultiConstraint::generateToys ( int N = 100 )

inline

Definition at line 172 of file Neg2LLMultiConstraint.h.

                                          {
             cout << "Smearing input values " << endl;
             RooDataSet* data_cov = _gauss_cov->generate(*_x, N);
             return data_cov;
     }

◆ getCovMatrix()

TMatrixDSym* MINT::Neg2LLMultiConstraint::getCovMatrix ( )

inline

Definition at line 178 of file Neg2LLMultiConstraint.h.

                                {
         return _cov;
     }

◆ getNewVal()

virtual double MINT::Neg2LLMultiConstraint::getNewVal ( )

inlinevirtual

Reimplemented from MINT::Minimisable.

Definition at line 127 of file Neg2LLMultiConstraint.h.

                                   { 
             parametersChanged();
             return getVal();
         }

◆ getNumberParams()

int MINT::Neg2LLMultiConstraint::getNumberParams ( )

inline

Definition at line 132 of file Neg2LLMultiConstraint.h.

                          {
         return _cov->GetNcols();
     }

◆ getVal()

virtual double MINT::Neg2LLMultiConstraint::getVal ( )

inlinevirtual

Implements MINT::Minimisable.

Definition at line 123 of file Neg2LLMultiConstraint.h.

                                {
         return -2. * log(_gauss_cov->getVal(*_x));
         }

◆ parametersChanged()

virtual void MINT::Neg2LLMultiConstraint::parametersChanged ( )

inlinevirtual

Reimplemented from MINT::Minimisable.

Definition at line 116 of file Neg2LLMultiConstraint.h.

                                         {
         for(int i=0; i < _x->getSize(); i++){
             ((RooRealVar*)_x->at(i))->setVal(_mps->getParPtr(_x->at(i)->GetName())->mean());
         }
     };

◆ smearInputValues()

void MINT::Neg2LLMultiConstraint::smearInputValues ( )

inline

Definition at line 136 of file Neg2LLMultiConstraint.h.

                            {
         
             cout << "Smearing input values " << endl;
 
             RooDataSet* data_cov = _gauss_cov->generate(*_x, 1);
             RooArgList* xvec_cov= (RooArgList*)data_cov->get(0);
 
             xvec_cov->Print();
         
             for(int i=0; i < xvec_cov->getSize(); i++){
                 _mps->getParPtr(xvec_cov->at(i)->GetName())->setCurrentFitVal(((RooRealVar*)xvec_cov->at(i))->getVal());
                 ((FitParameter*)_mps->getParPtr(xvec_cov->at(i)->GetName()))->setInit(((RooRealVar*)xvec_cov->at(i))->getVal());
                 cout << "Set parameter " << xvec_cov->at(i)->GetName() << " to " << _mps->getParPtr(xvec_cov->at(i)->GetName())->mean() << endl;   
             }
     }

◆ smearInputValuesChol()

double MINT::Neg2LLMultiConstraint::smearInputValuesChol	(	int	index = `0`,
		int	seed = `0`,
		int	offset = `1`
	)

inline

Definition at line 152 of file Neg2LLMultiConstraint.h.

                                                                             {
         
             // Use the same random seed for each chol index
             TRandom3 r(seed + offset);
             cout << "Smearing input values for parameter " << index << " using random seed = " << seed + offset << endl;
 
             double val = 0;
             for(int i = 0 ; i < _UT->GetNcols(); i++){
                 if(i != index)continue;
                 val = _mps->getParPtr(_x->at(i)->GetName())->mean();
                 for(int j = 0 ; j < _UT->GetNcols(); j++){
                     val += r.Gaus(0.,1.) * (*_UT)(i,j);
                 }   
                 _mps->getParPtr(_x->at(i)->GetName())->setCurrentFitVal(val);
                 ((FitParameter*)_mps->getParPtr(_x->at(i)->GetName()))->setInit(val);
                 cout << "Set parameter " << _x->at(i)->GetName() << " to " << val << endl;
             }
             return val;
     }

◆ split()

vector<string> MINT::Neg2LLMultiConstraint::split	(	string	s,
		string	delimiter
	)

inline

Definition at line 182 of file Neg2LLMultiConstraint.h.

                                                     {
         size_t pos = 0;
         string token;
         vector<string> list;
         while ((pos = s.find(delimiter)) != string::npos) {
             token = s.substr(0, pos);
             list.push_back(token);
             s.erase(0, pos + delimiter.length());
         }
         return list;
     }

◆ splitValues()

vector<double> MINT::Neg2LLMultiConstraint::splitValues	(	string	s,
		string	delimiter
	)

inline

Definition at line 194 of file Neg2LLMultiConstraint.h.

                                                           {
         size_t pos = 0;
         string token;
         vector<double> list;
         while ((pos = s.find(delimiter)) != string::npos) {
             token = s.substr(0, pos);
             list.push_back(atof(token.c_str()));
             s.erase(0, pos + delimiter.length());
         }
         return list;
     }

Member Data Documentation

◆ _cov

TMatrixTSym<double>* MINT::Neg2LLMultiConstraint::_cov

protected

Definition at line 27 of file Neg2LLMultiConstraint.h.

◆ _gauss_cov

RooMultiVarGaussian* MINT::Neg2LLMultiConstraint::_gauss_cov

protected

Definition at line 26 of file Neg2LLMultiConstraint.h.

◆ _mps

MinuitParameterSet* MINT::Neg2LLMultiConstraint::_mps

protected

Definition at line 28 of file Neg2LLMultiConstraint.h.

◆ _mu

RooArgList * MINT::Neg2LLMultiConstraint::_mu

protected

Definition at line 25 of file Neg2LLMultiConstraint.h.

◆ _UT

TMatrixD* MINT::Neg2LLMultiConstraint::_UT

protected

Definition at line 29 of file Neg2LLMultiConstraint.h.

◆ _x

RooArgList* MINT::Neg2LLMultiConstraint::_x

protected

Definition at line 25 of file Neg2LLMultiConstraint.h.

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

Mint/Neg2LLMultiConstraint.h

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Neg2LLMultiConstraint()

◆ ~Neg2LLMultiConstraint()

Member Function Documentation

◆ beginFit()

◆ endFit()

◆ generateToys()

◆ getCovMatrix()

◆ getNewVal()

◆ getNumberParams()

◆ getVal()

◆ parametersChanged()

◆ smearInputValues()

◆ smearInputValuesChol()

◆ split()

◆ splitValues()

Member Data Documentation

◆ _cov

◆ _gauss_cov

◆ _mps

◆ _mu

◆ _UT

◆ _x