diff --git a/.ipynb_checkpoints/raremodel-nb-checkpoint.ipynb b/.ipynb_checkpoints/raremodel-nb-checkpoint.ipynb index 2f54c56..18a9c1b 100644 --- a/.ipynb_checkpoints/raremodel-nb-checkpoint.ipynb +++ b/.ipynb_checkpoints/raremodel-nb-checkpoint.ipynb @@ -1861,19 +1861,7 @@ "cell_type": "code", "execution_count": 40, "metadata": {}, - "outputs": [ - { - "ename": "NameError", - "evalue": "name 'Ctt_list' is not defined", - "output_type": "error", - "traceback": [ - "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", - "\u001b[1;31mNameError\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 1\u001b[0m \u001b[1;32mwith\u001b[0m \u001b[0mopen\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"data/results/Ctt_list.pkl\"\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m\"wb\"\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 2\u001b[1;33m \u001b[0mpkl\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdump\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mCtt_list\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mpkl\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mHIGHEST_PROTOCOL\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 3\u001b[0m \u001b[1;32mwith\u001b[0m \u001b[0mopen\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"data/results/Ctt_error_list.pkl\"\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m\"wb\"\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 4\u001b[0m \u001b[0mpkl\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdump\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mCtt_error_list\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mpkl\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mHIGHEST_PROTOCOL\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", - "\u001b[1;31mNameError\u001b[0m: name 'Ctt_list' is not defined" - ] - } - ], + "outputs": [], "source": [ "# with open(\"data/results/Ctt_list.pkl\", \"wb\") as f:\n", "# pkl.dump(Ctt_list, f, pkl.HIGHEST_PROTOCOL)\n", @@ -1883,7 +1871,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 41, "metadata": {}, "outputs": [], "source": [ @@ -1895,7 +1883,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 42, "metadata": {}, "outputs": [], "source": [ @@ -1910,7 +1898,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 43, "metadata": {}, "outputs": [], "source": [ @@ -1919,7 +1907,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 44, "metadata": {}, "outputs": [], "source": [ @@ -1928,7 +1916,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 45, "metadata": {}, "outputs": [], "source": [ @@ -1937,7 +1925,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 46, "metadata": {}, "outputs": [], "source": [ @@ -1954,32 +1942,483 @@ { "cell_type": "code", "execution_count": 47, - "metadata": {}, + "metadata": { + "scrolled": false + }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ + "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: 0/11\n", "Current Ctt: 0.0\n", "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.977E+05 | Ncalls=1362 (1362 total) |\n", + "| EDM = 0.00848 (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", + "Step: 0/11\n", + "Current Ctt: 0.0\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.976E+05 | Ncalls=1372 (1372 total) |\n", + "| EDM = 0.00145 (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", + "Step: 0/11\n", + "Current Ctt: 0.0\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.977E+05 | Ncalls=1309 (1309 total) |\n", + "| EDM = 2.31E-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", + "Step: 0/11\n", + "Current Ctt: 0.0\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1172 (1172 total) |\n", + "| EDM = 0.000838 (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", + "Time Taken: 463\n", + "Predicted time left: 7.0 min\n", + "Step: 1/11\n", + "Current Ctt: 0.1543033499620919\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 3.026E+05 | Ncalls=505 (516 total) |\n", + "| EDM = 6.64E+08 (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", + "Step: 1/11\n", + "Current Ctt: 0.1543033499620919\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=1459 (1459 total) |\n", + "| EDM = 0.00013 (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", + "Step: 1/11\n", + "Current Ctt: 0.1543033499620919\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1300 (1300 total) |\n", + "| EDM = 1.54E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 813\n", + "Predicted time left: 11.0 min, 5.0 s\n", + "Step: 2/11\n", + "Current Ctt: 0.2182178902359924\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1331 (1331 total) |\n", + "| EDM = 6.73E-06 (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", + "Step: 2/11\n", + "Current Ctt: 0.2182178902359924\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1507 (1507 total) |\n", + "| EDM = 0.000231 (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", + "Time Taken: 1111\n", + "Predicted time left: 13.0 min, 28.0 s\n", + "Step: 3/11\n", + "Current Ctt: 0.2672612419124244\n", + "Ctt floating: True\n", "Toy 0/1\n" ] }, { - "ename": "RuntimeError", - "evalue": "exception was raised in user function\nUser function arguments:\n p4415_s = +0.738610\n psi2s_p = +1.596704\n p4415_p = +5.534902\n phi_s = +14.820266\n DDstar_s = +0.055634\n p3770_s = +3.056326\n p4160_s = +3.013799\n omega_s = +8.221504\n p3770_p = +3.227282\n phi_p = -1.956867\n Ctt = -0.074235\n DDstar_p = -3.658252\n rho_p = +3.220475\n Dbar_p = +1.482929\n p4040_p = -1.211697\n rho_s = +0.426552\n jpsi_p = -4.776017\n omega_p = +2.899561\n p4160_p = -1.589904\n p4040_s = +1.161008\n Dbar_s = +0.130791\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 76\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 77\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---> 78\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 79\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 80\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~\\.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 p4415_s = +0.738610\n psi2s_p = +1.596704\n p4415_p = +5.534902\n phi_s = +14.820266\n DDstar_s = +0.055634\n p3770_s = +3.056326\n p4160_s = +3.013799\n omega_s = +8.221504\n p3770_p = +3.227282\n phi_p = -1.956867\n Ctt = -0.074235\n DDstar_p = -3.658252\n rho_p = +3.220475\n Dbar_p = +1.482929\n p4040_p = -1.211697\n rho_s = +0.426552\n jpsi_p = -4.776017\n omega_p = +2.899561\n p4160_p = -1.589904\n p4040_s = +1.161008\n Dbar_s = +0.130791\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" + "name": "stdout", + "output_type": "stream", + "text": [ + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=1331 (1331 total) |\n", + "| EDM = 0.00103 (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", + "Step: 3/11\n", + "Current Ctt: 0.2672612419124244\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.98E+05 | Ncalls=1407 (1407 total) |\n", + "| EDM = 5.43E-06 (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", + "Step: 3/11\n", + "Current Ctt: 0.2672612419124244\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.989E+05 | Ncalls=566 (577 total) |\n", + "| EDM = 1.99E+08 (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", + "Step: 3/11\n", + "Current Ctt: 0.2672612419124244\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=366 (366 total) |\n", + "| EDM = 6.89E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 1656\n", + "Predicted time left: 17.0 min, 33.0 s\n", + "Step: 4/11\n", + "Current Ctt: 0.3086066999241838\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1364 (1364 total) |\n", + "| EDM = 2.06E-06 (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", + "Step: 4/11\n", + "Current Ctt: 0.3086066999241838\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1497 (1497 total) |\n", + "| EDM = 5.28E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 2044\n", + "Predicted time left: 18.0 min, 34.0 s\n", + "Step: 5/11\n", + "Current Ctt: 0.3450327796711771\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.98E+05 | Ncalls=1249 (1249 total) |\n", + "| EDM = 0.00017 (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", + "Step: 5/11\n", + "Current Ctt: 0.3450327796711771\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.976E+05 | Ncalls=1276 (1276 total) |\n", + "| EDM = 0.000344 (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", + "Time Taken: 2428\n", + "Predicted time left: 18.0 min, 23.0 s\n", + "Step: 6/11\n", + "Current Ctt: 0.37796447300922725\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1446 (1446 total) |\n", + "| EDM = 0.000166 (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", + "Step: 6/11\n", + "Current Ctt: 0.37796447300922725\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=792 (792 total) |\n", + "| EDM = 7.63E-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", + "Time Taken: 2836\n", + "Predicted time left: 17.0 min, 11.0 s\n", + "Step: 7/11\n", + "Current Ctt: 0.408248290463863\n", + "Ctt floating: True\n", + "Toy 0/1\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=1391 (1391 total) |\n", + "| EDM = 0.000185 (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", + "Step: 7/11\n", + "Current Ctt: 0.408248290463863\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=273 (273 total) |\n", + "| EDM = 7.9E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 3232\n", + "Predicted time left: 14.0 min, 41.0 s\n", + "Step: 8/11\n", + "Current Ctt: 0.4364357804719848\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=831 (831 total) |\n", + "| EDM = 5.14E-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", + "Step: 8/11\n", + "Current Ctt: 0.4364357804719848\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.98E+05 | Ncalls=1364 (1364 total) |\n", + "| EDM = 3.96E-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", + "Time Taken: 3683\n", + "Predicted time left: 11.0 min, 9.0 s\n", + "Step: 9/11\n", + "Current Ctt: 0.46291004988627577\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.98E+05 | Ncalls=1822 (1822 total) |\n", + "| EDM = 5.91E-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", + "Step: 9/11\n", + "Current Ctt: 0.46291004988627577\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=1290 (1290 total) |\n", + "| EDM = 1.25E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 4233\n", + "Predicted time left: 6.0 min, 24.0 s\n", + "Step: 10/11\n", + "Current Ctt: 0.4879500364742666\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=1696 (1696 total) |\n", + "| EDM = 9.5E-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", + "Step: 10/11\n", + "Current Ctt: 0.4879500364742666\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1317 (1317 total) |\n", + "| EDM = 5.84E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 4835\n", + "Predicted time left: \n" ] } ], @@ -1993,7 +2432,7 @@ "Ctt_list = []\n", "Ctt_error_list = []\n", "\n", - "nr_of_toys = 1\n", + "nr_of_toys = 5\n", "nevents = int(pdg[\"number_of_decays\"]*cut_BR)\n", "# nevents = pdg[\"number_of_decays\"]\n", "event_stack = 1000000\n", @@ -2028,7 +2467,7 @@ " \n", " Nll_list.append([])\n", "\n", - " while len(Nll_list[-1]) <= nr_of_toys:\n", + " while len(Nll_list[-1]) < nr_of_toys:\n", " \n", " print('Step: {0}/{1}'.format(__, ste))\n", " \n", @@ -2043,14 +2482,15 @@ " Ctt.set_value(Ctt_step)\n", " else:\n", " Ctt.set_value(0.0)\n", + " \n", + " if floaty:\n", + " sampler.resample(n=nevents)\n", + " s = sampler.unstack_x()\n", + " total_samp = zfit.run(s)\n", + " calls = 0\n", + " c = 1\n", "\n", - " sampler.resample(n=nevents)\n", - " s = sampler.unstack_x()\n", - " total_samp = zfit.run(s)\n", - " calls = 0\n", - " c = 1\n", - "\n", - " data = zfit.data.Data.from_numpy(array=total_samp[:int(nevents)], obs=obs_fit)\n", + " data = zfit.data.Data.from_numpy(array=total_samp[:int(nevents)], obs=obs_fit)\n", "\n", " ### Fit data\n", "\n", @@ -2069,12 +2509,18 @@ " params = result.params\n", "\n", " if result.converged:\n", - " Nll_list[-1].append(result.fmin)\n" + " Nll_list[-1].append(result.fmin)\n", + " \n", + " _t = int(time.time()-start)\n", + " \n", + " print('Time Taken: {}'.format(display_time(int(_t))))\n", + " \n", + " print('Predicted time left: {}'.format(display_time(int((_t/(__+1)*(ste-__-1))))))\n" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 50, "metadata": {}, "outputs": [], "source": [ @@ -2082,36 +2528,38 @@ "\n", "CLs_values = []\n", "\n", - "for i in range(len(Nll_list)/2):\n", + "for i in range(int(len(Nll_list)/2)):\n", " CLs_values.append([])\n", " for j in range(nr_of_toys):\n", - " CLs_values[i].append(Nll_list[i][j]-Nll_list[i+1][j])\n", + " CLs_values[i].append(Nll_list[2*i][j]-Nll_list[2*i+1][j])\n", "\n", "if not os.path.exists(dirName):\n", " os.mkdir(dirName)\n", " print(\"Directory \" , dirName , \" Created \")\n", "\n", - "with open(\"'data/CLs/CLs_Nll_list.pkl\", \"wb\") as f:\n", + "with open(\"data/CLs/{}-{}_{}s--CLs_Nll_list.pkl\".format(mi,ma,ste), \"wb\") as f:\n", " pkl.dump(Nll_list, f, pkl.HIGHEST_PROTOCOL)\n", "\n", - "with open(\"'data/CLs/CLs_list.pkl\", \"wb\") as f:\n", + "with open(\"data/CLs/{}-{}_{}s--CLs_list.pkl\".format(mi,ma,ste), \"wb\") as f:\n", " pkl.dump(CLs_values, f, pkl.HIGHEST_PROTOCOL)" ] }, { "cell_type": "code", - "execution_count": 45, + "execution_count": 51, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ - " 1.0E-03\n" + "[[-158.64372490125243], [46.42150753445458], [15.524937722249888], [189.72476294625085], [-24.962952912901528], [375.9000742339995], [43.24977700418094], [-46.182621752261184], [-114.10263936838601], [73.08447201456875], [83.91141392168356]]\n" ] } ], - "source": [] + "source": [ + "print(CLs_values)" + ] }, { "cell_type": "code", diff --git a/__pycache__/pdg_const.cpython-37.pyc b/__pycache__/pdg_const.cpython-37.pyc index fbbaa17..fbac169 100644 --- a/__pycache__/pdg_const.cpython-37.pyc +++ b/__pycache__/pdg_const.cpython-37.pyc Binary files differ diff --git a/data/CLs/0.0-0.001_11s--CLs_Nll_list.pkl b/data/CLs/0.0-0.001_11s--CLs_Nll_list.pkl new file mode 100644 index 0000000..b0ae92d --- /dev/null +++ b/data/CLs/0.0-0.001_11s--CLs_Nll_list.pkl Binary files differ diff --git a/data/CLs/0.0-0.001_11s--CLs_list.pkl b/data/CLs/0.0-0.001_11s--CLs_list.pkl new file mode 100644 index 0000000..57987b3 --- /dev/null +++ b/data/CLs/0.0-0.001_11s--CLs_list.pkl Binary files differ diff --git a/raremodel-nb.ipynb b/raremodel-nb.ipynb index 2f54c56..18a9c1b 100644 --- a/raremodel-nb.ipynb +++ b/raremodel-nb.ipynb @@ -1861,19 +1861,7 @@ "cell_type": "code", "execution_count": 40, "metadata": {}, - "outputs": [ - { - "ename": "NameError", - "evalue": "name 'Ctt_list' is not defined", - "output_type": "error", - "traceback": [ - "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", - "\u001b[1;31mNameError\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 1\u001b[0m \u001b[1;32mwith\u001b[0m \u001b[0mopen\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"data/results/Ctt_list.pkl\"\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m\"wb\"\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 2\u001b[1;33m \u001b[0mpkl\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdump\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mCtt_list\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mpkl\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mHIGHEST_PROTOCOL\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 3\u001b[0m \u001b[1;32mwith\u001b[0m \u001b[0mopen\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"data/results/Ctt_error_list.pkl\"\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m\"wb\"\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 4\u001b[0m \u001b[0mpkl\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdump\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mCtt_error_list\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mpkl\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mHIGHEST_PROTOCOL\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", - "\u001b[1;31mNameError\u001b[0m: name 'Ctt_list' is not defined" - ] - } - ], + "outputs": [], "source": [ "# with open(\"data/results/Ctt_list.pkl\", \"wb\") as f:\n", "# pkl.dump(Ctt_list, f, pkl.HIGHEST_PROTOCOL)\n", @@ -1883,7 +1871,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 41, "metadata": {}, "outputs": [], "source": [ @@ -1895,7 +1883,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 42, "metadata": {}, "outputs": [], "source": [ @@ -1910,7 +1898,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 43, "metadata": {}, "outputs": [], "source": [ @@ -1919,7 +1907,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 44, "metadata": {}, "outputs": [], "source": [ @@ -1928,7 +1916,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 45, "metadata": {}, "outputs": [], "source": [ @@ -1937,7 +1925,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 46, "metadata": {}, "outputs": [], "source": [ @@ -1954,32 +1942,483 @@ { "cell_type": "code", "execution_count": 47, - "metadata": {}, + "metadata": { + "scrolled": false + }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ + "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: 0/11\n", "Current Ctt: 0.0\n", "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.977E+05 | Ncalls=1362 (1362 total) |\n", + "| EDM = 0.00848 (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", + "Step: 0/11\n", + "Current Ctt: 0.0\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.976E+05 | Ncalls=1372 (1372 total) |\n", + "| EDM = 0.00145 (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", + "Step: 0/11\n", + "Current Ctt: 0.0\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.977E+05 | Ncalls=1309 (1309 total) |\n", + "| EDM = 2.31E-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", + "Step: 0/11\n", + "Current Ctt: 0.0\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1172 (1172 total) |\n", + "| EDM = 0.000838 (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", + "Time Taken: 463\n", + "Predicted time left: 7.0 min\n", + "Step: 1/11\n", + "Current Ctt: 0.1543033499620919\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 3.026E+05 | Ncalls=505 (516 total) |\n", + "| EDM = 6.64E+08 (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", + "Step: 1/11\n", + "Current Ctt: 0.1543033499620919\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=1459 (1459 total) |\n", + "| EDM = 0.00013 (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", + "Step: 1/11\n", + "Current Ctt: 0.1543033499620919\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1300 (1300 total) |\n", + "| EDM = 1.54E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 813\n", + "Predicted time left: 11.0 min, 5.0 s\n", + "Step: 2/11\n", + "Current Ctt: 0.2182178902359924\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1331 (1331 total) |\n", + "| EDM = 6.73E-06 (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", + "Step: 2/11\n", + "Current Ctt: 0.2182178902359924\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1507 (1507 total) |\n", + "| EDM = 0.000231 (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", + "Time Taken: 1111\n", + "Predicted time left: 13.0 min, 28.0 s\n", + "Step: 3/11\n", + "Current Ctt: 0.2672612419124244\n", + "Ctt floating: True\n", "Toy 0/1\n" ] }, { - "ename": "RuntimeError", - "evalue": "exception was raised in user function\nUser function arguments:\n p4415_s = +0.738610\n psi2s_p = +1.596704\n p4415_p = +5.534902\n phi_s = +14.820266\n DDstar_s = +0.055634\n p3770_s = +3.056326\n p4160_s = +3.013799\n omega_s = +8.221504\n p3770_p = +3.227282\n phi_p = -1.956867\n Ctt = -0.074235\n DDstar_p = -3.658252\n rho_p = +3.220475\n Dbar_p = +1.482929\n p4040_p = -1.211697\n rho_s = +0.426552\n jpsi_p = -4.776017\n omega_p = +2.899561\n p4160_p = -1.589904\n p4040_s = +1.161008\n Dbar_s = +0.130791\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 76\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 77\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---> 78\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 79\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 80\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~\\.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 p4415_s = +0.738610\n psi2s_p = +1.596704\n p4415_p = +5.534902\n phi_s = +14.820266\n DDstar_s = +0.055634\n p3770_s = +3.056326\n p4160_s = +3.013799\n omega_s = +8.221504\n p3770_p = +3.227282\n phi_p = -1.956867\n Ctt = -0.074235\n DDstar_p = -3.658252\n rho_p = +3.220475\n Dbar_p = +1.482929\n p4040_p = -1.211697\n rho_s = +0.426552\n jpsi_p = -4.776017\n omega_p = +2.899561\n p4160_p = -1.589904\n p4040_s = +1.161008\n Dbar_s = +0.130791\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" + "name": "stdout", + "output_type": "stream", + "text": [ + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=1331 (1331 total) |\n", + "| EDM = 0.00103 (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", + "Step: 3/11\n", + "Current Ctt: 0.2672612419124244\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.98E+05 | Ncalls=1407 (1407 total) |\n", + "| EDM = 5.43E-06 (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", + "Step: 3/11\n", + "Current Ctt: 0.2672612419124244\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.989E+05 | Ncalls=566 (577 total) |\n", + "| EDM = 1.99E+08 (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", + "Step: 3/11\n", + "Current Ctt: 0.2672612419124244\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=366 (366 total) |\n", + "| EDM = 6.89E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 1656\n", + "Predicted time left: 17.0 min, 33.0 s\n", + "Step: 4/11\n", + "Current Ctt: 0.3086066999241838\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1364 (1364 total) |\n", + "| EDM = 2.06E-06 (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", + "Step: 4/11\n", + "Current Ctt: 0.3086066999241838\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1497 (1497 total) |\n", + "| EDM = 5.28E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 2044\n", + "Predicted time left: 18.0 min, 34.0 s\n", + "Step: 5/11\n", + "Current Ctt: 0.3450327796711771\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.98E+05 | Ncalls=1249 (1249 total) |\n", + "| EDM = 0.00017 (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", + "Step: 5/11\n", + "Current Ctt: 0.3450327796711771\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.976E+05 | Ncalls=1276 (1276 total) |\n", + "| EDM = 0.000344 (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", + "Time Taken: 2428\n", + "Predicted time left: 18.0 min, 23.0 s\n", + "Step: 6/11\n", + "Current Ctt: 0.37796447300922725\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1446 (1446 total) |\n", + "| EDM = 0.000166 (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", + "Step: 6/11\n", + "Current Ctt: 0.37796447300922725\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=792 (792 total) |\n", + "| EDM = 7.63E-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", + "Time Taken: 2836\n", + "Predicted time left: 17.0 min, 11.0 s\n", + "Step: 7/11\n", + "Current Ctt: 0.408248290463863\n", + "Ctt floating: True\n", + "Toy 0/1\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=1391 (1391 total) |\n", + "| EDM = 0.000185 (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", + "Step: 7/11\n", + "Current Ctt: 0.408248290463863\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=273 (273 total) |\n", + "| EDM = 7.9E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 3232\n", + "Predicted time left: 14.0 min, 41.0 s\n", + "Step: 8/11\n", + "Current Ctt: 0.4364357804719848\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=831 (831 total) |\n", + "| EDM = 5.14E-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", + "Step: 8/11\n", + "Current Ctt: 0.4364357804719848\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.98E+05 | Ncalls=1364 (1364 total) |\n", + "| EDM = 3.96E-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", + "Time Taken: 3683\n", + "Predicted time left: 11.0 min, 9.0 s\n", + "Step: 9/11\n", + "Current Ctt: 0.46291004988627577\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.98E+05 | Ncalls=1822 (1822 total) |\n", + "| EDM = 5.91E-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", + "Step: 9/11\n", + "Current Ctt: 0.46291004988627577\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=1290 (1290 total) |\n", + "| EDM = 1.25E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 4233\n", + "Predicted time left: 6.0 min, 24.0 s\n", + "Step: 10/11\n", + "Current Ctt: 0.4879500364742666\n", + "Ctt floating: True\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.979E+05 | Ncalls=1696 (1696 total) |\n", + "| EDM = 9.5E-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", + "Step: 10/11\n", + "Current Ctt: 0.4879500364742666\n", + "Ctt floating: False\n", + "Toy 0/1\n", + "------------------------------------------------------------------\n", + "| FCN = 2.978E+05 | Ncalls=1317 (1317 total) |\n", + "| EDM = 5.84E-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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Time Taken: 4835\n", + "Predicted time left: \n" ] } ], @@ -1993,7 +2432,7 @@ "Ctt_list = []\n", "Ctt_error_list = []\n", "\n", - "nr_of_toys = 1\n", + "nr_of_toys = 5\n", "nevents = int(pdg[\"number_of_decays\"]*cut_BR)\n", "# nevents = pdg[\"number_of_decays\"]\n", "event_stack = 1000000\n", @@ -2028,7 +2467,7 @@ " \n", " Nll_list.append([])\n", "\n", - " while len(Nll_list[-1]) <= nr_of_toys:\n", + " while len(Nll_list[-1]) < nr_of_toys:\n", " \n", " print('Step: {0}/{1}'.format(__, ste))\n", " \n", @@ -2043,14 +2482,15 @@ " Ctt.set_value(Ctt_step)\n", " else:\n", " Ctt.set_value(0.0)\n", + " \n", + " if floaty:\n", + " sampler.resample(n=nevents)\n", + " s = sampler.unstack_x()\n", + " total_samp = zfit.run(s)\n", + " calls = 0\n", + " c = 1\n", "\n", - " sampler.resample(n=nevents)\n", - " s = sampler.unstack_x()\n", - " total_samp = zfit.run(s)\n", - " calls = 0\n", - " c = 1\n", - "\n", - " data = zfit.data.Data.from_numpy(array=total_samp[:int(nevents)], obs=obs_fit)\n", + " data = zfit.data.Data.from_numpy(array=total_samp[:int(nevents)], obs=obs_fit)\n", "\n", " ### Fit data\n", "\n", @@ -2069,12 +2509,18 @@ " params = result.params\n", "\n", " if result.converged:\n", - " Nll_list[-1].append(result.fmin)\n" + " Nll_list[-1].append(result.fmin)\n", + " \n", + " _t = int(time.time()-start)\n", + " \n", + " print('Time Taken: {}'.format(display_time(int(_t))))\n", + " \n", + " print('Predicted time left: {}'.format(display_time(int((_t/(__+1)*(ste-__-1))))))\n" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 50, "metadata": {}, "outputs": [], "source": [ @@ -2082,36 +2528,38 @@ "\n", "CLs_values = []\n", "\n", - "for i in range(len(Nll_list)/2):\n", + "for i in range(int(len(Nll_list)/2)):\n", " CLs_values.append([])\n", " for j in range(nr_of_toys):\n", - " CLs_values[i].append(Nll_list[i][j]-Nll_list[i+1][j])\n", + " CLs_values[i].append(Nll_list[2*i][j]-Nll_list[2*i+1][j])\n", "\n", "if not os.path.exists(dirName):\n", " os.mkdir(dirName)\n", " print(\"Directory \" , dirName , \" Created \")\n", "\n", - "with open(\"'data/CLs/CLs_Nll_list.pkl\", \"wb\") as f:\n", + "with open(\"data/CLs/{}-{}_{}s--CLs_Nll_list.pkl\".format(mi,ma,ste), \"wb\") as f:\n", " pkl.dump(Nll_list, f, pkl.HIGHEST_PROTOCOL)\n", "\n", - "with open(\"'data/CLs/CLs_list.pkl\", \"wb\") as f:\n", + "with open(\"data/CLs/{}-{}_{}s--CLs_list.pkl\".format(mi,ma,ste), \"wb\") as f:\n", " pkl.dump(CLs_values, f, pkl.HIGHEST_PROTOCOL)" ] }, { "cell_type": "code", - "execution_count": 45, + "execution_count": 51, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ - " 1.0E-03\n" + "[[-158.64372490125243], [46.42150753445458], [15.524937722249888], [189.72476294625085], [-24.962952912901528], [375.9000742339995], [43.24977700418094], [-46.182621752261184], [-114.10263936838601], [73.08447201456875], [83.91141392168356]]\n" ] } ], - "source": [] + "source": [ + "print(CLs_values)" + ] }, { "cell_type": "code",