diff --git a/__pycache__/pdg_const.cpython-37.pyc b/__pycache__/pdg_const.cpython-37.pyc index e91547c..8eb1fd0 100644 --- a/__pycache__/pdg_const.cpython-37.pyc +++ b/__pycache__/pdg_const.cpython-37.pyc Binary files differ diff --git a/data/zfit_toys/toy_0/0.pkl b/data/zfit_toys/toy_0/0.pkl index 57b3cea..4ff68ba 100644 --- a/data/zfit_toys/toy_0/0.pkl +++ b/data/zfit_toys/toy_0/0.pkl Binary files differ diff --git a/data/zfit_toys/toy_0/1.pkl b/data/zfit_toys/toy_0/1.pkl index ea1ce79..dc62767 100644 --- a/data/zfit_toys/toy_0/1.pkl +++ b/data/zfit_toys/toy_0/1.pkl Binary files differ diff --git a/data/zfit_toys/toy_0/2.pkl b/data/zfit_toys/toy_0/2.pkl index 370b3f4..a58bfe8 100644 --- a/data/zfit_toys/toy_0/2.pkl +++ b/data/zfit_toys/toy_0/2.pkl Binary files differ diff --git a/data/zfit_toys/toy_0/3.pkl b/data/zfit_toys/toy_0/3.pkl index aa74233..d0626f8 100644 --- a/data/zfit_toys/toy_0/3.pkl +++ b/data/zfit_toys/toy_0/3.pkl Binary files differ diff --git a/data/zfit_toys/toy_0/4.pkl b/data/zfit_toys/toy_0/4.pkl index 1f325cd..7ff11ac 100644 --- a/data/zfit_toys/toy_0/4.pkl +++ b/data/zfit_toys/toy_0/4.pkl Binary files differ diff --git a/data/zfit_toys/toy_0/5.pkl b/data/zfit_toys/toy_0/5.pkl index 6b4f129..c17eb4e 100644 --- a/data/zfit_toys/toy_0/5.pkl +++ b/data/zfit_toys/toy_0/5.pkl Binary files differ diff --git a/raremodel-nb.ipynb b/raremodel-nb.ipynb index cf17488..8b9d01e 100644 --- a/raremodel-nb.ipynb +++ b/raremodel-nb.ipynb @@ -598,7 +598,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 9, @@ -1323,7 +1323,7 @@ }, { "cell_type": "code", - "execution_count": 37, + "execution_count": 39, "metadata": { "scrolled": false }, @@ -1332,20 +1332,235 @@ "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", - "WARNING:tensorflow:From C:\\Users\\sa_li\\.conda\\envs\\rmd\\lib\\site-packages\\tensorflow_probability\\python\\distributions\\categorical.py:263: multinomial (from tensorflow.python.ops.random_ops) is deprecated and will be removed in a future version.\n", - "Instructions for updating:\n", - "Use tf.random.categorical instead.\n", "Toy 0: Generating data...\n", "Toy 0: Data generation finished\n", "Toy 0: Loading data...\n", "Toy 0: Loading data finished\n", "Toy 0: Fitting pdf...\n", "------------------------------------------------------------------\n", - "| FCN = 4.133E+04 | Ncalls=95 (106 total) |\n", - "| EDM = 12.2 (Goal: 5E-06) | up = 0.5 |\n", + "| FCN = 4.145E+04 | Ncalls=189 (189 total) |\n", + "| EDM = 0.000434 (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", + "Toy 0: Fitting finished\n", + "Function minimum: 41446.17176214419\n", + "---------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "---------------------------------------------------------------------------------------------\n", + "| 0 | psi2s_p | -4.4 | 0.4 | | |-6.28319 | 6.28319 | |\n", + "| 1 | Dbar_s | -0.7 | 1.8 | | | -1.464 | 1.464 | |\n", + "| 2 | Ctt | -0.5 | 1.0 | | | -0.5 | 0.5 | |\n", + "| 3 | jpsi_p | -1.9 | 0.7 | | |-6.28319 | 6.28319 | |\n", + "| 4 | jpsi_s | 0.980E4 | 0.011E4 | | | 9797.52 | 9996.48 | |\n", + "| 5 | psi2s_s | 1430 | 50 | | | 1358.64 | 1433.36 | |\n", + "---------------------------------------------------------------------------------------------\n", + "-------------------------------------------------------------\n", + "| | psi2s_p Dbar_s Ctt jpsi_p jpsi_s psi2s_s |\n", + "-------------------------------------------------------------\n", + "| psi2s_p | 1.000 0.896 0.021 0.977 -0.145 -0.023 |\n", + "| Dbar_s | 0.896 1.000 0.236 0.941 -0.223 0.004 |\n", + "| Ctt | 0.021 0.236 1.000 0.013 -0.053 -0.007 |\n", + "| jpsi_p | 0.977 0.941 0.013 1.000 -0.133 -0.004 |\n", + "| jpsi_s | -0.145 -0.223 -0.053 -0.133 1.000 -0.001 |\n", + "| psi2s_s | -0.023 0.004 -0.007 -0.004 -0.001 1.000 |\n", + "-------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 0.41840024605493187}), (, {'error': 1.8317143176053508}), (, {'error': 0.996289918615876}), (, {'error': 0.7473453412456892}), (, {'error': 113.9126570143626}), (, {'error': 47.11513567257623})])\n", + "Toy 1/10\n", + "Time taken: 2 min, 32 s\n", + "Projected time left: 22 min, 50 s\n", + "Toy 1: Generating data...\n", + "Toy 1: Data generation finished\n", + "Toy 1: Loading data...\n", + "Toy 1: Loading data finished\n", + "Toy 1: Fitting pdf...\n", + "------------------------------------------------------------------\n", + "| FCN = 4.145E+04 | Ncalls=174 (174 total) |\n", + "| EDM = 0.000246 (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", + "Toy 1: Fitting finished\n", + "Function minimum: 41446.21391501494\n", + "---------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "---------------------------------------------------------------------------------------------\n", + "| 0 | psi2s_p | 1.51 | 0.17 | | |-6.28319 | 6.28319 | |\n", + "| 1 | Dbar_s | -0.7 | 1.8 | | | -1.464 | 1.464 | |\n", + "| 2 | Ctt | -0.5 | 0.8 | | | -0.5 | 0.5 | |\n", + "| 3 | jpsi_p | 2.0 | 0.5 | | |-6.28319 | 6.28319 | |\n", + "| 4 | jpsi_s | 0.980E4 | 0.013E4 | | | 9797.52 | 9996.48 | |\n", + "| 5 | psi2s_s | 1430 | 50 | | | 1358.64 | 1433.36 | |\n", + "---------------------------------------------------------------------------------------------\n", + "-------------------------------------------------------------\n", + "| | psi2s_p Dbar_s Ctt jpsi_p jpsi_s psi2s_s |\n", + "-------------------------------------------------------------\n", + "| psi2s_p | 1.00 0.80 0.84 0.72 0.16 -0.27 |\n", + "| Dbar_s | 0.80 1.00 0.98 0.88 0.17 -0.22 |\n", + "| Ctt | 0.84 0.98 1.00 0.94 0.21 -0.23 |\n", + "| jpsi_p | 0.72 0.88 0.94 1.00 0.17 -0.22 |\n", + "| jpsi_s | 0.16 0.17 0.21 0.17 1.00 -0.05 |\n", + "| psi2s_s | -0.27 -0.22 -0.23 -0.22 -0.05 1.00 |\n", + "-------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 0.17027941762338372}), (, {'error': 1.7969490420133767}), (, {'error': 0.8149920548026126}), (, {'error': 0.4859306435087767}), (, {'error': 129.95009225722333}), (, {'error': 47.350217126767575})])\n", + "Toy 2/10\n", + "Time taken: 5 min, 16 s\n", + "Projected time left: 21 min, 7 s\n", + "Toy 2: Generating data...\n", + "Toy 2: Data generation finished\n", + "Toy 2: Loading data...\n", + "Toy 2: Loading data finished\n", + "Toy 2: Fitting pdf...\n", + "------------------------------------------------------------------\n", + "| FCN = 4.145E+04 | Ncalls=204 (204 total) |\n", + "| EDM = 0.00128 (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 | True | True | False |\n", + "------------------------------------------------------------------\n", + "Toy 2: Fitting finished\n", + "Function minimum: 41446.21446623224\n", + "---------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "---------------------------------------------------------------------------------------------\n", + "| 0 | psi2s_p | -4.78 | 0.18 | | |-6.28319 | 6.28319 | |\n", + "| 1 | Dbar_s | -0.7 | 1.8 | | | -1.464 | 1.464 | |\n", + "| 2 | Ctt | -0.5 | 0.7 | | | -0.5 | 0.5 | |\n", + "| 3 | jpsi_p | -4.3 | 0.5 | | |-6.28319 | 6.28319 | |\n", + "| 4 | jpsi_s | 0.980E4 | 0.013E4 | | | 9797.52 | 9996.48 | |\n", + "| 5 | psi2s_s | 1430 | 40 | | | 1358.64 | 1433.36 | |\n", + "---------------------------------------------------------------------------------------------\n", + "-------------------------------------------------------------\n", + "| | psi2s_p Dbar_s Ctt jpsi_p jpsi_s psi2s_s |\n", + "-------------------------------------------------------------\n", + "| psi2s_p | 1.00 0.82 0.86 0.76 0.20 0.32 |\n", + "| Dbar_s | 0.82 1.00 0.98 0.89 0.21 0.26 |\n", + "| Ctt | 0.86 0.98 1.00 0.95 0.24 0.27 |\n", + "| jpsi_p | 0.76 0.89 0.95 1.00 0.21 0.26 |\n", + "| jpsi_s | 0.20 0.21 0.24 0.21 1.00 0.07 |\n", + "| psi2s_s | 0.32 0.26 0.27 0.26 0.07 1.00 |\n", + "-------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 0.17884003769401335}), (, {'error': 1.7953075848497315}), (, {'error': 0.6507663138579277}), (, {'error': 0.49871222470060794}), (, {'error': 127.32075889710268}), (, {'error': 38.958773985027506})])\n", + "Toy 3/10\n", + "Time taken: 7 min, 58 s\n", + "Projected time left: 18 min, 37 s\n", + "Toy 3: Generating data...\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Toy 3: Data generation finished\n", + "Toy 3: Loading data...\n", + "Toy 3: Loading data finished\n", + "Toy 3: Fitting pdf...\n", + "------------------------------------------------------------------\n", + "| FCN = 4.145E+04 | Ncalls=233 (233 total) |\n", + "| EDM = 0.000301 (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", + "Toy 3: Fitting finished\n", + "Function minimum: 41446.213143894536\n", + "---------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "---------------------------------------------------------------------------------------------\n", + "| 0 | psi2s_p | -4.77 | 0.20 | | |-6.28319 | 6.28319 | |\n", + "| 1 | Dbar_s | -0.7 | 1.8 | | | -1.464 | 1.464 | |\n", + "| 2 | Ctt | -0.5 | 1.0 | | | -0.5 | 0.5 | |\n", + "| 3 | jpsi_p | -4.3 | 0.5 | | |-6.28319 | 6.28319 | |\n", + "| 4 | jpsi_s | 0.980E4 | 0.012E4 | | | 9797.52 | 9996.48 | |\n", + "| 5 | psi2s_s | 1430 | 50 | | | 1358.64 | 1433.36 | |\n", + "---------------------------------------------------------------------------------------------\n", + "-------------------------------------------------------------\n", + "| | psi2s_p Dbar_s Ctt jpsi_p jpsi_s psi2s_s |\n", + "-------------------------------------------------------------\n", + "| psi2s_p | 1.00 0.89 0.92 0.85 0.41 -0.32 |\n", + "| Dbar_s | 0.89 1.00 0.98 0.93 0.43 -0.26 |\n", + "| Ctt | 0.92 0.98 1.00 0.97 0.46 -0.27 |\n", + "| jpsi_p | 0.85 0.93 0.97 1.00 0.43 -0.26 |\n", + "| jpsi_s | 0.41 0.43 0.46 0.43 1.00 -0.12 |\n", + "| psi2s_s | -0.32 -0.26 -0.27 -0.26 -0.12 1.00 |\n", + "-------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 0.2028644492314955}), (, {'error': 1.8154514362409575}), (, {'error': 0.9807202234423524}), (, {'error': 0.5159537418531039}), (, {'error': 117.77974601317328}), (, {'error': 45.73863717851896})])\n", + "Toy 4/10\n", + "Time taken: 10 min, 46 s\n", + "Projected time left: 16 min, 9 s\n", + "Toy 4: Generating data...\n", + "Toy 4: Data generation finished\n", + "Toy 4: Loading data...\n", + "Toy 4: Loading data finished\n", + "Toy 4: Fitting pdf...\n", + "------------------------------------------------------------------\n", + "| FCN = 4.145E+04 | Ncalls=180 (180 total) |\n", + "| EDM = 0.000259 (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", + "Toy 4: Fitting finished\n", + "Function minimum: 41446.213230540496\n", + "---------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "---------------------------------------------------------------------------------------------\n", + "| 0 | psi2s_p | -4.77 | 0.21 | | |-6.28319 | 6.28319 | |\n", + "| 1 | Dbar_s | -0.7 | 1.8 | | | -1.464 | 1.464 | |\n", + "| 2 | Ctt | -0.5 | 0.7 | | | -0.5 | 0.5 | |\n", + "| 3 | jpsi_p | -4.3 | 0.5 | | |-6.28319 | 6.28319 | |\n", + "| 4 | jpsi_s | 0.980E4 | 0.012E4 | | | 9797.52 | 9996.48 | |\n", + "| 5 | psi2s_s | 1430 | 60 | | | 1358.64 | 1433.36 | |\n", + "---------------------------------------------------------------------------------------------\n", + "-------------------------------------------------------------\n", + "| | psi2s_p Dbar_s Ctt jpsi_p jpsi_s psi2s_s |\n", + "-------------------------------------------------------------\n", + "| psi2s_p | 1.00 0.89 0.92 0.87 0.45 -0.45 |\n", + "| Dbar_s | 0.89 1.00 0.98 0.93 0.47 -0.37 |\n", + "| Ctt | 0.92 0.98 1.00 0.97 0.50 -0.38 |\n", + "| jpsi_p | 0.87 0.93 0.97 1.00 0.47 -0.38 |\n", + "| jpsi_s | 0.45 0.47 0.50 0.47 1.00 -0.19 |\n", + "| psi2s_s | -0.45 -0.37 -0.38 -0.38 -0.19 1.00 |\n", + "-------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 0.21090063108815427}), (, {'error': 1.8277281246469832}), (, {'error': 0.7144740345016457}), (, {'error': 0.5130723061292268}), (, {'error': 120.0500270283901}), (, {'error': 59.20149657515435})])\n", + "Toy 5/10\n", + "Time taken: 13 min, 31 s\n", + "Projected time left: 13 min, 31 s\n", + "Toy 5: Generating data...\n", + "Toy 5: Data generation finished\n", + "Toy 5: Loading data...\n", + "Toy 5: Loading data finished\n", + "Toy 5: Fitting pdf...\n", + "------------------------------------------------------------------\n", + "| FCN = 4.145E+04 | Ncalls=304 (304 total) |\n", + "| EDM = 0.00152 (Goal: 5E-06) | up = 0.5 |\n", "------------------------------------------------------------------\n", "| Valid Min. | Valid Param. | Above EDM | Reached call limit |\n", "------------------------------------------------------------------\n", @@ -1355,31 +1570,215 @@ "------------------------------------------------------------------\n", "| False | True | False | False | True |\n", "------------------------------------------------------------------\n", - "Toy 0: Fitting finished\n", - "Function minimum: 41332.70666221551\n", + "Toy 5: Fitting finished\n", + "Function minimum: 41446.26376718663\n", "---------------------------------------------------------------------------------------------\n", "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", "---------------------------------------------------------------------------------------------\n", - "| 0 | Dbar_s | -0.6 | 1.6 | | | -1.464 | 1.464 | |\n", - "| 1 | jpsi_p | -1.19 | 0.27 | | |-6.28319 | 6.28319 | |\n", - "| 2 | psi2s_p | 2.7 | 1.1 | | |-6.28319 | 6.28319 | |\n", - "| 3 | Ctt | 0.12 | 0.80 | | | -0.5 | 0.5 | |\n", - "| 4 | psi2s_s | 1360 | 70 | | | 1358.64 | 1433.36 | |\n", - "| 5 | jpsi_s | 0.987E4 | 0.012E4 | | | 9797.52 | 9996.48 | |\n", + "| 0 | psi2s_p | 4.67 | 0.32 | | |-6.28319 | 6.28319 | |\n", + "| 1 | Dbar_s | -0.8 | 1.8 | | | -1.464 | 1.464 | |\n", + "| 2 | Ctt | -0.5 | 0.5 | | | -0.5 | 0.5 | |\n", + "| 3 | jpsi_p | -2.1 | 0.5 | | |-6.28319 | 6.28319 | |\n", + "| 4 | jpsi_s | 0.980E4 | 0.015E4 | | | 9797.52 | 9996.48 | |\n", + "| 5 | psi2s_s | 1430 | 50 | | | 1358.64 | 1433.36 | |\n", "---------------------------------------------------------------------------------------------\n", "-------------------------------------------------------------\n", - "| | Dbar_s jpsi_p psi2s_p Ctt psi2s_s jpsi_s |\n", + "| | psi2s_p Dbar_s Ctt jpsi_p jpsi_s psi2s_s |\n", "-------------------------------------------------------------\n", - "| Dbar_s | 1.00 -0.82 -0.96 -0.90 0.29 0.97 |\n", - "| jpsi_p | -0.82 1.00 0.91 0.61 -0.25 -0.76 |\n", - "| psi2s_p | -0.96 0.91 1.00 0.77 -0.28 -0.95 |\n", - "| Ctt | -0.90 0.61 0.77 1.00 -0.20 -0.81 |\n", - "| psi2s_s | 0.29 -0.25 -0.28 -0.20 1.00 0.29 |\n", - "| jpsi_s | 0.97 -0.76 -0.95 -0.81 0.29 1.00 |\n", + "| psi2s_p | 1.000 -0.915 -0.954 0.850 -0.143 0.107 |\n", + "| Dbar_s | -0.915 1.000 0.967 -0.848 0.117 -0.004 |\n", + "| Ctt | -0.954 0.967 1.000 -0.937 0.157 -0.019 |\n", + "| jpsi_p | 0.850 -0.848 -0.937 1.000 -0.126 0.030 |\n", + "| jpsi_s | -0.143 0.117 0.157 -0.126 1.000 -0.005 |\n", + "| psi2s_s | 0.107 -0.004 -0.019 0.030 -0.005 1.000 |\n", "-------------------------------------------------------------\n", - "Hesse errors: OrderedDict([(, {'error': 1.5923595875666954}), (, {'error': 0.27012452493223504}), (, {'error': 1.1017684775063215}), (, {'error': 0.8037764107824681}), (, {'error': 73.13977088861418}), (, {'error': 122.7425419080937})])\n", - "Toy 1/1\n", - "Time taken: 11 min, 47 s\n", + "Hesse errors: OrderedDict([(, {'error': 0.31586736018571226}), (, {'error': 1.7931181136417407}), (, {'error': 0.505275767707557}), (, {'error': 0.5046726165950495}), (, {'error': 147.58315697207308}), (, {'error': 45.95130250092461})])\n", + "Toy 6/10\n", + "Time taken: 16 min, 29 s\n", + "Projected time left: 10 min, 59 s\n", + "Toy 6: Generating data...\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Toy 6: Data generation finished\n", + "Toy 6: Loading data...\n", + "Toy 6: Loading data finished\n", + "Toy 6: Fitting pdf...\n", + "------------------------------------------------------------------\n", + "| FCN = 4.145E+04 | Ncalls=216 (216 total) |\n", + "| EDM = 0.000883 (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", + "Toy 6: Fitting finished\n", + "Function minimum: 41446.21393607266\n", + "---------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "---------------------------------------------------------------------------------------------\n", + "| 0 | psi2s_p | -4.77 | 0.19 | | |-6.28319 | 6.28319 | |\n", + "| 1 | Dbar_s | -0.7 | 1.8 | | | -1.464 | 1.464 | |\n", + "| 2 | Ctt | -0.5 | 0.5 | | | -0.5 | 0.5 | |\n", + "| 3 | jpsi_p | 2.0 | 0.5 | | |-6.28319 | 6.28319 | |\n", + "| 4 | jpsi_s | 0.980E4 | 0.010E4 | | | 9797.52 | 9996.48 | |\n", + "| 5 | psi2s_s | 1430 | 40 | | | 1358.64 | 1433.36 | |\n", + "---------------------------------------------------------------------------------------------\n", + "-------------------------------------------------------------\n", + "| | psi2s_p Dbar_s Ctt jpsi_p jpsi_s psi2s_s |\n", + "-------------------------------------------------------------\n", + "| psi2s_p | 1.00 0.87 0.90 0.83 0.49 0.16 |\n", + "| Dbar_s | 0.87 1.00 0.98 0.92 0.51 0.13 |\n", + "| Ctt | 0.90 0.98 1.00 0.96 0.56 0.13 |\n", + "| jpsi_p | 0.83 0.92 0.96 1.00 0.51 0.13 |\n", + "| jpsi_s | 0.49 0.51 0.56 0.51 1.00 0.07 |\n", + "| psi2s_s | 0.16 0.13 0.13 0.13 0.07 1.00 |\n", + "-------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 0.1917856359807728}), (, {'error': 1.8156013659952255}), (, {'error': 0.5117900295417959}), (, {'error': 0.49818141138885563}), (, {'error': 100.16204116151494}), (, {'error': 38.14017204277275})])\n", + "Toy 7/10\n", + "Time taken: 19 min, 23 s\n", + "Projected time left: 8 min, 18 s\n", + "Toy 7: Generating data...\n", + "Toy 7: Data generation finished\n", + "Toy 7: Loading data...\n", + "Toy 7: Loading data finished\n", + "Toy 7: Fitting pdf...\n", + "------------------------------------------------------------------\n", + "| FCN = 4.145E+04 | Ncalls=202 (202 total) |\n", + "| EDM = 0.000126 (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", + "Toy 7: Fitting finished\n", + "Function minimum: 41446.171485007675\n", + "---------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "---------------------------------------------------------------------------------------------\n", + "| 0 | psi2s_p | 1.9 | 0.4 | | |-6.28319 | 6.28319 | |\n", + "| 1 | Dbar_s | -0.7 | 1.8 | | | -1.464 | 1.464 | |\n", + "| 2 | Ctt | -0.5 | 1.0 | | | -0.5 | 0.5 | |\n", + "| 3 | jpsi_p | -1.9 | 0.7 | | |-6.28319 | 6.28319 | |\n", + "| 4 | jpsi_s | 0.980E4 | 0.012E4 | | | 9797.52 | 9996.48 | |\n", + "| 5 | psi2s_s | 1430 | 60 | | | 1358.64 | 1433.36 | |\n", + "---------------------------------------------------------------------------------------------\n", + "-------------------------------------------------------------\n", + "| | psi2s_p Dbar_s Ctt jpsi_p jpsi_s psi2s_s |\n", + "-------------------------------------------------------------\n", + "| psi2s_p | 1.000 0.861 0.035 0.970 -0.112 -0.115 |\n", + "| Dbar_s | 0.861 1.000 0.322 0.916 -0.184 0.022 |\n", + "| Ctt | 0.035 0.322 1.000 0.018 -0.058 -0.048 |\n", + "| jpsi_p | 0.970 0.916 0.018 1.000 -0.102 -0.013 |\n", + "| jpsi_s | -0.112 -0.184 -0.058 -0.102 1.000 -0.004 |\n", + "| psi2s_s | -0.115 0.022 -0.048 -0.013 -0.004 1.000 |\n", + "-------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 0.39988882358765476}), (, {'error': 1.8338950485749796}), (, {'error': 0.9813742089966919}), (, {'error': 0.709065375300536}), (, {'error': 116.26239022795016}), (, {'error': 64.86664680863748})])\n", + "Toy 8/10\n", + "Time taken: 22 min, 12 s\n", + "Projected time left: 5 min, 33 s\n", + "Toy 8: Generating data...\n", + "Toy 8: Data generation finished\n", + "Toy 8: Loading data...\n", + "Toy 8: Loading data finished\n", + "Toy 8: Fitting pdf...\n", + "------------------------------------------------------------------\n", + "| FCN = 4.145E+04 | Ncalls=185 (185 total) |\n", + "| EDM = 7.47E-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", + "Toy 8: Fitting finished\n", + "Function minimum: 41446.171426607005\n", + "---------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "---------------------------------------------------------------------------------------------\n", + "| 0 | psi2s_p | -4.4 | 0.4 | | |-6.28319 | 6.28319 | |\n", + "| 1 | Dbar_s | -0.7 | 1.8 | | | -1.464 | 1.464 | |\n", + "| 2 | Ctt | -0.5 | 1.0 | | | -0.5 | 0.5 | |\n", + "| 3 | jpsi_p | 4.4 | 0.7 | | |-6.28319 | 6.28319 | |\n", + "| 4 | jpsi_s | 0.980E4 | 0.012E4 | | | 9797.52 | 9996.48 | |\n", + "| 5 | psi2s_s | 1430 | 60 | | | 1358.64 | 1433.36 | |\n", + "---------------------------------------------------------------------------------------------\n", + "-------------------------------------------------------------\n", + "| | psi2s_p Dbar_s Ctt jpsi_p jpsi_s psi2s_s |\n", + "-------------------------------------------------------------\n", + "| psi2s_p | 1.000 0.850 0.027 0.966 -0.108 -0.122 |\n", + "| Dbar_s | 0.850 1.000 0.309 0.913 -0.184 0.021 |\n", + "| Ctt | 0.027 0.309 1.000 0.012 -0.053 -0.042 |\n", + "| jpsi_p | 0.966 0.913 0.012 1.000 -0.097 -0.019 |\n", + "| jpsi_s | -0.108 -0.184 -0.053 -0.097 1.000 -0.004 |\n", + "| psi2s_s | -0.122 0.021 -0.042 -0.019 -0.004 1.000 |\n", + "-------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 0.37581893654929166}), (, {'error': 1.8374920961421453}), (, {'error': 0.995543793163386}), (, {'error': 0.6592821643495412}), (, {'error': 118.1086035526505}), (, {'error': 62.73818368176751})])\n", + "Toy 9/10\n", + "Time taken: 25 min, 2 s\n", + "Projected time left: 2 min, 46 s\n", + "Toy 9: Generating data...\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Toy 9: Data generation finished\n", + "Toy 9: Loading data...\n", + "Toy 9: Loading data finished\n", + "Toy 9: Fitting pdf...\n", + "------------------------------------------------------------------\n", + "| FCN = 4.145E+04 | Ncalls=181 (181 total) |\n", + "| EDM = 0.000129 (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", + "Toy 9: Fitting finished\n", + "Function minimum: 41446.27791906745\n", + "---------------------------------------------------------------------------------------------\n", + "| | Name | Value | Hesse Err | Minos Err- | Minos Err+ | Limit- | Limit+ | Fixed |\n", + "---------------------------------------------------------------------------------------------\n", + "| 0 | psi2s_p | -1.94 | 0.26 | | |-6.28319 | 6.28319 | |\n", + "| 1 | Dbar_s | -0.7 | 0.8 | | | -1.464 | 1.464 | |\n", + "| 2 | Ctt | -0.5 | 0.8 | | | -0.5 | 0.5 | |\n", + "| 3 | jpsi_p | 1.8 | 0.4 | | |-6.28319 | 6.28319 | |\n", + "| 4 | jpsi_s | 0.980E4 | 0.014E4 | | | 9797.52 | 9996.48 | |\n", + "| 5 | psi2s_s | 1430 | 40 | | | 1358.64 | 1433.36 | |\n", + "---------------------------------------------------------------------------------------------\n", + "-------------------------------------------------------------\n", + "| | psi2s_p Dbar_s Ctt jpsi_p jpsi_s psi2s_s |\n", + "-------------------------------------------------------------\n", + "| psi2s_p | 1.000 -0.766 -0.254 0.902 0.074 0.130 |\n", + "| Dbar_s | -0.766 1.000 0.391 -0.817 -0.169 0.028 |\n", + "| Ctt | -0.254 0.391 1.000 -0.005 -0.040 -0.028 |\n", + "| jpsi_p | 0.902 -0.817 -0.005 1.000 0.041 0.019 |\n", + "| jpsi_s | 0.074 -0.169 -0.040 0.041 1.000 -0.004 |\n", + "| psi2s_s | 0.130 0.028 -0.028 0.019 -0.004 1.000 |\n", + "-------------------------------------------------------------\n", + "Hesse errors: OrderedDict([(, {'error': 0.2574778505665698}), (, {'error': 0.8041361051173593}), (, {'error': 0.800071724878314}), (, {'error': 0.40940340284378474}), (, {'error': 136.70046081025066}), (, {'error': 41.38771918974419})])\n", + "Toy 10/10\n", + "Time taken: 27 min, 51 s\n", "Projected time left: \n" ] } @@ -1390,7 +1789,7 @@ "Ctt_list = []\n", "Ctt_error_list = []\n", "\n", - "nr_of_toys = 1\n", + "nr_of_toys = 10\n", "nevents = int(pdg[\"number_of_decays\"])\n", "nevents = pdg[\"number_of_decays\"]\n", "event_stack = 1000000\n", @@ -1480,15 +1879,15 @@ }, { "cell_type": "code", - "execution_count": 38, + "execution_count": 40, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ - "Mean Ctt value = 0.11936585533648691\n", - "Mean Ctt error = 0.8037764107824681\n" + "Mean Ctt value = -0.4955320725102208\n", + "Mean Ctt error = 0.7951298069508159\n" ] } ],