333 const char charal[29] =
" .ABCDEFGHIJKLMNOPQRSTUVWXYZ";
422 fGraphicsMode =
kTRUE;
440 fGraphicsMode =
kTRUE;
447 gROOT->GetListOfSpecials()->Add(gMinuit);
473 gROOT->GetListOfSpecials()->Add(gMinuit);
481 Error(
"TMinuit",
"can not copy construct TMinuit");
493 if (
gROOT != 0 &&
gROOT->GetListOfSpecials() != 0)
gROOT->GetListOfSpecials()->Remove(
this);
494 if (gMinuit ==
this) gMinuit = 0;
575 for (
int i = 0; i <
fMaxpar; i++) {
656 mncont(pa1,pa2,npoints,xcoor,ycoor,npfound);
659 Warning(
"Contour",
"Cannot find more than 4 points, no TGraph returned");
660 fStatus= (npfound==0 ? 1 : npfound);
665 if (npfound!=npoints) {
667 Warning(
"Contour",
"Returning a TGraph with %d points only",npfound);
672 xcoor[npoints] = xcoor[0];
673 ycoor[npoints] = ycoor[0];
676 if ((h =
gROOT->GetPluginManager()->FindHandler(
"TMinuitGraph"))) {
693 mnparm( parNo, sname, initVal, initErr, lowerLimit, upperLimit, err);
827 mnexcm(
"FIX", tmp, 1, err );
841 mnpout( parNo, name, currentValue, currentError, bnd1, bnd2, err );
880 mnexcm(
"MIGRAD", tmp, 0, err );
894 mnexcm(
"RELEASE", tmp, 1, err );
906 mnexcm(
"SET ERRDEF", &up, 1, err );
972 mnexcm(
"SET PRINT", tmp, 1, err );
974 if (printLevel <=-1)
mnexcm(
"SET NOWarnings",tmp,0,err);
995 Printf(
" FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4.");
1016 Double_t awid,ah, al, sigfig, sigrnd, alb;
1021 if (al == ah) ah = al + 1;
1024 if (naa == -1)
goto L150;
1033 if (awid <= 1) --log_;
1036 if (sigfig > 2)
goto L40;
1040 if (sigfig > 2.5)
goto L50;
1044 if (sigfig > 5)
goto L60;
1055 if (bwid <= 0)
goto L10;
1059 if (alb < 0) --lwid;
1061 alb = ah / bwid + 1;
1063 if (alb < 0) --kwid;
1066 if (naa > 5)
goto L240;
1067 if (naa == -1)
return;
1069 if (naa > 1 || nb == 1)
return;
1074 if (nb << 1 != naa)
return;
1090 Int_t ndex, i, j,
m,
n, nparx;
1096 for (i = 1; i <=
fNpar; ++i) {
1098 for (j = 1; j <=
fNpar; ++j) {
1101 ndex = m*(m-1) / 2 + n;
1106 for (i = 1; i <=
fNpar; ++i) {denom +=
fGrd[i-1]*(
fXt[i-1] - pvec[i-1]); }
1112 ycalf = (f -
fApsi) / denom;
1131 for (i = 1; i <=
fMaxext; ++i) {
1154 static TString clabel =
"0123456789ABCDEFGHIJ";
1158 Double_t fcna[115], fcnb[115], contur[20];
1159 Double_t ylabel, fmn, fmx, xlo, ylo, xup, yup;
1160 Double_t devs, xsav, ysav, bwidx, bwidy, unext, ff, xb4;
1161 Int_t i, ngrid, ixmid, nparx, ix,
nx,
ny, ki1, ki2, ixzero, iy, ics;
1166 if (ke1 <= 0 || ke2 <= 0)
goto L1350;
1167 if (ke1 >
fNu || ke2 >
fNu)
goto L1350;
1170 if (ki1 <= 0 || ki2 <= 0)
goto L1350;
1171 if (ki1 == ki2)
goto L1350;
1181 if (devs <= 0) devs = 2;
1182 xlo =
fU[ke1-1] - devs*
fWerr[ki1-1];
1183 xup =
fU[ke1-1] + devs*fWerr[ki1-1];
1184 ylo =
fU[ke2-1] - devs*fWerr[ki2-1];
1185 yup =
fU[ke2-1] + devs*fWerr[ki2-1];
1197 if (nx < 11) nx = 11;
1198 if (ny < 11) ny = 11;
1199 if (nx >= 115) nx = 114;
1210 bwidx = (xup - xlo) /
Double_t(nx);
1211 bwidy = (yup - ylo) /
Double_t(ny);
1212 ixmid =
Int_t(((xsav - xlo)*
Double_t(nx) / (xup - xlo)) + 1);
1213 if (ixmid < 1) ixmid = 1;
1216 for (i = 1; i <= 20; ++i) { contur[i-1] =
fAmin +
fUp*(i-1)*(i-1); }
1217 contur[0] +=
fUp*.01;
1226 for (ix = 1; ix <= nx + 1; ++ix) {
1230 if (xb4 < 0 &&
fU[ke1-1] > 0) ixzero = ix - 1;
1235 Printf(
" Y-AXIS: PARAMETER %3d: %s",ke2,(
const char*)
fCpnam[ke2-1]);
1237 chzero[ixzero-1] =
'+';
1242 for (iy = 1; iy <=
ny; ++iy) {
1243 unext =
fU[ke2-1] - bwidy;
1248 chln[ixmid-1] =
'*';
1249 if (ixzero != 0) chln[ixzero-1] =
':';
1250 if (
fU[ke2-1] > ysav && unext < ysav) chln = chmid;
1251 if (
fU[ke2-1] > 0 && unext < 0) chln = chzero;
1253 ylabel =
fU[ke2-1] + bwidy*.5;
1255 for (ix = 1; ix <= nx + 1; ++ix) {
1256 fcna[ix-1] = fcnb[ix-1];
1262 for (ix = 1; ix <=
nx; ++ix) {
1269 for (ics = 1; ics <= 20; ++ics) {
1270 if (contur[ics-1] > fmn)
goto L240;
1274 if (contur[ics-1] < fmx) chln[ix-1] = clabel[ics-1];
1277 Printf(
" %12.4g %s",ylabel,(
const char*)chln);
1282 chln(ixmid-1,1) =
'I';
1284 Printf(
" %s",(
const char*)chln);
1289 Printf(
" %12.4g%s%12.4g",xlo,(
const char*)chln,xup);
1290 Printf(
" %s%12.4g",(
const char*)chln,xsav);
1292 Printf(
" %12.4g%s%12.4g%s%12.4g",xlo,(
const char*)chln,xsav,(
const char*)chln,xup);
1294 Printf(
" X-AXIS: PARAMETER %3d %s ONE COLUMN=%12.4g"
1295 ,ke1,(
const char*)
fCpnam[ke1-1],bwidx);
1303 Printf(
" INVALID PARAMETER NUMBER(S) REQUESTED. IGNORED.");
1333 Int_t ierr, ipos, i, llist, lenbuf, lnc;
1335 TString comand, crdbuf, ctemp;
1339 lenbuf = crdbuf.
Length();
1344 for (i = 1; i <=
TMath::Min(20,lenbuf); ++i) {
1345 if (crdbuf[i-1] ==
'\'')
break;
1346 if (crdbuf[i-1] ==
' ') {
1354 if (ipos > lenbuf) {
1355 Printf(
" BLANK COMMAND IGNORED.");
1361 if (crdbuf(ipos-1,3) ==
"PAR") {
1367 if (crdbuf(ipos-1,3) ==
"SET INP") {
1373 if (crdbuf(ipos-1,7) ==
"SET TIT") {
1379 if (crdbuf(ipos-1,7) ==
"SET COV") {
1385 ctemp = crdbuf(ipos-1,lenbuf-ipos+1);
1388 Printf(
" COMMAND CANNOT BE INTERPRETED");
1422 Double_t a1, a2, val2mi, val2pl, dc, sclfac, bigdis, sigsav;
1423 Int_t nall, iold,
line, mpar, ierr, inew, move,
next, i, j, nfcol, iercr;
1424 Int_t idist=0, npcol, kints, i2, i1, lr, nfcnco=0, ki1, ki2, ki3, ke3;
1425 Int_t nowpts, istrav, nfmxin, isw2, isw4;
1431 ldebug =
fIdbg[6] >= 1;
1432 if (ke1 <= 0 || ke2 <= 0)
goto L1350;
1433 if (ke1 >
fNu || ke2 >
fNu)
goto L1350;
1436 if (ki1 <= 0 || ki2 <= 0)
goto L1350;
1437 if (ki1 == ki2)
goto L1350;
1438 if (nptu < 4)
goto L1400;
1450 Printf(
" START MNCONTOUR CALCULATION OF %4d POINTS ON CONTOUR.",nptu);
1453 ki3 = 6 - ki1 - ki2;
1455 Printf(
" EACH POINT IS A MINIMUM WITH RESPECT TO PARAMETER %3d %s",ke3,(
const char*)
fCpnam[ke3-1]);
1457 Printf(
" EACH POINT IS A MINIMUM WITH RESPECT TO THE OTHER %3d VARIABLE PARAMETERS.",
fNpar - 2);
1464 mnmnot(ke1, ke2, val2pl, val2mi);
1466 xptu[0] =
fAlim[ke1-1];
1467 mnwarn(
"W",
"MNContour ",
"Contour squeezed by parameter limits.");
1469 if (
fErn[ki1-1] >= 0)
goto L1500;
1470 xptu[0] = u1min +
fErn[ki1-1];
1475 xptu[2] =
fBlim[ke1-1];
1476 mnwarn(
"W",
"MNContour ",
"Contour squeezed by parameter limits.");
1478 if (
fErp[ki1-1] <= 0)
goto L1500;
1479 xptu[2] = u1min +
fErp[ki1-1];
1482 scalx = 1 / (xptu[2] - xptu[0]);
1484 mnmnot(ke2, ke1, val2pl, val2mi);
1486 yptu[1] =
fAlim[ke2-1];
1487 mnwarn(
"W",
"MNContour ",
"Contour squeezed by parameter limits.");
1489 if (
fErn[ki2-1] >= 0)
goto L1500;
1490 yptu[1] = u2min +
fErn[ki2-1];
1494 yptu[3] =
fBlim[ke2-1];
1495 mnwarn(
"W",
"MNContour ",
"Contour squeezed by parameter limits.");
1497 if (
fErp[ki2-1] <= 0)
goto L1500;
1498 yptu[3] = u2min +
fErp[ki2-1];
1501 scaly = 1 / (yptu[3] - yptu[1]);
1505 Printf(
" Plot of four points found by MINOS");
1511 for (i = 2; i <= nall; ++i) {
1512 fXpt[i-1] = xptu[i-2];
1513 fYpt[i-1] = yptu[i-2];
1515 sprintf(
fChpt,
"%s",
" ABCD");
1529 for (i = 1; i <= mpar; ++i) {
fXt[i-1] =
fX[i-1]; }
1530 i__1 = mpar*(mpar + 1) / 2;
1532 for (i = 1; i <= mpar; ++i) {
1542 for (inew = next; inew <= nptu; ++inew) {
1545 for (iold = 1; iold <= inew - 1; ++iold) {
1547 if (i2 == inew) i2 = 1;
1548 d__1 = scalx*(xptu[iold-1] - xptu[i2-1]);
1549 d__2 = scaly*(yptu[iold-1] - yptu[i2-1]);
1550 dist = d__1*d__1 + d__2*d__2;
1551 if (dist > bigdis) {
1558 if (i2 == inew) i2 = 1;
1563 fXmidcr = a1*xptu[i1-1] + a2*xptu[i2-1];
1564 fYmidcr = a1*yptu[i1-1] + a2*yptu[i2-1];
1565 xdir = yptu[i2-1] - yptu[i1-1];
1566 ydir = xptu[i1-1] - xptu[i2-1];
1579 Printf(
" MNCONT CANNOT FIND NEXT POINT ON CONTOUR. ONLY %3d POINTS FOUND.",nowpts);
1583 mnwarn(
"W",
"MNContour ",
"Cannot find midpoint, try closer.");
1589 for (move = nowpts; move >= i1 + 1; --move) {
1590 xptu[move] = xptu[move-1];
1591 yptu[move] = yptu[move-1];
1601 if (nowpts < nptu)
fCstatu =
"INCOMPLETE";
1609 for (i = 2; i <= nall; ++i) {
1610 fXpt[i-1] = xptu[i-2];
1611 fYpt[i-1] = yptu[i-2];
1615 Printf(
" Y-AXIS: PARAMETER %3d %s",ke2,(
const char*)
fCpnam[ke2-1]);
1619 Printf(
" X-AXIS: PARAMETER %3d %s",ke1,(
const char*)
fCpnam[ke1-1]);
1623 npcol = (nowpts + 1) / 2;
1625 Printf(
"%5d POINTS ON CONTOUR. FMIN=%13.5e ERRDEF=%11.3g",nowpts,abest,
fUp);
1627 (
const char*)
fCpnam[ke2-1],
1628 (
const char*)
fCpnam[ke1-1],
1629 (
const char*)
fCpnam[ke2-1]);
1630 for (line = 1; line <= nfcol; ++
line) {
1632 Printf(
" %5d%13.5e%13.5e %5d%13.5e%13.5e",line,xptu[line-1],yptu[line-1],lr,xptu[lr-1],yptu[lr-1]);
1634 if (nfcol < npcol) {
1635 Printf(
" %5d%13.5e%13.5e",npcol,xptu[npcol-1],yptu[npcol-1]);
1642 i__1 = mpar*(mpar + 1) / 2;
1644 for (i = 1; i <= mpar; ++i) {
1663 Printf(
" INVALID PARAMETER NUMBERS.");
1666 Printf(
" LESS THAN FOUR POINTS REQUESTED.");
1672 Printf(
" MNCONT UNABLE TO FIND FOUR POINTS.");
1701 const char *cnumer =
"123456789-.0+";
1704 Int_t ifld, iend, lend, left, nreq, ipos, kcmnd, nextb, ic, ibegin, ltoadd;
1705 Int_t ielmnt, lelmnt[25], nelmnt;
1711 char *crdbuf = (
char*)cardbuf.
Data();
1718 for (ipos = nextb; ipos <= lend; ++ipos) {
1720 if (crdbuf[ipos-1] ==
' ')
continue;
1721 if (crdbuf[ipos-1] ==
',')
goto L250;
1727 for (ipos = ibegin + 1; ipos <= lend; ++ipos) {
1728 if (crdbuf[ipos-1] ==
' ')
goto L250;
1729 if (crdbuf[ipos-1] ==
',')
goto L250;
1735 if (iend >= ibegin) celmnt[ielmnt-1] = &crdbuf[ibegin-1];
1736 else celmnt[ielmnt-1] = cnull;
1737 lelmnt[ielmnt-1] = iend - ibegin + 1;
1738 if (lelmnt[ielmnt-1] > 19) {
1739 Printf(
" MINUIT WARNING: INPUT DATA WORD TOO LONG.");
1740 ctemp = cardbuf(ibegin-1,iend-ibegin+1);
1742 Printf(
" TRUNCATED TO:%s",celmnt[ielmnt-1]);
1743 lelmnt[ielmnt-1] = 19;
1745 if (ipos >= lend)
goto L300;
1746 if (ielmnt >= 25)
goto L300;
1748 for (ipos = iend + 1; ipos <= lend; ++ipos) {
1749 if (crdbuf[ipos-1] ==
' ')
continue;
1751 if (crdbuf[ipos-1] ==
',') nextb = ipos + 1;
1758 command[0] =
' '; command[1] = 0;
1762 if (ielmnt == 0)
goto L900;
1764 for (ielmnt = 1; ielmnt <= nelmnt; ++ielmnt) {
1765 if ( celmnt[ielmnt-1] == cnull)
goto L450;
1766 for (ic = 1; ic <= 13; ++ic) {
1767 if (*celmnt[ielmnt-1] == cnumer[ic-1])
goto L450;
1769 if (kcmnd >= maxcwd)
continue;
1770 left = maxcwd - kcmnd;
1771 ltoadd = lelmnt[ielmnt-1];
1772 if (ltoadd > left) ltoadd = left;
1773 strncpy(&command[kcmnd],celmnt[ielmnt-1],ltoadd);
1775 if (kcmnd == maxcwd)
continue;
1776 command[kcmnd] =
' ';
1786 for (ifld = ielmnt; ifld <= nelmnt; ++ifld) {
1789 nreq = nelmnt - ielmnt + 1;
1790 Printf(
" MINUIT WARNING IN MNCRCK: ");
1791 Printf(
" COMMAND HAS INPUT %5d NUMERIC FIELDS, BUT MINUIT CAN ACCEPT ONLY%3d",nreq,mxp);
1794 if (celmnt[ifld-1] == cnull) plist[llist-1] = 0;
1796 sscanf(celmnt[ifld-1],
"%lf",&plist[llist-1]);
1801 if (lnc <= 0) lnc = 1;
1820 Double_t alsb[3], flsb[3], bmin, bmax, zmid, sdev, zdir, zlim;
1821 Double_t coeff[3], aleft, aulim, fdist, adist, aminsv;
1823 Double_t determ, rt, smalla, aright, aim, tla, tlf, dfda,ecart;
1824 Int_t iout=0, i, ileft, ierev, maxlk, ibest, ik, it;
1825 Int_t noless, iworst=0, iright, itoohi, kex, ipt;
1830 ldebug =
fIdbg[6] >= 1;
1849 for (ik = 1; ik <= 2; ++ik) {
1855 if (
fKe2cr == 0)
continue;
1860 if (
fNvarl[kex-1] <= 1)
continue;
1861 if (zdir == 0)
continue;
1862 zlim =
fAlim[kex-1];
1863 if (zdir > 0) zlim =
fBlim[kex-1];
1864 aulim =
TMath::Min(aulim,(zlim - zmid) / zdir);
1872 mneval(anext, fnext, ierev);
1875 Printf(
" MNCROS: calls=%8d AIM=%10.5f F,A=%10.5f%10.5f",
fNfcn,aim,fnext,aopt);
1877 if (ierev > 0)
goto L900;
1878 if (
fLimset && fnext <= aim)
goto L930;
1880 fXpt[ipt-1] = anext;
1881 fYpt[ipt-1] = fnext;
1882 fChpt[ipt-1] = charal[ipt-1];
1887 if (
TMath::Abs(fnext - aim) < tlf)
goto L800;
1889 if (aopt < -.5)aopt = -.5;
1890 if (aopt > 1) aopt = 1;
1896 mneval(aopt, fnext, ierev);
1899 Printf(
" MNCROS: calls=%8d AIM=%10.5f F,A=%10.5f%10.5f",
fNfcn,aim,fnext,aopt);
1901 if (ierev > 0)
goto L900;
1902 if (
fLimset && fnext <= aim)
goto L930;
1905 fXpt[ipt-1] = alsb[1];
1906 fYpt[ipt-1] = fnext;
1907 fChpt[ipt-1] = charal[ipt-1];
1909 dfda = (flsb[1] - flsb[0]) / (alsb[1] - alsb[0]);
1911 if (dfda > 0)
goto L460;
1913 mnwarn(
"D",
"MNCROS ",
"Looking for slope of the right sign");
1915 for (it = 1; it <= maxlk; ++it) {
1924 mneval(aopt, fnext, ierev);
1927 Printf(
" MNCROS: calls=%8d AIM=%10.5f F,A=%10.5f%10.5f",
fNfcn,aim,fnext,aopt);
1929 if (ierev > 0)
goto L900;
1930 if (
fLimset && fnext <= aim)
goto L930;
1933 fXpt[ipt-1] = alsb[1];
1934 fYpt[ipt-1] = fnext;
1935 fChpt[ipt-1] = charal[ipt-1];
1937 dfda = (flsb[1] - flsb[0]) / (alsb[1] - alsb[0]);
1938 if (dfda > 0)
goto L450;
1940 mnwarn(
"W",
"MNCROS ",
"Cannot find slope of the right sign");
1945 aopt = alsb[1] + (aim - flsb[1]) / dfda;
1950 if (adist < tla && fdist < tlf)
goto L800;
1951 if (ipt >= 15)
goto L950;
1953 if (aopt < bmin) aopt = bmin;
1955 if (aopt > bmax) aopt = bmax;
1962 mneval(aopt, fnext, ierev);
1965 Printf(
" MNCROS: calls=%8d AIM=%10.5f F,A=%10.5f%10.5f",
fNfcn,aim,fnext,aopt);
1967 if (ierev > 0)
goto L900;
1968 if (
fLimset && fnext <= aim)
goto L930;
1971 fXpt[ipt-1] = alsb[2];
1972 fYpt[ipt-1] = fnext;
1973 fChpt[ipt-1] = charal[ipt-1];
1980 for (i = 1; i <= 3; ++i) {
1982 if (ecart > ecarmx) { ecarmx = ecart; iworst = i; }
1983 if (ecart < ecarmn) { ecarmn = ecart; ibest = i; }
1984 if (flsb[i-1] < aim) ++noless;
1987 if (noless == 1 || noless == 2)
goto L500;
1989 if (noless == 0 && ibest != 3)
goto L950;
1992 if (noless == 3 && ibest != 3) {
1998 alsb[iworst-1] = alsb[2];
1999 flsb[iworst-1] = flsb[2];
2000 dfda = (flsb[1] - flsb[0]) / (alsb[1] - alsb[0]);
2004 mnpfit(alsb, flsb, 3, coeff, sdev);
2005 if (coeff[2] <= 0) {
2006 mnwarn(
"D",
"MNCROS ",
"Curvature is negative near contour line.");
2008 determ = coeff[1]*coeff[1] - coeff[2]*4*(coeff[0] - aim);
2010 mnwarn(
"D",
"MNCROS ",
"Problem 2, impossible determinant");
2015 x1 = (-coeff[1] + rt) / (coeff[2]*2);
2016 x2 = (-coeff[1] - rt) / (coeff[2]*2);
2017 s1 = coeff[1] + x1*2*coeff[2];
2018 s2 = coeff[1] + x2*2*coeff[2];
2020 Printf(
" MNCONTour problem 1");
2034 if (ipt >= 15)
goto L950;
2043 for (i = 1; i <= 3; ++i) {
2045 if (ecart < ecarmn) { ecarmn = ecart; ibest = i; }
2046 if (ecart > ecarmx) { ecarmx = ecart; }
2047 if (flsb[i-1] > aim) {
2048 if (iright == 0) iright = i;
2049 else if (flsb[i-1] > flsb[iright-1]) iout = i;
2050 else { iout = iright; iright = i; }
2052 else if (ileft == 0) ileft = i;
2053 else if (flsb[i-1] < flsb[ileft-1]) iout = i;
2054 else { iout = ileft; ileft = i; }
2057 if (ecarmx >
TMath::Abs(flsb[iout-1] - aim)*10) {
2058 aopt = aopt*.5 + (alsb[iright-1] + alsb[ileft-1])*.25;
2062 if (slope*smalla > tlf) smalla = tlf / slope;
2063 aleft = alsb[ileft-1] + smalla;
2064 aright = alsb[iright-1] - smalla;
2066 if (aopt < aleft) aopt = aleft;
2067 if (aopt > aright) aopt = aright;
2068 if (aleft > aright) aopt = (aleft + aright)*.5;
2077 mneval(aopt, fnext, ierev);
2080 Printf(
" MNCROS: calls=%8d AIM=%10.5f F,A=%10.5f%10.5f",
fNfcn,aim,fnext,aopt);
2082 if (ierev > 0)
goto L900;
2083 if (
fLimset && fnext <= aim)
goto L930;
2086 fYpt[ipt-1] = fnext;
2087 fChpt[ipt-1] = charal[ipt-1];
2089 alsb[iout-1] = aopt;
2090 flsb[iout-1] = fnext;
2102 if (ierev == 1)
goto L940;
2119 for (i = 1; i <= ipt; ++i) {
2120 if (
fYpt[i-1] > aim + fUp) {
2128 if (
fXdircr < 0) chsign =
"NEGA";
2130 Printf(
" %sTIVE MINOS ERROR, PARAMETER %3d",chsign,
fKe1cr);
2133 Printf(
"POINTS LABELLED '+' WERE TOO HIGH TO PLOT.");
2136 Printf(
"RIGHTMOST POINT IS UP AGAINST LIMIT.");
2154 Int_t ndex, iext, i, j;
2157 Printf(
" FUNCTION MUST BE MINIMIZED BEFORE CALLING %s",(
const char*)
fCfrom);
2164 mnwarn(
"W",
fCfrom,
"NO ERROR MATRIX. WILL IMPROVISE.");
2165 for (i = 1; i <=
fNpar; ++i) {
2167 for (j = 1; j <= i-1; ++j) {
2172 if (
fG2[i-1] <= 0) {
2175 if (
fNvarl[iext-1] > 1) {
2180 fG2[i-1] =
fUp / (wint*wint);
2202 Double_t step, dfmin, stepb4, dd, df, fs1;
2203 Double_t tlrstp, tlrgrd, epspri, optstp, stpmax, stpmin, fs2, grbfor=0, d1d2, xtf;
2204 Int_t icyc, ncyc, iint, iext, i, nparx;
2208 ldebug =
fIdbg[2] >= 1;
2210 if (
fISW[2] == 1)
goto L100;
2222 Printf(
" FIRST DERIVATIVE DEBUG PRINTOUT. MNDERI");
2223 Printf(
" PAR DERIV STEP MINSTEP OPTSTEP D1-D2 2ND DRV");
2240 for (i = 1; i <=
fNpar; ++i) {
2247 for (icyc = 1; icyc <= ncyc; ++icyc) {
2253 if (
fGstep[i-1] < 0 && step > .5) step = .5;
2256 if (step > stpmax) step = stpmax;
2259 if (step < stpmin) step = stpmin;
2261 if (TMath::Abs((step - stepb4) / step) < tlrstp)
goto L50;
2267 fX[i-1] = xtf + step;
2271 fX[i-1] = xtf - step;
2275 fGrd[i-1] = (fs1 - fs2) / (step*2);
2276 fG2[i-1] = (fs1 + fs2 -
fAmin*2) / (step*step);
2279 d1d2 = (fs1 + fs2 -
fAmin*2) / step;
2280 Printf(
"%4d%11.3g%11.3g%10.2g%10.2g%10.2g%10.2g",i,
fGrd[i-1],step,stpmin,optstp,d1d2,
fG2[i-1]);
2283 if (TMath::Abs(grbfor -
fGrd[i-1]) / (TMath::Abs(
fGrd[i-1]) + dfmin/step) < tlrgrd)
2287 if (ncyc == 1)
goto L50;
2296 for (iint = 1; iint <=
fNpar; ++iint) {
2298 if (
fNvarl[iext-1] <= 1) {
2336 Int_t i, j, k,
l,
m=0, i0, i1, j1, m1, n1;
2340 a_offset = ndima + 1;
2348 for (i1 = 2; i1 <=
n; ++i1) {
2350 f = a[i + (i-1)*ndima];
2353 if (l < 1)
goto L25;
2355 for (k = 1; k <=
l; ++k) {
2356 d__1 = a[i + k*ndima];
2362 if (gl > 1e-35)
goto L30;
2370 if (f >= 0) gl = -gl;
2373 a[i + (i-1)*ndima] = f - gl;
2375 for (j = 1; j <=
l; ++j) {
2376 a[j + i*ndima] = a[i + j*ndima] /
h;
2378 for (k = 1; k <= j; ++k) { gl += a[j + k*ndima]*a[i + k*ndima]; }
2379 if (j >= l)
goto L47;
2381 for (k = j1; k <=
l; ++k) { gl += a[k + j*ndima]*a[i + k*ndima]; }
2383 work[n + j] = gl /
h;
2384 f += gl*a[j + i*ndima];
2387 for (j = 1; j <=
l; ++j) {
2389 gl = work[n + j] - hh*
f;
2391 for (k = 1; k <= j; ++k) {
2392 a[j + k*ndima] = a[j + k*ndima] - f*work[n + k] - gl*a[i + k*ndima];
2401 for (i = 1; i <=
n; ++i) {
2403 if (work[i] == 0 || l == 0)
goto L100;
2405 for (j = 1; j <=
l; ++j) {
2407 for (k = 1; k <=
l; ++k) { gl += a[i + k*ndima]*a[k + j*ndima]; }
2408 for (k = 1; k <=
l; ++k) { a[k + j*ndima] -= gl*a[k + i*ndima]; }
2411 work[i] = a[i + i*ndima];
2413 if (l == 0)
continue;
2415 for (j = 1; j <=
l; ++j) {
2422 for (i = 2; i <=
n; ++i) {
2424 work[i0] = work[i0 + 1];
2429 for (l = 1; l <=
n; ++
l) {
2433 for (m1 = l; m1 <=
n; ++m1) {
2439 if (m == l)
goto L205;
2442 if (j == mits)
return;
2444 pt = (work[l + 1] - work[
l]) / (work[n + l]*2);
2447 if (pt < 0) pr = pt -
r;
2449 h = work[
l] - work[n +
l] / pr;
2450 for (i = l; i <=
n; ++i) { work[i] -=
h; }
2457 for (i1 = l; i1 <= m1; ++i1) {
2464 c = pt / work[n + i];
2466 work[n + j] = s*work[n + i]*
r;
2471 c = work[n + i] /
pt;
2473 work[n + j] = s*pt*
r;
2477 pt = c*work[i] - s*gl;
2478 work[j] = h + s*(c*gl + s*work[i]);
2479 for (k = 1; k <=
n; ++k) {
2481 a[k + j*ndima] = s*a[k + i*ndima] + c*
h;
2482 a[k + i*ndima] = c*a[k + i*ndima] - s*
h;
2493 for (i = 1; i <= n1; ++i) {
2497 for (j = i1; j <=
n; ++j) {
2498 if (work[j] >= pt)
continue;
2503 if (k == i)
continue;
2507 for (j = 1; j <=
n; ++j) {
2508 pt = a[j + i*ndima];
2509 a[j + i*ndima] = a[j + k*ndima];
2510 a[j + k*ndima] =
pt;
2526 Int_t emat_dim1, emat_offset;
2530 Int_t i, j, k, npard, k2, kk, iz, nperln, kga, kgb;
2535 emat_offset = emat_dim1 + 1;
2536 emat -= emat_offset;
2539 if (
fISW[1] < 1)
return;
2541 Printf(
" EXTERNAL ERROR MATRIX. NDIM=%4d NPAR=%3d ERR DEF=%g",ndim,
fNpar,
fUp);
2548 Printf(
" USER-DIMENSIONED ARRAY EMAT NOT BIG ENOUGH. REDUCED MATRIX CALCULATED.");
2555 if (
fISW[4] >= 1 && npard > nperln) {
2556 Printf(
" ELEMENTS ABOVE DIAGONAL ARE NOT PRINTED.");
2559 for (i = 1; i <= npard; ++i) {
2562 for (j = 1; j <= i; ++j) {
2565 emat[i + j*emat_dim1] = dxdi*
fVhmat[kgb-1]*dxdj*
fUp;
2566 emat[j + i*emat_dim1] = emat[i + j*emat_dim1];
2571 for (i = 1; i <= npard; ++i) {
2573 if (npard >= nperln) iz = i;
2575 for (k = 1; nperln < 0 ? k >= iz : k <= iz; k += nperln) {
2576 k2 = k + nperln - 1;
2577 if (k2 > iz) k2 = iz;
2578 for (kk = k; kk <= k2; ++kk) {
2581 Printf(
"%s",(
const char*)ctemp);
2602 Int_t ndiag, iin, iex;
2606 if (iex >
fNu || iex <= 0)
goto L900;
2608 if (iin <= 0)
goto L900;
2611 eplus =
fErp[iin-1];
2612 if (eplus ==
fUndefi) eplus = 0;
2613 eminus =
fErn[iin-1];
2614 if (eminus ==
fUndefi) eminus = 0;
2616 ndiag = iin*(iin + 1) / 2;
2620 if (
fISW[1] < 2)
return;
2659 if (
fISW[0] >= 1) ierev = 1;
2660 if (
fISW[3] < 1) ierev = 2;
2693 static const char *cname[40] = {
2738 Double_t step, xptu[101], yptu[101],
f, rno;
2739 Int_t icol, kcol, ierr, iint, iext, lnow, nptu, i, iflag, ierrf;
2740 Int_t ilist, nparx, izero, nf, lk, it, iw, inonde, nsuper;
2741 Int_t it2, ke1, ke2, nowprt, kll, krl;
2742 TString chwhy, c26, cvblnk, cneway, comd;
2744 Bool_t lfreed, ltofix, lfixed;
2751 if (lk > 20) lk = 20;
2756 for (iw = 1; iw <=
fMaxpar; ++iw) {
2758 if (iw <= llist)
fWord7[iw-1] = plist[iw-1];
2765 if (lnow > 4) lnow = 4;
2768 for (i = 1; i <= lnow; ++i) {
2771 Printf(
"%s",(
const char*)ctemp);
2780 for (i = lnow + 1; i <= kll; ++i) {
2781 Printf(
"%12.4g",plist[i-1]);
2786 Printf(
" ERROR: ABOVE CALL TO MNEXCM TRIED TO PASS MORE THAN %d PARAMETERS.",
fMaxpar);
2804 for (i = 1; i <= nntot; ++i) {
2805 if (strncmp(ctemp.
Data(),cname[i-1],3) == 0)
goto L90;
2807 Printf(
"UNKNOWN COMMAND IGNORED:%s", comand.
Data());
2812 if (
fCword(0,4) ==
"MINO") i = 5;
2813 if (i != 6 && i != 7 && i != 8 && i != 23) {
2828 case 10:
goto L1000;
2829 case 11:
goto L1100;
2830 case 12:
goto L1200;
2831 case 13:
goto L1300;
2832 case 14:
goto L1400;
2833 case 15:
goto L1500;
2834 case 16:
goto L1600;
2835 case 17:
goto L1700;
2836 case 18:
goto L1800;
2837 case 19:
goto L1900;
2838 case 20:
goto L1900;
2839 case 21:
goto L1900;
2840 case 22:
goto L2200;
2841 case 23:
goto L2300;
2842 case 24:
goto L2400;
2843 case 25:
goto L1900;
2844 case 26:
goto L2600;
2845 case 27:
goto L3300;
2846 case 28:
goto L3300;
2847 case 29:
goto L3300;
2848 case 30:
goto L3300;
2849 case 31:
goto L3300;
2850 case 32:
goto L3300;
2851 case 33:
goto L3300;
2852 case 34:
goto L3400;
2853 case 35:
goto L3500;
2854 case 36:
goto L3600;
2855 case 37:
goto L3700;
2856 case 38:
goto L3800;
2857 case 39:
goto L3900;
2858 case 40:
goto L4000;
2867 if (
fISW[3] < 1) ierflg = 4;
2875 if (
fISW[3] >= 1)
return;
2877 if (
fISW[0] == 1)
return;
2878 if (
fCword(0,3) ==
"MIG")
return;
2883 if (
fISW[0] == 1)
return;
2886 if (
fISW[3] >= 1) ierflg = 0;
2902 if (fNfcn < nsuper)
goto L510;
2903 Printf(
" TOO MANY FUNCTION CALLS. MINOS GIVES UP");
2923 Printf(
"%s: NO PARAMETERS REQUESTED ",(
const char*)
fCword);
2926 for (ilist = 1; ilist <= llist; ++ilist) {
2927 iext =
Int_t(plist[ilist-1]);
2928 chwhy =
" IS UNDEFINED.";
2929 if (iext <= 0)
goto L930;
2930 if (iext >
fNu)
goto L930;
2931 if (
fNvarl[iext-1] < 0)
goto L930;
2932 chwhy =
" IS CONSTANT. ";
2933 if (
fNvarl[iext-1] == 0)
goto L930;
2936 chwhy =
" ALREADY FIXED.";
2937 if (iint == 0)
goto L930;
2939 if (ierr == 0) lfixed =
kTRUE;
2942 chwhy =
" ALREADY VARIABLE.";
2943 if (iint > 0)
goto L930;
2950 if (
fISW[4] >= 0)
Printf(
" PARAMETER %4d %s IGNORED.",iext,(
const char*)chwhy);
2952 if (lfreed || lfixed)
mnrset(0);
2965 if (it > 1 || it < 0)
goto L1005;
2976 Printf(
" IGNORED. UNKNOWN ARGUMENT:%4d",it);
2986 if (iext <= 0)
goto L1210;
2989 if (it2 <= 0)
goto L1250;
2995 Printf(
" PARAMETER %4d NOT VARIABLE.",iext);
3007 Printf(
"%s: NO PARAMETERS REQUESTED ",(
const char*)
fCword);
3013 mncntr(ke1-1, ke2-1, ierrf);
3014 if (ierrf > 0) ierflg = 3;
3045 }
else if (f <
fAmin) {
3049 if (
fISW[4] >= 0 && iflag <= 5 && nowprt == 1) {
3054 if (iflag > 5)
mnrset(1);
3065 if (
fISW[4] >= 0)
Printf(
" CALL TO USER FUNCTION WITH IFLAG = 3");
3070 if (
fCword(0,3) ==
"END") ierflg = 10;
3071 if (
fCword(0,3) ==
"RET") ierflg = 12;
3077 Printf(
" MINUIT MEMORY CLEARED. NO PARAMETERS NOW DEFINED.");
3083 for (icol = 5; icol <= lk; ++icol) {
3084 if (
fCword[icol-1] ==
' ')
continue;
3089 if (kcol == 0) comd =
"* ";
3090 else comd =
fCword(kcol-1,lk-kcol+1);
3098 if (ke1 == 0 &&
fNpar == 2) {
3103 if (nptu <= 0) nptu = 20;
3104 if (nptu > 101) nptu = 101;
3106 mncont(ke1-1, ke2-1, nptu, xptu, yptu, ierrf);
3107 if (ierrf < nptu) ierflg = 4;
3108 if (ierrf == -1) ierflg = 3;
3113 if (step <= 0) step = 2;
3116 for (i = 1; i <=
fNpar; ++i) {
3119 fX[i-1] += rno*step*
fWerr[i-1];
3127 Printf(
" BLANK COMMAND IGNORED.");
3133 Printf(
" THE *COVARIANCE* COMMAND IS OSBSOLETE. THE COVARIANCE MATRIX IS NOW SAVED IN A DIFFERENT FORMAT WITH THE *SAVE* COMMAND AND READ IN WITH:*SET COVARIANCE*");
3138 cneway =
"SET PRInt ";
3142 cneway =
"SET GRAd ";
3146 cneway =
"SHOW COVar";
3150 cneway =
"SET ERRdef";
3154 cneway =
"SET LIMits";
3161 Printf(
" OBSOLETE COMMAND:%s PLEASE USE:%s",(
const char*)
fCword
3162 ,(
const char*)cneway);
3164 if (
fCword ==
"SAVE ")
goto L1500;
3182 for (iint = 1; iint <=
fNpar; ++iint) {
3185 pint[iint-1] = pinti;
3199 Int_t kold, nold, ndex, knew, iext, i, j,
m,
n, lc, ik;
3203 Int_t iint = iint1+1;
3204 if (iint >
fNpar || iint <= 0) {
3206 Printf(
" MINUIT ERROR. ARGUMENT TO MNFIXP=%4d",iint);
3212 Printf(
" MINUIT CANNOT FIX PARAMETER %4d MAXIMUM NUMBER THAT CAN BE FIXED IS %d",iext,
fMaxpar);
3232 for (ik = iext + 1; ik <=
fNu; ++ik) {
3246 if (
fISW[1] <= 0)
return;
3248 if (
fNpar <= 0)
return;
3249 for (i = 1; i <= nold; ++i) {
3252 ndex = m*(m-1) / 2 + n;
3258 for (i = 1; i <= nold; ++i) {
3259 for (j = 1; j <= i; ++j) {
3261 if (j == iint || i == iint)
continue;
3286 Double_t grdv, xv, dirinv, g2v, gstepv, xtv;
3287 Int_t i, ipsav, ka, lc, ik,
iq, ir, is;
3290 Printf(
" CALL TO MNFREE IGNORED. ARGUMENT GREATER THAN ONE");
3293 Printf(
" CALL TO MNFREE IGNORED. THERE ARE NO FIXED PARAMETERS");
3295 if (k == 1 || k == 0)
goto L40;
3299 if (
fNiofex[ka-1] == 0)
goto L15;
3300 Printf(
" IGNORED. PARAMETER SPECIFIED IS ALREADY VARIABLE.");
3303 if (
fNpfix < 1)
goto L21;
3304 for (ik = 1; ik <=
fNpfix; ++ik) {
if (
fIpfix[ik-1] == ka)
goto L24; }
3306 Printf(
" PARAMETER %4d NOT FIXED. CANNOT BE RELEASED.",ka);
3309 if (ik ==
fNpfix)
goto L40;
3319 for (i = ik + 1; i <=
fNpfix; ++i) {
3337 if (
fNpfix < 1)
goto L300;
3340 for (ik =
fNu; ik >= ir; --ik) {
3346 fX[lc-1] =
fX[lc-2];
3356 if (is == 0) is =
fNpar;
3360 fX[is-1] =
fXs[iq-1];
3371 Printf(
" PARAMETER %4d %s RESTORED TO VARIABLE.",ir,
3372 (
const char*)
fCpnam[ir-1]);
3374 if (k == 0)
goto L40;
3394 Int_t i, nparx, lc, istsav;
3400 if (
fWord7[0] > 0)
goto L2000;
3414 Printf(
" CHECK OF GRADIENT CALCULATION IN FCN");
3415 Printf(
" PARAMETER G(IN FCN) G(MINUIT) DG(MINUIT) AGREEMENT");
3418 for (lc = 1; lc <=
fNpar; ++lc) {
3428 if (cwd !=
"GOOD")
fISW[2] = 0;
3429 Printf(
" %5d %10s%12.4e%12.4e%12.4e %s",i
3430 ,(
const char*)
fCpnam[i-1]
3434 Printf(
" AGREEMENT=NONE MEANS FCN DID NOT CALCULATE THE DERIVATIVE");
3437 Printf(
" MINUIT DOES NOT ACCEPT DERIVATIVE CALCULATIONS BY FCN");
3438 Printf(
" TO FORCE ACCEPTANCE, ENTER *SET GRAD 1*");
3476 if( comd.
Length() == 0 || comd[0] ==
'*' || comd[0] ==
'?' || comd[0] == 0 || comd==
"HELP" ) {
3477 Printf(
" ==>List of MINUIT Interactive commands:");
3478 Printf(
" CLEar Reset all parameter names and values undefined");
3479 Printf(
" CONtour Make contour map of the user function");
3480 Printf(
" EXIT Exit from Interactive Minuit");
3481 Printf(
" FIX Cause parameter(s) to remain constant");
3482 Printf(
" HESse Calculate the Hessian or error matrix.");
3483 Printf(
" IMPROVE Search for a new minimum around current minimum");
3484 Printf(
" MIGrad Minimize by the method of Migrad");
3485 Printf(
" MINImize MIGRAD + SIMPLEX method if Migrad fails");
3486 Printf(
" MINOs Exact (non-linear) parameter error analysis");
3487 Printf(
" MNContour Calculate one MINOS function contour");
3488 Printf(
" PARameter Define or redefine new parameters and values");
3489 Printf(
" RELease Make previously FIXed parameters variable again");
3490 Printf(
" REStore Release last parameter fixed");
3491 Printf(
" SAVe Save current parameter values on a file");
3492 Printf(
" SCAn Scan the user function by varying parameters");
3493 Printf(
" SEEk Minimize by the method of Monte Carlo");
3494 Printf(
" SET Set various MINUIT constants or conditions");
3495 Printf(
" SHOw Show values of current constants or conditions");
3496 Printf(
" SIMplex Minimize by the method of Simplex");
3505 if( !strncmp(comd.
Data(),
"CLE",3) ) {
3507 Printf(
" Resets all parameter names and values to undefined.");
3508 Printf(
" Must normally be followed by a PARameters command or ");
3509 Printf(
" equivalent, in order to define parameter values.");
3517 if( !strncmp(comd.
Data(),
"CON",3) ) {
3518 Printf(
" ***>CONTOUR <par1> <par2> [devs] [ngrid]");
3519 Printf(
" Instructs Minuit to trace contour lines of the user function");
3520 Printf(
" with respect to the two parameters whose external numbers");
3521 Printf(
" are <par1> and <par2>.");
3522 Printf(
" Other variable parameters of the function, if any, will have");
3523 Printf(
" their values fixed at the current values during the contour");
3524 Printf(
" tracing. The optional parameter [devs] (default value 2.)");
3525 Printf(
" gives the number of standard deviations in each parameter");
3526 Printf(
" which should lie entirely within the plotting area.");
3527 Printf(
" Optional parameter [ngrid] (default value 25 unless page");
3528 Printf(
" size is too small) determines the resolution of the plot,");
3529 Printf(
" i.e. the number of rows and columns of the grid at which the");
3530 Printf(
" function will be evaluated. [See also MNContour.]");
3538 if( !strncmp(comd.
Data(),
"END",3) ) {
3540 Printf(
" Signals the end of a data block (i.e., the end of a fit),");
3541 Printf(
" and implies that execution should continue, because another");
3542 Printf(
" Data Block follows. A Data Block is a set of Minuit data");
3543 Printf(
" consisting of");
3545 Printf(
" (2) One or more Parameter Definitions,");
3546 Printf(
" (3) A blank line, and");
3547 Printf(
" (4) A set of Minuit Commands.");
3548 Printf(
" The END command is used when more than one Data Block is to");
3549 Printf(
" be used with the same FCN function. It first causes Minuit");
3550 Printf(
" to issue a CALL FCN with IFLAG=3, in order to allow FCN to");
3551 Printf(
" perform any calculations associated with the final fitted");
3552 Printf(
" parameter values, unless a CALL FCN 3 command has already");
3553 Printf(
" been executed at the current FCN value.");
3561 if( !strncmp(comd.
Data(),
"EXI",3) ) {
3563 Printf(
" Signals the end of execution.");
3564 Printf(
" The EXIT command first causes Minuit to issue a CALL FCN");
3565 Printf(
" with IFLAG=3, to allow FCN to perform any calculations");
3566 Printf(
" associated with the final fitted parameter values, unless a");
3567 Printf(
" CALL FCN 3 command has already been executed.");
3575 if( !strncmp(comd.
Data(),
"FIX",3) ) {
3576 Printf(
" ***>FIX} <parno> [parno] ... [parno]");
3577 Printf(
" Causes parameter(s) <parno> to be removed from the list of");
3578 Printf(
" variable parameters, and their value(s) will remain constant");
3579 Printf(
" during subsequent minimizations, etc., until another command");
3580 Printf(
" changes their value(s) or status.");
3588 if( !strncmp(comd.
Data(),
"HES",3) ) {
3589 Printf(
" ***>HESse [maxcalls]");
3590 Printf(
" Calculate, by finite differences, the Hessian or error matrix.");
3591 Printf(
" That is, it calculates the full matrix of second derivatives");
3592 Printf(
" of the function with respect to the currently variable");
3593 Printf(
" parameters, and inverts it, printing out the resulting error");
3594 Printf(
" matrix. The optional argument [maxcalls] specifies the");
3595 Printf(
" (approximate) maximum number of function calls after which");
3596 Printf(
" the calculation will be stopped.");
3604 if( !strncmp(comd.
Data(),
"IMP",3) ) {
3605 Printf(
" ***>IMPROVE [maxcalls]");
3606 Printf(
" If a previous minimization has converged, and the current");
3607 Printf(
" values of the parameters therefore correspond to a local");
3608 Printf(
" minimum of the function, this command requests a search for");
3609 Printf(
" additional distinct local minima.");
3610 Printf(
" The optional argument [maxcalls] specifies the (approximate");
3611 Printf(
" maximum number of function calls after which the calculation");
3612 Printf(
" will be stopped.");
3620 if( !strncmp(comd.
Data(),
"MIG",3) ) {
3621 Printf(
" ***>MIGrad [maxcalls] [tolerance]");
3622 Printf(
" Causes minimization of the function by the method of Migrad,");
3623 Printf(
" the most efficient and complete single method, recommended");
3624 Printf(
" for general functions (see also MINImize).");
3625 Printf(
" The minimization produces as a by-product the error matrix");
3626 Printf(
" of the parameters, which is usually reliable unless warning");
3627 Printf(
" messages are produced.");
3628 Printf(
" The optional argument [maxcalls] specifies the (approximate)");
3629 Printf(
" maximum number of function calls after which the calculation");
3630 Printf(
" will be stopped even if it has not yet converged.");
3631 Printf(
" The optional argument [tolerance] specifies required tolerance");
3632 Printf(
" on the function value at the minimum.");
3633 Printf(
" The default tolerance is 0.1, and the minimization will stop");
3634 Printf(
" when the estimated vertical distance to the minimum (EDM) is");
3635 Printf(
" less than 0.001*[tolerance]*UP (see [SET ERRordef]).");
3643 if( !strncmp(comd.
Data(),
"MINI",4) ) {
3644 Printf(
" ***>MINImize [maxcalls] [tolerance]");
3645 Printf(
" Causes minimization of the function by the method of Migrad,");
3646 Printf(
" as does the MIGrad command, but switches to the SIMplex method");
3647 Printf(
" if Migrad fails to converge. Arguments are as for MIGrad.");
3648 Printf(
" Note that command requires four characters to be unambiguous.");
3656 if( !strncmp(comd.
Data(),
"MIN0",4) ) {
3657 Printf(
" ***>MINOs [maxcalls] [parno] [parno] ...");
3658 Printf(
" Causes a Minos error analysis to be performed on the parameters");
3659 Printf(
" whose numbers [parno] are specified. If none are specified,");
3660 Printf(
" Minos errors are calculated for all variable parameters.");
3661 Printf(
" Minos errors may be expensive to calculate, but are very");
3662 Printf(
" reliable since they take account of non-linearities in the");
3663 Printf(
" problem as well as parameter correlations, and are in general");
3665 Printf(
" The optional argument [maxcalls] specifies the (approximate)");
3666 Printf(
" maximum number of function calls per parameter requested,");
3667 Printf(
" after which the calculation will stop for that parameter.");
3675 if( !strncmp(comd.
Data(),
"MNC",3) ) {
3676 Printf(
" ***>MNContour <par1> <par2> [npts]");
3677 Printf(
" Calculates one function contour of FCN with respect to");
3678 Printf(
" parameters par1 and par2, with FCN minimized always with");
3679 Printf(
" respect to all other NPAR-2 variable parameters (if any).");
3680 Printf(
" Minuit will try to find npts points on the contour (default 20)");
3681 Printf(
" If only two parameters are variable at the time, it is not");
3682 Printf(
" necessary to specify their numbers. To calculate more than");
3683 Printf(
" one contour, it is necessary to SET ERRordef to the appropriate");
3684 Printf(
" value and issue the MNContour command for each contour.");
3692 if( !strncmp(comd.
Data(),
"PAR",3) ) {
3693 Printf(
" ***>PARameters");
3694 Printf(
" followed by one or more parameter definitions.");
3695 Printf(
" Parameter definitions are of the form:");
3696 Printf(
" <number> ''name'' <value> <step> [lolim] [uplim] ");
3698 Printf(
" 3 ''K width'' 1.2 0.1");
3699 Printf(
" the last definition is followed by a blank line or a zero.");
3707 if( !strncmp(comd.
Data(),
"REL",3) ) {
3708 Printf(
" ***>RELease <parno> [parno] ... [parno]");
3709 Printf(
" If <parno> is the number of a previously variable parameter");
3710 Printf(
" which has been fixed by a command: FIX <parno>, then that");
3711 Printf(
" parameter will return to variable status. Otherwise a warning");
3712 Printf(
" message is printed and the command is ignored.");
3713 Printf(
" Note that this command operates only on parameters which were");
3714 Printf(
" at one time variable and have been FIXed. It cannot make");
3715 Printf(
" constant parameters variable; that must be done by redefining");
3716 Printf(
" the parameter with a PARameters command.");
3724 if( !strncmp(comd.
Data(),
"RES",3) ) {
3725 Printf(
" ***>REStore [code]");
3726 Printf(
" If no [code] is specified, this command restores all previously");
3727 Printf(
" FIXed parameters to variable status. If [code]=1, then only");
3728 Printf(
" the last parameter FIXed is restored to variable status.");
3729 Printf(
" If code is neither zero nor one, the command is ignored.");
3737 if( !strncmp(comd.
Data(),
"RET",3) ) {
3739 Printf(
" Signals the end of a data block, and instructs Minuit to return");
3740 Printf(
" to the program which called it. The RETurn command first");
3741 Printf(
" causes Minuit to CALL FCN with IFLAG=3, in order to allow FCN");
3742 Printf(
" to perform any calculations associated with the final fitted");
3743 Printf(
" parameter values, unless a CALL FCN 3 command has already been");
3744 Printf(
" executed at the current FCN value.");
3752 if( !strncmp(comd.
Data(),
"SAV",3) ) {
3754 Printf(
" Causes the current parameter values to be saved on a file in");
3755 Printf(
" such a format that they can be read in again as Minuit");
3756 Printf(
" parameter definitions. If the covariance matrix exists, it is");
3757 Printf(
" also output in such a format. The unit number is by default 7,");
3758 Printf(
" or that specified by the user in their call to MINTIO or");
3759 Printf(
" MNINIT. The user is responsible for opening the file previous");
3760 Printf(
" to issuing the [SAVe] command (except where this can be done");
3761 Printf(
" interactively).");
3769 if( !strncmp(comd.
Data(),
"SCA",3) ) {
3770 Printf(
" ***>SCAn [parno] [numpts] [from] [to]");
3771 Printf(
" Scans the value of the user function by varying parameter");
3772 Printf(
" number [parno], leaving all other parameters fixed at the");
3773 Printf(
" current value. If [parno] is not specified, all variable");
3774 Printf(
" parameters are scanned in sequence.");
3775 Printf(
" The number of points [numpts] in the scan is 40 by default,");
3776 Printf(
" and cannot exceed 100. The range of the scan is by default");
3777 Printf(
" 2 standard deviations on each side of the current best value,");
3778 Printf(
" but can be specified as from [from] to [to].");
3779 Printf(
" After each scan, if a new minimum is found, the best parameter");
3780 Printf(
" values are retained as start values for future scans or");
3781 Printf(
" minimizations. The curve resulting from each scan is plotted");
3782 Printf(
" on the output unit in order to show the approximate behaviour");
3783 Printf(
" of the function.");
3784 Printf(
" This command is not intended for minimization, but is sometimes");
3785 Printf(
" useful for debugging the user function or finding a");
3786 Printf(
" reasonable starting point.");
3794 if( !strncmp(comd.
Data(),
"SEE",3) ) {
3795 Printf(
" ***>SEEk [maxcalls] [devs]");
3796 Printf(
" Causes a Monte Carlo minimization of the function, by choosing");
3797 Printf(
" random values of the variable parameters, chosen uniformly");
3798 Printf(
" over a hypercube centered at the current best value.");
3799 Printf(
" The region size is by default 3 standard deviations on each");
3800 Printf(
" side, but can be changed by specifying the value of [devs].");
3808 if( !strncmp(comd.
Data(),
"SET",3) ) {
3809 Printf(
" ***>SET <option_name>");
3811 Printf(
" Informs Minuit that it is running in batch mode.");
3814 Printf(
" SET EPSmachine <accuracy>");
3815 Printf(
" Informs Minuit that the relative floating point arithmetic");
3816 Printf(
" precision is <accuracy>. Minuit determines the nominal");
3817 Printf(
" precision itself, but the SET EPSmachine command can be");
3818 Printf(
" used to override Minuit own determination, when the user");
3819 Printf(
" knows that the FCN function value is not calculated to");
3820 Printf(
" the nominal machine accuracy. Typical values of <accuracy>");
3821 Printf(
" are between 10**-5 and 10**-14.");
3824 Printf(
" SET ERRordef <up>");
3825 Printf(
" Sets the value of UP (default value= 1.), defining parameter");
3826 Printf(
" errors. Minuit defines parameter errors as the change");
3827 Printf(
" in parameter value required to change the function value");
3828 Printf(
" by UP. Normally, for chisquared fits UP=1, and for negative");
3829 Printf(
" log likelihood, UP=0.5.");
3832 Printf(
" SET GRAdient [force]");
3833 Printf(
" Informs Minuit that the user function is prepared to");
3834 Printf(
" calculate its own first derivatives and return their values");
3835 Printf(
" in the array GRAD when IFLAG=2 (see specs of FCN).");
3836 Printf(
" If [force] is not specified, Minuit will calculate");
3837 Printf(
" the FCN derivatives by finite differences at the current");
3838 Printf(
" point and compare with the user calculation at that point,");
3839 Printf(
" accepting the user values only if they agree.");
3840 Printf(
" If [force]=1, Minuit does not do its own derivative");
3841 Printf(
" calculation, and uses the derivatives calculated in FCN.");
3844 Printf(
" SET INPut [unitno] [filename]");
3845 Printf(
" Causes Minuit, in data-driven mode only, to read subsequent");
3846 Printf(
" commands (or parameter definitions) from a different input");
3847 Printf(
" file. If no [unitno] is specified, reading reverts to the");
3848 Printf(
" previous input file, assuming that there was one.");
3849 Printf(
" If [unitno] is specified, and that unit has not been opened,");
3850 Printf(
" then Minuit attempts to open the file [filename]} if a");
3851 Printf(
" name is specified. If running in interactive mode and");
3852 Printf(
" [filename] is not specified and [unitno] is not opened,");
3853 Printf(
" Minuit prompts the user to enter a file name.");
3854 Printf(
" If the word REWIND is added to the command (note:no blanks");
3855 Printf(
" between INPUT and REWIND), the file is rewound before");
3856 Printf(
" reading. Note that this command is implemented in standard");
3857 Printf(
" Fortran 77 and the results may depend on the system;");
3858 Printf(
" for example, if a filename is given under VM/CMS, it must");
3859 Printf(
" be preceded by a slash.");
3862 Printf(
" SET INTeractive");
3863 Printf(
" Informs Minuit that it is running interactively.");
3866 Printf(
" SET LIMits [parno] [lolim] [uplim]");
3867 Printf(
" Allows the user to change the limits on one or all");
3868 Printf(
" parameters. If no arguments are specified, all limits are");
3869 Printf(
" removed from all parameters. If [parno] alone is specified,");
3870 Printf(
" limits are removed from parameter [parno].");
3871 Printf(
" If all arguments are specified, then parameter [parno] will");
3872 Printf(
" be bounded between [lolim] and [uplim].");
3873 Printf(
" Limits can be specified in either order, Minuit will take");
3874 Printf(
" the smaller as [lolim] and the larger as [uplim].");
3875 Printf(
" However, if [lolim] is equal to [uplim], an error condition");
3879 Printf(
" SET LINesperpage");
3880 Printf(
" Sets the number of lines for one page of output.");
3881 Printf(
" Default value is 24 for interactive mode");
3884 Printf(
" SET NOGradient");
3885 Printf(
" The inverse of SET GRAdient, instructs Minuit not to");
3886 Printf(
" use the first derivatives calculated by the user in FCN.");
3889 Printf(
" SET NOWarnings");
3890 Printf(
" Supresses Minuit warning messages.");
3893 Printf(
" SET OUTputfile <unitno>");
3894 Printf(
" Instructs Minuit to write further output to unit <unitno>.");
3897 Printf(
" SET PAGethrow <integer>");
3898 Printf(
" Sets the carriage control character for ``new page'' to");
3899 Printf(
" <integer>. Thus the value 1 produces a new page, and 0");
3900 Printf(
" produces a blank line, on some devices (see TOPofpage)");
3904 Printf(
" SET PARameter <parno> <value>");
3905 Printf(
" Sets the value of parameter <parno> to <value>.");
3906 Printf(
" The parameter in question may be variable, fixed, or");
3907 Printf(
" constant, but must be defined.");
3910 Printf(
" SET PRIntout <level>");
3911 Printf(
" Sets the print level, determining how much output will be");
3912 Printf(
" produced. Allowed values and their meanings are displayed");
3913 Printf(
" after a SHOw PRInt command, and are currently <level>=:");
3914 Printf(
" [-1] no output except from SHOW commands");
3915 Printf(
" [0] minimum output");
3916 Printf(
" [1] default value, normal output");
3917 Printf(
" [2] additional output giving intermediate results.");
3918 Printf(
" [3] maximum output, showing progress of minimizations.");
3919 Printf(
" Note: See also the SET WARnings command.");
3922 Printf(
" SET RANdomgenerator <seed>");
3923 Printf(
" Sets the seed of the random number generator used in SEEk.");
3924 Printf(
" This can be any integer between 10000 and 900000000, for");
3925 Printf(
" example one which was output from a SHOw RANdom command of");
3926 Printf(
" a previous run.");
3929 Printf(
" SET STRategy <level>");
3930 Printf(
" Sets the strategy to be used in calculating first and second");
3931 Printf(
" derivatives and in certain minimization methods.");
3932 Printf(
" In general, low values of <level> mean fewer function calls");
3933 Printf(
" and high values mean more reliable minimization.");
3934 Printf(
" Currently allowed values are 0, 1 (default), and 2.");
3938 Printf(
" Informs Minuit that the next input line is to be considered");
3939 Printf(
" the (new) title for this task or sub-task. This is for");
3940 Printf(
" the convenience of the user in reading their output.");
3944 Printf(
" Instructs Minuit to output warning messages when suspicious");
3945 Printf(
" conditions arise which may indicate unreliable results.");
3946 Printf(
" This is the default.");
3949 Printf(
" SET WIDthpage");
3950 Printf(
" Informs Minuit of the output page width.");
3951 Printf(
" Default values are 80 for interactive jobs");
3959 if( !strncmp(comd.
Data(),
"SHO",3) ) {
3960 Printf(
" ***>SHOw <option_name>");
3961 Printf(
" All SET XXXX commands have a corresponding SHOw XXXX command.");
3962 Printf(
" In addition, the SHOw commands listed starting here have no");
3963 Printf(
" corresponding SET command for obvious reasons.");
3966 Printf(
" SHOw CORrelations");
3967 Printf(
" Calculates and prints the parameter correlations from the");
3968 Printf(
" error matrix.");
3971 Printf(
" SHOw COVariance");
3972 Printf(
" Prints the (external) covariance (error) matrix.");
3975 Printf(
" SHOw EIGenvalues");
3976 Printf(
" Calculates and prints the eigenvalues of the covariance");
3980 Printf(
" SHOw FCNvalue");
3981 Printf(
" Prints the current value of FCN.");
3989 if( !strncmp(comd.
Data(),
"SIM",3) ) {
3990 Printf(
" ***>SIMplex [maxcalls] [tolerance]");
3991 Printf(
" Performs a function minimization using the simplex method of");
3992 Printf(
" Nelder and Mead. Minimization terminates either when the");
3993 Printf(
" function has been called (approximately) [maxcalls] times,");
3994 Printf(
" or when the estimated vertical distance to minimum (EDM) is");
3995 Printf(
" less than [tolerance].");
3996 Printf(
" The default value of [tolerance] is 0.1*UP(see SET ERRordef).");
4004 if( !strncmp(comd.
Data(),
"STA",3) ) {
4013 if( !strncmp(comd.
Data(),
"STO",3) ) {
4015 Printf(
" Same as EXIT.");
4023 if( !strncmp(comd.
Data(),
"TOP",3) ) {
4024 Printf(
" ***>TOPofpage");
4025 Printf(
" Causes Minuit to write the character specified in a");
4026 Printf(
" SET PAGethrow command (default = 1) to column 1 of the output");
4027 Printf(
" file, which may or may not position your output medium to");
4028 Printf(
" the top of a page depending on the device and system.");
4032 Printf(
" Unknown MINUIT command. Type HELP for list of commands.");
4051 Double_t dmin_, dxdi, elem, wint, tlrg2,
d, dlast, ztemp, g2bfor;
4052 Double_t df, aimsag, fs1, tlrstp, fs2, stpinm, g2i, sag=0, xtf, xti, xtj;
4053 Int_t icyc, ncyc, ndex, idrv, iext, npar2, i, j, ifail, npard, nparx,
id, multpy;
4056 ldebug =
fIdbg[3] >= 1;
4073 if (
fISW[4] >= 2 || ldebug) {
4074 Printf(
" START COVARIANCE MATRIX CALCULATION.");
4090 Printf(
" PAR D GSTEP D G2 GRD SAG ");
4100 for (i = 1; i <= npar2; ++i) {
fVhmat[i-1] = 0; }
4104 for (
id = 1;
id <= npard; ++
id) {
4105 i =
id +
fNpar - npard;
4107 if (
fG2[i-1] == 0) {
4108 mnwarn(
"W",
"HESSE",
Form(
"Second derivative enters zero, param %d",iext));
4110 if (
fNvarl[iext-1] > 1) {
4115 fG2[i-1] =
fUp / (wint*wint);
4123 for (icyc = 1; icyc <= ncyc; ++icyc) {
4125 for (multpy = 1; multpy <= 5; ++multpy) {
4135 sag = (fs1 + fs2 -
fAmin*2)*.5;
4136 if (sag != 0)
goto L30;
4138 if (d >= .5)
goto L26;
4140 if (d > .5) d = .51;
4151 fG2[i-1] = sag*2 / (d*
d);
4152 fGrd[i-1] = (fs1 - fs2) / (d*2);
4154 Printf(
"%4d%2d%12.5g%12.5g%12.5g%12.5g%12.5g",i,idrv,
fGstep[i-1],d,
fG2[i-1],
fGrd[i-1],sag);
4165 if (d < dmin_) d = dmin_;
4167 if (TMath::Abs((d - dlast) / d) < tlrstp ||
4168 TMath::Abs((
fG2[i-1] - g2bfor) /
fG2[i-1]) < tlrg2) {
4179 ndex = i*(i + 1) / 2;
4190 if (
fNpar == 1)
goto L214;
4191 for (i = 1; i <=
fNpar; ++i) {
4192 for (j = 1; j <= i-1; ++j) {
4196 fX[j-1] = xtj + fDirin[j-1];
4202 fDirin[i-1]*fDirin[j-1]);
4203 ndex = i*(i-1) / 2 + j;
4211 for (i = 1; i <=
fNpar; ++i) {
4212 for (j = 1; j <= i; ++j) {
4213 ndex = i*(i-1) / 2 + j;
4220 mnwarn(
"W",
"HESSE",
"Matrix inversion fails.");
4226 for (i = 1; i <=
fNpar; ++i) {
4229 for (j = 1; j <= i-1; ++j) {
4241 Printf(
" COVARIANCE MATRIX CALCULATED SUCCESSFULLY");
4250 Printf(
" MNHESS FAILS AND WILL RETURN DIAGONAL MATRIX. ");
4252 for (i = 1; i <=
fNpar; ++i) {
4254 for (j = 1; j <= i-1; ++j) {
4260 if (g2i <= 0) g2i = 1;
4261 fVhmat[ndex-1] = 2 / g2i;
4277 Double_t dmin_,
d, dfmin, dgmin=0, change, chgold, grdold=0, epspri;
4278 Double_t fs1, optstp, fs2, grdnew=0, sag, xtf;
4279 Int_t icyc, ncyc=0, idrv, i, nparx;
4282 ldebug =
fIdbg[5] >= 1;
4290 for (i = 1; i <=
fNpar; ++i) {
4296 if (d > optstp) d = optstp;
4297 if (d < dmin_) d = dmin_;
4300 for (icyc = 1; icyc <= ncyc; ++icyc) {
4309 sag = (fs1 + fs2 -
fAmin*2)*.5;
4311 grdnew = (fs1 - fs2) / (d*2);
4314 Printf(
"%4d%2d%12.5g%12.5g%12.5g%12.5g%12.5g",i,idrv,
fGstep[i-1],d,
fG2[i-1],grdnew,sag);
4316 if (grdnew == 0)
goto L60;
4317 change =
TMath::Abs((grdold - grdnew) / grdnew);
4318 if (change > chgold && icyc > 1)
goto L60;
4324 if (change < .05)
goto L60;
4325 if (TMath::Abs(grdold - grdnew) < dgmin)
goto L60;
4327 mnwarn(
"D",
"MNHES1",
"Step size too small for 1st drv.");
4358 Double_t amax, ycalf, ystar, ystst;
4360 Int_t npfn, ndex, loop=0, i, j, ifail, iseed=0;
4361 Int_t jhold, nloop, nparx, nparp1, jh, jl, iswtr;
4363 if (
fNpar <= 0)
return;
4370 if (nloop <= 0) nloop =
fNpar + 4;
4377 for (i = 1; i <=
fNpar; ++i) {
4380 for (j = 1; j <= i; ++j) {
4381 ndex = i*(i-1) / 2 + j;
4387 if (ifail >= 1)
goto L280;
4389 for (i = 1; i <=
fNpar; ++i) {
4391 for (j = 1; j <= i; ++j) {
4399 for (i = 1; i <=
fNpar; ++i) {
4407 Printf(
"START ATTEMPT NO.%2d TO FIND NEW MINIMUM",loop);
4417 for (i = 1; i <=
fNpar; ++i) {
4420 fX[i-1] = xi -
fDirin[i-1]*(rnum - .5);
4426 }
else if (
fIMPRy[i-1] > amax) {
4439 if (
fAmin < 0)
goto L95;
4440 if (
fISW[1] <= 2)
goto L280;
4442 if (sig2 < ep &&
fEDM < ep)
goto L100;
4446 for (i = 1; i <=
fNpar; ++i) {
4454 if (ystar >=
fAmin)
goto L70;
4456 for (i = 1; i <=
fNpar; ++i) {
4461 if (ystst <
fIMPRy[jl-1])
goto L67;
4469 if (ystar >=
fIMPRy[jh-1])
goto L73;
4472 if (jhold != jh)
goto L50;
4475 for (i = 1; i <=
fNpar; ++i) {
4480 if (ystst >
fIMPRy[jh-1])
goto L30;
4482 if (ystst <
fAmin)
goto L67;
4488 Printf(
" AN IMPROVEMENT ON THE PREVIOUS MINIMUM HAS BEEN FOUND");
4495 for (i = 1; i <=
fNpar; ++i) {
4505 for (i = 1; i <=
fNpar; ++i) {
4523 Printf(
" IMPROVE HAS FOUND A TRULY NEW MINIMUM");
4524 Printf(
" *************************************");
4530 Printf(
" COVARIANCE MATRIX WAS NOT POSITIVE-DEFINITE");
4536 for (i = 1; i <=
fNpar; ++i) {
4545 Printf(
" IMPROVE HAS RETURNED TO REGION OF ORIGINAL MINIMUM");
4549 if (
fISW[1] < 2)
goto L380;
4550 if (loop < nloop &&
fISW[0] < 1)
goto L20;
4567 for (j = 0; j <
fNpar; ++j) {
4588 Double_t piby2, epstry, epsbak, distnn;
4606 fCovmes[0] =
"NO ERROR MATRIX ";
4607 fCovmes[1] =
"ERR MATRIX APPROXIMATE";
4608 fCovmes[2] =
"ERR MATRIX NOT POS-DEF";
4609 fCovmes[3] =
"ERROR MATRIX ACCURATE ";
4627 for (idb = 0; idb <= 10; ++
idb) {
fIdbg[
idb] = 0; }
4647 for (i = 1; i <= 100; ++i) {
4651 if (epsbak < epstry)
goto L35;
4655 Printf(
" MNINIT UNABLE TO DETERMINE ARITHMETIC PRECISION. WILL ASSUME:%g",
fEpsmac);
4685 if (i2 >
fMaxext || i2 < 0)
goto L900;
4686 if (i2 > 0)
goto L30;
4690 for (inu = 1; inu <=
fNu; ++inu) {
4691 if (
fNvarl[inu-1] <= 0)
continue;
4692 if (
fNvarl[inu-1] == 1 && newcod == 1)
continue;
4697 Printf(
" LIMITS NOT CHANGED FOR FIXED PARAMETER:%4d",inu);
4704 Printf(
" LIMITS REMOVED FROM PARAMETER :%3d",inu);
4708 snew =
fGstep[kint-1]*dxdi;
4716 Printf(
" PARAMETER %3d LIMITS SET TO %15.5g%15.5g",inu,
fAlim[inu-1],
fBlim[inu-1]);
4727 Printf(
" PARAMETER %3d IS NOT VARIABLE.", i2);
4733 Printf(
" REQUEST TO CHANGE LIMITS ON FIXED PARAMETER:%3d",i2);
4734 for (ifx = 1; ifx <=
fNpfix; ++ifx) {
4735 if (i2 ==
fIpfix[ifx-1])
goto L92;
4737 Printf(
" MINUIT BUG IN MNLIMS. SEE F. JAMES");
4745 Printf(
" LIMITS REMOVED FROM PARAMETER %2d",i2);
4754 fGrd[kint-1] *= dxdi;
4758 Printf(
" NO LIMITS SPECIFIED. PARAMETER %3d IS ALREADY UNLIMITED. NO CHANGE.",i2);
4767 Printf(
" PARAMETER %3d LIMITS SET TO %15.5g%15.5g",i2,
fAlim[i2-1],
fBlim[i2-1]);
4770 if (kint <= 0)
fGsteps[ifx-1] = -.1;
4771 else fGstep[kint-1] = -.1;
4774 if (
fCstatu !=
"NO CHANGE ") {
4799 Double_t xpq[12], ypq[12], slam, sdev, coeff[3], denom, flast;
4800 Double_t fvals[3], xvals[3],
f1, fvmin, xvmin, ratio,
f2,
f3 = 0., fvmax;
4801 Double_t toler8, toler9, overal, undral, slamin, slamax, slopem;
4802 Int_t i, nparx=0, nvmax=0, nxypt, kk, ipt;
4809 l65 = 0; l70 = 0; l80 = 0;
4810 ldebug =
fIdbg[1] >= 1;
4819 Printf(
" MNLINE start point not consistent, F values, parameters=");
4820 for (kk = 1; kk <=
fNpar; ++kk) {
4829 chpq[0] = charal[0];
4834 for (i = 1; i <=
fNpar; ++i) {
4835 if (step[i-1] != 0) {
4837 if (slamin == 0) slamin = ratio;
4838 if (ratio < slamin) slamin = ratio;
4840 fX[i-1] = start[i-1] + step[i-1];
4842 if (slamin == 0) slamin =
fEpsmac;
4849 chpq[nxypt-1] = charal[nxypt-1];
4864 denom = (flast - fstart - slope*slam)*2 / (slam*slam);
4866 if (denom != 0) slam = -slope / denom;
4867 if (slam < 0) slam = slamax;
4868 if (slam > slamax) slam = slamax;
4869 if (slam < toler8) slam = toler8;
4870 if (slam < slamin) {
4874 if (
TMath::Abs(slam - 1) < toler8 && f1 < fstart) {
4878 if (
TMath::Abs(slam - 1) < toler8) slam = toler8 + 1;
4883 for (i = 1; i <=
fNpar; ++i) {
fX[i-1] = start[i-1] + slam*step[i-1]; }
4888 chpq[nxypt-1] = charal[nxypt-1];
4889 xpq[nxypt-1] = slam;
4895 if (fstart == fvmin) {
4897 toler8 = toler*slam;
4898 overal = slam - toler8;
4901 }
while (fstart == fvmin);
4903 if (!l65 && !l70 && !l80) {
4907 xvals[1] = xpq[nxypt-2];
4908 fvals[1] = ypq[nxypt-2];
4909 xvals[2] = xpq[nxypt-1];
4910 fvals[2] = ypq[nxypt-1];
4914 mnpfit(xvals, fvals, 3, coeff, sdev);
4915 if (coeff[2] <= 0) {
4916 slopem = coeff[2]*2*xvmin + coeff[1];
4917 if (slopem <= 0) slam = xvmin + slamax;
4918 else slam = xvmin - slamax;
4920 slam = -coeff[1] / (coeff[2]*2);
4921 if (slam > xvmin + slamax) slam = xvmin + slamax;
4922 if (slam < xvmin - slamax) slam = xvmin - slamax;
4927 else if (slam < undral)
4934 for (ipt = 1; ipt <= 3; ++ipt) {
4935 if (
TMath::Abs(slam - xvals[ipt-1]) < toler9) {
4946 for (i = 1; i <=
fNpar; ++i) {
fX[i-1] = start[i-1] + slam*step[i-1]; }
4950 chpq[nxypt-1] = charal[nxypt-1];
4951 xpq[nxypt-1] = slam;
4956 if (fvals[1] > fvmax) {
4960 if (fvals[2] > fvmax) {
4970 if (slam > xvmin) overal =
TMath::Min(overal,slam - toler8);
4971 if (slam < xvmin) undral =
TMath::Max(undral,slam + toler8);
4972 slam = (slam + xvmin)*.5;
4974 }
while (f3 >= fvmax);
4977 if (l65 || l70)
break;
4979 xvals[nvmax-1] = slam;
4980 fvals[nvmax-1] =
f3;
4985 if (slam > xvmin) overal =
TMath::Min(overal,slam - toler8);
4986 if (slam < xvmin) undral =
TMath::Max(undral,slam + toler8);
4988 }
while (nxypt < 12);
4994 cmess =
" LINE SEARCH HAS EXHAUSTED THE LIMIT OF FUNCTION CALLS ";
4996 Printf(
" MNLINE DEBUG: steps=");
4997 for (kk = 1; kk <=
fNpar; ++kk) {
4998 Printf(
" %12.4g",step[kk-1]);
5003 if (l70) cmess =
" LINE SEARCH HAS ATTAINED TOLERANCE ";
5004 if (l80) cmess =
" STEP SIZE AT ARITHMETICALLY ALLOWED MINIMUM";
5007 for (i = 1; i <=
fNpar; ++i) {
5008 fDirin[i-1] = step[i-1]*xvmin;
5009 fX[i-1] = start[i-1] +
fDirin[i-1];
5013 mnwarn(
"D",
"MNLINE",
" LINE MINIMUM IN BACKWARDS DIRECTION");
5015 if (fvmin == fstart) {
5016 mnwarn(
"D",
"MNLINE",
" LINE SEARCH FINDS NO IMPROVEMENT ");
5019 Printf(
" AFTER %3d POINTS,%s",nxypt,(
const char*)cmess);
5034 Int_t ndex, i, j,
m,
n, ncoef, nparm,
id, it, ix;
5035 Int_t nsofar, ndi, ndj, iso, isw2, isw5;
5044 Printf(
" MNMATU: NPAR=0");
5058 if (
fNpar <= 1)
return;
5064 Printf(
" PARAMETER CORRELATION COEFFICIENTS ");
5065 ctemp =
" NO. GLOBAL";
5066 for (
id = 1;
id <= nparm; ++
id) {
5069 Printf(
"%s",(
const char*)ctemp);
5070 for (i = 1; i <=
fNpar; ++i) {
5072 ndi = i*(i + 1) / 2;
5073 for (j = 1; j <=
fNpar; ++j) {
5076 ndex = m*(m-1) / 2 + n;
5077 ndj = j*(j + 1) / 2;
5082 for (it = 1; it <= nparm; ++it) {
5085 Printf(
"%s",(
const char*)ctemp);
5086 if (i <= nparm)
continue;
5088 for (iso = 1; iso <= 10; ++iso) {
5091 for (it = nsofar + 1; it <= nparm; ++it) {
5094 Printf(
"%s",(
const char*)ctemp);
5095 if (i <= nparm)
break;
5114 Double_t gdel, gami, vlen, dsum, gssq, vsum,
d;
5115 Double_t fzero, fs, ri, delgam, rhotol;
5117 Int_t npfn, ndex, iext, i, j,
m,
n, npsdf, nparx;
5118 Int_t iswtr, lined2, kk, nfcnmg, nrstrt,
iter;
5122 if (
fNpar <= 0)
return;
5137 rhotol =
fApsi*.001;
5139 Printf(
" START MIGRAD MINIMIZATION. STRATEGY %2d. CONVERGENCE WHEN EDM .LT.%9.2e",
fIstrat,rhotol);
5154 if (
fISW[1] >= 1)
goto L10;
5162 if (
fISW[1] >= 1)
goto L10;
5164 for (i = 1; i <=
fNpar; ++i) {
5171 for (i = 1; i <=
fNpar; ++i) {
5172 if (
fG2[i-1] > 0)
continue;
5178 mnwarn(
"D",
"MNMIGR",
"Negative G2 line search");
5181 Printf(
" Negative G2 line search, param %3d %13.3g%13.3g",iext,fs,
fAmin);
5187 for (i = 1; i <=
fNpar; ++i) {
5189 for (j = 1; j <= i-1; ++j) {
5194 if (
fG2[i-1] <= 0)
fG2[i-1] = 1;
5199 Printf(
" DEBUG MNMIGR, STARTING MATRIX DIAGONAL, VHMAT=");
5200 for (kk = 1; kk <=
Int_t(vlen); ++kk) {
5214 for (i = 1; i <=
fNpar; ++i) {
5218 for (j = 1; j <= i-1; ++j) {
5227 mnwarn(
"W",
"MIGRAD",
"STARTING MATRIX NOT POS-DEFINITE.");
5237 if (iswtr >= 2)
mnmatu(0);
5243 for (i = 1; i <=
fNpar; ++i) {
5246 for (j = 1; j <=
fNpar; ++j) {
5249 ndex = m*(m-1) / 2 + n;
5250 ri +=
fVhmat[ndex-1]*fMIGRgs[j-1];
5256 mnwarn(
"D",
"MIGRAD",
" FIRST DERIVATIVES OF FCN ARE ALL ZERO");
5261 mnwarn(
"D",
"MIGRAD",
" NEWTON STEP NOT DESCENT.");
5262 if (npsdf == 1)
goto L1;
5269 if (
fAmin == fs)
goto L200;
5286 for (i = 1; i <=
fNpar; ++i) {
5289 for (j = 1; j <=
fNpar; ++j) {
5292 ndex = m*(m-1) / 2 + n;
5298 gdgssq += gami*gami;
5300 delgam +=
fDirin[i-1]*gami;
5305 if (
fEDM < 0 || gvg <= 0) {
5306 mnwarn(
"D",
"MIGRAD",
"NOT POS-DEF. EDM OR GVG NEGATIVE.");
5308 if (npsdf == 1)
goto L230;
5315 if (iswtr >= 3 || (iswtr == 2 && iter % 10 == 1)) {
5320 mnwarn(
"D",
"MIGRAD",
"NO CHANGE IN FIRST DERIVATIVES OVER LAST STEP");
5323 mnwarn(
"D",
"MIGRAD",
"FIRST DERIVATIVES INCREASING ALONG SEARCH LINE");
5329 for (kk = 1; kk <= 10; ++kk) {
5335 for (i = 1; i <=
fNpar; ++i) {
5336 for (j = 1; j <= i; ++j) {
5337 if(delgam == 0 || gvg == 0) d = 0;
5340 ndex = i*(i-1) / 2 + j;
5347 if (iswtr >= 3 || ldebug) {
5348 Printf(
" RELATIVE CHANGE IN COV. MATRIX=%5.1f per cent",
fDcovar*100);
5352 for (kk = 1; kk <= 10; ++kk) {
5356 if (delgam <= gvg)
goto L135;
5357 for (i = 1; i <=
fNpar; ++i) {
5360 for (i = 1; i <=
fNpar; ++i) {
5361 for (j = 1; j <= i; ++j) {
5362 ndex = i*(i-1) / 2 + j;
5368 if (
fEDM < rhotol*.1)
goto L300;
5370 for (i = 1; i <=
fNpar; ++i) {
5384 Printf(
" CALL LIMIT EXCEEDED IN MIGRAD.");
5391 Printf(
" MIGRAD FAILS TO FIND IMPROVEMENT");
5394 if (
fEDM < rhotol)
goto L300;
5397 Printf(
" MACHINE ACCURACY LIMITS FURTHER IMPROVEMENT.");
5403 Printf(
" MIGRAD FAILS WITH STRATEGY=0. WILL TRY WITH STRATEGY=1.");
5411 Printf(
" MIGRAD TERMINATED WITHOUT CONVERGENCE.");
5419 Printf(
" MIGRAD MINIMIZATION HAS CONVERGED.");
5424 Printf(
" MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.");
5429 if (
fEDM > rhotol)
goto L10;
5440 if (iswtr >= 0)
mnprin(3, fAmin);
5441 if (iswtr >= 1)
mnmatu(1);
5456 Int_t nbad, ilax, ilax2, ngood, nfcnmi, iin, knt;
5458 if (
fNpar <= 0)
goto L700;
5463 for (knt = 1; knt <=
fNpar; ++knt) {
5467 if (knt >= 7)
break;
5469 if (ilax == 0)
break;
5470 if (ilax > 0 && ilax <=
fNu) {
5471 if (
fNiofex[ilax-1] > 0)
goto L565;
5473 Printf(
" PARAMETER NUMBER %3d NOT A VARIABLE. IGNORED.",ilax);
5479 mnmnot(ilax, ilax2, val2pl, val2mi);
5483 if (
fErp[iin-1] > 0) ++ngood;
5485 if (
fErn[iin-1] < 0) ++ngood;
5493 if (ngood == 0 && nbad == 0)
goto L700;
5494 if (ngood > 0 && nbad == 0)
fCstatu =
"SUCCESSFUL";
5495 if (ngood == 0 && nbad > 0)
fCstatu =
"FAILURE ";
5496 if (ngood > 0 && nbad > 0)
fCstatu =
"PROBLEMS ";
5506 Printf(
" NEW MINIMUM FOUND. GO BACK TO MINIMIZATION STEP.");
5507 Printf(
" =================================================");
5517 Printf(
" THERE ARE NO MINOS ERRORS TO CALCULATE.");
5535 Double_t abest, xunit, dc, ut, sigsav, du1;
5537 Int_t marc, isig, mpar, ndex, imax, indx, ierr, i, j;
5538 Int_t iercr, it, istrav, nfmxin, nlimit, isw2, isw4;
5552 for (i = 1; i <= mpar; ++i) {
fXt[i-1] =
fX[i-1]; }
5553 i__1 = mpar*(mpar + 1) / 2;
5555 for (i = 1; i <= mpar; ++i) {
5564 if (
fNvarl[ilax-1] == 1) {
5566 fBlim[ilax-1] = ut + fMNOTw[it-1]*100;
5568 ndex = it*(it + 1) / 2;
5571 for (i = 1; i <= mpar; ++i) {
5572 if (i == it)
continue;
5581 Printf(
" MINUIT ERROR. CANNOT FIX PARAMETER %4d INTERNAL %3d",ilax,it);
5587 for (isig = 1; isig <= 2; ++isig) {
5597 Printf(
" DETERMINATION OF %sTIVE MINOS ERROR FOR PARAMETER %d %s"
5598 ,(
const char*)csig,ilax
5599 ,(
const char*)
fCpnam[ilax-1]);
5602 mnwarn(
"D",
"MINOS",
"NO COVARIANCE MATRIX.");
5604 nlimit =
fNfcn + nfmxin;
5607 fU[ilax-1] = ut + sig*du1;
5610 delu =
fU[ilax-1] - ut;
5613 fac = delu /
fMNOTw[it-1];
5614 for (i = 1; i <=
fNpar; ++i) {
5618 Printf(
" PARAMETER %4d SET TO%11.3e + %10.3e = %12.3e",ilax,ut,delu,
fU[ilax-1]);
5629 if (abest -
fAmin >
fUp*.01)
goto L650;
5630 if (iercr == 1)
goto L440;
5631 if (iercr == 2)
goto L450;
5632 if (iercr == 3)
goto L460;
5636 Printf(
" THE %4sTIVE MINOS ERROR OF PARAMETER %3d %10s, IS %12.4e"
5637 ,(
const char*)csig,ilax
5638 ,(
const char*)
fCpnam[ilax-1],eros);
5644 Printf(
" THE %4sTIVE MINOS ERROR OF PARAMETER %3d, %s EXCEEDS ITS LIMIT."
5645 ,(
const char*)csig,ilax
5646 ,(
const char*)
fCpnam[ilax-1]);
5652 Printf(
" THE %4sTIVE MINOS ERROR %4d REQUIRES MORE THAN %5d FUNCTION CALLS."
5653 ,(
const char*)csig,ilax,nfmxin);
5659 Printf(
" %4sTIVE MINOS ERROR NOT CALCULATED FOR PARAMETER %d"
5660 ,(
const char*)csig,ilax);
5666 Printf(
" **************************************************************************");
5670 if (ilax2 > 0 && ilax2 <=
fNu) val2mi =
fU[ilax2-1];
5673 if (ilax2 > 0 && ilax2 <=
fNu) val2pl =
fU[ilax2-1];
5680 i__1 = mpar*(mpar + 1) / 2;
5682 for (i = 1; i <= mpar; ++i) {
5731 Double_t vplu, a_small, gsmin, pinti, vminu, danger, sav, sav2;
5732 Int_t ierr,
kint, in, ix, ktofix, lastin, kinfix, nvl;
5740 Printf(
" MINUIT USER ERROR. PARAMETER NUMBER IS %3d ALLOWED RANGE IS ONE TO %4d",k,fMaxext);
5745 if (
fNvarl[k-1] < 0)
goto L50;
5748 for (ix = 1; ix <=
fNpfix; ++ix) {
5749 if (
fIpfix[ix-1] == k) ktofix = k;
5752 mnwarn(
"W",
"PARAM DEF",
"REDEFINING A FIXED PARAMETER.");
5754 Printf(
" CANNOT RELEASE. MAX NPAR EXCEEDED.");
5765 Printf(
" PARAMETER DEFINITIONS:");
5766 Printf(
" NO. NAME VALUE STEP SIZE LIMITS");
5769 if (wk > 0)
goto L122;
5772 Printf(
" %5d %-10s %13.5e constant",k,(
const char*)cnamk,uk);
5777 if (a == 0 && b == 0) {
5781 Printf(
" %5d %-10s %13.5e%13.5e no limits",k,(
const char*)cnamk,uk,wk);
5788 Printf(
" %5d %-10s %13.5e%13.5e %13.5e%13.5e",k,(
const char*)cnamk,uk,wk,a,b);
5794 Printf(
" MINUIT USER ERROR. TOO MANY VARIABLE PARAMETERS.");
5797 if (nvl == 1)
goto L200;
5799 Printf(
" USER ERROR IN MINUIT PARAMETER");
5801 Printf(
" UPPER AND LOWER LIMITS EQUAL.");
5808 mnwarn(
"W",
"PARAM DEF",
"PARAMETER LIMITS WERE REVERSED.");
5815 danger = (b - uk)*(uk - a);
5817 mnwarn(
"W",
"PARAM DEF",
"STARTING VALUE OUTSIDE LIMITS.");
5820 mnwarn(
"W",
"PARAM DEF",
"STARTING VALUE IS AT LIMIT.");
5838 for (ix = 1; ix <= k-1; ++ix) {
if (
fNiofex[ix-1] > 0) ++lastin; }
5840 if (kint ==
fNpar)
goto L280;
5843 for (in =
fNpar; in >= lastin + 1; --in) {
5857 for (in = lastin + 1; in <=
kint; ++in) {
5880 mnpint(sav, ix-1, pinti);
5882 fXt[in-1] =
fX[in-1];
5885 mnpint(sav2, ix-1, pinti);
5886 vplu = pinti -
fX[in-1];
5888 mnpint(sav2, ix-1, pinti);
5889 vminu = pinti - fX[in-1];
5908 if (kinfix > 0)
mnfixp(kinfix-1, ierr);
5909 if (ierr > 0)
goto L800;
5934 Double_t a=0, b=0, fk=0, uk=0, wk=0, xk=0;
5935 Int_t ierr, kapo1, kapo2;
5936 Int_t k, llist, ibegin, lenbuf, istart, lnc, icy;
5937 TString cnamk, comand, celmnt, ctemp;
5940 lenbuf = strlen((
const char*)crdbuf);
5942 kapo1 = strspn((
const char*)crdbuf,
"'");
5943 if (kapo1 == 0)
goto L150;
5944 kapo2 = strspn((
const char*)crdbuf + kapo1,
"'");
5945 if (kapo2 == 0)
goto L150;
5949 for (istart = 1; istart <= kapo1-1; ++istart) {
5950 if (crdbuf(istart-1,1) !=
' ')
goto L120;
5955 celmnt = crdbuf(istart-1, kapo1-istart);
5956 if (scanf((
const char*)celmnt,&fk)) {;}
5958 if (k <= 0)
goto L210;
5959 cnamk =
"PARAM " + celmnt;
5960 if (kapo2 - kapo1 > 1) {
5961 cnamk = crdbuf(kapo1, kapo2-1-kapo1);
5964 for (icy = kapo2 + 1; icy <= lenbuf; ++icy) {
5965 if (crdbuf(icy-1,1) ==
',')
goto L139;
5966 if (crdbuf(icy-1,1) !=
' ')
goto L140;
5977 ctemp = crdbuf(ibegin-1,lenbuf-ibegin);
5979 if (ierr > 0)
goto L180;
5990 if (scanf((
const char*)crdbuf,&xk,stmp,&uk,&wk,&a,&b)) {;}
5993 if (k == 0)
goto L210;
5996 mnparm(k-1, cnamk, uk, wk, a, b, ierr);
6024 Double_t a,
f, s,
t,
y,
s2,
x2,
x3,
x4, y2, cz[3], xm,
xy, x2y;
6034 for (i = 1; i <= 3; ++i) { cz[i-1] = 0; }
6036 if (npar2p < 3)
goto L10;
6040 for (i = 1; i <= npar2p; ++i) { xm += parx2p[i]; }
6049 for (i = 1; i <= npar2p; ++i) {
6061 a = (f*x4 - x2*
x2)*x2 - f*(x3*x3);
6062 if (a == 0)
goto L10;
6063 cz[2] = (x2*(f*x2y - x2*
y) - f*x3*xy) /
a;
6064 cz[1] = (xy - x3*cz[2]) / x2;
6065 cz[0] = (y - x2*cz[2]) / f;
6066 if (npar2p == 3)
goto L6;
6067 sdev2p = y2 - (cz[0]*y + cz[1]*xy + cz[2]*x2y);
6068 if (sdev2p < 0) sdev2p = 0;
6071 cz[0] += xm*(xm*cz[2] - cz[1]);
6072 cz[1] -= xm*2*cz[2];
6074 for (i = 1; i <= 3; ++i) { coef2p[i] = cz[i-1]; }
6097 yy = (pexti - alimi)*2 / (blimi - alimi) - 1;
6102 chbuf2 =
" IS AT ITS LOWER ALLOWED LIMIT.";
6105 chbuf2 =
" IS AT ITS UPPER ALLOWED LIMIT.";
6108 pexti = alimi + (blimi - alimi)*.5*(
TMath::Sin(a) + 1);
6110 if (yy2 > 1) chbuf2 =
" BROUGHT BACK INSIDE LIMITS.";
6140 if ((h =
gROOT->GetPluginManager()->FindHandler(
"TMinuitGraph"))) {
6153 Double_t bwidx, bwidy, xbest, ybest, ax, ay, bx, by;
6154 Double_t xvalus[12], any, dxx, dyy;
6155 Int_t iten, i, j, k, maxnx, maxny, iquit, ni, linodd;
6156 Int_t nxbest, nybest, km1, ibk, isp1,
nx,
ny, ks, ix;
6165 if (maxnx < 10) maxnx = 10;
6167 if (maxny < 10) maxny = 10;
6168 if (nxypt <= 1)
return;
6174 for (i = 1; i <= km1; ++i) {
6177 for (j = 1; j <= ni; ++j) {
6178 if (ypt[j-1] > ypt[j])
continue;
6190 if (iquit == 0)
break;
6195 for (i = 1; i <= nxypt; ++i) {
6196 if (xpt[i-1] > xmax) xmax = xpt[i-1];
6197 if (xpt[i-1] < xmin) xmin = xpt[i-1];
6199 dxx = (xmax -
xmin)*.001;
6202 mnbins(xmin, xmax, maxnx, xmin, xmax, nx, bwidx);
6204 ymin = ypt[nxypt-1];
6205 if (ymax == ymin) ymax = ymin + 1;
6206 dyy = (ymax -
ymin)*.001;
6209 mnbins(ymin, ymax, maxny, ymin, ymax, ny, bwidy);
6212 if (chbest ==
' ')
goto L50;
6213 xbest = (xmax +
xmin)*.5;
6214 ybest = (ymax +
ymin)*.5;
6222 for (i = 1; i <= nxypt; ++i) {
6223 xpt[i-1] = ax*xpt[i-1] + bx;
6224 ypt[i-1] = any - ay*ypt[i-1] - by;
6226 nxbest =
Int_t((ax*xbest + bx));
6227 nybest =
Int_t((any - ay*ybest - by));
6234 for (i = 1; i <=
ny; ++i) {
6235 for (ibk = 1; ibk <=
nx; ++ibk) { cline[ibk-1] =
' '; }
6240 if (nx>0) cline[nx-1] =
'.';
6241 cline[nxbest-1] =
'.';
6242 if (i != 1 && i != nybest && i != ny)
goto L320;
6243 for (j = 1; j <=
nx; ++j) { cline[j-1] =
'.'; }
6246 if (isp1 > nxypt)
goto L350;
6248 for (k = isp1; k <= nxypt; ++k) {
6249 ks =
Int_t(ypt[k-1]);
6250 if (ks > i)
goto L345;
6251 ix =
Int_t(xpt[k-1]);
6252 if (cline[ix-1] ==
'.')
goto L340;
6253 if (cline[ix-1] ==
' ')
goto L340;
6254 if (cline[ix-1] == chpt[k-1])
continue;
6261 cline[ix-1] = chpt[k-1];
6268 if (linodd == 1 || i == ny)
goto L380;
6271 Printf(
" %s",(
const char*)ctemp);
6275 Printf(
" %14.7g ..%s",yprt,(
const char*)ctemp);
6281 for (ibk = 1; ibk <=
nx; ++ibk) {
6283 if (ibk % 10 == 1) cline[ibk-1] =
'/';
6287 for (ibk = 1; ibk <= 12; ++ibk) {
6288 xvalus[ibk-1] = xmin +
Double_t(ibk-1)*10*bwidx;
6290 iten = (nx + 9) / 10;
6291 for (ibk = 1; ibk <= iten && ibk <= 12; ++ibk) {
6292 snprintf(cline + (ibk-1)*10, 11-(ibk == 12),
"%#9.3g ", xvalus[ibk-1]);
6296 if (overpr) chmess =
" Overprint character is &";
6297 Printf(
" ONE COLUMN=%13.7g%s",bwidx,(
const char*)chmess);
6320 Int_t iint, iext, nvl;
6322 Int_t iuext = iuext1 + 1;
6326 if (iuext == 0)
goto L100;
6330 if (iint >
fNpar)
goto L100;
6336 if (iext >
fNu)
goto L100;
6342 if (nvl < 0)
goto L100;
6345 if (iint > 0) err =
fWerr[iint-1];
6347 xlolim =
fAlim[iext-1];
6348 xuplim =
fBlim[iext-1];
6354 chnam =
"undefined";
6383 Int_t nadd, i, k,
l,
m, ikode, ic, nc, ntrail, lbl;
6385 TString colhdl[6], colhdu[6], cx2, cx3, cheval;
6388 Printf(
" THERE ARE CURRENTLY NO PARAMETERS DEFINED");
6394 ikode =
fISW[1] + 1;
6395 if (ikode > 3) ikode = 3;
6398 for (k = 1; k <= 6; ++k) {
6399 colhdu[k-1] =
"UNDEFINED";
6400 colhdl[k-1] =
"COLUMN HEAD";
6407 if (fval ==
fUndefi) cheval =
" unknown ";
6408 else cheval.
Form(
"%g",fval);
6414 Printf(
" FCN=%s FROM %8s STATUS=%10s %6d CALLS %9d TOTAL"
6415 ,(
const char*)cheval
6419 if (m == 0 || m == 2 ||
fDcovar == 0) {
6420 Printf(
" EDM=%s STRATEGY=%2d %s"
6426 Printf(
" EDM=%s STRATEGY=%2d ERROR MATRIX UNCERTAINTY %5.1f per cent"
6427 ,(
const char*)chedm,
fIstrat,dc);
6430 if (ikode == 0)
return;
6433 for (i = 1; i <=
fNu; ++i) {
6434 if (
fNvarl[i-1] < 0)
continue;
6435 for (ic = 10; ic >= 1; --ic) {
6436 if (
fCpnam[i-1](ic-1,1) !=
" ")
goto L16;
6441 if (lbl < ntrail) ntrail = lbl;
6443 nadd = ntrail / 2 + 1;
6446 colhdl[0] =
" ERROR ";
6447 colhdu[1] =
" PHYSICAL";
6448 colhdu[2] =
" LIMITS ";
6449 colhdl[1] =
" NEGATIVE ";
6450 colhdl[2] =
" POSITIVE ";
6454 colhdl[0] =
" ERROR ";
6455 colhdu[1] =
" INTERNAL ";
6456 colhdl[1] =
" STEP SIZE ";
6457 colhdu[2] =
" INTERNAL ";
6458 colhdl[2] =
" VALUE ";
6462 colhdl[0] =
" ERROR ";
6463 colhdu[1] =
" STEP ";
6464 colhdl[1] =
" SIZE ";
6465 colhdu[2] =
" FIRST ";
6466 colhdl[2] =
" DERIVATIVE ";
6469 colhdu[0] =
" PARABOLIC ";
6470 colhdl[0] =
" ERROR ";
6471 colhdu[1] =
" MINOS ";
6472 colhdu[2] =
"ERRORS ";
6473 colhdl[1] =
" NEGATIVE ";
6474 colhdl[2] =
" POSITIVE ";
6478 if (
fISW[1] < 3) colhdu[0] =
" APPROXIMATE ";
6479 if (
fISW[1] < 1) colhdu[0] =
" CURRENT GUESS";
6481 Printf(
" EXT PARAMETER %-14s%-14s%-14s",(
const char*)colhdu[0]
6482 ,(
const char*)colhdu[1]
6483 ,(
const char*)colhdu[2]);
6484 Printf(
" NO. NAME VALUE %-14s%-14s%-14s",(
const char*)colhdl[0]
6485 ,(
const char*)colhdl[1]
6486 ,(
const char*)colhdl[2]);
6488 for (i = 1; i <=
fNu; ++i) {
6489 if (
fNvarl[i-1] < 0)
continue;
6491 cnambf = cblank(0,nadd) +
fCpnam[i-1];
6492 if (l == 0)
goto L55;
6495 cx2 =
"PLEASE GET X..";
6496 cx3 =
"PLEASE GET X..";
6499 Printf(
"%4d %-11s%14.5e%14.5e",i,(
const char*)cnambf,
fU[i-1],x1);
6514 cx3 =
"** at limit **";
6519 if (x2 == 0) cx2 =
" ";
6520 if (x2 ==
fUndefi) cx2 =
" at limit ";
6522 if (x3 == 0) cx3 =
" ";
6523 if (x3 ==
fUndefi) cx3 =
" at limit ";
6525 if (cx2 ==
"PLEASE GET X..") cx2.
Form(
"%14.5e",x2);
6526 if (cx3 ==
"PLEASE GET X..") cx3.
Form(
"%14.5e",x3);
6527 Printf(
"%4d %-11s%14.5e%14.5e%-14s%-14s",i
6528 ,(
const char*)cnambf,
fU[i-1],x1
6529 ,(
const char*)cx2,(
const char*)cx3);
6532 if (
fNvarl[i-1] <= 1 || ikode == 3)
continue;
6534 Printf(
" WARNING - - ABOVE PARAMETER IS AT LIMIT.");
6540 colhdu[0] =
" constant ";
6541 if (
fNvarl[i-1] > 0) colhdu[0] =
" fixed ";
6542 if (
fNvarl[i-1] == 4 && ikode == 1) {
6543 Printf(
"%4d %-11s%14.5e%-14s%14.5e%14.5e",i
6544 ,(
const char*)cnambf,
fU[i-1]
6545 ,(
const char*)colhdu[0],
fAlim[i-1],
fBlim[i-1]);
6547 Printf(
"%4d %-11s%14.5e%s",i
6548 ,(
const char*)cnambf,
fU[i-1],(
const char*)colhdu[0]);
6567 Double_t dgmin, padd, pmin, pmax, dg, epspdf, epsmin;
6568 Int_t ndex, i, j, ndexd, ip, ifault;
6575 for (i = 1; i <=
fNpar; ++i) {
6576 ndex = i*(i + 1) / 2;
6577 if (
fVhmat[ndex-1] <= 0) {
6583 dg = epspdf + 1 - dgmin;
6589 for (i = 1; i <=
fNpar; ++i) {
6593 if (
fVhmat[ndexd-1]==0) {
6598 for (j = 1; j <= i; ++j) {
6607 for (ip = 2; ip <=
fNpar; ++ip) {
6613 Printf(
" EIGENVALUES OF SECOND-DERIVATIVE MATRIX:");
6615 for (ip = 1; ip <=
fNpar; ++ip) {
6620 if (pmin > epspdf*pmax)
return;
6622 padd = pmax*.001 - pmin;
6623 for (ip = 1; ip <=
fNpar; ++ip) {
6624 ndex = ip*(ip + 1) / 2;
6625 fVhmat[ndex-1] *= padd + 1;
6628 mnwarn(
"W",
fCfrom,
Form(
"MATRIX FORCED POS-DEF BY ADDING %f TO DIAGONAL.",padd));
6649 for (i = 1; i <=
fNpar; ++i) {
fX[i-1] = pnew[i-1]; }
6657 for (j = 2; j <= nparp1; ++j) {
if (y[j-1] > y[jh-1]) jh = j; }
6658 fEDM = y[jh-1] - y[jl-1];
6659 if (
fEDM <= 0)
goto L45;
6660 for (i = 1; i <=
fNpar; ++i) {
6663 for (j = 2; j <= nparp1; ++j) {
6667 fDirin[i-1] = pbig - plit;
6672 Printf(
" FUNCTION VALUE DOES NOT SEEM TO DEPEND ON ANY OF THE %d VARIABLE PARAMETERS.",
fNpar);
6673 Printf(
" VERIFY THAT STEP SIZES ARE BIG ENOUGH AND CHECK FCN LOGIC.");
6674 Printf(
" *******************************************************************************");
6675 Printf(
" *******************************************************************************");
6692 static Int_t iseed = 12345;
6696 if (val == 3)
goto L100;
6699 iseed = (iseed - k*53668)*40014 - k*12211;
6700 if (iseed < 0) iseed += 2147483563;
6701 val =
Double_t(iseed*4.656613e-10);
6732 for (i = 1; i <=
fNpar; ++i) {
6754 Printf(
"mnsave is dummy in TMinuit");
6769 Double_t step, uhigh, xhreq, xlreq, ubest, fnext, unext, xh, xl;
6771 Int_t nxypt, nccall, iparwd;
6776 if (ncall <= 1) ncall = 41;
6777 if (ncall > 98) ncall = 98;
6783 if (iparwd > 0)
goto L200;
6788 if (ipar >
fNu)
goto L900;
6790 if (iint <= 0)
goto L100;
6802 if (
fNvarl[ipar-1] > 1)
goto L300;
6805 if (xlreq == xhreq)
goto L250;
6807 step = (xhreq - xlreq) /
Double_t(ncall-1);
6810 xl = ubest -
fWerr[iint-1];
6811 xh = ubest + fWerr[iint-1];
6812 mnbins(xl, xh, ncall, unext, uhigh, nbins, step);
6817 if (xlreq == xhreq)
goto L350;
6822 if (xl >= xh)
goto L700;
6824 step = (xh - xl) /
Double_t(ncall-1);
6827 unext =
fAlim[ipar-1];
6831 for (icall = 1; icall <= nccall; ++icall) {
6836 fXpt[nxypt-1] = unext;
6837 fYpt[nxypt-1] = fnext;
6838 fChpt[nxypt-1] =
'*';
6839 if (fnext <
fAmin) {
6852 Printf(
"%dSCAN OF PARAMETER NO. %d, %s"
6857 Printf(
" REQUESTED RANGE OUTSIDE LIMITS FOR PARAMETER %d",ipar);
6859 if (iparwd <= 0)
goto L100;
6880 Double_t dxdi, rnum, ftry, rnum1, rnum2, alpha;
6882 Int_t ipar, iext, j, ifail, iseed=0, nparx, istep, ib, mxfail, mxstep;
6885 if (mxfail <= 0) mxfail =
fNpar*20 + 100;
6889 if (alpha <= 0) alpha = 3;
6891 Printf(
" MNSEEK: MONTE CARLO MINIMIZATION USING METROPOLIS ALGORITHM");
6892 Printf(
" TO STOP AFTER %6d SUCCESSIVE FAILURES, OR %7d STEPS",mxfail,mxstep);
6893 Printf(
" MAXIMUM STEP SIZE IS %9.3f ERROR BARS.",alpha);
6905 for (ipar = 1; ipar <=
fNpar; ++ipar) {
6908 if (
fNvarl[iext-1] > 1) {
6911 if (dxdi == 0) dxdi = 1;
6912 fDirin[ipar-1] = alpha*2*fWerr[ipar-1] / dxdi;
6914 fDirin[ipar-1] = 6.2831859999999997;
6921 for (istep = 1; istep <= mxstep; ++istep) {
6922 if (ifail >= mxfail)
break;
6923 for (ipar = 1; ipar <=
fNpar; ++ipar) {
6953 Printf(
" MNSEEK: %5d SUCCESSIVE UNSUCCESSFUL TRIALS.",ifail);
6977 static const char *cname[30] = {
7009 static Int_t nname = 25;
7010 static Int_t nntot = 30;
7012 "-1: NO OUTPUT EXCEPT FROM SHOW ",
7013 " 0: REDUCED OUTPUT ",
7014 " 1: NORMAL OUTPUT ",
7015 " 2: EXTRA OUTPUT FOR PROBLEM CASES",
7016 " 3: MAXIMUM OUTPUT "};
7019 " 0: MINIMIZE THE NUMBER OF CALLS TO FUNCTION",
7020 " 1: TRY TO BALANCE SPEED AGAINST RELIABILITY",
7021 " 2: MAKE SURE MINIMUM TRUE, ERRORS CORRECT "};
7024 "REPORT ALL EXCEPTIONAL CONDITIONS ",
7025 "MNLINE: LINE SEARCH MINIMIZATION ",
7026 "MNDERI: FIRST DERIVATIVE CALCULATIONS ",
7027 "MNHESS: SECOND DERIVATIVE CALCULATIONS ",
7028 "MNMIGR: COVARIANCE MATRIX UPDATES ",
7029 "MNHES1: FIRST DERIVATIVE UNCERTAINTIES ",
7030 "MNCONT: MNCONTOUR PLOT (MNCROS SEARCH) "};
7037 Int_t iset, iprm, i, jseed, kname, iseed, iunit,
id,
ii, kk;
7038 Int_t ikseed, idbopt, igrain=0, iswsav, isw2;
7039 TString cfname, cmode, ckind, cwarn, copt, ctemp, ctemp2;
7042 for (i = 1; i <= nntot; ++i) {
7046 if (strstr(ctemp2.
Data(),ckind.
Data()))
goto L5;
7054 if ( ctemp2.
Contains(
"HEL"))
goto L2000;
7055 if ( ctemp2.
Contains(
"SHO"))
goto L1000;
7056 if (!ctemp2.
Contains(
"SET"))
goto L1900;
7060 if (kname <= 0)
goto L1900;
7062 switch ((
int)kname) {
7080 case 18:
goto L3000;
7082 case 20:
goto L3000;
7087 case 25:
goto L3000;
7088 case 26:
goto L1900;
7098 if (iprm >
fNu)
goto L25;
7099 if (iprm <= 0)
goto L25;
7100 if (
fNvarl[iprm-1] < 0)
goto L25;
7112 Printf(
" UNDEFINED PARAMETER NUMBER. IGNORED.");
7143 for (i = 1; i <=
fNpar; ++i) {
7171 mnwarn(
"W",
"SHO",
"SHO");
7179 Printf(
" MINUIT RANDOM NUMBER SEED SET TO %d",jseed);
7191 if (
fISW[4] > 0)
goto L1172;
7209 if (
fISW[4] >= 0)
goto L1220;
7214 if (
fISW[4] >= 0)
goto L1100;
7219 if (
fISW[4] >= 0)
goto L1100;
7230 if (idbopt > 6)
goto L288;
7232 fIdbg[idbopt] = iset;
7233 if (iset == 1)
fIdbg[0] = 1;
7236 for (
id = 0;
id <= 6; ++
id) {
fIdbg[
id] = iset; }
7239 mnwarn(
"D",
"SHO",
"SHO");
7242 Printf(
" UNKNOWN DEBUG OPTION %d REQUESTED. IGNORED",idbopt);
7253 if (kname <= 0)
goto L1900;
7255 switch ((
int)kname) {
7265 case 10:
goto L1100;
7266 case 11:
goto L1110;
7267 case 12:
goto L1120;
7268 case 13:
goto L1130;
7269 case 14:
goto L1130;
7270 case 15:
goto L1150;
7271 case 16:
goto L1160;
7272 case 17:
goto L1170;
7273 case 18:
goto L1180;
7274 case 19:
goto L1190;
7275 case 20:
goto L1200;
7276 case 21:
goto L1210;
7277 case 22:
goto L1220;
7278 case 23:
goto L1100;
7279 case 24:
goto L1100;
7280 case 25:
goto L1250;
7281 case 26:
goto L1900;
7282 case 27:
goto L1270;
7283 case 28:
goto L1270;
7284 case 29:
goto L1290;
7285 case 30:
goto L1300;
7315 Printf(
" ALLOWED PRINT LEVELS ARE:");
7321 Printf(
" CURRENT PRINTOUT LEVEL IS %s",cprlev[
fISW[4]+1].
Data());
7326 Printf(
" NOGRAD IS SET. DERIVATIVES NOT COMPUTED IN FCN.");
7328 Printf(
" GRAD IS SET. USER COMPUTES DERIVATIVES IN FCN.");
7333 Printf(
" ERRORS CORRESPOND TO FUNCTION CHANGE OF %g",
fUp);
7357 cmode =
"BATCH MODE ";
7358 if (
fISW[5] == 1) cmode =
"INTERACTIVE MODE";
7359 if (! lname) cfname =
"unknown";
7360 Printf(
" INPUT NOW BEING READ IN %s FROM UNIT NO. %d FILENAME: %s"
7361 ,(
const char*)cmode,
fIsysrd,(
const char*)cfname);
7373 cwarn =
"SUPPRESSED";
7374 if (
fLwarn) cwarn =
"REPORTED ";
7375 Printf(
"%s",(
const char*)cwarn);
7383 Printf(
" MINUIT RNDM SEED IS CURRENTLY=%d",ikseed);
7390 Printf(
" TITLE OF CURRENT TASK IS:%s",(
const char*)
fCtitl);
7394 Printf(
" ALLOWED STRATEGIES ARE:");
7399 Printf(
" NOW USING STRATEGY %s",(
const char*)cstrat[
fIstrat]);
7416 Printf(
" NO PAGE THROWS IN MINUIT OUTPUT");
7421 for (ii = 1; ii <=
fNpar; ++
ii) {
7422 if (
fErp[ii-1] > 0 ||
fErn[ii-1] < 0)
goto L1204;
7424 Printf(
" THERE ARE NO MINOS ERRORS CURRENTLY VALID.");
7431 Printf(
" FLOATING-POINT NUMBERS ASSUMED ACCURATE TO %g",
fEpsmac);
7439 Printf(
" THIS IS MINUIT VERSION:%s",(
const char*)
fCvrsn);
7443 for (
id = 0;
id <= 6; ++
id) {
7445 if (
fIdbg[
id] >= 1) copt =
"ON ";
7446 Printf(
" DEBUG OPTION %3d IS %3s :%s"
7447 ,
id,(
const char*)copt,(
const char*)cdbopt[
id]);
7462 Printf(
" THE COMMAND:%10s IS UNKNOWN.",(
const char*)
fCword);
7469 if (strcmp(ctemp2.
Data(),
"SHO")) ckind =
"SHOW";
7471 Printf(
" THE FORMAT OF THE %4s COMMAND IS:",(
const char*)ckind);
7472 Printf(
" %s xxx [numerical arguments if any]",(
const char*)ckind);
7473 Printf(
" WHERE xxx MAY BE ONE OF THE FOLLOWING:");
7474 for (kk = 1; kk <= nname; ++kk) {
7475 Printf(
" %s",cname[kk-1]);
7481 Printf(
" ABOVE COMMAND IS ILLEGAL. IGNORED");
7503 Double_t dmin_, dxdi, yrho,
f, ynpp1, aming, ypbar;
7504 Double_t bestx, ystar, y1, y2, ystst, pb, wg;
7505 Double_t absmin, rho, sig2, rho1, rho2;
7506 Int_t npfn, i, j, k, jhold, ncycl, nparx;
7507 Int_t nparp1, kg, jh, nf, jl, ns;
7509 if (
fNpar <= 0)
return;
7518 rho2 = rho1 + alpha*
gamma;
7521 Printf(
" START SIMPLEX MINIMIZATION. CONVERGENCE WHEN EDM .LT. %g",
fEpsi);
7523 for (i = 1; i <=
fNpar; ++i) {
7536 for (i = 1; i <=
fNpar; ++i) {
7547 if (f <= aming)
goto L6;
7549 if (kg == 1)
goto L8;
7553 if (nf < 3)
goto L4;
7563 if (ns < 6)
goto L4;
7567 if (aming < absmin) jl = i;
7568 if (aming < absmin) absmin = aming;
7588 for (i = 1; i <=
fNpar; ++i) {
7596 if (ystar >=
fAmin)
goto L70;
7598 for (i = 1; i <=
fNpar; ++i) {
7604 y1 = (ystar -
fSIMPy[jh-1])*rho2;
7605 y2 = (ystst -
fSIMPy[jh-1])*rho1;
7606 rho = (rho2*y1 - rho1*y2)*.5 / (y1 - y2);
7607 if (rho < rhomin)
goto L66;
7608 if (rho > rhomax) rho = rhomax;
7609 for (i = 1; i <=
fNpar; ++i) {
7614 if (yrho <
fSIMPy[jl-1] && yrho < ystst)
goto L65;
7615 if (ystst <
fSIMPy[jl-1])
goto L67;
7616 if (yrho >
fSIMPy[jl-1])
goto L66;
7622 if (ystst <
fSIMPy[jl-1])
goto L67;
7629 if (
fISW[4] < 2)
goto L50;
7630 if (
fISW[4] >= 3 || ncycl % 10 == 0) {
7636 if (ystar >=
fSIMPy[jh-1])
goto L73;
7639 if (jhold != jh)
goto L50;
7642 for (i = 1; i <=
fNpar; ++i) {
7647 if (ystst >
fSIMPy[jh-1])
goto L1;
7649 if (ystst <
fAmin)
goto L67;
7655 Printf(
" SIMPLEX MINIMIZATION HAS CONVERGED.");
7661 Printf(
" SIMPLEX TERMINATES WITHOUT CONVERGENCE.");
7667 for (i = 1; i <=
fNpar; ++i) {
7737 static TString cpt =
" ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz1234567890./;:[]$%*_!@#&+()";
7740 l = strlen((
const char*)cfname);
7741 for (i = 1; i <=
l; ++i) {
7742 for (ic = 1; ic <= 80; ++ic) {
7743 if (cfname[i-1] == cpt[ic-1])
goto L100;
7767 Int_t i, j, k, kp1, km1;
7775 if (n < 1)
goto L100;
7778 for (i = 1; i <=
n; ++i) {
7780 if (si <= 0)
goto L100;
7783 for (i = 1; i <=
n; ++i) {
7784 for (j = 1; j <=
n; ++j) {
7789 for (i = 1; i <=
n; ++i) {
7792 if (a[k + k*l] != 0)
fVERTq[k-1] = 1 / a[k + k*
l];
7798 if (km1 < 0)
goto L100;
7799 else if (km1 == 0)
goto L50;
7802 for (j = 1; j <= km1; ++j) {
7808 if (k - n < 0)
goto L51;
7809 else if (k - n == 0)
goto L60;
7812 for (j = kp1; j <=
n; ++j) {
7819 for (j = 1; j <=
n; ++j) {
7824 for (j = 1; j <=
n; ++j) {
7825 for (k = 1; k <= j; ++k) {
7827 a[j + k*
l] = a[k + j*
l];
7855 const Int_t kMAXMES = 10;
7856 Int_t ityp, i, ic, nm;
7859 if (corg(0,3) !=
"SHO" || cmes(0,3) !=
"SHO") {
7865 Printf(
" MINUIT WARNING IN %s",(
const char*)corg);
7866 Printf(
" ============== %s",(
const char*)cmes);
7872 Printf(
" MINUIT DEBUG FOR %s",(
const char*)corg);
7873 Printf(
" =============== %s ",(
const char*)cmes);
7898 englsh =
" WAS SUPPRESSED. ";
7899 if (
fNwrmes[ityp-1] > 1) englsh =
"S WERE SUPPRESSED.";
7901 ,(
const char*)ctyp,(
const char*)englsh);
7905 Printf(
" ONLY THE MOST RECENT 10 WILL BE LISTED BELOW.");
7909 Printf(
" CALLS ORIGIN MESSAGE");
7910 for (i = 1; i <= nm; ++i) {
7912 if (ic > kMAXMES) ic = 1;
7929 Double_t denom, ba, al, dx, du1, du2;
7930 Int_t ndex, ierr, i, j, k,
l, ndiag, k1, iin;
7934 for (l = 1; l <=
fNpar; ++
l) {
7935 ndex = l*(l + 1) / 2;
7940 ba =
fBlim[i-1] - al;
7943 if (dx > 1) du1 = ba;
7951 for (i = 1; i <=
fNpar; ++i) {
7954 for (j = 1; j <= i; ++j) {
7962 for (iin = 1; iin <=
fNpar; ++iin) {
7963 ndiag = iin*(iin + 1) / 2;
7965 if (denom <= 1 && denom >= 0)
fGlobcc[iin-1] = 0;
virtual void mnmnos()
*-*-*-*-*-*-*-*-*-*-*Performs a MINOS error analysis*-*-*-*-*-*-*-*-*-*-*-* *-* =====================...
virtual void mnsimp()
*-*-*-*-*Minimization using the simplex method of Nelder and Mead*-*-*-*-* *-* ======================...
virtual Int_t GetNumFreePars() const
returns the number of currently free parameters
virtual void mnstat(Double_t &fmin, Double_t &fedm, Double_t &errdef, Int_t &npari, Int_t &nparx, Int_t &istat)
*-*-*-*-*Returns concerning the current status of the minimization*-*-*-*-* *-* =====================...
double dist(Rotation3D const &r1, Rotation3D const &r2)
virtual void mneval(Double_t anext, Double_t &fnext, Int_t &ierev)
*-*-*-*-*-*-*Evaluates the function being analyzed by MNCROS*-*-*-*-*-*-*-* *-* =====================...
virtual void mnprin(Int_t inkode, Double_t fval)
*-*-*-*Prints the values of the parameters at the time of the call*-*-*-*-* *-* =====================...
virtual Int_t SetPrintLevel(Int_t printLevel=0)
set Minuit print level printlevel = -1 quiet (also suppresse all warnings) = 0 normal = 1 verbose ...
virtual void mnexin(Double_t *pint)
*-*-*-*-*Transforms the external parameter values U to internal values*-*-* *-* =====================...
virtual void mnhes1()
*-*-*-*Calculate first derivatives (GRD) and uncertainties (DGRD)*-*-*-*-*-* *-* ====================...
void fcn(Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag)
virtual void mnvert(Double_t *a, Int_t l, Int_t m, Int_t n, Int_t &ifail)
*-*-*-*-*-*-*-*-*-*-*-*Inverts a symmetric matrix*-*-*-*-*-*-*-*-*-*-*-*-* *-* ======================...
Implementation in C++ of the Minuit package written by F.
ClassImp(TSeqCollection) Int_t TSeqCollection TIter next(this)
Return index of object in collection.
virtual void mnemat(Double_t *emat, Int_t ndim)
Calculates the external error matrix from the internal matrix.
virtual void mnwarn(const char *copt, const char *corg, const char *cmes)
*-*-*-*-*-*-*-*-*-*-*-*Prints Warning messages*-*-*-*-*-*-*-*-*-*-*-*-*-* *-* =======================...
virtual void mnline(Double_t *start, Double_t fstart, Double_t *step, Double_t slope, Double_t toler)
*-*-*-*-*-*-*-*-*-*Perform a line search from position START*-*-*-*-*-*-*-* *-* =====================...
void ToUpper()
Change string to upper case.
virtual void mnseek()
*-*-*-*Performs a rough (but global) minimization by monte carlo search*-* *-* ======================...
Int_t LoadPlugin()
Load the plugin library for this handler.
Short_t Min(Short_t a, Short_t b)
virtual void mnderi()
*-*-*-*-*-*-*-*Calculates the first derivatives of FCN (GRD)*-*-*-*-*-*-*-* *-* =====================...
virtual void mncros(Double_t &aopt, Int_t &iercr)
*-*-*-*-*-*-*-*-*-*-*Find point where MNEVAL=AMIN+UP*-*-*-*-*-*-*-*-*-*-*-* *-* =====================...
virtual void mngrad()
*-*-*-*-*-*-*-*-*-*Interprets the SET GRAD command*-*-*-*-*-*-*-*-*-*-*-*-* *-* =====================...
virtual void mnbins(Double_t a1, Double_t a2, Int_t naa, Double_t &bl, Double_t &bh, Int_t &nb, Double_t &bwid)
*-*-*-*-*-*-*-*-*-*-*Compute reasonable histogram intervals*-*-*-*-*-*-*-*-* *-* ====================...
virtual void mnrset(Int_t iopt)
*-*-*-*-*-*-*-*Resets function value and errors to UNDEFINED*-*-*-*-*-*-*-* *-* =====================...
Long_t ExecPlugin(int nargs, const T &...params)
virtual void mntiny(volatile Double_t epsp1, Double_t &epsbak)
*-*-*-*-*-*-*-*To find the machine precision*-*-*-*-*-*-*-*-*-*-*-*-*-*-* *-* =======================...
LongDouble_t Power(LongDouble_t x, LongDouble_t y)
double beta(double x, double y)
Calculates the beta function.
virtual Int_t GetNumPars() const
returns the total number of parameters that have been defined as fixed or free.
TString fOrigin[kMAXWARN]
virtual Int_t FixParameter(Int_t parNo)
fix a parameter
virtual Int_t Eval(Int_t npar, Double_t *grad, Double_t &fval, Double_t *par, Int_t flag)
Evaluate the minimisation function Input parameters: npar: number of currently variable parameters pa...
virtual void mneig(Double_t *a, Int_t ndima, Int_t n, Int_t mits, Double_t *work, Double_t precis, Int_t &ifault)
*-*-*-*-*-*-*-*-*-*-*-*Compute matrix eigen values*-*-*-*-*-*-*-*-*-*-*-*-* *-* =====================...
void SetParamPtrs(void *paramArr, Int_t nparam=-1)
ParamArr is an array containing the function argument values.
const char * Data() const
virtual void mnmigr()
*-*-*-*-*-*-*-*-*Performs a local function minimization*-*-*-*-*-*-*-*-*-* *-* ======================...
virtual void SetFCN(void *fcn)
*-*-*-*-*-*-*To set the address of the minimization function*-*-*-*-*-*-*-* *-* =====================...
virtual void mnfree(Int_t k)
*-*-*-*Restores one or more fixed parameter(s) to variable status*-*-*-*-*-* *-* ====================...
static const double x2[5]
static TString Format(const char *fmt,...)
Static method which formats a string using a printf style format descriptor and return a TString...
The TNamed class is the base class for all named ROOT classes.
virtual void mncrck(TString crdbuf, Int_t maxcwd, TString &comand, Int_t &lnc, Int_t mxp, Double_t *plist, Int_t &llist, Int_t &ierr, Int_t isyswr)
*-*-*-*-*-*-*-*-*-*-*-*Cracks the free-format input*-*-*-*-*-*-*-*-*-*-*-*-* *-* ====================...
ClassImp(TMinuit) TMinuit
*-*-*-*-*-*-*-*-*-*-*Minuit normal constructor*-*-*-*-*-*-*-*-*-*-*-*-*-*-* *-* =====================...
virtual void mnrn15(Double_t &val, Int_t &inseed)
*-*-*-*-*-*-*This is a super-portable random number generator*-*-*-*-*-*-* *-* ======================...
std::map< std::string, std::string >::const_iterator iter
virtual TObject * Contour(Int_t npoints=10, Int_t pa1=0, Int_t pa2=1)
Creates a TGraph object describing the n-sigma contour of a TMinuit fit.
Double_t Log10(Double_t x)
virtual ~TMinuit()
*-*-*-*-*-*-*-*-*-*-*Minuit default destructor*-*-*-*-*-*-*-*-*-*-*-*-*-*-* *-* =====================...
static const double x4[22]
static Vc_ALWAYS_INLINE Vector< T > abs(const Vector< T > &x)
virtual void DeleteArrays()
*-*-*-*-*-*-*-*-*-*-*-*Delete internal Minuit arrays*-*-*-*-*-*-*-*-* *-* ===========================...
std::vector< std::vector< double > > Data
virtual void mncuve()
*-*-*-*-*-*-*-*Makes sure that the current point is a local minimum*-*-*-*-* *-* ====================...
virtual TObject * Clone(const char *newname="") const
Make a clone of an object using the Streamer facility.
virtual void BuildArrays(Int_t maxpar=15)
*-*-*-*-*-*-*Create internal Minuit arrays for the maxpar parameters*-*-*-* *-* =====================...
Method or function calling interface.
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
virtual void mnmnot(Int_t ilax, Int_t ilax2, Double_t &val2pl, Double_t &val2mi)
*-*-*-*-*-*Performs a MINOS error analysis on one parameter*-*-*-*-*-*-*-*-* *-* ====================...
virtual Int_t DefineParameter(Int_t parNo, const char *name, Double_t initVal, Double_t initErr, Double_t lowerLimit, Double_t upperLimit)
Define a parameter.
virtual void mncont(Int_t ke1, Int_t ke2, Int_t nptu, Double_t *xptu, Double_t *yptu, Int_t &ierrf)
*-*-*-*-*-*-*Find points along a contour where FCN is minimum*-*-*-*-*-*-* *-* ======================...
TString * fCpnam
Character to be plotted at the X,Y contour positions.
virtual void mncler()
*-*-*-*-*-*-*-*-*-*-*Resets the parameter list to UNDEFINED*-*-*-*-*-*-*-* *-* ======================...
virtual void mnparm(Int_t k, TString cnamj, Double_t uk, Double_t wk, Double_t a, Double_t b, Int_t &ierflg)
*-*-*-*-*-*-*-*-*Implements one parameter definition*-*-*-*-*-*-*-*-*-*-*-* *-* =====================...
virtual Int_t GetParameter(Int_t parNo, Double_t ¤tValue, Double_t ¤tError) const
return parameter value and error
virtual void mnfixp(Int_t iint, Int_t &ierr)
*-*-*-*-*-*-*Removes parameter IINT from the internal parameter list*-*-* *-* =======================...
virtual void mncalf(Double_t *pvec, Double_t &ycalf)
*-*-*-*-*-*-*-*-*-*Transform FCN to find further minima*-*-*-*-*-*-*-*-*-* *-* ======================...
void Form(const char *fmt,...)
Formats a string using a printf style format descriptor.
char * Form(const char *fmt,...)
virtual Int_t SetErrorDef(Double_t up)
To get the n-sigma contour the error def parameter "up" has to set to n^2.
TMethodCall * GetMethodCall() const
TMethodCall * fMethodCall
virtual TObject * Clone(const char *newname="") const
Make a clone of an object using the Streamer facility.
virtual void mnerrs(Int_t number, Double_t &eplus, Double_t &eminus, Double_t &eparab, Double_t &gcc)
*-*-*-*-*-*-*-*-*-*Utility routine to get MINOS errors*-*-*-*-*-*-*-*-*-*-* *-* =====================...
virtual void mnhess()
*-*-*-*-*-*Calculates the full second-derivative matrix of FCN*-*-*-*-*-*-*-* *-* ===================...
void InitWithPrototype(TClass *cl, const char *method, const char *proto, Bool_t objectIsConst=kFALSE, ROOT::EFunctionMatchMode mode=ROOT::kConversionMatch)
Initialize the method invocation environment.
virtual Int_t GetNumFixedPars() const
returns the number of currently fixed parameters
virtual void mnpsdf()
*-*-*-*-*-*Calculates the eigenvalues of v to see if positive-def*-*-*-*-* *-* ======================...
virtual void mndxdi(Double_t pint, Int_t ipar, Double_t &dxdi)
*-*-*-*Calculates the transformation factor between ext/internal values*-* *-* ======================...
virtual void mnamin()
*-*-*-*-*-*-*-*-*-*-*-*-*Initialize AMIN*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-* *-* =============== *-*C Ca...
static const double x1[5]
virtual void mnscan()
*-*-*-*-*Scans the values of FCN as a function of one parameter*-*-*-*-*-* *-* ======================...
TString fWarmes[kMAXWARN]
virtual void mnhelp(TString comd)
*-*-*-*-*-*-*-*HELP routine for MINUIT interactive commands*-*-*-*-*-*-*-*-* *-* ====================...
void InteractiveFCNm(Int_t &npar, Double_t *gin, Double_t &f, Double_t *u, Int_t flag)
*-*-*-*-*-*-*Static function called when SetFCN is called in interactive mode *-* ===================...
virtual Int_t Release(Int_t parNo)
release a parameter
virtual void mnpfit(Double_t *parx2p, Double_t *pary2p, Int_t npar2p, Double_t *coef2p, Double_t &sdev2p)
*-*-*-*-*-*-*-*-*-*To fit a parabola to npar2p points*-*-*-*-*-*-*-*-*-*-* *-* ======================...
virtual void mnlims()
*-*-*-*Interprets the SET LIM command, to reset the parameter limits*-*-*-* *-* =====================...
virtual void mnexcm(const char *comand, Double_t *plist, Int_t llist, Int_t &ierflg)
*-*-*-*-*-*Interprets a command and takes appropriate action*-*-*-*-*-*-*-* *-* =====================...
virtual void mncntr(Int_t ke1, Int_t ke2, Int_t &ierrf)
*-*-*-*-*Print function contours in two variables, on line printer*-*-*-*-* *-* =====================...
virtual void SetMaxIterations(Int_t maxiter=500)
virtual void mncomd(const char *crdbin, Int_t &icondn)
*-*-*-*-*-*-*-*-*-*-*Reads a command string and executes*-*-*-*-*-*-*-*-*-* *-* =====================...
virtual Int_t Command(const char *command)
execute a Minuit command Equivalent to MNEXCM except that the command is given as a character string...
virtual void mninit(Int_t i1, Int_t i2, Int_t i3)
*-*-*-*-*-*Main initialization member function for MINUIT*-*-*-*-*-*-*-*-* *-* ======================...
Mother of all ROOT objects.
virtual Int_t Migrad()
invokes the MIGRAD minimizer
virtual void mnpars(TString &crdbuf, Int_t &icondn)
*-*-*-*-*-*-*-*Implements one parameter definition*-*-*-*-*-*-*-*-*-*-*-*-* *-* =========== =========...
virtual void mnplot(Double_t *xpt, Double_t *ypt, char *chpt, Int_t nxypt, Int_t npagwd, Int_t npagln)
*-*-*-*Plots points in array xypt onto one page with labelled axes*-*-*-*-* *-* =====================...
virtual void mnpint(Double_t &pexti, Int_t i, Double_t &pinti)
*-*-*-*-*-*-*Calculates the internal parameter value PINTI*-*-*-*-*-*-*-* *-* =======================...
virtual void mnsave()
*-*-*-*Writes current parameter values and step sizes onto file ISYSSA*-*-* *-* =====================...
virtual void mninex(Double_t *pint)
*-*-*-*-*Transforms from internal coordinates (PINT) to external (U)*-*-*-* *-* =====================...
void Execute(const char *, const char *, int *=0)
Execute method on this object with the given parameter string, e.g.
Short_t Max(Short_t a, Short_t b)
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
virtual void mnmatu(Int_t kode)
*-*-*-*-*-*-*-*Prints the covariance matrix v when KODE=1*-*-*-*-*-*-*-*-* *-* ======================...
virtual void mnset()
*-*-*-*-*Interprets the commands that start with SET and SHOW*-*-*-*-*-*-* *-* ======================...
ClassImp(TSlaveInfo) Int_t TSlaveInfo const TSlaveInfo * si
Used to sort slaveinfos by ordinal.
Bool_t mnunpt(TString &cfname)
*-*-*-*-*-*Returns .TRUE.
virtual void mnrazz(Double_t ynew, Double_t *pnew, Double_t *y, Int_t &jh, Int_t &jl)
*-*-*-*-*Called only by MNSIMP (and MNIMPR) to add a new point*-*-*-*-*-*-* *-* =====================...
virtual const char * Getp2f2funcname(void *) const
Double_t Sqrt(Double_t x)
R__EXTERN TInterpreter * gCling
virtual void mnpout(Int_t iuext, TString &chnam, Double_t &val, Double_t &err, Double_t &xlolim, Double_t &xuplim, Int_t &iuint) const
*-*-*-*Provides the user with information concerning the current status*-*-* *-* ====================...
virtual void mnimpr()
*-*-*-*-*-*-*Attempts to improve on a good local minimum*-*-*-*-*-*-*-*-*-* *-* =====================...
void Resize(Ssiz_t n)
Resize the string. Truncate or add blanks as necessary.
static const double x3[11]
virtual void mnwerr()
*-*-*-*-*-*-*-*Calculates the WERR, external parameter errors*-*-*-*-*-*-* *-* ======================...
void(* fFCN)(Int_t &npar, Double_t *gin, Double_t &f, Double_t *u, Int_t flag)
virtual void Warning(const char *method, const char *msgfmt,...) const
Issue warning message.