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);