{ "cells": [ { "cell_type": "markdown", "id": "3bc82063-09d6-48b1-93d0-4c7263c7b9da", "metadata": { "tags": [] }, "source": [ "# nll-fast" ] }, { "cell_type": "code", "execution_count": 1, "id": "b3184bee-8948-4b61-94b5-1a41deccb0f7", "metadata": { "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.2.10.6\n", "/home/apn/git/nll-fast/nll-fast\n" ] } ], "source": [ "import hepi\n", "print(hepi.__version__)\n", "import smpl\n", "import numpy as np\n", "import hepi.util as util\n", "import matplotlib.pyplot as plt\n", "from hepi.run import nllfast as nll\n", "# set the folder where the resummino binary can be found either in either ./{,bin,bin/build}/resummino\n", "nll.set_path(\"nll-fast-3.1\")\n", "# By default hepi will run with nice -n5 to prevent overloading the system if more scans than cores are running\n", "#rs.set_pre(\"\") disables any prefixing with nice\n", "print (nll.get_path())" ] }, { "cell_type": "code", "execution_count": 4, "id": "a922de4c-d019-4b67-bc73-5b062987efce", "metadata": { "tags": [] }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/apn/git/HEPi/hepi/util.py:101: UserWarning: LHAPDF python binding not installed? Make sure you set PYTHONPATH correctly (i.e. correct python version) if you want to compute PDF uncertainties.\n", " warnings.warn(\n" ] }, { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "f523ad0a59c04f4c8fb87155391f3960", "version_major": 2, "version_minor": 0 }, "text/plain": [ "Checking input: 0%| | 0/16 [00:00" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%%time\n", "params = [\n", " \"mastercode_with_gm2.in\", # baseline slha file in the relative ./output folder by default unless set_output_dir was used \n", "]\n", "pss = [ \n", " (1000021,-1000002), # Final state particles for resummino to run\n", " ]\n", "\n", "for pa,pb in pss:\n", " for param in params:\n", " # All the inputs Order, CMS in GeV, particle 1, particle 2, slha, pdf_lo, pdf_nlo,mu_f, mu_r\n", " i = hepi.Input(hepi.Order.NLO_PLUS_NLL,13000,pa,pb,param,\"cteq6l1\",\"cteq66\",1., 1.,id=4)\n", " li = [i] # li is our list of inputs that we want resummino to run\n", " li = hepi.mass_scan([i],pa, np.linspace(1000,2000,16)) # we scan the slepton mass from 100 to 1000 at 15 equidistant points \n", " rs_dl = nll.run(li,skip=False,n_jobs=1) # run resummino, skipping if the result already exists.\n", " # rs_dl is a panda table(dataframe) with all inputs and result \n", " _,axs = hepi.init_double_plot()\n", " # hepi has some useful plotting routines, but the results from rs_dl can easily be accessed\n", " # Now plot the mass of PDG id pa from the results at LO ,NLO and aNNLO+NLL\n", " hepi.mass_plot(rs_dl,[\"LO\",\"NLO\",\"NLO_PLUS_NLL\"],pa,axes=axs[0],tight=False)\n", " # Plot K factors vs LO and aNNLO+NLL/NLO\n", " hepi.mass_plot(rs_dl,[\"K_LO\",\"K_NLO\",\"K_NLO_PLUS_NLL\",\"NLO_PLUS_NLL_OVER_NLO\"],pa, yaxis=None, axes=axs[1],logy=False,tight=False)\n", " # construct a title from the inputs\n", " hepi.title(li[0],axs[0],scenario=\"mastercode\")" ] }, { "cell_type": "code", "execution_count": 5, "id": "595f203c-4255-4b39-8c0c-e36815aa2b37", "metadata": { "tags": [] }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
LONLONLO_PLUS_NLLaNNLO_PLUS_NNLLK_LOK_NLOK_NLO_PLUS_NLLNLO_PLUS_NLL_OVER_NLOK_aNNLO_PLUS_NNLLaNNLO_PLUS_NNLL_OVER_NLO...idmodelmumass_1000021runnernf_final_state_innf_pdfnf_squark_massnf_gluino_massnf_deg
0(2.8099999999999996+/-0)e-05(7.71+/-0)e-05(8.4+/-0)e-05None1.0+/-02.743772241992883+/-02.9893238434163707+/-01.0894941634241246+/-0NoneNone...42535.3290501000.000000NLLfastRunner-?sgcteq4070.65811000.00000010
1(2.3+/-0)e-05(6.390000000000001+/-0)e-05(6.9799999999999995+/-0)e-05None1.0+/-02.7782608695652176+/-03.0347826086956524+/-01.0923317683881064+/-0NoneNone...42568.6623831066.666667NLLfastRunner-?sgcteq4070.65811066.66666710
2(1.88+/-0)e-05(5.289999999999999+/-0)e-05(5.8+/-0)e-05None1.0+/-02.8138297872340425+/-03.0851063829787235+/-01.0964083175803403+/-0NoneNone...42601.9957171133.333333NLLfastRunner-?sgcteq4070.65811133.33333310
3(1.54+/-0)e-05(4.39+/-0)e-05(4.819999999999999+/-0)e-05None1.0+/-02.85064935064935+/-03.1298701298701297+/-01.0979498861047836+/-0NoneNone...42635.3290501200.000000NLLfastRunner-?sgcteq4070.65811200.00000010
4(1.27+/-0)e-05(3.6599999999999997+/-0)e-05(4.02+/-0)e-05None1.0+/-02.8818897637795278+/-03.1653543307086616+/-01.0983606557377048+/-0NoneNone...42668.6623831266.666667NLLfastRunner-?sgcteq4070.65811266.66666710
5(1.0399999999999998+/-0)e-05(3.0399999999999996+/-0)e-05(3.3499999999999996+/-0)e-05None1.0+/-02.9230769230769234+/-03.2211538461538467+/-01.1019736842105263+/-0NoneNone...42701.9957171333.333333NLLfastRunner-?sgcteq4070.65811333.33333310
6(8.54+/-0)e-06(2.53+/-0)e-05(2.7999999999999994+/-0)e-05None1.0+/-02.9625292740046842+/-03.2786885245901636+/-01.1067193675889326+/-0NoneNone...42735.3290501400.000000NLLfastRunner-?sgcteq4070.65811400.00000010
7(7.040000000000001+/-0)e-06(2.11+/-0)e-05(2.34+/-0)e-05None1.0+/-02.997159090909091+/-03.3238636363636362+/-01.1090047393364928+/-0NoneNone...42768.6623831466.666667NLLfastRunner-?sgcteq4070.65811466.66666710
8(5.79+/-0)e-06(1.76+/-0)e-05(1.9499999999999997+/-0)e-05None1.0+/-03.03972366148532+/-03.3678756476683938+/-01.1079545454545454+/-0NoneNone...42801.9957171533.333333NLLfastRunner-?sgcteq4070.65811533.33333310
9(4.760000000000001+/-0)e-06(1.4599999999999997+/-0)e-05(1.63+/-0)e-05None1.0+/-03.0672268907563023+/-03.4243697478991595+/-01.1164383561643836+/-0NoneNone...42835.3290501600.000000NLLfastRunner-?sgcteq4070.65811600.00000010
10(3.9299999999999997+/-0)e-06(1.23+/-0)e-05(1.3699999999999999+/-0)e-05None1.0+/-03.1297709923664128+/-03.486005089058524+/-01.113821138211382+/-0NoneNone...42868.6623831666.666667NLLfastRunner-?sgcteq4070.65811666.66666710
11(3.2400000000000007+/-0)e-06(1.02+/-0)e-05(1.15+/-0)e-05None1.0+/-03.148148148148148+/-03.5493827160493825+/-01.1274509803921569+/-0NoneNone...42901.9957171733.333333NLLfastRunner-?sgcteq4070.65811733.33333310
12(2.67+/-0)e-06(8.56+/-0)e-06(9.629999999999999+/-0)e-06None1.0+/-03.2059925093632957+/-03.606741573033708+/-01.125+/-0NoneNone...42935.3290501800.000000NLLfastRunner-?sgcteq4070.65811800.00000010
13(2.21+/-0)e-06(7.160000000000001+/-0)e-06(8.090000000000002+/-0)e-06None1.0+/-03.2398190045248874+/-03.6606334841628962+/-01.1298882681564244+/-0NoneNone...42968.6623831866.666667NLLfastRunner-?sgcteq4070.65811866.66666710
14(1.82+/-0)e-06(5.9799999999999995+/-0)e-06(6.779999999999999+/-0)e-06None1.0+/-03.2857142857142856+/-03.725274725274725+/-01.1337792642140467+/-0NoneNone...43001.9957171933.333333NLLfastRunner-?sgcteq4070.65811933.33333310
15(1.5+/-0)e-06(5.000000000000001+/-0)e-06(5.68+/-0)e-06None1.0+/-03.3333333333333335+/-03.7866666666666666+/-01.136+/-0NoneNone...43035.3290502000.000000NLLfastRunner-?sgcteq4070.65812000.00000010
\n", "

16 rows × 43 columns

\n", "
" ], "text/plain": [ " LO NLO \n", "0 (2.8099999999999996+/-0)e-05 (7.71+/-0)e-05 \\\n", "1 (2.3+/-0)e-05 (6.390000000000001+/-0)e-05 \n", "2 (1.88+/-0)e-05 (5.289999999999999+/-0)e-05 \n", "3 (1.54+/-0)e-05 (4.39+/-0)e-05 \n", "4 (1.27+/-0)e-05 (3.6599999999999997+/-0)e-05 \n", "5 (1.0399999999999998+/-0)e-05 (3.0399999999999996+/-0)e-05 \n", "6 (8.54+/-0)e-06 (2.53+/-0)e-05 \n", "7 (7.040000000000001+/-0)e-06 (2.11+/-0)e-05 \n", "8 (5.79+/-0)e-06 (1.76+/-0)e-05 \n", "9 (4.760000000000001+/-0)e-06 (1.4599999999999997+/-0)e-05 \n", "10 (3.9299999999999997+/-0)e-06 (1.23+/-0)e-05 \n", "11 (3.2400000000000007+/-0)e-06 (1.02+/-0)e-05 \n", "12 (2.67+/-0)e-06 (8.56+/-0)e-06 \n", "13 (2.21+/-0)e-06 (7.160000000000001+/-0)e-06 \n", "14 (1.82+/-0)e-06 (5.9799999999999995+/-0)e-06 \n", "15 (1.5+/-0)e-06 (5.000000000000001+/-0)e-06 \n", "\n", " NLO_PLUS_NLL aNNLO_PLUS_NNLL K_LO \n", "0 (8.4+/-0)e-05 None 1.0+/-0 \\\n", "1 (6.9799999999999995+/-0)e-05 None 1.0+/-0 \n", "2 (5.8+/-0)e-05 None 1.0+/-0 \n", "3 (4.819999999999999+/-0)e-05 None 1.0+/-0 \n", "4 (4.02+/-0)e-05 None 1.0+/-0 \n", "5 (3.3499999999999996+/-0)e-05 None 1.0+/-0 \n", "6 (2.7999999999999994+/-0)e-05 None 1.0+/-0 \n", "7 (2.34+/-0)e-05 None 1.0+/-0 \n", "8 (1.9499999999999997+/-0)e-05 None 1.0+/-0 \n", "9 (1.63+/-0)e-05 None 1.0+/-0 \n", "10 (1.3699999999999999+/-0)e-05 None 1.0+/-0 \n", "11 (1.15+/-0)e-05 None 1.0+/-0 \n", "12 (9.629999999999999+/-0)e-06 None 1.0+/-0 \n", "13 (8.090000000000002+/-0)e-06 None 1.0+/-0 \n", "14 (6.779999999999999+/-0)e-06 None 1.0+/-0 \n", "15 (5.68+/-0)e-06 None 1.0+/-0 \n", "\n", " K_NLO K_NLO_PLUS_NLL NLO_PLUS_NLL_OVER_NLO \n", "0 2.743772241992883+/-0 2.9893238434163707+/-0 1.0894941634241246+/-0 \\\n", "1 2.7782608695652176+/-0 3.0347826086956524+/-0 1.0923317683881064+/-0 \n", "2 2.8138297872340425+/-0 3.0851063829787235+/-0 1.0964083175803403+/-0 \n", "3 2.85064935064935+/-0 3.1298701298701297+/-0 1.0979498861047836+/-0 \n", "4 2.8818897637795278+/-0 3.1653543307086616+/-0 1.0983606557377048+/-0 \n", "5 2.9230769230769234+/-0 3.2211538461538467+/-0 1.1019736842105263+/-0 \n", "6 2.9625292740046842+/-0 3.2786885245901636+/-0 1.1067193675889326+/-0 \n", "7 2.997159090909091+/-0 3.3238636363636362+/-0 1.1090047393364928+/-0 \n", "8 3.03972366148532+/-0 3.3678756476683938+/-0 1.1079545454545454+/-0 \n", "9 3.0672268907563023+/-0 3.4243697478991595+/-0 1.1164383561643836+/-0 \n", "10 3.1297709923664128+/-0 3.486005089058524+/-0 1.113821138211382+/-0 \n", "11 3.148148148148148+/-0 3.5493827160493825+/-0 1.1274509803921569+/-0 \n", "12 3.2059925093632957+/-0 3.606741573033708+/-0 1.125+/-0 \n", "13 3.2398190045248874+/-0 3.6606334841628962+/-0 1.1298882681564244+/-0 \n", "14 3.2857142857142856+/-0 3.725274725274725+/-0 1.1337792642140467+/-0 \n", "15 3.3333333333333335+/-0 3.7866666666666666+/-0 1.136+/-0 \n", "\n", " K_aNNLO_PLUS_NNLL aNNLO_PLUS_NNLL_OVER_NLO ... id model mu \n", "0 None None ... 4 2535.329050 \\\n", "1 None None ... 4 2568.662383 \n", "2 None None ... 4 2601.995717 \n", "3 None None ... 4 2635.329050 \n", "4 None None ... 4 2668.662383 \n", "5 None None ... 4 2701.995717 \n", "6 None None ... 4 2735.329050 \n", "7 None None ... 4 2768.662383 \n", "8 None None ... 4 2801.995717 \n", "9 None None ... 4 2835.329050 \n", "10 None None ... 4 2868.662383 \n", "11 None None ... 4 2901.995717 \n", "12 None None ... 4 2935.329050 \n", "13 None None ... 4 2968.662383 \n", "14 None None ... 4 3001.995717 \n", "15 None None ... 4 3035.329050 \n", "\n", " mass_1000021 runner nf_final_state_in nf_pdf nf_squark_mass \n", "0 1000.000000 NLLfastRunner-? sg cteq 4070.6581 \\\n", "1 1066.666667 NLLfastRunner-? sg cteq 4070.6581 \n", "2 1133.333333 NLLfastRunner-? sg cteq 4070.6581 \n", "3 1200.000000 NLLfastRunner-? sg cteq 4070.6581 \n", "4 1266.666667 NLLfastRunner-? sg cteq 4070.6581 \n", "5 1333.333333 NLLfastRunner-? sg cteq 4070.6581 \n", "6 1400.000000 NLLfastRunner-? sg cteq 4070.6581 \n", "7 1466.666667 NLLfastRunner-? sg cteq 4070.6581 \n", "8 1533.333333 NLLfastRunner-? sg cteq 4070.6581 \n", "9 1600.000000 NLLfastRunner-? sg cteq 4070.6581 \n", "10 1666.666667 NLLfastRunner-? sg cteq 4070.6581 \n", "11 1733.333333 NLLfastRunner-? sg cteq 4070.6581 \n", "12 1800.000000 NLLfastRunner-? sg cteq 4070.6581 \n", "13 1866.666667 NLLfastRunner-? sg cteq 4070.6581 \n", "14 1933.333333 NLLfastRunner-? sg cteq 4070.6581 \n", "15 2000.000000 NLLfastRunner-? sg cteq 4070.6581 \n", "\n", " nf_gluino_mass nf_deg \n", "0 1000.000000 10 \n", "1 1066.666667 10 \n", "2 1133.333333 10 \n", "3 1200.000000 10 \n", "4 1266.666667 10 \n", "5 1333.333333 10 \n", "6 1400.000000 10 \n", "7 1466.666667 10 \n", "8 1533.333333 10 \n", "9 1600.000000 10 \n", "10 1666.666667 10 \n", "11 1733.333333 10 \n", "12 1800.000000 10 \n", "13 1866.666667 10 \n", "14 1933.333333 10 \n", "15 2000.000000 10 \n", "\n", "[16 rows x 43 columns]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "rs_dl" ] }, { "cell_type": "code", "execution_count": null, "id": "fff65bb6-66ce-405b-8345-f4e1627d065f", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.3" } }, "nbformat": 4, "nbformat_minor": 5 }