Dear Experts,
when I run these instructions:
RooAbsReal* NLL;
NLL = (*TotalPDF)->createNLL(*dataSet,Extended(true),ExternalConstraints(*vecConstr));
RooMinuit RooMin(NLL);
RooPlot myFrameNLL = RooMin.contour(*GetVar(*TotalPDF,“AfbS”),*GetVar(*TotalPDF,“FlS”),1.0,2.0,3.0);
in particular the last one: RooPlot* myFrameNLL = RooMin.contour(*GetVar(*TotalPDF,“AfbS”),*GetVar(*TotalPDF,“FlS”),1.0,2.0,3.0);
the program crash badly with the following error message:
START MNCONTOUR CALCULATION OF 50 POINTS ON CONTOUR.
EACH POINT IS A MINIMUM WITH RESPECT TO THE OTHER 10 VARIABLE PARAMETERS.
POSITIVE MINOS ERROR NOT CALCULATED FOR PARAMETER 1
NEGATIVE MINOS ERROR NOT CALCULATED FOR PARAMETER 1
MNCONT UNABLE TO FIND FOUR POINTS.
Warning in TMinuit::Contour: Cannot find more than 4 points, no TGraph returned
** 10 **SET ERRDEF 2
START MNCONTOUR CALCULATION OF 50 POINTS ON CONTOUR.
EACH POINT IS A MINIMUM WITH RESPECT TO THE OTHER 10 VARIABLE PARAMETERS.
POSITIVE MINOS ERROR NOT CALCULATED FOR PARAMETER 1
NEGATIVE MINOS ERROR NOT CALCULATED FOR PARAMETER 1
MNCONT UNABLE TO FIND FOUR POINTS.
Warning in TMinuit::Contour: Cannot find more than 4 points, no TGraph returned
** 11 **SET ERRDEF 4.5
START MNCONTOUR CALCULATION OF 50 POINTS ON CONTOUR.
EACH POINT IS A MINIMUM WITH RESPECT TO THE OTHER 10 VARIABLE PARAMETERS.
MINUIT WARNING IN MIGRAD
============== STARTING MATRIX NOT POS-DEFINITE.
POSITIVE MINOS ERROR NOT CALCULATED FOR PARAMETER 1
MIGRAD FAILS WITH STRATEGY=0. WILL TRY WITH STRATEGY=1.
NEGATIVE MINOS ERROR NOT CALCULATED FOR PARAMETER 1
MNCONT UNABLE TO FIND FOUR POINTS.
Warning in TMinuit::Contour: Cannot find more than 4 points, no TGraph returned
** 12 **SET ERRDEF 0.5
*** Break *** segmentation violation
Generating stack trace…
0x0000000108c92769 in IterativeMass2AnglesFitq2Bins(RooDataSet*, bool, double, double, double, RooRealVar*, RooRealVar*, RooRealVar*, int, unsigned int, std::__1::vector<TH1D*, std::__1::allocator<TH1D*> >*, std::__1::vect
0x0000000108c9cf12 in main (in ExtractYield) + 22770
0x00007fff91fe45fd in start (in libdyld.dylib) + 1
0x0000000000000005 in
ExtractYield(18319,0x7fff74ece310) malloc: *** error for object 0x7f961754d990: pointer being freed was not allocated
*** set a breakpoint in malloc_error_break to debug
Abort trap: 6
what can I do to avoid the crash ?
Many thanks,
- Mauro.