root [0] .x fit_3D.C Warning: Automatic variable ncount is allocated FILE:fit_3D.C LINE:176 start reading 3D histos into arrays reading Dlnu Stat histos into arrays reading Dstar Stat histos into arrays selecting input files making binned arrays Starting timer PARAMETER DEFINITIONS: NO. NAME VALUE STEP SIZE LIMITS 1 dummy 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 2 D slope 1.40000e+00 1.00000e-02 0.00000e+00 2.00000e+00 3 Br_BpD 3.00000e-02 1.00000e-03 0.00000e+00 1.00000e-01 4 Br_B0D 2.10000e-02 1.00000e-03 0.00000e+00 1.00000e-01 5 Br_BpD* 6.00000e-02 1.00000e-03 0.00000e+00 1.00000e-01 6 Br_B0D* 4.50000e-02 1.00000e-03 0.00000e+00 1.00000e-01 7 BpDDPiAll 1.80000e-02 1.00000e-03 0.00000e+00 1.00000e-01 8 f_BpDD12 7.67000e-01 1.00000e-02 0.00000e+00 2.00000e+00 9 f_BpDPiDD12 9.11000e-01 1.00000e-02 0.00000e+00 2.00000e+00 10 f_BpDSPiDD12 9.92000e-01 1.00000e-02 0.00000e+00 2.00000e+00 11 f_BpDPiDD0 4.00000e-01 1.00000e-02 0.00000e+00 1.00000e+01 12 f_BpDSPiDDP1 1.00000e-01 1.00000e-02 0.00000e+00 1.00000e+01 13 UncorrDL 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 14 UncorrCL 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 15 CascL 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 16 f_DchD0Bkg 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 17 f_+0 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 18 c_0 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 19 c_+ 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 20 D* slope 1.14500e+00 1.00000e-02 0.00000e+00 2.00000e+00 21 D* R1 1.39600e+00 1.00000e-02 0.00000e+00 2.00000e+00 22 D* R2 8.85000e-01 1.00000e-02 0.00000e+00 2.00000e+00 start fitting ********** ** 1 **CALL FCN 1 ********** rho2 = 1.4 oldDlnuBp = 0.268702 oldDlnuB0 = 0.26542 newDlnuBp = 0.159356 newDlnuB0 = 0.157885 RNBp = 1.68617 RNB0 = 1.68109 rho2Ds = 1.145 oldDstarBp = 0.865978 oldDstarB0 = 0.860858 newDstarBp = 0.824051 newDstarB0 = 0.818848 RNDstarBp = 1.05088 RNDstarB0 = 1.0513 Number of bins into chisq = 267 from D0 = 138 from Dch = 129 chisq_D0 = 459.65 chisq_Dch = 212.387 chisq = 673.126 FCN=673.126 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 dummy 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 2 D slope 1.40000e+00 1.00000e-02 0.00000e+00 2.00000e+00 3 Br_BpD 3.00000e-02 1.00000e-03 0.00000e+00 1.00000e-01 4 Br_B0D 2.10000e-02 1.00000e-03 0.00000e+00 1.00000e-01 5 Br_BpD* 6.00000e-02 1.00000e-03 0.00000e+00 1.00000e-01 6 Br_B0D* 4.50000e-02 1.00000e-03 0.00000e+00 1.00000e-01 7 BpDDPiAll 1.80000e-02 1.00000e-03 0.00000e+00 1.00000e-01 8 f_BpDD12 7.67000e-01 1.00000e-02 0.00000e+00 2.00000e+00 9 f_BpDPiDD12 9.11000e-01 1.00000e-02 0.00000e+00 2.00000e+00 10 f_BpDSPiDD12 9.92000e-01 1.00000e-02 0.00000e+00 2.00000e+00 11 f_BpDPiDD0 4.00000e-01 1.00000e-02 0.00000e+00 1.00000e+01 12 f_BpDSPiDDP1 1.00000e-01 1.00000e-02 0.00000e+00 1.00000e+01 13 UncorrDL 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 14 UncorrCL 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 15 CascL 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 16 f_DchD0Bkg 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 17 f_+0 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 18 c_0 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 19 c_+ 1.00000e+00 1.00000e-02 0.00000e+00 2.00000e+00 20 D* slope 1.14500e+00 1.00000e-02 0.00000e+00 2.00000e+00 21 D* R1 1.39600e+00 1.00000e-02 0.00000e+00 2.00000e+00 22 D* R2 8.85000e-01 1.00000e-02 0.00000e+00 2.00000e+00 ********** ** 2 **SET PRINT 0 ********** ********** ** 3 **FIX 1 ********** ********** ** 4 **FIX 4 ********** ********** ** 5 **FIX 8 ********** ********** ** 6 **FIX 9 ********** ********** ** 7 **FIX 10 ********** ********** ** 8 **FIX 11 ********** ********** ** 9 **FIX 12 ********** ********** ** 10 **FIX 18 ********** ********** ** 11 **FIX 19 ********** ********** ** 12 **MIGRAD ********** MIGRAD MINIMIZATION HAS CONVERGED. MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX. FCN=415.674 FROM MIGRAD STATUS=CONVERGED 608 CALLS 609 TOTAL EDM=1.92368e-05 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 dummy 1.00000e+00 fixed 2 D slope 1.43855e+00 5.82728e-02 4.92415e-04 -1.37202e-02 3 Br_BpD 2.64007e-02 6.14085e-04 8.32388e-05 2.55440e-01 4 Br_B0D 2.10000e-02 fixed 5 Br_BpD* 6.49857e-02 2.17009e-03 1.26916e-04 1.04558e-02 6 Br_B0D* 4.71940e-02 2.43180e-03 1.61602e-04 8.52579e-03 7 BpDDPiAll 2.12912e-02 1.04109e-03 1.36078e-04 5.76851e-02 8 f_BpDD12 7.67000e-01 fixed 9 f_BpDPiDD12 9.11000e-01 fixed 10 f_BpDSPiDD12 9.92000e-01 fixed 11 f_BpDPiDD0 4.00000e-01 fixed 12 f_BpDSPiDDP1 1.00000e-01 fixed 13 UncorrDL 6.79853e-01 4.38716e-02 3.14008e-04 3.93360e-02 14 UncorrCL 9.89314e-01 6.85844e-02 4.59460e-04 3.94529e-02 15 CascL 8.97480e-01 9.89461e-02 6.52221e-04 2.19144e-02 16 f_DchD0Bkg 1.05464e+00 6.97769e-02 5.71320e-04 2.69478e-02 17 f_+0 1.05329e+00 1.98349e-02 6.01842e-05 2.08834e-03 18 c_0 1.00000e+00 fixed 19 c_+ 1.00000e+00 fixed 20 D* slope 1.30677e+00 6.82556e-02 2.52250e-04 -2.40246e-02 21 D* R1 1.07779e+00 7.00046e-02 3.17982e-04 -2.95226e-02 22 D* R2 8.40443e-01 7.61448e-02 2.28982e-04 -9.71426e-02 ********** ** 13 **SHOW COVARIANCE ********** EXTERNAL ERROR MATRIX. NDIM= 22 NPAR= 13 ERR DEF=1 ELEMENTS ABOVE DIAGONAL ARE NOT PRINTED. 3.400e-03 1.593e-05 3.771e-07 1.525e-06 3.670e-07 4.713e-06 -6.205e-05 -8.525e-07 -2.429e-06 5.918e-06 1.383e-05 1.492e-07 1.864e-07 -1.160e-06 1.084e-06 1.348e-04 1.075e-06 1.163e-06 -9.303e-06 -8.042e-07 1.926e-03 -8.050e-04 -4.766e-06 1.337e-05 2.880e-05 -1.960e-07 -6.417e-04 4.711e-03 3.192e-04 7.513e-06 1.831e-05 1.835e-05 -4.914e-05 7.977e-05 -2.580e-03 9.823e-03 -2.995e-04 -6.228e-06 2.655e-05 -3.063e-05 -1.572e-06 -5.925e-04 -8.537e-04 -9.313e-04 4.877e-03 1.377e-05 -2.992e-06 -3.289e-05 4.326e-06 -3.134e-06 -1.282e-04 -2.613e-04 -2.217e-04 1.111e-04 3.935e-04 -1.389e-03 -6.534e-07 7.622e-06 3.551e-05 -1.807e-05 -3.648e-04 1.096e-03 -2.643e-04 -1.577e-05 1.583e-04 4.667e-03 -6.651e-04 3.039e-06 -2.648e-06 1.693e-05 -1.661e-05 3.536e-04 1.505e-03 -6.154e-04 -4.135e-04 1.655e-04 3.698e-03 -6.651e-04 3.039e-06 -2.648e-06 1.693e-05 -1.661e-05 3.536e-04 1.505e-03 -6.154e-04 -4.135e-04 1.655e-04 3.698e-03 4.909e-03 1.274e-03 4.065e-06 -1.629e-05 -4.193e-05 1.898e-05 -2.478e-04 -1.369e-03 8.241e-04 1.483e-04 -6.871e-05 -4.677e-03 1.274e-03 4.065e-06 -1.629e-05 -4.193e-05 1.898e-05 -2.478e-04 -1.369e-03 8.241e-04 1.483e-04 -6.871e-05 -4.677e-03 -4.468e-03 5.810e-03 PARAMETER CORRELATION COEFFICIENTS NO. GLOBAL 2 3 5 6 7 13 14 15 16 17 20 21 22 2 0.64847 1.000 0.445 0.012 -0.437 0.228 0.053 -0.201 0.055 -0.074 0.012 -0.349 -0.163 0.287 3 0.80043 0.445 1.000 0.275 -0.571 0.233 0.040 -0.113 0.123 -0.145 -0.246 -0.016 0.071 0.087 5 0.96007 0.012 0.275 1.000 -0.460 0.082 0.012 0.090 0.085 0.175 -0.764 0.051 -0.017 -0.098 6 0.94303 -0.437 -0.571 -0.460 1.000 -0.458 -0.087 0.173 0.076 -0.180 0.090 0.214 0.099 -0.226 7 0.84375 0.228 0.233 0.082 -0.458 1.000 -0.018 -0.003 -0.476 -0.022 -0.152 -0.254 -0.228 0.239 13 0.73318 0.053 0.040 0.012 -0.087 -0.018 1.000 -0.213 0.018 -0.193 -0.147 -0.122 0.115 -0.074 14 0.74098 -0.201 -0.113 0.090 0.173 -0.003 -0.213 1.000 -0.379 -0.178 -0.192 0.234 0.313 -0.262 15 0.75412 0.055 0.123 0.085 0.076 -0.476 0.018 -0.379 1.000 -0.135 -0.113 -0.039 -0.089 0.109 16 0.57288 -0.074 -0.145 0.175 -0.180 -0.022 -0.193 -0.178 -0.135 1.000 0.080 -0.003 -0.085 0.028 17 0.95268 0.012 -0.246 -0.764 0.090 -0.152 -0.147 -0.192 -0.113 0.080 1.000 0.117 0.119 -0.045 20 0.93576 -0.349 -0.016 0.051 0.214 -0.254 -0.122 0.234 -0.039 -0.003 0.117 1.000 0.773 -0.898 21 0.89101 -0.163 0.071 -0.017 0.099 -0.228 0.115 0.313 -0.089 -0.085 0.119 0.773 1.000 -0.837 22 0.95470 0.287 0.087 -0.098 -0.226 0.239 -0.074 -0.262 0.109 0.028 -0.045 -0.898 -0.837 1.000 ********** ** 14 **CALL FCN 2 ********** rho2 = 1.43855 oldDlnuBp = 0.268702 oldDlnuB0 = 0.26542 newDlnuBp = 0.153524 newDlnuB0 = 0.152143 RNBp = 1.75023 RNB0 = 1.74454 rho2Ds = 1.30677 oldDstarBp = 0.865978 oldDstarB0 = 0.860858 newDstarBp = 0.749881 newDstarB0 = 0.74539 RNDstarBp = 1.15482 RNDstarB0 = 1.15491 Number of bins into chisq = 267 from D0 = 138 from Dch = 129 chisq_D0 = 272.151 chisq_Dch = 141.901 chisq = 415.674 ######## Now printout information ######## f_DD2 = 2.1 SlopeDstar = 1.30677 +- 0.0682556 (5.22323 %) R1Dstar = 1.07779 +- 0.0700046 (6.49521 %) R2Dstar = 0.840443 +- 0.0761448 (9.06008 %) Slope = 1.43855 +- 0.0582728 (4.05081 %) BpDlnu = 0.0264007 +- 0.000614085 B0Dlnu = 0.0246505 +- 0 BpDstar = 0.0649857 +- 0.00217009 B0Dstar = 0.047194 +- 0.0024318 BpDDPi = 0.0212912 +- 0.00104109 B0DDPi = 0.0198797 fDD12 = 0.767 +- 0 fDPiDD12 = 0.59 +- 0 fDSPiDD12 = 0.65 +- 0 fDPiDD0 = 0.4 +- 0 fDSPiDDP1 = 0.1 +- 0 UncorDL = 0.679853 +- 0.0438716 UncorCL = 0.989314 +- 0.0685844 CascL = 0.89748 +- 0.0989461 fDchD0Bkg = 1.05464 +- 0.0697769 fp0 = 1.05329 +- 0.0198349 C0 = 1 +- 0 Cp = 1 +- 0 TotBp = 0.112678 BpDlnu = 23.4303 % BpDstar = 57.674 % BpDDPi = 18.8957 % TotB0 = 0.0917243 B0Dlnu = 26.8746 % B0Dstar = 51.4521 % B0DDPi = 21.6733 % ************************************ ** Slopes ** ************************************ D* slope = 1.30677 +- 0.0682556 (5.22323 %) D slope = 1.43855 +- 0.0582728 (4.05081 %) ************************************ ** Normalised Branching Fractions ** ************************************ B(Bp->Dlnu) = 0.0256328 +- 0.000596223 (2.32602 %) B(Bp->D*lnu) = 0.0630954 +- 0.00210697 (3.33934 %) B(Bp->D**/D(*)Pilnu) = 0.0206719 +- 0.00101081 (4.88978 %) Bp sum = 0.1094 ************************* B(B0->Dlnu) = 0.0274517 +- 0 (0 %) B(B0->D*lnu) = 0.052557 +- 0.00270815 (5.15278 %) B(B0->D**/D(*)Pilnu) = 0.0221388 +- 0 (0 %) B0 sum = 0.102148 ************************************ ** Ratios of branching fractions ** ************************************ f_BpDD12 = 0.767 +- 0 (0 %) f_BpDPiDD12 = 0.59 +- 0 (0 %) f_BpDSPiDD12 = 0.65 +- 0 (0 %) f_BpDPiDD0 = 0.4 +- 0 (0 %) f_BpDSPiDDP1 = 0.1 +- 0 (0 %) ************************************ ** Others ** ************************************ UncorDL background = 0.679853 +- 0.0438716 (6.4531 %) UncorCL background = 0.989314 +- 0.0685844 (6.93253 %) CascL background = 0.89748 +- 0.0989461 (11.0249 %) Dch/D0 background ratio = 1.05464 +- 0.0697769 (6.61615 %) ************************* f_+0 = 1.05329 +- 0.0198349 (1.88313 %) c_0 = 1 +- 0 (0 %) c_+ = 1 +- 0 (0 %) ************************* Time at the end of job = 743.46 seconds