Integration error for RooMultiVarGaussian

Hello,

I am trying to apply a constraint using Hfitter and a RooMultiVarGaussian defined below. I am not exactly sure how Hfitter is handling the multivariate Gaussian, but when calling the Hfitter method that is using the Gaussian, I receive the error:

[#0] ERROR:Integration -- RooNumIntFactory::createIntegrator: No integration method has been defined for an open ended 2-dimensional integral
[#0] ERROR:Integration -- RooRealIntegral::signal_constraint_Int[mean_gaussian,mean_gaussian_obs,sigma_gaussian,sigma_gaussian_obs]: failed to create valid integrator.
[#0] ERROR:Integration -- RooRealIntegral::signal_constraint_Int[mean_gaussian,mean_gaussian_obs,sigma_gaussian,sigma_gaussian_obs]:evaluate: cannot initialize numerical integrator

Does anyone have suggestions for how to debug this error?

Best regards

RooWorkspace* w = new RooWorkspace("w");
w->factory("mean_gaussian_obs[41.6495]");
w->factory("sigma_gaussian_obs[8.43191]");
w->factory("mean_gaussian[41.6495, 0, 100]");
w->factory("sigma_gaussian[8.43191, 0, 100]");
w->defineSet("signal_shape_parameters", "sigma_gaussian,mean_gaussian");
w->defineSet("signal_shape_parameters_obs", "sigma_gaussian_obs,mean_gaussian_obs");
RooMultiVarGaussian* signal_constraint = new RooMultiVarGaussian("signal_constraint", "", *w->set("signal_shape_parameters"), *w->set("signal_shape_parameters_obs"), cov);

I fixed your example by adding the missing covariance matrix [ 1 0, 0 1] in my case. Now, nothing happens. What has to be done to trigger the problem?

RooWorkspace* w = new RooWorkspace("w");
w->factory("mean_gaussian_obs[41.6495]");
w->factory("sigma_gaussian_obs[8.43191]");
w->factory("mean_gaussian[41.6495, 0, 100]");
w->factory("sigma_gaussian[8.43191, 0, 100]");
w->defineSet("signal_shape_parameters", "sigma_gaussian,mean_gaussian");
w->defineSet("signal_shape_parameters_obs", "sigma_gaussian_obs,mean_gaussian_obs");

TMatrixDSym matr(2);
matr = 1.;
matr[0][1] = 0.;
matr[1][0] = 0.;
matr.Print();

RooMultiVarGaussian* signal_constraint = new RooMultiVarGaussian("signal_constraint", "", *w->set("signal_shape_parameters"), *w->set("signal_shape_parameters_obs"), matr);

This topic was automatically closed 14 days after the last reply. New replies are no longer allowed.