GetConfidenceIntervals with TMinuit?

From: Martin Simard <msimard_at_slac.stanford.edu>
Date: Wed, 27 Jun 2012 12:07:01 -0700


Hi,

I want to use the GetConfidenceIntervals function using TMinuit. Below is a short example of what I'm trying to do. I know where is my error, but i don't know how to get it work. Is it possible to obtain the confidence interval with TMinuit? Do I have to compute it manually?

Thanks,
Martin


TMinuit* minuit = new TMinuit;
minuit->SetFCN(&Chi2);
...
minuit->Migrad();

TGraph gConf = new TGraph(50); // the graph which holds the contours TVirtualFitter* tvf = minuit->GetFitter(); //This doesn't compile... tvf->GetConfidenceIntervals(gConf,0.95); // 95% confidence interval

//Instead of the 3 lines above I also try the following but I get a
//Break - Segmentation when running.

TVirtualFitter* tvf = TVirtualFitter::GetFitter(); TGraph gConf = new TGraph(50);
(TVirtualFitter::GetFitter())->GetConfidenceIntervals(gConf,0.95); Received on Wed Jun 27 2012 - 21:07:10 CEST

This archive was generated by hypermail 2.2.0 : Thu Jun 28 2012 - 05:50:02 CEST