root [2] .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 = 420.796 chisq_Dch = 306.469 chisq = 728.354 FCN=728.354 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 6 ********** ********** ** 6 **FIX 8 ********** ********** ** 7 **FIX 9 ********** ********** ** 8 **FIX 10 ********** ********** ** 9 **FIX 11 ********** ********** ** 10 **FIX 12 ********** ********** ** 11 **FIX 18 ********** ********** ** 12 **FIX 19 ********** ********** ** 13 **MIGRAD ********** MIGRAD MINIMIZATION HAS CONVERGED. MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX. FCN=428.47 FROM MIGRAD STATUS=CONVERGED 619 CALLS 620 TOTAL EDM=5.32462e-08 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.37333e+00 6.18831e-02 5.29986e-04 -1.39566e-03 3 Br_BpD 2.52901e-02 5.23545e-04 8.42610e-05 -4.45113e-03 4 Br_B0D 2.10000e-02 fixed 5 Br_BpD* 5.89254e-02 1.15179e-03 7.48617e-05 -2.02887e-02 6 Br_B0D* 4.50000e-02 fixed 7 BpDDPiAll 2.00305e-02 9.64475e-04 1.39034e-04 2.61196e-03 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.70852e-01 4.36553e-02 3.17988e-04 -3.58868e-03 14 UncorrCL 1.00494e+00 6.74145e-02 4.63242e-04 1.05627e-03 15 CascL 9.01679e-01 9.84228e-02 6.58408e-04 2.17691e-03 16 f_DchD0Bkg 1.00225e+00 6.67247e-02 5.74053e-04 3.14362e-03 17 f_+0 1.08400e+00 1.78168e-02 6.56229e-05 -2.38827e-02 18 c_0 1.00000e+00 fixed 19 c_+ 1.00000e+00 fixed 20 D* slope 1.33600e+00 6.86967e-02 2.51852e-04 -2.76023e-03 21 D* R1 1.09586e+00 7.20870e-02 3.28438e-04 -1.37943e-03 22 D* R2 8.11332e-01 7.78484e-02 2.32140e-04 -5.39756e-03 ********** ** 14 **SHOW COVARIANCE ********** EXTERNAL ERROR MATRIX. NDIM= 22 NPAR= 12 ERR DEF=1 ELEMENTS ABOVE DIAGONAL ARE NOT PRINTED. 3.835e-03 1.195e-05 2.741e-07 -3.046e-05 -1.806e-07 1.327e-06 8.831e-06 3.650e-08 -4.128e-07 9.304e-07 1.171e-04 4.309e-07 -4.334e-06 -1.536e-06 1.907e-03 -8.316e-04 -3.479e-06 1.924e-05 1.518e-06 -6.283e-04 4.552e-03 3.800e-04 7.747e-06 1.698e-05 -4.776e-05 6.092e-05 -2.581e-03 9.719e-03 -6.330e-04 -1.049e-05 2.894e-06 -7.393e-06 -6.137e-04 -6.002e-04 -8.775e-04 4.459e-03 1.871e-04 -2.469e-07 -1.516e-05 -9.954e-08 -1.027e-04 -2.999e-04 -2.147e-04 2.223e-04 3.175e-04 -1.450e-03 2.242e-06 2.084e-05 -1.546e-05 -3.115e-04 1.030e-03 -2.569e-04 1.409e-04 8.069e-05 4.728e-03 -6.949e-04 5.224e-06 6.872e-06 -1.529e-05 4.064e-04 1.444e-03 -6.115e-04 -2.482e-04 1.166e-04 3.868e-03 5.206e-03 1.361e-03 1.121e-06 -2.793e-05 1.667e-05 -2.937e-04 -1.309e-03 8.338e-04 -3.474e-05 -1.725e-06 -4.850e-03 -4.742e-03 1.361e-03 1.121e-06 -2.793e-05 1.667e-05 -2.937e-04 -1.309e-03 8.338e-04 -3.474e-05 -1.725e-06 -4.850e-03 -4.742e-03 6.073e-03 PARAMETER CORRELATION COEFFICIENTS NO. GLOBAL 2 3 5 7 13 14 15 16 17 20 21 22 2 0.61997 1.000 0.368 -0.427 0.148 0.043 -0.199 0.062 -0.153 0.170 -0.341 -0.156 0.282 3 0.72253 0.368 1.000 -0.300 0.072 0.019 -0.098 0.150 -0.300 -0.026 0.062 0.138 0.027 5 0.94877 -0.427 -0.300 1.000 -0.372 -0.086 0.248 0.149 0.038 -0.738 0.263 0.083 -0.311 7 0.82158 0.148 0.072 -0.372 1.000 -0.036 0.023 -0.502 -0.115 -0.006 -0.233 -0.220 0.222 13 0.73368 0.043 0.019 -0.086 -0.036 1.000 -0.213 0.014 -0.210 -0.132 -0.104 0.129 -0.086 14 0.73456 -0.199 -0.098 0.248 0.023 -0.213 1.000 -0.388 -0.133 -0.249 0.222 0.297 -0.249 15 0.75407 0.062 0.150 0.149 -0.502 0.014 -0.388 1.000 -0.133 -0.122 -0.038 -0.086 0.109 16 0.52744 -0.153 -0.300 0.038 -0.115 -0.210 -0.133 -0.133 1.000 0.187 0.031 -0.052 -0.007 17 0.93191 0.170 -0.026 -0.738 -0.006 -0.132 -0.249 -0.122 0.187 1.000 0.066 0.091 -0.001 20 0.94009 -0.341 0.062 0.263 -0.233 -0.104 0.222 -0.038 0.031 0.066 1.000 0.780 -0.905 21 0.89415 -0.156 0.138 0.083 -0.220 0.129 0.297 -0.086 -0.052 0.091 0.780 1.000 -0.843 22 0.95731 0.282 0.027 -0.311 0.222 -0.086 -0.249 0.109 -0.007 -0.001 -0.905 -0.843 1.000 ********** ** 15 **CALL FCN 2 ********** rho2 = 1.37333 oldDlnuBp = 0.268702 oldDlnuB0 = 0.26542 newDlnuBp = 0.163465 newDlnuB0 = 0.16193 RNBp = 1.64379 RNB0 = 1.6391 rho2Ds = 1.336 oldDstarBp = 0.865978 oldDstarB0 = 0.860858 newDstarBp = 0.749016 newDstarB0 = 0.744551 RNDstarBp = 1.15615 RNDstarB0 = 1.15621 Number of bins into chisq = 267 from D0 = 138 from Dch = 129 chisq_D0 = 276.733 chisq_Dch = 144.932 chisq = 428.47 ######## Now printout information ######## f_DD2 = 2.1 SlopeDstar = 1.336 +- 0.0686967 (5.14197 %) R1Dstar = 1.09586 +- 0.072087 (6.57813 %) R2Dstar = 0.811332 +- 0.0778484 (9.59513 %) Slope = 1.37333 +- 0.0618831 (4.50605 %) BpDlnu = 0.0252901 +- 0.000523545 B0Dlnu = 0.0236136 +- 0 BpDstar = 0.0589254 +- 0.00115179 B0Dstar = 0.055019 +- 0 BpDDPi = 0.0200305 +- 0.000964475 B0DDPi = 0.0187026 fDD12 = 0.767 +- 0 fDPiDD12 = 0.59 +- 0 fDSPiDD12 = 0.65 +- 0 fDPiDD0 = 0.4 +- 0 fDSPiDDP1 = 0.1 +- 0 UncorDL = 0.670852 +- 0.0436553 UncorCL = 1.00494 +- 0.0674145 CascL = 0.901679 +- 0.0984228 fDchD0Bkg = 1.00225 +- 0.0667247 fp0 = 1.084 +- 0.0178168 C0 = 1 +- 0 Cp = 1 +- 0 TotBp = 0.104246 BpDlnu = 24.2601 % BpDstar = 56.5253 % BpDDPi = 19.2146 % TotB0 = 0.0973352 B0Dlnu = 24.2601 % B0Dstar = 56.5253 % B0DDPi = 19.2146 % ************************************ ** Slopes ** ************************************ D* slope = 1.336 +- 0.0686967 (5.14197 %) D slope = 1.37333 +- 0.0618831 (4.50605 %) ************************************ ** Normalised Branching Fractions ** ************************************ B(Bp->Dlnu) = 0.0265405 +- 0.000549429 (2.07015 %) B(Bp->D*lnu) = 0.0618387 +- 0.00120873 (1.95465 %) B(Bp->D**/D(*)Pilnu) = 0.0210208 +- 0.00101216 (4.81504 %) Bp sum = 0.1094 ************************* B(B0->Dlnu) = 0.024781 +- 0 (0 %) B(B0->D*lnu) = 0.0577392 +- 0 (0 %) B(B0->D**/D(*)Pilnu) = 0.0196273 +- 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.670852 +- 0.0436553 (6.50745 %) UncorCL background = 1.00494 +- 0.0674145 (6.70832 %) CascL background = 0.901679 +- 0.0984228 (10.9155 %) Dch/D0 background ratio = 1.00225 +- 0.0667247 (6.65749 %) ************************* f_+0 = 1.084 +- 0.0178168 (1.64362 %) c_0 = 1 +- 0 (0 %) c_+ = 1 +- 0 (0 %) ************************* Time at the end of job = 767.09 seconds