From 500b941f078aa5772f813f6f8659b3dc40fd1afc Mon Sep 17 00:00:00 2001 From: rkansal47 Date: Mon, 22 Jul 2024 07:55:16 -0700 Subject: [PATCH] notebook --- .../scale_factors/top_reweighting.ipynb | 544 ++---------------- 1 file changed, 41 insertions(+), 503 deletions(-) diff --git a/src/HHbbVV/scale_factors/top_reweighting.ipynb b/src/HHbbVV/scale_factors/top_reweighting.ipynb index 94fd1a04..8de28262 100644 --- a/src/HHbbVV/scale_factors/top_reweighting.ipynb +++ b/src/HHbbVV/scale_factors/top_reweighting.ipynb @@ -1045,20 +1045,9 @@ }, { "cell_type": "code", - "execution_count": 59, - "metadata": {}, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "2\n", - "2\n", - "2\n", - "2\n" - ] - } - ], + "execution_count": null, + "metadata": {}, + "outputs": [], "source": [ "pt_extrap_lookups_dict = {\"params\": [], \"errs\": [], \"sys_up_params\": [], \"sys_down_params\": []}\n", "max_fparams = 4 # max order (+1) of functions\n", @@ -1128,7 +1117,7 @@ }, { "cell_type": "code", - "execution_count": 51, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -1156,7 +1145,7 @@ }, { "cell_type": "code", - "execution_count": 53, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -1188,7 +1177,7 @@ }, { "cell_type": "code", - "execution_count": 92, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -1198,24 +1187,36 @@ }, { "cell_type": "code", - "execution_count": 93, + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "type(lds.layout) is ak._ext.ListArray64" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "lp_hist" + ] + }, + { + "cell_type": "code", + "execution_count": null, "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "Hist(\n", - " Variable([10, 65, 110, 175, 240, 300, 99999], name='subjet_pt', label='Subjet pT [GeV]'),\n", - " Variable([-1, -0.55, -0.1, 0.35, 0.8, 1.25, 1.7, 2.15, 2.6, 3.05, 3.5, 3.95, 4.4, 4.85, 5.3, 5.75, 6.2, 6.65, 7.1, 7.55, 8], name='logD', label='ln(0.8/Delta)'),\n", - " Variable([-5, -4.40461, -3.80922, -3.21384, -2.61845, -2.02306, -1.42767, -0.832286, -0.236898, 0.35849, 0.953878, 1.54927, 2.14465, 2.74004, 3.33543, 3.93082, 4.5262, 5.12159, 5.71698, 6.31237, 6.90776], name='logkt', label='ln(kT/GeV)'),\n", - " storage=Double()) # Sum: 435622.09045410156 (435925.44845581055 with flow)" - ] - }, - "execution_count": 93, - "metadata": {}, - "output_type": "execute_result" - } - ], + "outputs": [], + "source": [ + "sum([lp_hist, []])" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], "source": [ "lp_hist.fill(\n", " subjet_pt=flat_subjet_pt,\n", @@ -1227,458 +1228,9 @@ }, { "cell_type": "code", - "execution_count": 97, - "metadata": {}, - "outputs": [ - { - "data": { - "text/html": [ - "\n", - "
\n", - "
\n", - "\n", - "\n", - "-1\n", - "\n", - "\n", - "8\n", - "\n", - "\n", - "-5\n", - "\n", - "\n", - "6.91\n", - "\n", - "\n", - "ln(0.8/Delta)\n", - "\n", - "\n", - "ln(kT/GeV)\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "
\n", - "
\n", - "Variable([-1, -0.55, -0.1, 0.35, 0.8, 1.25, 1.7, 2.15, 2.6, 3.05, 3.5, 3.95, 4.4, 4.85, 5.3, 5.75, 6.2, 6.65, 7.1, 7.55, 8], name='logD', label='ln(0.8/Delta)')
\n", - "Variable([-5, -4.40461, -3.80922, -3.21384, -2.61845, -2.02306, -1.42767, -0.832286, -0.236898, 0.35849, 0.953878, 1.54927, 2.14465, 2.74004, 3.33543, 3.93082, 4.5262, 5.12159, 5.71698, 6.31237, 6.90776], name='logkt', label='ln(kT/GeV)')
\n", - "
\n", - "Double() Σ=39436.54022216797\n", - "\n", - "
\n", - "
\n", - "" - ], - "text/plain": [ - "Hist(\n", - " Variable([-1, -0.55, -0.1, 0.35, 0.8, 1.25, 1.7, 2.15, 2.6, 3.05, 3.5, 3.95, 4.4, 4.85, 5.3, 5.75, 6.2, 6.65, 7.1, 7.55, 8], name='logD', label='ln(0.8/Delta)'),\n", - " Variable([-5, -4.40461, -3.80922, -3.21384, -2.61845, -2.02306, -1.42767, -0.832286, -0.236898, 0.35849, 0.953878, 1.54927, 2.14465, 2.74004, 3.33543, 3.93082, 4.5262, 5.12159, 5.71698, 6.31237, 6.90776], name='logkt', label='ln(kT/GeV)'),\n", - " storage=Double()) # Sum: 39436.54022216797" - ] - }, - "execution_count": 97, - "metadata": {}, - "output_type": "execute_result" - } - ], + "execution_count": null, + "metadata": {}, + "outputs": [], "source": [ "lp_hist[0, ...]" ] @@ -1701,21 +1253,7 @@ "cell_type": "code", "execution_count": null, "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "Hist(\n", - " Variable([10, 65, 110, 175, 240, 300, 99999], name='subjet_pt', label='Subjet pT [GeV]'),\n", - " Variable([-1, -0.55, -0.1, 0.35, 0.8, 1.25, 1.7, 2.15, 2.6, 3.05, 3.5, 3.95, 4.4, 4.85, 5.3, 5.75, 6.2, 6.65, 7.1, 7.55, 8], name='logD', label='ln(0.8/Delta)'),\n", - " Variable([-5, -4.40461, -3.80922, -3.21384, -2.61845, -2.02306, -1.42767, -0.832286, -0.236898, 0.35849, 0.953878, 1.54927, 2.14465, 2.74004, 3.33543, 3.93082, 4.5262, 5.12159, 5.71698, 6.31237, 6.90776], name='logkt', label='ln(kT/GeV)'),\n", - " storage=Double()) # Sum: 438022.09045410156 (438628.80645751953 with flow)" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "sum([lp_hist, lp_hist2])" ] @@ -1729,7 +1267,7 @@ }, { "cell_type": "code", - "execution_count": 61, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -1745,7 +1283,7 @@ }, { "cell_type": "code", - "execution_count": 62, + "execution_count": null, "metadata": {}, "outputs": [], "source": [