diff --git a/.ipynb_checkpoints/raremodel-nb-checkpoint.ipynb b/.ipynb_checkpoints/raremodel-nb-checkpoint.ipynb index 05cd12c..6f23704 100644 --- a/.ipynb_checkpoints/raremodel-nb-checkpoint.ipynb +++ b/.ipynb_checkpoints/raremodel-nb-checkpoint.ipynb @@ -795,18 +795,7 @@ "cell_type": "code", "execution_count": 11, "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "" - ] - }, - "execution_count": 11, - "metadata": {}, - "output_type": "execute_result" - } - ], + "outputs": [], "source": [ "# total_f_fit.normalization(obs_fit)" ] @@ -1590,65 +1579,74 @@ }, { "cell_type": "code", - "execution_count": 38, + "execution_count": 54, "metadata": {}, "outputs": [], "source": [ - "def reset_param_values(): \n", - " jpsi_m.set_value(jpsi_mass)\n", - " jpsi_s.set_value(jpsi_scale)\n", - " jpsi_p.set_value(jpsi_phase)\n", - " jpsi_w.set_value(jpsi_width)\n", - " psi2s_m.set_value(psi2s_mass)\n", - " psi2s_s.set_value(psi2s_scale)\n", - " psi2s_p.set_value(psi2s_phase)\n", - " psi2s_w.set_value(psi2s_width)\n", - " p3770_m.set_value(p3770_mass)\n", - " p3770_s.set_value(p3770_scale)\n", - " p3770_p.set_value(p3770_phase)\n", - " p3770_w.set_value(p3770_width)\n", - " p4040_m.set_value(p4040_mass)\n", - " p4040_s.set_value(p4040_scale)\n", - " p4040_p.set_value(p4040_phase)\n", - " p4040_w.set_value(p4040_width)\n", - " p4160_m.set_value(p4160_mass)\n", - " p4160_s.set_value(p4160_scale)\n", - " p4160_p.set_value(p4160_phase)\n", - " p4160_w.set_value(p4160_width)\n", - " p4415_m.set_value(p4415_mass)\n", - " p4415_s.set_value(p4415_scale)\n", - " p4415_p.set_value(p4415_phase)\n", - " p4415_w.set_value(p4415_width)\n", - " rho_m.set_value(rho_mass)\n", - " rho_s.set_value(rho_scale)\n", - " rho_p.set_value(rho_phase)\n", - " rho_w.set_value(rho_width)\n", - " omega_m.set_value(omega_mass)\n", - " omega_s.set_value(omega_scale)\n", - " omega_p.set_value(omega_phase)\n", - " omega_w.set_value(omega_width)\n", - " phi_m.set_value(phi_mass)\n", - " phi_s.set_value(phi_scale)\n", - " phi_p.set_value(phi_phase)\n", - " phi_w.set_value(phi_width)\n", - " Dstar_m.set_value(Dstar_mass)\n", - " DDstar_s.set_value(0.0)\n", - " DDstar_p.set_value(0.0)\n", - " D_m.set_value(D_mass)\n", - " Dbar_m.set_value(Dbar_mass)\n", - " Dbar_s.set_value(0.0)\n", - " Dbar_p.set_value(0.0)\n", - " tau_m.set_value(pdg['tau_M'])\n", - " Ctt.set_value(0.0)\n", - " b0_0.set_value(0.292)\n", - " b0_1.set_value(0.281)\n", - " b0_2.set_value(0.150)\n", - " bplus_0.set_value(0.466)\n", - " bplus_1.set_value(-0.885)\n", - " bplus_2.set_value(-0.213)\n", - " bT_0.set_value(0.460)\n", - " bT_1.set_value(-1.089)\n", - " bT_2.set_value(-1.114)\n" + "param_values_dic = {\n", + " 'jpsi_m': jpsi_mass,\n", + " 'jpsi_s': jpsi_scale,\n", + " 'jpsi_p': jpsi_phase,\n", + " 'jpsi_w': jpsi_width,\n", + " 'psi2s_m': psi2s_mass,\n", + " 'psi2s_s': psi2s_scale,\n", + " 'psi2s_p': psi2s_phase,\n", + " 'psi2s_w': psi2s_width,\n", + " 'p3770_m': p3770_mass,\n", + " 'p3770_s': p3770_scale,\n", + " 'p3770_p': p3770_phase,\n", + " 'p3770_w': p3770_width,\n", + " 'p4040_m': p4040_mass,\n", + " 'p4040_s': p4040_scale,\n", + " 'p4040_p': p4040_phase,\n", + " 'p4040_w': p4040_width,\n", + " 'p4160_m': p4160_mass,\n", + " 'p4160_s': p4160_scale,\n", + " 'p4160_p': p4160_phase,\n", + " 'p4160_w': p4160_width,\n", + " 'p4415_m': p4415_mass,\n", + " 'p4415_s': p4415_scale,\n", + " 'p4415_p': p4415_phase,\n", + " 'p4415_w': p4415_width,\n", + " 'rho_m': rho_mass,\n", + " 'rho_s': rho_scale,\n", + " 'rho_p': rho_phase,\n", + " 'rho_w': rho_width,\n", + " 'omega_m': omega_mass,\n", + " 'omega_s': omega_scale,\n", + " 'omega_p': omega_phase,\n", + " 'omega_w': omega_width,\n", + " 'phi_m': phi_mass,\n", + " 'phi_s': phi_scale,\n", + " 'phi_p': phi_phase,\n", + " 'phi_w': phi_width,\n", + " 'Dstar_m': Dstar_mass,\n", + " 'DDstar_s': 0.0,\n", + " 'DDstar_p': 0.0,\n", + " 'D_m': D_mass,\n", + " 'Dbar_m': Dbar_mass,\n", + " 'Dbar_s': 0.0,\n", + " 'Dbar_p': 0.0,\n", + " 'tau_m': pdg['tau_M'],\n", + " 'Ctt': 0.0,\n", + " 'b0_0': 0.292,\n", + " 'b0_1': 0.281,\n", + " 'b0_2': 0.150,\n", + " 'bplus_0': 0.466,\n", + " 'bplus_1': -0.885,\n", + " 'bplus_2': -0.213,\n", + " 'bT_0': 0.460,\n", + " 'bT_1': -1.089,\n", + " 'bT_2': -1.114}\n", + "\n", + "\n", + "def reset_param_values(variation = 0.05):\n", + " for param in total_f_fit.get_dependents():\n", + " if param.floating:\n", + " param.set_value(param_values_dic[param.name] + random.uniform(-1, 1)*param_values_dic[param.name]* variation)\n", + "# print(param.name)\n", + "# for param in totalf.get_dependents():\n", + "# param.set_value()" ] }, { @@ -1953,7 +1951,7 @@ }, { "cell_type": "code", - "execution_count": 47, + "execution_count": 48, "metadata": { "scrolled": false }, @@ -1964,29 +1962,360 @@ "text": [ "[0.15430335 0.53452248 0.74001287 0.89973541 1.03509834 1.15470054\n", " 1.26302735 1.36277029 1.45569489 1.5430335 ]\n", - "WARNING:tensorflow:From C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\zfit\\core\\sample.py:163: to_int32 (from tensorflow.python.ops.math_ops) is deprecated and will be removed in a future version.\n", - "Instructions for updating:\n", - "Use tf.cast instead.\n", - "Step: 1/10\n", - "Current Ctt: 0.1543033499620919\n", + "Time Taken: 2 s\n", + "Predicted time left: 8 s\n", + "Time Taken: 2 s\n", + "Predicted time left: 8 s\n", + "Time Taken: 2 s\n", + "Predicted time left: 4 s\n", + "Time Taken: 2 s\n", + "Predicted time left: 4 s\n", + "Time Taken: 2 s\n", + "Predicted time left: 3 s\n", + "Time Taken: 2 s\n", + "Predicted time left: 3 s\n", + "Time Taken: 2 s\n", + "Predicted time left: 2 s\n", + "Time Taken: 2 s\n", + "Predicted time left: 2 s\n", + "Time Taken: 2 s\n", + "Predicted time left: 1 \n", + "Time Taken: 2 s\n", + "Predicted time left: 1 \n", + "Step: 6/10\n", + "Current Ctt: 1.1547005383792517\n", "Ctt floating: True\n", - "Toy 0/1 - Fit 0/5\n" + "Toy 0/1 - Fit 0/5\n", + "------------------------------------------------------------------\n", + "| FCN = 2.967E+05 | Ncalls=1243 (1243 total) |\n", + "| EDM = 0.000701 (Goal: 5E-06) | up = 0.5 |\n", + "------------------------------------------------------------------\n", + "| Valid Min. | Valid Param. | Above EDM | Reached call limit |\n", + "------------------------------------------------------------------\n", + "| True | True | False | False |\n", + "------------------------------------------------------------------\n", + "| Hesse failed | Has cov. | Accurate | Pos. def. | Forced |\n", + "------------------------------------------------------------------\n", + "| False | True | True | True | False |\n", + "------------------------------------------------------------------\n", + "Function minimum: 296683.1850103811\n", + "----------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "----------------------------------------------------------------------------------------------\n", + "| 0 | omega_s | 9 | 4 | | | 4.19232 | 9.40768 | |\n", + "| 1 | bplus_0 | 0.512 | 0.022 | | | -2 | 2 | |\n", + "| 2 | p3770_p | 2.78 | 0.11 | | |-6.28319 | 6.28319 | |\n", + "| 3 | Dbar_s | -0.30 | 0.04 | | | -0.3 | 0.3 | |\n", + "| 4 | p4160_p | 3.82 | 0.10 | | |-6.28319 | 6.28319 | |\n", + "| 5 | DDstar_p | 5.2 | 0.4 | | |-6.28319 | 6.28319 | |\n", + "| 6 | Dbar_p | 5.6 | 0.4 | | |-6.28319 | 6.28319 | |\n", + "| 7 | p4040_s | 1.20 | 0.18 | | |0.00501244| 2.01499 | |\n", + "| 8 | rho_p | 5.95 | 0.30 | | |-6.28319 | 6.28319 | |\n", + "| 9 | bplus_1 | -0.96 | 0.05 | | | -2 | 2 | |\n", + "| 10| jpsi_p | 1.681 | 0.025 | | |-6.28319 | 6.28319 | |\n", + "| 11| omega_p | -5.33 | 0.22 | | |-6.28319 | 6.28319 | |\n", + "| 12| phi_p | -5.66 | 0.25 | | |-6.28319 | 6.28319 | |\n", + "| 13| p4040_p | 3.29 | 0.15 | | |-6.28319 | 6.28319 | |\n", + "| 14| p4415_p | -2.47 | 0.14 | | |-6.28319 | 6.28319 | |\n", + "| 15| Ctt | 1.61 | 0.14 | | | -1.7 | 1.7 | |\n", + "| 16| bplus_2 | -0.37 | 0.08 | | | -2 | 2 | |\n", + "| 17| DDstar_s | -0.30 | 0.05 | | | -0.3 | 0.3 | |\n", + "| 18| p4415_s | 1.55 | 0.20 | | |0.126447 | 2.35355 | |\n", + "| 19| p4160_s | 2.36 | 0.17 | | | 0.71676 | 3.68324 | |\n", + "| 20| psi2s_p | 2.02 | 0.03 | | |-6.28319 | 6.28319 | |\n", + "| 21| rho_s | 1.4 | 0.4 | | |0.0253049| 2.0747 | |\n", + "| 22| p3770_s | 2.50 | 0.25 | | |0.918861 | 4.08114 | |\n", + "| 23| phi_s | 19.8 | 1.2 | | | 14.8182 | 23.5818 | |\n", + "----------------------------------------------------------------------------------------------\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "| | omega_s bplus_0 p3770_p Dbar_s p4160_p DDstar_p Dbar_p p4040_s rho_p bplus_1 jpsi_p omega_p phi_p p4040_p p4415_p Ctt bplus_2 DDstar_s p4415_s p4160_s psi2s_p rho_s p3770_s phi_s |\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "| omega_s | 1.000 -0.001 -0.000 -0.000 0.000 0.000 0.000 -0.000 -0.002 0.001 0.000 -0.031 0.001 0.000 0.000 0.000 0.001 0.000 -0.000 0.000 0.000 -0.007 0.000 -0.000 |\n", + "| bplus_0 | -0.001 1.000 0.069 0.004 0.052 -0.039 -0.035 0.004 0.030 -0.824 0.054 0.122 0.004 0.054 0.042 0.011 -0.077 0.003 0.014 0.037 0.015 0.218 0.023 -0.016 |\n", + "| p3770_p | -0.000 0.069 1.000 -0.058 0.288 -0.295 -0.411 -0.058 0.037 0.167 0.250 0.025 -0.006 0.185 0.126 -0.231 -0.085 -0.046 -0.011 -0.030 -0.081 0.039 -0.284 -0.018 |\n", + "| Dbar_s | -0.000 0.004 -0.058 1.000 -0.066 -0.003 0.026 0.027 0.005 0.012 -0.072 0.002 -0.002 -0.055 -0.043 -0.017 0.016 -0.006 0.017 0.009 -0.042 0.004 -0.015 -0.002 |\n", + "| p4160_p | 0.000 0.052 0.288 -0.066 1.000 -0.334 -0.182 -0.505 0.047 0.085 0.162 0.019 -0.012 0.090 0.317 -0.278 0.044 -0.104 -0.018 -0.208 -0.075 0.027 -0.142 -0.025 |\n", + "| DDstar_p | 0.000 -0.039 -0.295 -0.003 -0.334 1.000 -0.185 0.121 -0.017 -0.371 -0.160 -0.034 -0.007 -0.171 -0.102 0.084 0.284 0.052 0.020 0.065 -0.060 -0.065 -0.029 -0.008 |\n", + "| Dbar_p | 0.000 -0.035 -0.411 0.026 -0.182 -0.185 1.000 0.026 -0.046 -0.138 -0.379 -0.020 0.004 -0.084 -0.046 0.225 0.037 -0.013 -0.052 0.037 -0.146 -0.027 0.075 0.018 |\n", + "| p4040_s | -0.000 0.004 -0.058 0.027 -0.505 0.121 0.026 1.000 0.002 0.031 -0.015 0.001 -0.001 -0.212 -0.252 0.278 -0.074 0.008 0.130 -0.024 0.086 0.003 0.171 -0.001 |\n", + "| rho_p | -0.002 0.030 0.037 0.005 0.047 -0.017 -0.046 0.002 1.000 0.009 0.042 0.048 -0.237 0.054 0.063 0.080 0.239 0.002 -0.012 0.040 0.008 0.250 0.023 -0.170 |\n", + "| bplus_1 | 0.001 -0.824 0.167 0.012 0.085 -0.371 -0.138 0.031 0.009 1.000 0.183 -0.077 -0.010 0.075 0.001 -0.135 -0.255 0.013 0.143 0.046 0.086 -0.133 0.024 0.005 |\n", + "| jpsi_p | 0.000 0.054 0.250 -0.072 0.162 -0.160 -0.379 -0.015 0.042 0.183 1.000 0.019 -0.031 0.060 0.026 -0.370 -0.192 -0.104 0.011 -0.057 0.106 0.029 0.045 -0.050 |\n", + "| omega_p | -0.031 0.122 0.025 0.002 0.019 -0.034 -0.020 0.001 0.048 -0.077 0.019 1.000 0.087 0.016 0.012 0.004 -0.007 0.001 0.008 0.009 0.009 0.664 0.005 0.043 |\n", + "| phi_p | 0.001 0.004 -0.006 -0.002 -0.012 -0.007 0.004 -0.001 -0.237 -0.010 -0.031 0.087 1.000 -0.017 -0.019 -0.016 -0.075 -0.002 0.004 -0.013 -0.004 0.019 -0.009 0.729 |\n", + "| p4040_p | 0.000 0.054 0.185 -0.055 0.090 -0.171 -0.084 -0.212 0.054 0.075 0.060 0.016 -0.017 1.000 0.090 -0.109 0.013 -0.075 0.110 0.350 -0.170 0.023 -0.130 -0.032 |\n", + "| p4415_p | 0.000 0.042 0.126 -0.043 0.317 -0.102 -0.046 -0.252 0.063 0.001 0.026 0.012 -0.019 0.090 1.000 -0.040 0.148 -0.091 -0.182 -0.209 -0.080 0.014 -0.071 -0.036 |\n", + "| Ctt | 0.000 0.011 -0.231 -0.017 -0.278 0.084 0.225 0.278 0.080 -0.135 -0.370 0.004 -0.016 -0.109 -0.040 1.000 0.539 0.015 0.138 0.234 0.199 -0.005 0.052 -0.036 |\n", + "| bplus_2 | 0.001 -0.077 -0.085 0.016 0.044 0.284 0.037 -0.074 0.239 -0.255 -0.192 -0.007 -0.075 0.013 0.148 0.539 1.000 -0.000 -0.241 -0.021 -0.064 -0.041 -0.026 -0.131 |\n", + "| DDstar_s | 0.000 0.003 -0.046 -0.006 -0.104 0.052 -0.013 0.008 0.002 0.013 -0.104 0.001 -0.002 -0.075 -0.091 0.015 -0.000 1.000 0.015 -0.009 -0.050 0.003 -0.021 -0.001 |\n", + "| p4415_s | -0.000 0.014 -0.011 0.017 -0.018 0.020 -0.052 0.130 -0.012 0.143 0.011 0.008 0.004 0.110 -0.182 0.138 -0.241 0.015 1.000 0.324 0.051 0.018 0.065 0.008 |\n", + "| p4160_s | 0.000 0.037 -0.030 0.009 -0.208 0.065 0.037 -0.024 0.040 0.046 -0.057 0.009 -0.013 0.350 -0.209 0.234 -0.021 -0.009 0.324 1.000 -0.046 0.013 0.048 -0.023 |\n", + "| psi2s_p | 0.000 0.015 -0.081 -0.042 -0.075 -0.060 -0.146 0.086 0.008 0.086 0.106 0.009 -0.004 -0.170 -0.080 0.199 -0.064 -0.050 0.051 -0.046 1.000 0.016 -0.158 -0.006 |\n", + "| rho_s | -0.007 0.218 0.039 0.004 0.027 -0.065 -0.027 0.003 0.250 -0.133 0.029 0.664 0.019 0.023 0.014 -0.005 -0.041 0.003 0.018 0.013 0.016 1.000 0.007 -0.021 |\n", + "| p3770_s | 0.000 0.023 -0.284 -0.015 -0.142 -0.029 0.075 0.171 0.023 0.024 0.045 0.005 -0.009 -0.130 -0.071 0.052 -0.026 -0.021 0.065 0.048 -0.158 0.007 1.000 -0.015 |\n", + "| phi_s | -0.000 -0.016 -0.018 -0.002 -0.025 -0.008 0.018 -0.001 -0.170 0.005 -0.050 0.043 0.729 -0.032 -0.036 -0.036 -0.131 -0.001 0.008 -0.023 -0.006 -0.021 -0.015 1.000 |\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 3.525680379073008}), (, {'error': 0.02207235527349316}), (, {'error': 0.11343595775266557}), (, {'error': 0.03618310665932811}), (, {'error': 0.09949384100148517}), (, {'error': 0.3808783318160245}), (, {'error': 0.3535462990363385}), (, {'error': 0.1750004739296166}), (, {'error': 0.3007303858774737}), (, {'error': 0.04619582464412908}), (, {'error': 0.024960652656695004}), (, {'error': 0.22130711307575224}), (, {'error': 0.2461879166128753}), (, {'error': 0.14614016890582704}), (, {'error': 0.14435100437895199}), (, {'error': 0.1375287156270839}), (, {'error': 0.08217225999599143}), (, {'error': 0.04911032929419584}), (, {'error': 0.19713982348564196}), (, {'error': 0.17244193473434555}), (, {'error': 0.032470914854760835}), (, {'error': 0.352537831138418}), (, {'error': 0.2497690977962057}), (, {'error': 1.1844440265383138})])\n", + "Step: 6/10\n", + "Current Ctt: 1.1547005383792517\n", + "Ctt floating: True\n", + "Toy 0/1 - Fit 1/5\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "------------------------------------------------------------------\n", + "| FCN = 2.967E+05 | Ncalls=1255 (1255 total) |\n", + "| EDM = 0.0112 (Goal: 5E-06) | up = 0.5 |\n", + "------------------------------------------------------------------\n", + "| Valid Min. | Valid Param. | Above EDM | Reached call limit |\n", + "------------------------------------------------------------------\n", + "| False | True | True | False |\n", + "------------------------------------------------------------------\n", + "| Hesse failed | Has cov. | Accurate | Pos. def. | Forced |\n", + "------------------------------------------------------------------\n", + "| False | True | False | False | True |\n", + "------------------------------------------------------------------\n", + "Function minimum: 296679.8728610612\n", + "----------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "----------------------------------------------------------------------------------------------\n", + "| 0 | omega_s | 9 | 4 | | | 4.19232 | 9.40768 | |\n", + "| 1 | bplus_0 | 0.511 | 0.017 | | | -2 | 2 | |\n", + "| 2 | p3770_p | -3.34 | 0.17 | | |-6.28319 | 6.28319 | |\n", + "| 3 | Dbar_s | -0.30 | 0.46 | | | -0.3 | 0.3 | |\n", + "| 4 | p4160_p | 3.85 | 0.17 | | |-6.28319 | 6.28319 | |\n", + "| 5 | DDstar_p | -5.0 | 1.1 | | |-6.28319 | 6.28319 | |\n", + "| 6 | Dbar_p | 0.8 | 0.7 | | |-6.28319 | 6.28319 | |\n", + "| 7 | p4040_s | 1.13 | 0.19 | | |0.00501244| 2.01499 | |\n", + "| 8 | rho_p | -0.38 | 0.31 | | |-6.28319 | 6.28319 | |\n", + "| 9 | bplus_1 | -0.99 | 0.05 | | | -2 | 2 | |\n", + "| 10| jpsi_p | 4.69 | 0.04 | | |-6.28319 | 6.28319 | |\n", + "| 11| omega_p | 0.93 | 0.21 | | |-6.28319 | 6.28319 | |\n", + "| 12| phi_p | 0.61 | 0.23 | | |-6.28319 | 6.28319 | |\n", + "| 13| p4040_p | 3.46 | 0.21 | | |-6.28319 | 6.28319 | |\n", + "| 14| p4415_p | 3.80 | 0.18 | | |-6.28319 | 6.28319 | |\n", + "| 15| Ctt | 1.39 | 0.18 | | | -1.7 | 1.7 | |\n", + "| 16| bplus_2 | -0.32 | 0.11 | | | -2 | 2 | |\n", + "| 17| DDstar_s | -0.30 | 0.57 | | | -0.3 | 0.3 | |\n", + "| 18| p4415_s | 1.49 | 0.20 | | |0.126447 | 2.35355 | |\n", + "| 19| p4160_s | 2.38 | 0.17 | | | 0.71676 | 3.68324 | |\n", + "| 20| psi2s_p | 2.09 | 0.04 | | |-6.28319 | 6.28319 | |\n", + "| 21| rho_s | 1.3 | 0.3 | | |0.0253049| 2.0747 | |\n", + "| 22| p3770_s | 2.72 | 0.26 | | |0.918861 | 4.08114 | |\n", + "| 23| phi_s | 19.9 | 1.1 | | | 14.8182 | 23.5818 | |\n", + "----------------------------------------------------------------------------------------------\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "| | omega_s bplus_0 p3770_p Dbar_s p4160_p DDstar_p Dbar_p p4040_s rho_p bplus_1 jpsi_p omega_p phi_p p4040_p p4415_p Ctt bplus_2 DDstar_s p4415_s p4160_s psi2s_p rho_s p3770_s phi_s |\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "| omega_s | 1.000 -0.004 -0.004 0.003 -0.003 -0.005 0.003 0.002 -0.007 -0.003 -0.002 -0.122 0.003 -0.002 -0.001 0.000 0.006 -0.005 0.001 0.002 -0.001 -0.021 -0.001 -0.002 |\n", + "| bplus_0 | -0.004 1.000 0.402 -0.215 0.353 0.517 -0.204 -0.205 0.046 0.047 0.245 0.136 0.012 0.241 0.140 -0.005 -0.421 0.474 -0.051 -0.109 0.069 0.267 0.113 0.005 |\n", + "| p3770_p | -0.004 0.402 1.000 -0.384 0.738 0.778 -0.612 -0.415 -0.001 0.595 0.584 0.096 0.032 0.582 0.419 0.054 -0.523 0.756 -0.245 -0.233 0.214 0.193 0.093 0.040 |\n", + "| Dbar_s | 0.003 -0.215 -0.384 1.000 -0.449 -0.611 0.622 0.313 0.036 -0.351 -0.215 -0.051 -0.027 -0.396 -0.246 -0.324 0.393 -0.488 0.206 0.159 -0.182 -0.109 -0.242 -0.043 |\n", + "| p4160_p | -0.003 0.353 0.738 -0.449 1.000 0.732 -0.615 -0.622 0.001 0.529 0.609 0.087 0.030 0.596 0.552 0.016 -0.450 0.780 -0.234 -0.300 0.266 0.174 0.203 0.037 |\n", + "| DDstar_p | -0.005 0.517 0.778 -0.611 0.732 1.000 -0.673 -0.500 0.009 0.783 0.547 0.124 0.034 0.533 0.314 0.161 -0.675 0.827 -0.257 -0.313 0.219 0.252 0.243 0.048 |\n", + "| Dbar_p | 0.003 -0.204 -0.612 0.622 -0.615 -0.673 1.000 0.391 0.051 -0.346 -0.438 -0.050 -0.036 -0.586 -0.373 -0.421 0.423 -0.409 0.255 0.181 -0.398 -0.106 -0.332 -0.053 |\n", + "| p4040_s | 0.002 -0.205 -0.415 0.313 -0.622 -0.500 0.391 1.000 0.025 -0.328 -0.296 -0.050 -0.024 -0.417 -0.347 0.127 0.313 -0.417 0.256 0.161 -0.095 -0.105 -0.018 -0.037 |\n", + "| rho_p | -0.007 0.046 -0.001 0.036 0.001 0.009 0.051 0.025 1.000 0.030 -0.076 0.058 -0.245 0.005 0.027 0.058 0.192 0.016 0.007 0.045 -0.027 0.263 0.010 -0.163 |\n", + "| bplus_1 | -0.003 0.047 0.595 -0.351 0.529 0.783 -0.346 -0.328 0.030 1.000 0.422 0.078 0.017 0.357 0.194 -0.016 -0.618 0.715 -0.091 -0.201 0.129 0.165 0.159 0.031 |\n", + "| jpsi_p | -0.002 0.245 0.584 -0.215 0.609 0.547 -0.438 -0.296 -0.076 0.422 1.000 0.067 0.039 0.456 0.394 0.203 -0.328 0.685 -0.195 -0.171 0.285 0.140 0.220 0.060 |\n", + "| omega_p | -0.122 0.136 0.096 -0.051 0.087 0.124 -0.050 -0.050 0.058 0.078 0.067 1.000 0.077 0.059 0.038 0.009 -0.080 0.118 -0.017 -0.028 0.019 0.639 0.026 0.038 |\n", + "| phi_p | 0.003 0.012 0.032 -0.027 0.030 0.034 -0.036 -0.024 -0.245 0.017 0.039 0.077 1.000 0.019 0.007 -0.006 -0.080 0.035 -0.010 -0.023 0.016 0.017 0.005 0.689 |\n", + "| p4040_p | -0.002 0.241 0.582 -0.396 0.596 0.533 -0.586 -0.417 0.005 0.357 0.456 0.059 0.019 1.000 0.422 0.088 -0.326 0.554 -0.131 0.096 0.178 0.116 0.154 0.022 |\n", + "| p4415_p | -0.001 0.140 0.419 -0.246 0.552 0.314 -0.373 -0.347 0.027 0.194 0.394 0.038 0.007 0.422 1.000 0.073 -0.108 0.439 -0.251 -0.225 0.190 0.070 0.138 0.003 |\n", + "| Ctt | 0.000 -0.005 0.054 -0.324 0.016 0.161 -0.421 0.127 0.058 -0.016 0.203 0.009 -0.006 0.088 0.073 1.000 0.288 0.011 0.074 0.169 0.320 0.007 0.143 -0.018 |\n", + "| bplus_2 | 0.006 -0.421 -0.523 0.393 -0.450 -0.675 0.423 0.313 0.192 -0.618 -0.328 -0.080 -0.080 -0.326 -0.108 0.288 1.000 -0.569 0.042 0.221 -0.155 -0.183 -0.161 -0.131 |\n", + "| DDstar_s | -0.005 0.474 0.756 -0.488 0.780 0.827 -0.409 -0.417 0.016 0.715 0.685 0.118 0.035 0.554 0.439 0.011 -0.569 1.000 -0.220 -0.244 0.263 0.235 0.229 0.044 |\n", + "| p4415_s | 0.001 -0.051 -0.245 0.206 -0.234 -0.257 0.255 0.256 0.007 -0.091 -0.195 -0.017 -0.010 -0.131 -0.251 0.074 0.042 -0.220 1.000 0.358 -0.078 -0.034 -0.043 -0.015 |\n", + "| p4160_s | 0.002 -0.109 -0.233 0.159 -0.300 -0.313 0.181 0.161 0.045 -0.201 -0.171 -0.028 -0.023 0.096 -0.225 0.169 0.221 -0.244 0.358 1.000 -0.085 -0.063 -0.023 -0.039 |\n", + "| psi2s_p | -0.001 0.069 0.214 -0.182 0.266 0.219 -0.398 -0.095 -0.027 0.129 0.285 0.019 0.016 0.178 0.190 0.320 -0.155 0.263 -0.078 -0.085 1.000 0.041 0.014 0.023 |\n", + "| rho_s | -0.021 0.267 0.193 -0.109 0.174 0.252 -0.106 -0.105 0.263 0.165 0.140 0.639 0.017 0.116 0.070 0.007 -0.183 0.235 -0.034 -0.063 0.041 1.000 0.051 -0.015 |\n", + "| p3770_s | -0.001 0.113 0.093 -0.242 0.203 0.243 -0.332 -0.018 0.010 0.159 0.220 0.026 0.005 0.154 0.138 0.143 -0.161 0.229 -0.043 -0.023 0.014 0.051 1.000 0.004 |\n", + "| phi_s | -0.002 0.005 0.040 -0.043 0.037 0.048 -0.053 -0.037 -0.163 0.031 0.060 0.038 0.689 0.022 0.003 -0.018 -0.131 0.044 -0.015 -0.039 0.023 -0.015 0.004 1.000 |\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 3.6828714770844417}), (, {'error': 0.017213810473602598}), (, {'error': 0.1669872439856892}), (, {'error': 0.4587401965716338}), (, {'error': 0.16725710271948824}), (, {'error': 1.1070258307403715}), (, {'error': 0.6673753412101515}), (, {'error': 0.19273514016275872}), (, {'error': 0.31222268239435724}), (, {'error': 0.04914242208045405}), (, {'error': 0.03592824962744956}), (, {'error': 0.21419759833574847}), (, {'error': 0.22875082760152043}), (, {'error': 0.2097373545461556}), (, {'error': 0.17561131353278725}), (, {'error': 0.1765648280036929}), (, {'error': 0.10601597248571426}), (, {'error': 0.5671660697555896}), (, {'error': 0.19702036437940162}), (, {'error': 0.17344156266937172}), (, {'error': 0.03612138672594423}), (, {'error': 0.34707836809237436}), (, {'error': 0.25712641295716976}), (, {'error': 1.1117746399794743})])\n", + "Step: 6/10\n", + "Current Ctt: 1.1547005383792517\n", + "Ctt floating: True\n", + "Toy 0/1 - Fit 1/5\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "------------------------------------------------------------------\n", + "| FCN = 2.967E+05 | Ncalls=1475 (1475 total) |\n", + "| EDM = 0.00166 (Goal: 5E-06) | up = 0.5 |\n", + "------------------------------------------------------------------\n", + "| Valid Min. | Valid Param. | Above EDM | Reached call limit |\n", + "------------------------------------------------------------------\n", + "| False | True | True | False |\n", + "------------------------------------------------------------------\n", + "| Hesse failed | Has cov. | Accurate | Pos. def. | Forced |\n", + "------------------------------------------------------------------\n", + "| False | True | False | False | True |\n", + "------------------------------------------------------------------\n", + "Function minimum: 296679.8773376704\n", + "----------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "----------------------------------------------------------------------------------------------\n", + "| 0 | omega_s | 9 | 4 | | | 4.19232 | 9.40768 | |\n", + "| 1 | bplus_0 | 0.511 | 0.022 | | | -2 | 2 | |\n", + "| 2 | p3770_p | 2.94 | 0.12 | | |-6.28319 | 6.28319 | |\n", + "| 3 | Dbar_s | 0.30 | 0.38 | | | -0.3 | 0.3 | |\n", + "| 4 | p4160_p | -2.44 | 0.11 | | |-6.28319 | 6.28319 | |\n", + "| 5 | DDstar_p | 4.4 | 0.7 | | |-6.28319 | 6.28319 | |\n", + "| 6 | Dbar_p | 4.0 | 0.6 | | |-6.28319 | 6.28319 | |\n", + "| 7 | p4040_s | 1.13 | 0.18 | | |0.00501244| 2.01499 | |\n", + "| 8 | rho_p | 5.87 | 0.31 | | |-6.28319 | 6.28319 | |\n", + "| 9 | bplus_1 | -0.99 | 0.05 | | | -2 | 2 | |\n", + "| 10| jpsi_p | -1.587 | 0.027 | | |-6.28319 | 6.28319 | |\n", + "| 11| omega_p | -5.39 | 0.25 | | |-6.28319 | 6.28319 | |\n", + "| 12| phi_p | -0.24 | 0.29 | | |-6.28319 | 6.28319 | |\n", + "| 13| p4040_p | 3.46 | 0.18 | | |-6.28319 | 6.28319 | |\n", + "| 14| p4415_p | 3.80 | 0.17 | | |-6.28319 | 6.28319 | |\n", + "| 15| Ctt | 1.38 | 0.18 | | | -1.7 | 1.7 | |\n", + "| 16| bplus_2 | -0.31 | 0.10 | | | -2 | 2 | |\n", + "| 17| DDstar_s | 0.30 | 0.40 | | | -0.3 | 0.3 | |\n", + "| 18| p4415_s | 1.49 | 0.20 | | |0.126447 | 2.35355 | |\n", + "| 19| p4160_s | 2.38 | 0.17 | | | 0.71676 | 3.68324 | |\n", + "| 20| psi2s_p | 2.09 | 0.04 | | |-6.28319 | 6.28319 | |\n", + "| 21| rho_s | 1.3 | 0.4 | | |0.0253049| 2.0747 | |\n", + "| 22| p3770_s | 2.72 | 0.25 | | |0.918861 | 4.08114 | |\n", + "| 23| phi_s | 16.9 | 1.2 | | | 14.8182 | 23.5818 | |\n", + "----------------------------------------------------------------------------------------------\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "| | omega_s bplus_0 p3770_p Dbar_s p4160_p DDstar_p Dbar_p p4040_s rho_p bplus_1 jpsi_p omega_p phi_p p4040_p p4415_p Ctt bplus_2 DDstar_s p4415_s p4160_s psi2s_p rho_s p3770_s phi_s |\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "| omega_s | 1.000 -0.022 -0.004 0.000 0.001 -0.011 0.005 0.002 -0.034 0.008 0.005 -0.486 0.015 0.002 0.005 0.002 0.018 -0.004 -0.002 0.003 0.002 -0.121 0.001 -0.001 |\n", + "| bplus_0 | -0.022 1.000 0.029 0.009 0.000 0.058 0.020 0.012 0.036 -0.684 -0.081 0.122 0.002 0.002 -0.006 0.005 -0.064 0.020 0.025 0.029 -0.032 0.225 0.013 -0.016 |\n", + "| p3770_p | -0.004 0.029 1.000 -0.004 0.408 0.475 -0.587 -0.199 -0.024 0.145 0.175 0.013 -0.004 0.335 0.168 0.069 -0.208 0.209 -0.121 -0.085 0.049 0.030 -0.111 0.005 |\n", + "| Dbar_s | 0.000 0.009 -0.004 1.000 -0.042 -0.143 0.206 0.044 0.016 0.000 0.061 0.003 0.001 -0.062 -0.018 -0.143 0.042 -0.031 0.038 0.011 -0.029 0.004 -0.057 -0.005 |\n", + "| p4160_p | 0.001 0.000 0.408 -0.042 1.000 0.227 -0.551 -0.536 -0.018 -0.063 0.243 -0.005 -0.001 0.354 0.414 -0.006 0.003 0.193 -0.122 -0.183 0.130 -0.010 0.040 0.002 |\n", + "| DDstar_p | -0.011 0.058 0.475 -0.143 0.227 1.000 -0.555 -0.279 0.003 0.561 0.015 0.052 -0.003 0.115 -0.141 0.107 -0.513 0.268 -0.048 -0.194 -0.016 0.111 0.058 0.009 |\n", + "| Dbar_p | 0.005 0.020 -0.587 0.206 -0.551 -0.555 1.000 0.255 0.047 -0.093 -0.397 -0.003 -0.003 -0.487 -0.273 -0.381 0.192 -0.063 0.153 0.074 -0.378 -0.011 -0.245 -0.020 |\n", + "| p4040_s | 0.002 0.012 -0.199 0.044 -0.536 -0.279 0.255 1.000 0.033 -0.059 -0.048 -0.001 0.003 -0.254 -0.217 0.209 0.108 -0.081 0.184 0.080 0.019 -0.008 0.103 -0.009 |\n", + "| rho_p | -0.034 0.036 -0.024 0.016 -0.018 0.003 0.047 0.033 1.000 0.006 -0.121 0.089 0.143 -0.002 0.021 0.083 0.218 0.004 0.007 0.051 -0.028 0.257 0.012 0.146 |\n", + "| bplus_1 | 0.008 -0.684 0.145 0.000 -0.063 0.561 -0.093 -0.059 0.006 1.000 -0.129 -0.043 0.003 -0.076 -0.227 -0.091 -0.437 0.197 0.127 -0.057 -0.096 -0.066 -0.011 0.018 |\n", + "| jpsi_p | 0.005 -0.081 0.175 0.061 0.243 0.015 -0.397 -0.048 -0.121 -0.129 1.000 -0.024 -0.033 0.204 0.208 0.387 0.132 0.151 -0.103 -0.006 0.209 -0.032 0.134 -0.001 |\n", + "| omega_p | -0.486 0.122 0.013 0.003 -0.005 0.052 -0.003 -0.001 0.089 -0.043 -0.024 1.000 0.033 -0.006 -0.015 0.008 -0.021 0.020 0.011 0.002 -0.012 0.639 0.001 0.018 |\n", + "| phi_p | 0.015 0.002 -0.004 0.001 -0.001 -0.003 -0.003 0.003 0.143 0.003 -0.033 0.033 1.000 0.002 0.004 0.003 0.020 -0.000 0.000 0.005 -0.002 0.026 0.003 0.753 |\n", + "| p4040_p | 0.002 0.002 0.335 -0.062 0.354 0.115 -0.487 -0.254 -0.002 -0.076 0.204 -0.006 0.002 1.000 0.275 0.081 0.020 0.096 -0.015 0.289 0.069 -0.014 0.028 -0.001 |\n", + "| p4415_p | 0.005 -0.006 0.168 -0.018 0.414 -0.141 -0.273 -0.217 0.021 -0.227 0.208 -0.015 0.004 0.275 1.000 0.104 0.245 0.057 -0.216 -0.138 0.120 -0.038 0.048 -0.009 |\n", + "| Ctt | 0.002 0.005 0.069 -0.143 -0.006 0.107 -0.381 0.209 0.083 -0.091 0.387 0.008 0.003 0.081 0.104 1.000 0.471 -0.043 0.108 0.224 0.369 -0.001 0.123 -0.023 |\n", + "| bplus_2 | 0.018 -0.064 -0.208 0.042 0.003 -0.513 0.192 0.108 0.218 -0.437 0.132 -0.021 0.020 0.020 0.245 0.471 1.000 -0.154 -0.149 0.117 0.045 -0.075 -0.013 -0.055 |\n", + "| DDstar_s | -0.004 0.020 0.209 -0.031 0.193 0.268 -0.063 -0.081 0.004 0.197 0.151 0.020 -0.000 0.096 0.057 -0.043 -0.154 1.000 -0.014 -0.046 0.030 0.042 0.028 0.003 |\n", + "| p4415_s | -0.002 0.025 -0.121 0.038 -0.122 -0.048 0.153 0.184 0.007 0.127 -0.103 0.011 0.000 -0.015 -0.216 0.108 -0.149 -0.014 1.000 0.327 -0.029 0.023 0.021 -0.001 |\n", + "| p4160_s | 0.003 0.029 -0.085 0.011 -0.183 -0.194 0.074 0.080 0.051 -0.057 -0.006 0.002 0.005 0.289 -0.138 0.224 0.117 -0.046 0.327 1.000 -0.011 -0.004 0.045 -0.014 |\n", + "| psi2s_p | 0.002 -0.032 0.049 -0.029 0.130 -0.016 -0.378 0.019 -0.028 -0.096 0.209 -0.012 -0.002 0.069 0.120 0.369 0.045 0.030 -0.029 -0.011 1.000 -0.023 -0.038 0.004 |\n", + "| rho_s | -0.121 0.225 0.030 0.004 -0.010 0.111 -0.011 -0.008 0.257 -0.066 -0.032 0.639 0.026 -0.014 -0.038 -0.001 -0.075 0.042 0.023 -0.004 -0.023 1.000 -0.001 0.024 |\n", + "| p3770_s | 0.001 0.013 -0.111 -0.057 0.040 0.058 -0.245 0.103 0.012 -0.011 0.134 0.001 0.003 0.028 0.048 0.123 -0.013 0.028 0.021 0.045 -0.038 -0.001 1.000 -0.003 |\n", + "| phi_s | -0.001 -0.016 0.005 -0.005 0.002 0.009 -0.020 -0.009 0.146 0.018 -0.001 0.018 0.753 -0.001 -0.009 -0.023 -0.055 0.003 -0.001 -0.014 0.004 0.024 -0.003 1.000 |\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 3.9326813053812453}), (, {'error': 0.02168326819870492}), (, {'error': 0.11567469113834594}), (, {'error': 0.3764346664705315}), (, {'error': 0.11228630591590139}), (, {'error': 0.6687629557879688}), (, {'error': 0.5907017389159401}), (, {'error': 0.18083047386916495}), (, {'error': 0.31202550578599464}), (, {'error': 0.0536773058697807}), (, {'error': 0.027379017276564}), (, {'error': 0.25064367393786924}), (, {'error': 0.29421263114623963}), (, {'error': 0.1807335483053949}), (, {'error': 0.16618570579769454}), (, {'error': 0.1847024331965177}), (, {'error': 0.09957544026395992}), (, {'error': 0.39729640383541787}), (, {'error': 0.19704992509546404}), (, {'error': 0.1728843095092969}), (, {'error': 0.03621310620909046}), (, {'error': 0.3539774140551096}), (, {'error': 0.2530895447662429}), (, {'error': 1.2313732424326842})])\n", + "Step: 6/10\n", + "Current Ctt: 1.1547005383792517\n", + "Ctt floating: True\n", + "Toy 0/1 - Fit 1/5\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "------------------------------------------------------------------\n", + "| FCN = 2.967E+05 | Ncalls=1345 (1345 total) |\n", + "| EDM = 5.67E-05 (Goal: 5E-06) | up = 0.5 |\n", + "------------------------------------------------------------------\n", + "| Valid Min. | Valid Param. | Above EDM | Reached call limit |\n", + "------------------------------------------------------------------\n", + "| True | True | False | False |\n", + "------------------------------------------------------------------\n", + "| Hesse failed | Has cov. | Accurate | Pos. def. | Forced |\n", + "------------------------------------------------------------------\n", + "| False | True | False | False | True |\n", + "------------------------------------------------------------------\n", + "Function minimum: 296721.3183498728\n", + "----------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "----------------------------------------------------------------------------------------------\n", + "| 0 | omega_s | 7.3 | 1.4 | | | 4.19232 | 9.40768 | |\n", + "| 1 | bplus_0 | -0.52 | 0.04 | | | -2 | 2 | |\n", + "| 2 | p3770_p | 2.93 | 0.12 | | |-6.28319 | 6.28319 | |\n", + "| 3 | Dbar_s | 0.300 | 0.022 | | | -0.3 | 0.3 | |\n", + "| 4 | p4160_p | -2.41 | 0.11 | | |-6.28319 | 6.28319 | |\n", + "| 5 | DDstar_p | -4.1 | 0.3 | | |-6.28319 | 6.28319 | |\n", + "| 6 | Dbar_p | -3.6 | 0.3 | | |-6.28319 | 6.28319 | |\n", + "| 7 | p4040_s | 0.98 | 0.18 | | |0.00501244| 2.01499 | |\n", + "| 8 | rho_p | 5.85 | 0.27 | | |-6.28319 | 6.28319 | |\n", + "| 9 | bplus_1 | 1.01 | 0.08 | | | -2 | 2 | |\n", + "| 10| jpsi_p | -4.615 | 0.025 | | |-6.28319 | 6.28319 | |\n", + "| 11| omega_p | 0.4 | 0.3 | | |-6.28319 | 6.28319 | |\n", + "| 12| phi_p | -6.28 | 0.17 | | |-6.28319 | 6.28319 | |\n", + "| 13| p4040_p | -3.15 | 0.17 | | |-6.28319 | 6.28319 | |\n", + "| 14| p4415_p | -6.28 | 0.16 | | |-6.28319 | 6.28319 | |\n", + "| 15| Ctt | 1.24 | 0.16 | | | -1.7 | 1.7 | |\n", + "| 16| bplus_2 | 0.27 | 0.09 | | | -2 | 2 | |\n", + "| 17| DDstar_s | 0.300 | 0.026 | | | -0.3 | 0.3 | |\n", + "| 18| p4415_s | 0.126 | 0.017 | | |0.126447 | 2.35355 | |\n", + "| 19| p4160_s | 1.90 | 0.17 | | | 0.71676 | 3.68324 | |\n", + "| 20| psi2s_p | 2.00 | 0.03 | | |-6.28319 | 6.28319 | |\n", + "| 21| rho_s | 1.4 | 0.3 | | |0.0253049| 2.0747 | |\n", + "| 22| p3770_s | 2.42 | 0.25 | | |0.918861 | 4.08114 | |\n", + "| 23| phi_s | 18.5 | 0.9 | | | 14.8182 | 23.5818 | |\n", + "----------------------------------------------------------------------------------------------\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "| | omega_s bplus_0 p3770_p Dbar_s p4160_p DDstar_p Dbar_p p4040_s rho_p bplus_1 jpsi_p omega_p phi_p p4040_p p4415_p Ctt bplus_2 DDstar_s p4415_s p4160_s psi2s_p rho_s p3770_s phi_s |\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "| omega_s | 1.000 0.158 -0.014 -0.002 -0.011 0.037 0.014 0.001 0.083 -0.146 -0.047 0.731 0.008 -0.012 0.001 0.015 -0.063 -0.002 0.000 -0.002 -0.011 -0.318 -0.005 0.079 |\n", + "| bplus_0 | 0.158 1.000 -0.040 -0.005 -0.042 0.131 0.084 -0.005 0.006 -0.951 -0.187 0.050 -0.001 -0.071 0.003 0.051 -0.221 -0.005 -0.000 -0.036 -0.038 -0.326 -0.035 0.120 |\n", + "| p3770_p | -0.014 -0.040 1.000 -0.022 0.263 -0.241 -0.407 -0.113 0.012 -0.069 0.257 -0.004 -0.001 0.161 -0.001 -0.241 0.132 -0.017 0.005 -0.043 -0.105 0.024 -0.261 -0.020 |\n", + "| Dbar_s | -0.002 -0.005 -0.022 1.000 -0.028 -0.002 0.013 0.014 0.004 -0.004 -0.029 -0.001 0.000 -0.027 0.001 -0.022 -0.005 -0.002 -0.001 -0.001 -0.024 0.003 -0.014 -0.002 |\n", + "| p4160_p | -0.011 -0.042 0.263 -0.028 1.000 -0.275 -0.182 -0.460 0.017 -0.012 0.184 -0.003 -0.000 -0.072 -0.008 -0.316 0.042 -0.042 0.008 -0.207 -0.097 0.020 -0.114 -0.018 |\n", + "| DDstar_p | 0.037 0.131 -0.241 -0.002 -0.275 1.000 -0.196 0.113 -0.007 0.087 -0.188 0.011 0.000 -0.116 0.000 0.053 -0.259 0.029 -0.004 0.049 -0.072 -0.081 -0.052 0.024 |\n", + "| Dbar_p | 0.014 0.084 -0.407 0.013 -0.182 -0.196 1.000 0.068 -0.028 0.002 -0.399 0.000 -0.001 -0.060 0.002 0.179 -0.102 -0.006 -0.001 0.079 -0.136 -0.035 0.038 0.027 |\n", + "| p4040_s | 0.001 -0.005 -0.113 0.014 -0.460 0.113 0.068 1.000 0.016 0.002 -0.075 0.001 0.001 -0.196 0.004 0.344 -0.044 0.005 -0.008 -0.134 0.135 0.000 0.171 -0.007 |\n", + "| rho_p | 0.083 0.006 0.012 0.004 0.017 -0.007 -0.028 0.016 1.000 -0.017 0.021 -0.017 0.041 0.051 -0.001 0.110 -0.250 0.002 0.002 0.063 0.025 0.162 0.030 0.051 |\n", + "| bplus_1 | -0.146 -0.951 -0.069 -0.004 -0.012 0.087 0.002 0.002 -0.017 1.000 0.041 -0.047 -0.000 0.017 -0.000 0.034 0.033 -0.004 0.001 0.016 -0.011 0.291 0.001 -0.112 |\n", + "| jpsi_p | -0.047 -0.187 0.257 -0.029 0.184 -0.188 -0.399 -0.075 0.021 0.041 1.000 -0.018 -0.004 0.044 0.001 -0.392 0.294 -0.049 0.002 -0.103 0.070 0.076 0.053 -0.053 |\n", + "| omega_p | 0.731 0.050 -0.004 -0.001 -0.003 0.011 0.000 0.001 -0.017 -0.047 -0.018 1.000 0.013 -0.001 0.000 0.011 -0.036 -0.001 0.000 0.003 -0.003 0.019 0.000 0.045 |\n", + "| phi_p | 0.008 -0.001 -0.001 0.000 -0.000 0.000 -0.001 0.001 0.041 -0.000 -0.004 0.013 1.000 0.001 -0.000 0.003 -0.008 0.000 0.000 0.002 -0.000 0.003 0.001 0.268 |\n", + "| p4040_p | -0.012 -0.071 0.161 -0.027 -0.072 -0.116 -0.060 -0.196 0.051 0.017 0.044 -0.001 0.001 1.000 -0.004 -0.063 -0.017 -0.035 -0.002 0.340 -0.178 0.023 -0.090 -0.037 |\n", + "| p4415_p | 0.001 0.003 -0.001 0.001 -0.008 0.000 0.002 0.004 -0.001 -0.000 0.001 0.000 -0.000 -0.004 1.000 -0.004 -0.001 0.001 -0.011 -0.003 0.001 -0.001 0.000 0.001 |\n", + "| Ctt | 0.015 0.051 -0.241 -0.022 -0.316 0.053 0.179 0.344 0.110 0.034 -0.392 0.011 0.003 -0.063 -0.004 1.000 -0.668 0.007 -0.004 0.260 0.282 -0.039 0.031 -0.047 |\n", + "| bplus_2 | -0.063 -0.221 0.132 -0.005 0.042 -0.259 -0.102 -0.044 -0.250 0.033 0.294 -0.036 -0.008 -0.017 -0.001 -0.668 1.000 -0.000 -0.006 -0.125 -0.027 0.120 0.007 0.081 |\n", + "| DDstar_s | -0.002 -0.005 -0.017 -0.002 -0.042 0.029 -0.006 0.005 0.002 -0.004 -0.049 -0.001 0.000 -0.035 0.001 0.007 -0.000 1.000 -0.001 -0.008 -0.024 0.003 -0.012 -0.002 |\n", + "| p4415_s | 0.000 -0.000 0.005 -0.001 0.008 -0.004 -0.001 -0.008 0.002 0.001 0.002 0.000 0.000 -0.002 -0.011 -0.004 -0.006 -0.001 1.000 -0.011 -0.003 -0.000 -0.002 -0.001 |\n", + "| p4160_s | -0.002 -0.036 -0.043 -0.001 -0.207 0.049 0.079 -0.134 0.063 0.016 -0.103 0.003 0.002 0.340 -0.003 0.260 -0.125 -0.008 -0.011 1.000 -0.046 0.005 0.025 -0.035 |\n", + "| psi2s_p | -0.011 -0.038 -0.105 -0.024 -0.097 -0.072 -0.136 0.135 0.025 -0.011 0.070 -0.003 -0.000 -0.178 0.001 0.282 -0.027 -0.024 -0.003 -0.046 1.000 0.019 -0.193 -0.021 |\n", + "| rho_s | -0.318 -0.326 0.024 0.003 0.020 -0.081 -0.035 0.000 0.162 0.291 0.076 0.019 0.003 0.023 -0.001 -0.039 0.120 0.003 -0.000 0.005 0.019 1.000 0.010 -0.055 |\n", + "| p3770_s | -0.005 -0.035 -0.261 -0.014 -0.114 -0.052 0.038 0.171 0.030 0.001 0.053 0.000 0.001 -0.090 0.000 0.031 0.007 -0.012 -0.002 0.025 -0.193 0.010 1.000 -0.021 |\n", + "| phi_s | 0.079 0.120 -0.020 -0.002 -0.018 0.024 0.027 -0.007 0.051 -0.112 -0.053 0.045 0.268 -0.037 0.001 -0.047 0.081 -0.002 -0.001 -0.035 -0.021 -0.055 -0.021 1.000 |\n", + "--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 1.4493227649327043}), (, {'error': 0.04023703321701233}), (, {'error': 0.11988744429727483}), (, {'error': 0.021741030731362687}), (, {'error': 0.11330152974230501}), (, {'error': 0.32996699930309137}), (, {'error': 0.3194463581349356}), (, {'error': 0.17917956236492494}), (, {'error': 0.2680040701168531}), (, {'error': 0.07696747357347622}), (, {'error': 0.0249974764536276}), (, {'error': 0.33609626135950776}), (, {'error': 0.1680879865000633}), (, {'error': 0.17126572009411434}), (, {'error': 0.15524039825176805}), (, {'error': 0.15586721762207345}), (, {'error': 0.09372682486883277}), (, {'error': 0.026443324962821146}), (, {'error': 0.017171498127401916}), (, {'error': 0.16786715955856713}), (, {'error': 0.03253517173554954}), (, {'error': 0.32491772133145025}), (, {'error': 0.25177083236717657}), (, {'error': 0.8976335198871119})])\n", + "Step: 6/10\n", + "Current Ctt: 1.1547005383792517\n", + "Ctt floating: True\n", + "Toy 0/1 - Fit 2/5\n" ] }, { "ename": "RuntimeError", - "evalue": "exception was raised in user function\nUser function arguments:\n bplus_0 = -0.627540\n p4160_s = +3.373292\n DDstar_s = +0.299698\n p3770_s = +3.780920\n jpsi_p = -5.063578\n p3770_p = +2.196764\n bplus_1 = +1.062520\n DDstar_p = +4.628551\n omega_s = +5.708895\n phi_s = +16.739100\n Ctt = +3.195816\n p4415_s = +0.911410\n psi2s_p = +2.482877\n Dbar_p = +0.982372\n p4040_p = -2.632544\n p4160_p = +2.741969\n p4415_p = -2.983702\n phi_p = -4.941643\n rho_p = +1.587532\n rho_s = +0.149543\n bplus_2 = -0.691416\n Dbar_s = +0.150965\n omega_p = -0.692164\n p4040_s = +1.222875\nOriginal python exception in user function:\nKeyboardInterrupt: \n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\zfit\\minimizers\\minimizer_minuit.py\", line 101, in func\n loss_evaluated = self.sess.run(loss_val)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 929, in run\n run_metadata_ptr)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1152, in _run\n feed_dict_tensor, options, run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1328, in _do_run\n run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1334, in _do_call\n return fn(*args)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1319, in _run_fn\n options, feed_dict, fetch_list, target_list, run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1407, in _call_tf_sessionrun\n run_metadata)\n", + "evalue": "exception was raised in user function\nUser function arguments:\n omega_s = +6.558876\n bplus_0 = -0.537708\n p3770_p = +2.922976\n Dbar_s = +0.299117\n p4160_p = -2.459129\n DDstar_p = +4.501845\n Dbar_p = +2.685834\n p4040_s = +1.254795\n rho_p = -0.384056\n bplus_1 = +0.937284\n jpsi_p = -4.628444\n omega_p = +0.275809\n phi_p = +6.126912\n p4040_p = -3.035592\n p4415_p = +3.714495\n Ctt = +1.579321\n bplus_2 = +0.437803\n DDstar_s = -0.246898\n p4415_s = +1.615964\n p4160_s = +2.362930\n psi2s_p = +2.020719\n rho_s = +1.998466\n p3770_s = +2.569549\n phi_s = +17.305781\nOriginal python exception in user function:\nKeyboardInterrupt: \n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\zfit\\minimizers\\minimizer_minuit.py\", line 101, in func\n loss_evaluated = self.sess.run(loss_val)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 929, in run\n run_metadata_ptr)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1152, in _run\n feed_dict_tensor, options, run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1328, in _do_run\n run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1334, in _do_call\n return fn(*args)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1319, in _run_fn\n options, feed_dict, fetch_list, target_list, run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1407, in _call_tf_sessionrun\n run_metadata)\n", "output_type": "error", "traceback": [ "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[1;31mRuntimeError\u001b[0m Traceback (most recent call last)", - "\u001b[1;32m\u001b[0m in \u001b[0;36m\u001b[1;34m\u001b[0m\n\u001b[0;32m 101\u001b[0m \u001b[0mminimizer\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mzfit\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mminimize\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mMinuitMinimizer\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mverbosity\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;36m5\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 102\u001b[0m \u001b[1;31m# minimizer._use_tfgrad = False\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 103\u001b[1;33m \u001b[0mresult\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mminimizer\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mminimize\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mnll\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 104\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 105\u001b[0m \u001b[1;31m# print(\"Function minimum:\", result.fmin)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m\u001b[0m in \u001b[0;36m\u001b[1;34m\u001b[0m\n\u001b[0;32m 104\u001b[0m \u001b[0mminimizer\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mzfit\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mminimize\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mMinuitMinimizer\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mverbosity\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;36m5\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 105\u001b[0m \u001b[1;31m# minimizer._use_tfgrad = False\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 106\u001b[1;33m \u001b[0mresult\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mminimizer\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mminimize\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mnll\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 107\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 108\u001b[0m \u001b[0mprint\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"Function minimum:\"\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mresult\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfmin\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", "\u001b[1;32m~\\.conda\\envs\\rmd\\lib\\site-packages\\zfit\\minimizers\\baseminimizer.py\u001b[0m in \u001b[0;36mminimize\u001b[1;34m(self, loss, params)\u001b[0m\n\u001b[0;32m 205\u001b[0m \u001b[0mtuple\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mstack\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0menter_context\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mparam\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mset_sess\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msess\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mparam\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mparams\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 206\u001b[0m \u001b[1;32mtry\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 207\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_hook_minimize\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mloss\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mloss\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mparams\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mparams\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 208\u001b[0m \u001b[1;32mexcept\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mFailMinimizeNaN\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mRuntimeError\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0merror\u001b[0m\u001b[1;33m:\u001b[0m \u001b[1;31m# iminuit raises RuntimeError if user raises Error\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 209\u001b[0m \u001b[0mfail_result\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mstrategy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfit_result\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", "\u001b[1;32m~\\.conda\\envs\\rmd\\lib\\site-packages\\zfit\\minimizers\\baseminimizer.py\u001b[0m in \u001b[0;36m_hook_minimize\u001b[1;34m(self, loss, params)\u001b[0m\n\u001b[0;32m 214\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 215\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0m_hook_minimize\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mloss\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mparams\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 216\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_call_minimize\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mloss\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mloss\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mparams\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mparams\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 217\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 218\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0m_call_minimize\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mloss\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mparams\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", "\u001b[1;32m~\\.conda\\envs\\rmd\\lib\\site-packages\\zfit\\minimizers\\baseminimizer.py\u001b[0m in \u001b[0;36m_call_minimize\u001b[1;34m(self, loss, params)\u001b[0m\n\u001b[0;32m 218\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0m_call_minimize\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mloss\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mparams\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 219\u001b[0m \u001b[1;32mtry\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 220\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_minimize\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mloss\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mloss\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mparams\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mparams\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 221\u001b[0m \u001b[1;32mexcept\u001b[0m \u001b[0mNotImplementedError\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0merror\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 222\u001b[0m \u001b[1;32mtry\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", "\u001b[1;32m~\\.conda\\envs\\rmd\\lib\\site-packages\\zfit\\minimizers\\minimizer_minuit.py\u001b[0m in \u001b[0;36m_minimize\u001b[1;34m(self, loss, params)\u001b[0m\n\u001b[0;32m 136\u001b[0m minimizer_setter)\n\u001b[0;32m 137\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_minuit_minimizer\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mminimizer\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 138\u001b[1;33m \u001b[0mresult\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mminimizer\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmigrad\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m**\u001b[0m\u001b[0mminimize_options\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 139\u001b[0m \u001b[0mparams_result\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mp_dict\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mp_dict\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mresult\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 140\u001b[0m \u001b[0mresult_vals\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mres\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;34m\"value\"\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mres\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mparams_result\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", "\u001b[1;32miminuit\\_libiminuit.pyx\u001b[0m in \u001b[0;36miminuit._libiminuit.Minuit.migrad\u001b[1;34m()\u001b[0m\n", - "\u001b[1;31mRuntimeError\u001b[0m: exception was raised in user function\nUser function arguments:\n bplus_0 = -0.627540\n p4160_s = +3.373292\n DDstar_s = +0.299698\n p3770_s = +3.780920\n jpsi_p = -5.063578\n p3770_p = +2.196764\n bplus_1 = +1.062520\n DDstar_p = +4.628551\n omega_s = +5.708895\n phi_s = +16.739100\n Ctt = +3.195816\n p4415_s = +0.911410\n psi2s_p = +2.482877\n Dbar_p = +0.982372\n p4040_p = -2.632544\n p4160_p = +2.741969\n p4415_p = -2.983702\n phi_p = -4.941643\n rho_p = +1.587532\n rho_s = +0.149543\n bplus_2 = -0.691416\n Dbar_s = +0.150965\n omega_p = -0.692164\n p4040_s = +1.222875\nOriginal python exception in user function:\nKeyboardInterrupt: \n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\zfit\\minimizers\\minimizer_minuit.py\", line 101, in func\n loss_evaluated = self.sess.run(loss_val)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 929, in run\n run_metadata_ptr)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1152, in _run\n feed_dict_tensor, options, run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1328, in _do_run\n run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1334, in _do_call\n return fn(*args)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1319, in _run_fn\n options, feed_dict, fetch_list, target_list, run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1407, in _call_tf_sessionrun\n run_metadata)\n" + "\u001b[1;31mRuntimeError\u001b[0m: exception was raised in user function\nUser function arguments:\n omega_s = +6.558876\n bplus_0 = -0.537708\n p3770_p = +2.922976\n Dbar_s = +0.299117\n p4160_p = -2.459129\n DDstar_p = +4.501845\n Dbar_p = +2.685834\n p4040_s = +1.254795\n rho_p = -0.384056\n bplus_1 = +0.937284\n jpsi_p = -4.628444\n omega_p = +0.275809\n phi_p = +6.126912\n p4040_p = -3.035592\n p4415_p = +3.714495\n Ctt = +1.579321\n bplus_2 = +0.437803\n DDstar_s = -0.246898\n p4415_s = +1.615964\n p4160_s = +2.362930\n psi2s_p = +2.020719\n rho_s = +1.998466\n p3770_s = +2.569549\n phi_s = +17.305781\nOriginal python exception in user function:\nKeyboardInterrupt: \n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\zfit\\minimizers\\minimizer_minuit.py\", line 101, in func\n loss_evaluated = self.sess.run(loss_val)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 929, in run\n run_metadata_ptr)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1152, in _run\n feed_dict_tensor, options, run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1328, in _do_run\n run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1334, in _do_call\n return fn(*args)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1319, in _run_fn\n options, feed_dict, fetch_list, target_list, run_metadata)\n File \"C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow\\python\\client\\session.py\", line 1407, in _call_tf_sessionrun\n run_metadata)\n" ] } ], @@ -2033,7 +2362,8 @@ "\n", "Nll_list = []\n", "\n", - "sampler = total_f.create_sampler(n=nevents, limits = obs_fit)\n", + "sampler = total_f.create_sampler(n=nevents)\n", + "sampler.set_data_range(obs_fit)\n", "\n", "__ = 0\n", "\n", @@ -2077,19 +2407,20 @@ " \n", " if newset or len(Nll_list[-1])%bo5_set == 0:\n", " sampler.resample(n=nevents)\n", - " data = sampler.unstack_x()\n", + " data = sampler\n", "# total_samp = zfit.run(s)\n", " calls = 0\n", " c = 1\n", " newset = False\n", + " \n", "\n", "\n", "# data = zfit.data.Data.from_numpy(array=total_samp[:int(nevents)], obs=obs_fit)\n", "\n", " ### Fit data\n", "\n", - " for param in total_f_fit.get_dependents():\n", - " param.randomize()\n", + "# for param in total_f_fit.get_dependents():\n", + "# param.randomize()\n", "\n", " nll = zfit.loss.UnbinnedNLL(model=total_f_fit, data=data, constraints = constraints)\n", "\n", @@ -2119,7 +2450,7 @@ "\n", " print('Toy {0}/{1}'.format(len(Nll_list[-1]), nr_of_toys))\n", "\n", - " reset_param_values()\n", + " reset_param_values(variation = 0.0)\n", "\n", " if floaty:\n", " Ctt.set_value(Ctt_step)\n", @@ -2138,8 +2469,7 @@ "\n", " ### Fit data\n", "\n", - " for param in total_f_fit.get_dependents():\n", - " param.randomize()\n", + " reset_param_values()\n", "\n", " nll = zfit.loss.UnbinnedNLL(model=total_f_fit, data=data, constraints = constraints)\n", "\n", @@ -2328,6 +2658,44 @@ }, { "cell_type": "code", + "execution_count": 51, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "omega_s\n", + "bplus_0\n", + "Dbar_s\n", + "p3770_p\n", + "p4160_p\n", + "DDstar_p\n", + "Dbar_p\n", + "p4040_s\n", + "rho_p\n", + "bplus_1\n", + "jpsi_p\n", + "omega_p\n", + "phi_p\n", + "p4040_p\n", + "p4415_p\n", + "Ctt\n", + "bplus_2\n", + "DDstar_s\n", + "p4415_s\n", + "p4160_s\n", + "psi2s_p\n", + "rho_s\n", + "p3770_s\n", + "phi_s\n" + ] + } + ], + "source": [] + }, + { + "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], diff --git a/raremodel-nb.ipynb b/raremodel-nb.ipynb index 1c48d41..6f23704 100644 --- a/raremodel-nb.ipynb +++ b/raremodel-nb.ipynb @@ -1579,66 +1579,72 @@ }, { "cell_type": "code", - "execution_count": 38, + "execution_count": 54, "metadata": {}, "outputs": [], "source": [ - "def reset_param_values(): \n", - " jpsi_m.set_value(jpsi_mass)\n", - " jpsi_s.set_value(jpsi_scale)\n", - " jpsi_p.set_value(jpsi_phase)\n", - " jpsi_w.set_value(jpsi_width)\n", - " psi2s_m.set_value(psi2s_mass)\n", - " psi2s_s.set_value(psi2s_scale)\n", - " psi2s_p.set_value(psi2s_phase)\n", - " psi2s_w.set_value(psi2s_width)\n", - " p3770_m.set_value(p3770_mass)\n", - " p3770_s.set_value(p3770_scale)\n", - " p3770_p.set_value(p3770_phase)\n", - " p3770_w.set_value(p3770_width)\n", - " p4040_m.set_value(p4040_mass)\n", - " p4040_s.set_value(p4040_scale)\n", - " p4040_p.set_value(p4040_phase)\n", - " p4040_w.set_value(p4040_width)\n", - " p4160_m.set_value(p4160_mass)\n", - " p4160_s.set_value(p4160_scale)\n", - " p4160_p.set_value(p4160_phase)\n", - " p4160_w.set_value(p4160_width)\n", - " p4415_m.set_value(p4415_mass)\n", - " p4415_s.set_value(p4415_scale)\n", - " p4415_p.set_value(p4415_phase)\n", - " p4415_w.set_value(p4415_width)\n", - " rho_m.set_value(rho_mass)\n", - " rho_s.set_value(rho_scale)\n", - " rho_p.set_value(rho_phase)\n", - " rho_w.set_value(rho_width)\n", - " omega_m.set_value(omega_mass)\n", - " omega_s.set_value(omega_scale)\n", - " omega_p.set_value(omega_phase)\n", - " omega_w.set_value(omega_width)\n", - " phi_m.set_value(phi_mass)\n", - " phi_s.set_value(phi_scale)\n", - " phi_p.set_value(phi_phase)\n", - " phi_w.set_value(phi_width)\n", - " Dstar_m.set_value(Dstar_mass)\n", - " DDstar_s.set_value(0.0)\n", - " DDstar_p.set_value(0.0)\n", - " D_m.set_value(D_mass)\n", - " Dbar_m.set_value(Dbar_mass)\n", - " Dbar_s.set_value(0.0)\n", - " Dbar_p.set_value(0.0)\n", - " tau_m.set_value(pdg['tau_M'])\n", - " Ctt.set_value(0.0)\n", - " b0_0.set_value(0.292)\n", - " b0_1.set_value(0.281)\n", - " b0_2.set_value(0.150)\n", - " bplus_0.set_value(0.466)\n", - " bplus_1.set_value(-0.885)\n", - " bplus_2.set_value(-0.213)\n", - " bT_0.set_value(0.460)\n", - " bT_1.set_value(-1.089)\n", - " bT_2.set_value(-1.114)\n", + "param_values_dic = {\n", + " 'jpsi_m': jpsi_mass,\n", + " 'jpsi_s': jpsi_scale,\n", + " 'jpsi_p': jpsi_phase,\n", + " 'jpsi_w': jpsi_width,\n", + " 'psi2s_m': psi2s_mass,\n", + " 'psi2s_s': psi2s_scale,\n", + " 'psi2s_p': psi2s_phase,\n", + " 'psi2s_w': psi2s_width,\n", + " 'p3770_m': p3770_mass,\n", + " 'p3770_s': p3770_scale,\n", + " 'p3770_p': p3770_phase,\n", + " 'p3770_w': p3770_width,\n", + " 'p4040_m': p4040_mass,\n", + " 'p4040_s': p4040_scale,\n", + " 'p4040_p': p4040_phase,\n", + " 'p4040_w': p4040_width,\n", + " 'p4160_m': p4160_mass,\n", + " 'p4160_s': p4160_scale,\n", + " 'p4160_p': p4160_phase,\n", + " 'p4160_w': p4160_width,\n", + " 'p4415_m': p4415_mass,\n", + " 'p4415_s': p4415_scale,\n", + " 'p4415_p': p4415_phase,\n", + " 'p4415_w': p4415_width,\n", + " 'rho_m': rho_mass,\n", + " 'rho_s': rho_scale,\n", + " 'rho_p': rho_phase,\n", + " 'rho_w': rho_width,\n", + " 'omega_m': omega_mass,\n", + " 'omega_s': omega_scale,\n", + " 'omega_p': omega_phase,\n", + " 'omega_w': omega_width,\n", + " 'phi_m': phi_mass,\n", + " 'phi_s': phi_scale,\n", + " 'phi_p': phi_phase,\n", + " 'phi_w': phi_width,\n", + " 'Dstar_m': Dstar_mass,\n", + " 'DDstar_s': 0.0,\n", + " 'DDstar_p': 0.0,\n", + " 'D_m': D_mass,\n", + " 'Dbar_m': Dbar_mass,\n", + " 'Dbar_s': 0.0,\n", + " 'Dbar_p': 0.0,\n", + " 'tau_m': pdg['tau_M'],\n", + " 'Ctt': 0.0,\n", + " 'b0_0': 0.292,\n", + " 'b0_1': 0.281,\n", + " 'b0_2': 0.150,\n", + " 'bplus_0': 0.466,\n", + " 'bplus_1': -0.885,\n", + " 'bplus_2': -0.213,\n", + " 'bT_0': 0.460,\n", + " 'bT_1': -1.089,\n", + " 'bT_2': -1.114}\n", "\n", + "\n", + "def reset_param_values(variation = 0.05):\n", + " for param in total_f_fit.get_dependents():\n", + " if param.floating:\n", + " param.set_value(param_values_dic[param.name] + random.uniform(-1, 1)*param_values_dic[param.name]* variation)\n", + "# print(param.name)\n", "# for param in totalf.get_dependents():\n", "# param.set_value()" ] @@ -2413,8 +2419,8 @@ "\n", " ### Fit data\n", "\n", - " for param in total_f_fit.get_dependents():\n", - " param.randomize()\n", + "# for param in total_f_fit.get_dependents():\n", + "# param.randomize()\n", "\n", " nll = zfit.loss.UnbinnedNLL(model=total_f_fit, data=data, constraints = constraints)\n", "\n", @@ -2444,7 +2450,7 @@ "\n", " print('Toy {0}/{1}'.format(len(Nll_list[-1]), nr_of_toys))\n", "\n", - " reset_param_values()\n", + " reset_param_values(variation = 0.0)\n", "\n", " if floaty:\n", " Ctt.set_value(Ctt_step)\n", @@ -2463,8 +2469,7 @@ "\n", " ### Fit data\n", "\n", - " for param in total_f_fit.get_dependents():\n", - " param.randomize()\n", + " reset_param_values()\n", "\n", " nll = zfit.loss.UnbinnedNLL(model=total_f_fit, data=data, constraints = constraints)\n", "\n", @@ -2653,6 +2658,44 @@ }, { "cell_type": "code", + "execution_count": 51, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "omega_s\n", + "bplus_0\n", + "Dbar_s\n", + "p3770_p\n", + "p4160_p\n", + "DDstar_p\n", + "Dbar_p\n", + "p4040_s\n", + "rho_p\n", + "bplus_1\n", + "jpsi_p\n", + "omega_p\n", + "phi_p\n", + "p4040_p\n", + "p4415_p\n", + "Ctt\n", + "bplus_2\n", + "DDstar_s\n", + "p4415_s\n", + "p4160_s\n", + "psi2s_p\n", + "rho_s\n", + "p3770_s\n", + "phi_s\n" + ] + } + ], + "source": [] + }, + { + "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [],