root [1] .x fit_3DMC.C
Warning: Automatic variable ncount is allocated FILE:fit_3DMC.C LINE:39
start
reading histos into arrays
reading table : integFF_Dlnu_newSlope.txt
Starting timer
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
1 FF slope 1.35000e+00 2.00000e-01 -1.00000e+00 4.00000e+00
2 B_BpD 2.10000e-02 2.00000e-03 0.00000e+00 1.00000e+00
3 B_B0D 2.10000e-02 2.00000e-03 0.00000e+00 1.00000e+00
4 B_BpDS 5.60000e-02 5.00000e-03 0.00000e+00 1.00000e+00
5 B_B0DS 5.60000e-02 5.00000e-03 0.00000e+00 1.00000e+00
6 B_BpDD 1.50000e-02 5.00000e-03 0.00000e+00 1.00000e+00
7 f_BDD0 3.57000e-01 5.00000e-02 0.00000e+00 1.00000e+00
8 f_BDDP1 6.61000e-01 5.00000e-02 0.00000e+00 1.00000e+00
9 B_BpDPi 1.20000e-02 5.00000e-03 0.00000e+00 1.00000e+00
10 DiffB(D0) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
11 DBCascL(D0) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
12 CascL(D0) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
13 LMisID(D0) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
14 RealD(D0) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
15 DiffB(Dch) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
16 DBCascL(Dch) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
17 CascL(Dch) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
18 LMisID(Dch) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
19 RealD(Dch) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
20 f_+0 1.00600e+00 6.70000e-02 0.00000e+00 2.00000e+00
21 c_+0 1.00000e+00 6.80000e-02 0.00000e+00 2.00000e+00
**********
** 1 **CALL FCN 1
**********
c[0] = 1
c[1] = 1
c[2] = 1
c[3] = 1
c[4] = 1
c[5] = 1
c[6] = 1.006
c[7] = 1
c[8] = 1.006
c[9] = 1
c[10] = 1.006
c[11] = 1
c[12] = 1.006
c[13] = 1
c[14] = 1.006
c[15] = 1
c[16] = 1.00554
c[17] = 0.999547
c[18] = 1.00595
c[19] = 0.999947
c[20] = 1.00638
c[21] = 1.00038
c[22] = 1.00595
c[23] = 0.999947
c[24] = 1.006
c[25] = 1
c[26] = 1.006
c[27] = 1
c[28] = 1.006
c[29] = 1
c[30] = 1.006
c[31] = 1.006
c[32] = 1.006
c[33] = 1.006
c[34] = 1
c[35] = 1
c[36] = 1
c[37] = 1
c[38] = 1
c[39] = 1
c[40] = 1.006
c[41] = 1
c[42] = 1.006
c[43] = 1
c[44] = 1.006
c[45] = 1
c[46] = 1.006
c[47] = 1
c[48] = 1.006
c[49] = 1
c[50] = 1.00554
c[51] = 0.999547
c[52] = 1.00595
c[53] = 0.999947
c[54] = 1.00638
c[55] = 1.00038
c[56] = 1.00595
c[57] = 0.999947
c[58] = 1.006
c[59] = 1
c[60] = 1.006
c[61] = 1
c[62] = 1
c[63] = 1
c[64] = 1
c[65] = 1
c[66] = 1
Number of bins into chisq = 179
from D0 = 92
from Dch = 87
chisq_D0 = 102.237
chisq_Dch = 78.0128
chisq = 180.25
FCN=180.25 FROM CALl fcn STATUS=RESET 1 CALLS 1 TOTAL
EDM= unknown STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS PHYSICAL LIMITS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 FF slope 1.35000e+00 2.00000e-01 -1.00000e+00 4.00000e+00
2 B_BpD 2.10000e-02 2.00000e-03 0.00000e+00 1.00000e+00
3 B_B0D 2.10000e-02 2.00000e-03 0.00000e+00 1.00000e+00
4 B_BpDS 5.60000e-02 5.00000e-03 0.00000e+00 1.00000e+00
5 B_B0DS 5.60000e-02 5.00000e-03 0.00000e+00 1.00000e+00
6 B_BpDD 1.50000e-02 5.00000e-03 0.00000e+00 1.00000e+00
7 f_BDD0 3.57000e-01 5.00000e-02 0.00000e+00 1.00000e+00
8 f_BDDP1 6.61000e-01 5.00000e-02 0.00000e+00 1.00000e+00
9 B_BpDPi 1.20000e-02 5.00000e-03 0.00000e+00 1.00000e+00
10 DiffB(D0) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
11 DBCascL(D0) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
12 CascL(D0) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
13 LMisID(D0) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
14 RealD(D0) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
15 DiffB(Dch) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
16 DBCascL(Dch) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
17 CascL(Dch) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
18 LMisID(Dch) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
19 RealD(Dch) 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00
20 f_+0 1.00600e+00 6.70000e-02 0.00000e+00 2.00000e+00
21 c_+0 1.00000e+00 6.80000e-02 0.00000e+00 2.00000e+00
**********
** 2 **SET PRINT 0
**********
**********
** 3 **FIX 1
**********
**********
** 4 **FIX 3
**********
**********
** 5 **FIX 7
**********
**********
** 6 **FIX 8
**********
**********
** 7 **MIGRAD
**********
MIGRAD MINIMIZATION HAS CONVERGED.
FCN=151.217 FROM MIGRAD STATUS=CONVERGED 789 CALLS 790 TOTAL
EDM=1.06932e-06 STRATEGY= 1 ERROR MATRIX UNCERTAINTY 6.8 per cent
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 FF slope 1.35000e+00 fixed
2 B_BpD 2.13759e-02 1.01144e-03 -2.85458e-06 -1.66414e-02
3 B_B0D 2.10000e-02 fixed
4 B_BpDS 5.19661e-02 2.97987e-03 9.19515e-06 1.17008e-02
5 B_B0DS 6.29753e-02 4.21376e-03 -1.41130e-05 5.89150e-02
6 B_BpDD 1.57565e-02 1.00877e-03 1.52889e-05 -3.96468e-02
7 f_BDD0 3.57000e-01 fixed
8 f_BDDP1 6.61000e-01 fixed
9 B_BpDPi 1.10301e-02 1.49426e-03 -6.99449e-06 -1.16619e-01
10 DiffB(D0) 9.91279e-01 4.72823e-02 -3.92437e-05 -1.23911e-02
11 DBCascL(D0) 1.06017e+00 7.51007e-02 2.13912e-04 -2.84771e-03
12 CascL(D0) 1.05010e+00 1.41229e-01 -4.25385e-04 -5.01153e-03
13 LMisID(D0) 8.47052e-01 2.07092e-01 -1.79700e-04 -2.27451e-04
14 RealD(D0) 1.13349e+00 5.89146e-02 -7.24383e-05 1.36470e-02
15 DiffB(Dch) 8.68782e-01 1.16420e-01 1.23098e-04 -3.49654e-03
16 DBCascL(Dch) 1.14900e+00 1.46640e-01 -6.41508e-04 -6.85817e-03
17 CascL(Dch) 1.13887e+00 3.00987e-01 1.42480e-03 -2.11235e-03
18 LMisID(Dch) 1.00571e+00 3.64116e-01 -1.50750e-03 8.33304e-04
19 RealD(Dch) 1.01577e+00 1.29870e-01 3.56058e-04 -3.54994e-03
20 f_+0 1.00310e+00 4.50150e-02 4.24074e-06 -2.60273e-02
21 c_+0 9.89816e-01 4.64378e-02 6.06905e-05 1.73890e-02
**********
** 8 **SHOW COVARIANCE
**********
EXTERNAL ERROR MATRIX. NDIM= 21 NPAR= 17 ERR DEF=1
ELEMENTS ABOVE DIAGONAL ARE NOT PRINTED.
1.023e-06
-1.604e-07 8.880e-06
2.812e-06 -7.852e-06 1.776e-05
3.380e-07 -1.906e-07 7.367e-07 1.018e-06
3.508e-07 3.683e-07 4.321e-07 -4.533e-07 2.233e-06
1.386e-05 -5.546e-06 4.576e-05 9.901e-06 -1.725e-05 2.237e-03
2.117e-05 2.021e-05 4.953e-05 2.751e-05 -1.956e-05 -3.312e-04 5.651e-03
8.180e-06 1.895e-05 6.374e-05 -5.840e-05 -4.915e-05 7.824e-04 -2.223e-03 2.008e-02
1.541e-05 8.027e-05 1.543e-05 -1.636e-06 3.664e-05 -2.457e-04 -3.334e-03 -7.345e-03 4.353e-02
-2.579e-06 3.983e-05 -1.600e-05 -6.888e-06 1.828e-05 2.295e-04 -4.515e-04 -5.662e-04 -9.806e-04 3.475e-03
6.670e-05 3.851e-07 1.726e-04 3.034e-05 -5.163e-06 1.299e-03 1.535e-03 1.351e-03 1.132e-03 -4.618e-04 1.362e-02
4.811e-05 -6.576e-05 2.684e-04 3.743e-05 -3.244e-05 1.154e-03 2.455e-03 4.522e-04 5.410e-04 3.024e-04 3.965e-04
4.811e-05 -6.576e-05 2.684e-04 3.743e-05 -3.244e-05 1.154e-03 2.455e-03 4.522e-04 5.410e-04 3.024e-04 3.965e-04 2.166e-02
1.268e-04 4.255e-05 4.096e-04 1.869e-05 -5.939e-05 3.052e-03 3.933e-03 8.869e-03 2.826e-03 2.640e-04 1.403e-02
1.268e-04 4.255e-05 4.096e-04 1.869e-05 -5.939e-05 3.052e-03 3.933e-03 8.869e-03 2.826e-03 2.640e-04 1.403e-02 -3.518e-03 9.353e-02
7.720e-05 -8.362e-05 1.185e-04 7.854e-06 8.314e-05 6.259e-04 -1.863e-03 4.247e-04 3.600e-04 -2.529e-03 -6.078e-03
7.720e-05 -8.362e-05 1.185e-04 7.854e-06 8.314e-05 6.259e-04 -1.863e-03 4.247e-04 3.600e-04 -2.529e-03 -6.078e-03 5.047e-03 -2.934e-02 1.389e-01
2.123e-05 -1.181e-04 1.867e-04 -2.574e-06 3.620e-05 -2.165e-05 -3.201e-04 -5.291e-04 -7.531e-05 -1.293e-04 2.314e-03
2.123e-05 -1.181e-04 1.867e-04 -2.574e-06 3.620e-05 -2.165e-05 -3.201e-04 -5.291e-04 -7.531e-05 -1.293e-04 2.314e-03 6.831e-04 2.343e-03 -9.369e-03 1.696e-02
-3.779e-05 -4.397e-05 -9.125e-05 -1.382e-05 -1.640e-05 -6.712e-04 -1.037e-03 -3.179e-04 -2.090e-03 -4.867e-04 -2.506e-03
-3.779e-05 -4.397e-05 -9.125e-05 -1.382e-05 -1.640e-05 -6.712e-04 -1.037e-03 -3.179e-04 -2.090e-03 -4.867e-04 -2.506e-03 -2.224e-03 -5.662e-03 -1.096e-03 -4.103e-04 2.028e-03
-4.165e-05 4.957e-05 -1.652e-04 -1.415e-05 -1.505e-05 -5.560e-04 -7.480e-04 -2.167e-04 -4.376e-04 1.933e-04 -3.041e-03
-4.165e-05 4.957e-05 -1.652e-04 -1.415e-05 -1.505e-05 -5.560e-04 -7.480e-04 -2.167e-04 -4.376e-04 1.933e-04 -3.041e-03 -2.820e-03 -5.479e-03 -3.746e-03 -2.248e-03 1.421e-03 2.158e-03
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 2 4 5 6 9 10 11 12 13 14 15 16 17 18 19 20
2 0.96117 1.000 -0.053 0.660 0.331 0.232 0.290 0.278 0.057 0.073 -0.043 0.565 0.323 0.410 0.205 0.161 -0.830
4 0.98128 -0.053 1.000 -0.625 -0.063 0.083 -0.039 0.090 0.045 0.129 0.227 0.001 -0.150 0.047 -0.075 -0.304 -0.328
5 0.98304 0.660 -0.625 1.000 0.173 0.069 0.230 0.156 0.107 0.018 -0.064 0.351 0.433 0.318 0.075 0.340 -0.481
6 0.87876 0.331 -0.063 0.173 1.000 -0.301 0.208 0.363 -0.409 -0.008 -0.116 0.258 0.252 0.061 0.021 -0.020 -0.304
9 0.88815 0.232 0.083 0.069 -0.301 1.000 -0.244 -0.174 -0.232 0.118 0.208 -0.030 -0.148 -0.130 0.149 0.186 -0.244
10 0.71274 0.290 -0.039 0.230 0.208 -0.244 1.000 -0.093 0.117 -0.025 0.082 0.235 0.166 0.211 0.036 -0.004 -0.315
11 0.75879 0.278 0.090 0.156 0.363 -0.174 -0.093 1.000 -0.209 -0.213 -0.102 0.175 0.222 0.171 -0.067 -0.033 -0.306
12 0.81752 0.057 0.045 0.107 -0.409 -0.232 0.117 -0.209 1.000 -0.248 -0.068 0.082 0.022 0.205 0.008 -0.029 -0.050
13 0.76178 0.073 0.129 0.018 -0.008 0.118 -0.025 -0.213 -0.248 1.000 -0.080 0.047 0.018 0.044 0.005 -0.003 -0.222
14 0.62063 -0.043 0.227 -0.064 -0.116 0.208 0.082 -0.102 -0.068 -0.080 1.000 -0.067 0.035 0.015 -0.115 -0.017 -0.183
15 0.78431 0.565 0.001 0.351 0.258 -0.030 0.235 0.175 0.082 0.047 -0.067 1.000 0.023 0.393 -0.140 0.152 -0.477
16 0.74371 0.323 -0.150 0.433 0.252 -0.148 0.166 0.222 0.022 0.018 0.035 0.023 1.000 -0.078 0.092 0.036 -0.336
17 0.75323 0.410 0.047 0.318 0.061 -0.130 0.211 0.171 0.205 0.044 0.015 0.393 -0.078 1.000 -0.257 0.059 -0.411
18 0.73486 0.205 -0.075 0.075 0.021 0.149 0.036 -0.067 0.008 0.005 -0.115 -0.140 0.092 -0.257 1.000 -0.193 -0.065
19 0.69462 0.161 -0.304 0.340 -0.020 0.186 -0.004 -0.033 -0.029 -0.003 -0.017 0.152 0.036 0.059 -0.193 1.000 -0.070
20 0.98912 -0.830 -0.328 -0.481 -0.304 -0.244 -0.315 -0.306 -0.050 -0.222 -0.183 -0.477 -0.336 -0.411 -0.065 -0.070 1.000
21 0.98054 -0.887 0.358 -0.844 -0.302 -0.217 -0.253 -0.214 -0.033 -0.045 0.071 -0.561 -0.412 -0.386 -0.216 -0.371 0.679
1.000
**********
** 9 **CALL FCN 2
**********
c[0] = 1
c[1] = 1
c[2] = 1
c[3] = 1
c[4] = 1
c[5] = 1.13349
c[6] = 1.0031
c[7] = 1
c[8] = 0.994353
c[9] = 0.991279
c[10] = 1.06346
c[11] = 1.06017
c[12] = 1.05336
c[13] = 1.0501
c[14] = 0.849679
c[15] = 0.847052
c[16] = 1.05322
c[17] = 1.04996
c[18] = 1.05364
c[19] = 1.05038
c[20] = 1.05409
c[21] = 1.05083
c[22] = 1.05364
c[23] = 1.05038
c[24] = 0.922025
c[25] = 0.919174
c[26] = 0.922025
c[27] = 0.919174
c[28] = 0.930845
c[29] = 1.12456
c[30] = 1.02106
c[31] = 1.02106
c[32] = 1.02106
c[33] = 1.02106
c[34] = 1
c[35] = 1
c[36] = 0.989816
c[37] = 0.989816
c[38] = 0.989816
c[39] = 1.00542
c[40] = 0.992886
c[41] = 0.989816
c[42] = 0.862601
c[43] = 0.859934
c[44] = 1.14082
c[45] = 1.13729
c[46] = 1.13077
c[47] = 1.12727
c[48] = 0.99856
c[49] = 0.995473
c[50] = 1.04249
c[51] = 1.03927
c[52] = 1.04291
c[53] = 1.03968
c[54] = 1.04336
c[55] = 1.04013
c[56] = 1.04291
c[57] = 1.03968
c[58] = 0.912635
c[59] = 0.909813
c[60] = 0.912635
c[61] = 0.909813
c[62] = 1.11311
c[63] = 1.00753
c[64] = 1.00753
c[65] = 1.00753
c[66] = 1.00753
chisq_D0 = 86.2991
chisq_Dch = 64.894
chisq = 151.217
BFp1 = 0.0213759
BFp1_err = 0.00101144
BF01 = 0.021
BF01_err = 0
BFp2 = 0.0519661
BFp2_err = 0.00297987
BF02 = 0.0629753
BF02_err = 0.00421376
BF3T = 0.0157565
BF3T_err = 0.00100877
FDD0 = 0.357
FDD0_err = 0
FDDP1 = 0.661
FDDP1_err = 0
BF4T = 0.0110301
BF4T_err = 0.00149426
f_BDD2 = 0.660714
f_BDPi = 3
****************************
** B+ Branching Fractions **
****************************
TotFracBp = 0.100129
inclBpF = 0.104
B(B->Dlnu) = 0.0222024 +- 0.00105055 (4.7317 %)
B(B->D*lnu) = 0.0539753 +- 0.00309508 (5.73425 %)
B(B->D**lnu) = 0.0163657 +- 0.00104778 (6.40226 %)
B(B->D(*)Pilnu) = 0.0114566 +- 0.00155203 (13.5471 %)
B(B->D0*lnu) = 0.00218111
B(B->D1lnu) = 0.00610955
B(B->D1'lnu) = 0.00403841
B(B->D2*lnu) = 0.00403667
B(B->DPiplnu) = 0.00572828
B(B->DPi0lnu) = 0.00286414
B(B->D*Piplnu) = 0.00190943
B(B->D*Pi0lnu) = 0.000954713
****************************
** B0 Branching Fractions **
****************************
TotFracB0 = 0.111138
inclB0F = 0.104
B(B->Dlnu) = 0.020003 +- 0.000946483 (4.7317 %)
B(B->D*lnu) = 0.0589307 +- 0.00394313 (6.69113 %)
B(B->D**lnu) = 0.0147446 +- 0.000943986 (6.40226 %)
B(B->D(*)Pilnu) = 0.0103217 +- 0.00139829 (13.5471 %)
B(B->D0*lnu) = 0.00196505
B(B->D1lnu) = 0.00550435
B(B->D1'lnu) = 0.00363837
B(B->D2*lnu) = 0.0036368
B(B->DPiplnu) = 0.00516084
B(B->DPi0lnu) = 0.00258042
B(B->D*Piplnu) = 0.00172028
B(B->D*Pi0lnu) = 0.00086014
****************************
** B->Dlnu slope **
****************************
slope = 1.35 +- 0 (0 %)
****************************
** Background **
****************************
BkgD0(DiffB) = 0.991279 +- 0.0472823 (4.76983 %)
BkgD0(DBCascL) = 1.06017 +- 0.0751007 (7.08385 %)
BkgD0(CascL) = 1.0501 +- 0.141229 (13.4491 %)
BkgD0(LMisID) = 0.847052 +- 0.207092 (24.4486 %)
BkgD0(RealD) = 1.13349 +- 0.0589146 (5.19761 %)
BkgDch(DiffB) = 0.868782 +- 0.11642 (13.4003 %)
BkgDch(DBCascL) = 1.149 +- 0.14664 (12.7624 %)
BkgDch(CascL) = 1.13887 +- 0.300987 (26.4286 %)
BkgDch(LMisID) = 1.00571 +- 0.364116 (36.2047 %)
BkgDch(RealD) = 1.01577 +- 0.12987 (12.7854 %)
****************************
** Constants **
****************************
f_+0 = f+-/f00 = 1.0031 +- 0.045015 (4.48758 %)
c_+0 = D0/D+ = 0.989816 +- 0.0464378 (4.69156 %)
Time at the end of job = 242.21 seconds
root [2]