diff --git a/codigo/.ipynb_checkpoints/Data_Analysis-checkpoint.ipynb b/codigo/.ipynb_checkpoints/Data_Analysis-checkpoint.ipynb new file mode 100644 index 0000000000000000000000000000000000000000..764dd86cb2028cbf5c6049cedf7596a46822cf40 --- /dev/null +++ b/codigo/.ipynb_checkpoints/Data_Analysis-checkpoint.ipynb @@ -0,0 +1,1785 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Z -> ll (decaimiento del boson Z a dos leptones) \n", + "> Proveniente de https://atlas-opendata.web.cern.ch/atlas-opendata/samples/2020/2lep/MC/mc_361106.Zee.1largeRjet1lep.root" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Importamos los paquetes necesarios" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Welcome to JupyROOT 6.18/04\n" + ] + } + ], + "source": [ + "import ROOT\n", + "from ROOT import TMath\n", + "import time" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Usamos pandas para leer los archivos .root" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [], + "source": [ + "import pandas as pd\n", + "from root_pandas import read_root" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [], + "source": [ + "%jsroot on" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Z boson invariant mass (reconstruction)" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": {}, + "outputs": [], + "source": [ + "file2 = ROOT.TFile.Open(\"/home/student/ejercicios-clase-08-datos/data-used/mc_361106.Zee.1largeRjet1lep.root\")\n", + "# file2 = ROOT.TFile.Open(\"/https://atlas-opendata.web.cern.ch/atlas-opendata/samples/2020/2lep/MC/mc_361106.Zee.1largeRjet1lep.root\")\n", + "# desde la web" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Visualizamos el archivo aprovechando que no es tan grande, podemos visualizarlo como si fuera un archivo csv" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [], + "source": [ + "df = read_root('/home/student/ejercicios-clase-08-datos/data-used/mc_361106.Zee.1largeRjet1lep.root')" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "scrolled": true + }, + "outputs": [ + { + "data": { + "text/html": [ + "<div>\n", + "<style scoped>\n", + " .dataframe tbody tr th:only-of-type {\n", + " vertical-align: middle;\n", + " }\n", + "\n", + " .dataframe tbody tr th {\n", + " vertical-align: top;\n", + " }\n", + "\n", + " .dataframe thead th {\n", + " text-align: right;\n", + " }\n", + "</style>\n", + "<table border=\"1\" class=\"dataframe\">\n", + " <thead>\n", + " <tr style=\"text-align: right;\">\n", + " <th></th>\n", + " <th>runNumber</th>\n", + " <th>eventNumber</th>\n", + " <th>channelNumber</th>\n", + " <th>mcWeight</th>\n", + " <th>scaleFactor_PILEUP</th>\n", + " <th>scaleFactor_ELE</th>\n", + " <th>scaleFactor_MUON</th>\n", + " <th>scaleFactor_PHOTON</th>\n", + " <th>scaleFactor_TAU</th>\n", + " <th>scaleFactor_BTAG</th>\n", + " <th>...</th>\n", + " <th>largeRjet_pt</th>\n", + " <th>largeRjet_eta</th>\n", + " <th>largeRjet_phi</th>\n", + " <th>largeRjet_E</th>\n", + " <th>largeRjet_m</th>\n", + " <th>largeRjet_truthMatched</th>\n", + " <th>largeRjet_D2</th>\n", + " <th>largeRjet_tau32</th>\n", + " <th>largeRjet_pt_syst</th>\n", + " <th>tau_charge</th>\n", + " </tr>\n", + " </thead>\n", + " <tbody>\n", + " <tr>\n", + " <th>0</th>\n", + " <td>284500</td>\n", + " <td>15440</td>\n", + " <td>361106</td>\n", + " <td>1940.819946</td>\n", + " <td>1.264594</td>\n", + " <td>0.941115</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999384</td>\n", + " <td>...</td>\n", + " <td>[336389.22, 326111.62]</td>\n", + " <td>[-1.6073604, 0.34273493]</td>\n", + " <td>[-3.1196432, -0.05282484]</td>\n", + " <td>[878163.0, 352508.03]</td>\n", + " <td>[95666.125, 70169.38]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[0.06860469, 1.7374916]</td>\n", + " <td>[0.5651315, 0.6696122]</td>\n", + " <td>[17716.031, 4227.1562]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>1</th>\n", + " <td>284500</td>\n", + " <td>81085</td>\n", + " <td>361106</td>\n", + " <td>1941.910034</td>\n", + " <td>1.320948</td>\n", + " <td>0.992704</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.997777</td>\n", + " <td>...</td>\n", + " <td>[424948.16, 295925.38, 288773.56]</td>\n", + " <td>[-0.50103956, 0.00353557, 0.7989666]</td>\n", + " <td>[0.39738846, 3.0375211, -2.1750822]</td>\n", + " <td>[494285.72, 314167.06, 393288.84]</td>\n", + " <td>[120339.15, 105489.414, 75616.88]</td>\n", + " <td>[0.0, 0.0, 0.0]</td>\n", + " <td>[1.4906838, 0.1404522, 0.62883687]</td>\n", + " <td>[0.64305073, 0.19809924, 0.69501674]</td>\n", + " <td>[22269.969, 14174.656, 24763.469]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>2</th>\n", + " <td>284500</td>\n", + " <td>89512</td>\n", + " <td>361106</td>\n", + " <td>1940.760010</td>\n", + " <td>1.251332</td>\n", + " <td>0.993904</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999102</td>\n", + " <td>...</td>\n", + " <td>[255974.6]</td>\n", + " <td>[1.6731088]</td>\n", + " <td>[0.80495465]</td>\n", + " <td>[706647.8]</td>\n", + " <td>[29652.18]</td>\n", + " <td>[0.0]</td>\n", + " <td>[3.0340846]</td>\n", + " <td>[0.8060808]</td>\n", + " <td>[0.0]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>3</th>\n", + " <td>284500</td>\n", + " <td>92009</td>\n", + " <td>361106</td>\n", + " <td>1941.709961</td>\n", + " <td>1.320948</td>\n", + " <td>0.957973</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.000000</td>\n", + " <td>...</td>\n", + " <td>[357909.62]</td>\n", + " <td>[-0.23516123]</td>\n", + " <td>[-2.2556815]</td>\n", + " <td>[377280.16]</td>\n", + " <td>[83818.125]</td>\n", + " <td>[1.0]</td>\n", + " <td>[0.13190588]</td>\n", + " <td>[0.18124953]</td>\n", + " <td>[1438.0938]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>4</th>\n", + " <td>284500</td>\n", + " <td>102088</td>\n", + " <td>361106</td>\n", + " <td>1941.239990</td>\n", + " <td>1.320948</td>\n", + " <td>0.948454</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.997393</td>\n", + " <td>...</td>\n", + " <td>[339715.34, 337192.75]</td>\n", + " <td>[0.22991252, -1.1238769]</td>\n", + " <td>[-1.2795317, 1.6654552]</td>\n", + " <td>[362021.66, 574320.1]</td>\n", + " <td>[97183.07, 30168.135]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[0.042992596, 2.0432703]</td>\n", + " <td>[0.075145245, 0.71639764]</td>\n", + " <td>[13062.9375, 5456.125]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>...</th>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " </tr>\n", + " <tr>\n", + " <th>53648</th>\n", + " <td>284500</td>\n", + " <td>5740044</td>\n", + " <td>361106</td>\n", + " <td>1941.760010</td>\n", + " <td>0.934886</td>\n", + " <td>0.986970</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999438</td>\n", + " <td>...</td>\n", + " <td>[272509.03, 264178.5]</td>\n", + " <td>[-0.11495946, -0.88012326]</td>\n", + " <td>[-2.6213317, 0.5289878]</td>\n", + " <td>[274337.56, 373833.25]</td>\n", + " <td>[3766.6348, 20426.285]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[0.0, 2.0461884]</td>\n", + " <td>[-1.0, 0.7430232]</td>\n", + " <td>[17955.438, 0.0]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>53649</th>\n", + " <td>284500</td>\n", + " <td>5752321</td>\n", + " <td>361106</td>\n", + " <td>1941.680054</td>\n", + " <td>0.596562</td>\n", + " <td>0.889910</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999002</td>\n", + " <td>...</td>\n", + " <td>[272032.44]</td>\n", + " <td>[1.6677606]</td>\n", + " <td>[-2.4962265]</td>\n", + " <td>[746732.94]</td>\n", + " <td>[14772.887]</td>\n", + " <td>[0.0]</td>\n", + " <td>[3.902746]</td>\n", + " <td>[0.58085346]</td>\n", + " <td>[21161.61]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>53650</th>\n", + " <td>284500</td>\n", + " <td>5757742</td>\n", + " <td>361106</td>\n", + " <td>1941.000000</td>\n", + " <td>1.251332</td>\n", + " <td>0.984580</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.998199</td>\n", + " <td>...</td>\n", + " <td>[262801.75]</td>\n", + " <td>[-0.25520322]</td>\n", + " <td>[-1.0126089]</td>\n", + " <td>[280588.1]</td>\n", + " <td>[71192.08]</td>\n", + " <td>[0.0]</td>\n", + " <td>[0.036484536]</td>\n", + " <td>[0.39032233]</td>\n", + " <td>[7288.8906]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>53651</th>\n", + " <td>284500</td>\n", + " <td>5745755</td>\n", + " <td>361106</td>\n", + " <td>1940.719971</td>\n", + " <td>1.196105</td>\n", + " <td>0.964722</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999469</td>\n", + " <td>...</td>\n", + " <td>[266810.53, 260362.53]</td>\n", + " <td>[-0.47211164, 0.4879255]</td>\n", + " <td>[-1.3692597, 1.7747238]</td>\n", + " <td>[298000.4, 312144.3]</td>\n", + " <td>[23127.982, 110384.94]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[1.6352673, 0.0088554425]</td>\n", + " <td>[0.65562195, -1.0]</td>\n", + " <td>[1491.4062, 8128.9688]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>53652</th>\n", + " <td>284500</td>\n", + " <td>5760130</td>\n", + " <td>361106</td>\n", + " <td>1940.819946</td>\n", + " <td>0.506784</td>\n", + " <td>1.006096</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.000000</td>\n", + " <td>...</td>\n", + " <td>[326995.22, 289713.94]</td>\n", + " <td>[0.7875121, 0.041406196]</td>\n", + " <td>[-1.3435338, 1.6638559]</td>\n", + " <td>[433746.4, 289980.38]</td>\n", + " <td>[1893.2017, 3234.1577]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[77.1253, 34.656742]</td>\n", + " <td>[-1.0, -1.0]</td>\n", + " <td>[788.78125, 122.65625]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " </tbody>\n", + "</table>\n", + "<p>53653 rows × 81 columns</p>\n", + "</div>" + ], + "text/plain": [ + " runNumber eventNumber channelNumber mcWeight scaleFactor_PILEUP \\\n", + "0 284500 15440 361106 1940.819946 1.264594 \n", + "1 284500 81085 361106 1941.910034 1.320948 \n", + "2 284500 89512 361106 1940.760010 1.251332 \n", + "3 284500 92009 361106 1941.709961 1.320948 \n", + "4 284500 102088 361106 1941.239990 1.320948 \n", + "... ... ... ... ... ... \n", + "53648 284500 5740044 361106 1941.760010 0.934886 \n", + "53649 284500 5752321 361106 1941.680054 0.596562 \n", + "53650 284500 5757742 361106 1941.000000 1.251332 \n", + "53651 284500 5745755 361106 1940.719971 1.196105 \n", + "53652 284500 5760130 361106 1940.819946 0.506784 \n", + "\n", + " scaleFactor_ELE scaleFactor_MUON scaleFactor_PHOTON scaleFactor_TAU \\\n", + "0 0.941115 1.0 1.0 1.0 \n", + "1 0.992704 1.0 1.0 1.0 \n", + "2 0.993904 1.0 1.0 1.0 \n", + "3 0.957973 1.0 1.0 1.0 \n", + "4 0.948454 1.0 1.0 1.0 \n", + "... ... ... ... ... \n", + "53648 0.986970 1.0 1.0 1.0 \n", + "53649 0.889910 1.0 1.0 1.0 \n", + "53650 0.984580 1.0 1.0 1.0 \n", + "53651 0.964722 1.0 1.0 1.0 \n", + "53652 1.006096 1.0 1.0 1.0 \n", + "\n", + " scaleFactor_BTAG ... largeRjet_pt \\\n", + "0 0.999384 ... [336389.22, 326111.62] \n", + "1 0.997777 ... [424948.16, 295925.38, 288773.56] \n", + "2 0.999102 ... [255974.6] \n", + "3 1.000000 ... [357909.62] \n", + "4 0.997393 ... [339715.34, 337192.75] \n", + "... ... ... ... \n", + "53648 0.999438 ... [272509.03, 264178.5] \n", + "53649 0.999002 ... [272032.44] \n", + "53650 0.998199 ... [262801.75] \n", + "53651 0.999469 ... [266810.53, 260362.53] \n", + "53652 1.000000 ... [326995.22, 289713.94] \n", + "\n", + " largeRjet_eta \\\n", + "0 [-1.6073604, 0.34273493] \n", + "1 [-0.50103956, 0.00353557, 0.7989666] \n", + "2 [1.6731088] \n", + "3 [-0.23516123] \n", + "4 [0.22991252, -1.1238769] \n", + "... ... \n", + "53648 [-0.11495946, -0.88012326] \n", + "53649 [1.6677606] \n", + "53650 [-0.25520322] \n", + "53651 [-0.47211164, 0.4879255] \n", + "53652 [0.7875121, 0.041406196] \n", + "\n", + " largeRjet_phi largeRjet_E \\\n", + "0 [-3.1196432, -0.05282484] [878163.0, 352508.03] \n", + "1 [0.39738846, 3.0375211, -2.1750822] [494285.72, 314167.06, 393288.84] \n", + "2 [0.80495465] [706647.8] \n", + "3 [-2.2556815] [377280.16] \n", + "4 [-1.2795317, 1.6654552] [362021.66, 574320.1] \n", + "... ... ... \n", + "53648 [-2.6213317, 0.5289878] [274337.56, 373833.25] \n", + "53649 [-2.4962265] [746732.94] \n", + "53650 [-1.0126089] [280588.1] \n", + "53651 [-1.3692597, 1.7747238] [298000.4, 312144.3] \n", + "53652 [-1.3435338, 1.6638559] [433746.4, 289980.38] \n", + "\n", + " largeRjet_m largeRjet_truthMatched \\\n", + "0 [95666.125, 70169.38] [0.0, 0.0] \n", + "1 [120339.15, 105489.414, 75616.88] [0.0, 0.0, 0.0] \n", + "2 [29652.18] [0.0] \n", + "3 [83818.125] [1.0] \n", + "4 [97183.07, 30168.135] [0.0, 0.0] \n", + "... ... ... \n", + "53648 [3766.6348, 20426.285] [0.0, 0.0] \n", + "53649 [14772.887] [0.0] \n", + "53650 [71192.08] [0.0] \n", + "53651 [23127.982, 110384.94] [0.0, 0.0] \n", + "53652 [1893.2017, 3234.1577] [0.0, 0.0] \n", + "\n", + " largeRjet_D2 \\\n", + "0 [0.06860469, 1.7374916] \n", + "1 [1.4906838, 0.1404522, 0.62883687] \n", + "2 [3.0340846] \n", + "3 [0.13190588] \n", + "4 [0.042992596, 2.0432703] \n", + "... ... \n", + "53648 [0.0, 2.0461884] \n", + "53649 [3.902746] \n", + "53650 [0.036484536] \n", + "53651 [1.6352673, 0.0088554425] \n", + "53652 [77.1253, 34.656742] \n", + "\n", + " largeRjet_tau32 \\\n", + "0 [0.5651315, 0.6696122] \n", + "1 [0.64305073, 0.19809924, 0.69501674] \n", + "2 [0.8060808] \n", + "3 [0.18124953] \n", + "4 [0.075145245, 0.71639764] \n", + "... ... \n", + "53648 [-1.0, 0.7430232] \n", + "53649 [0.58085346] \n", + "53650 [0.39032233] \n", + "53651 [0.65562195, -1.0] \n", + "53652 [-1.0, -1.0] \n", + "\n", + " largeRjet_pt_syst tau_charge \n", + "0 [17716.031, 4227.1562] [] \n", + "1 [22269.969, 14174.656, 24763.469] [] \n", + "2 [0.0] [] \n", + "3 [1438.0938] [] \n", + "4 [13062.9375, 5456.125] [] \n", + "... ... ... \n", + "53648 [17955.438, 0.0] [] \n", + "53649 [21161.61] [] \n", + "53650 [7288.8906] [] \n", + "53651 [1491.4062, 8128.9688] [] \n", + "53652 [788.78125, 122.65625] [] \n", + "\n", + "[53653 rows x 81 columns]" + ] + }, + "execution_count": 6, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "df" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Vemos que tiene 81 columnas, para poder verlas todas con la finalidad de ver que variables podemos usar, usamos la siguiente lÃnea." + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": {}, + "outputs": [], + "source": [ + "pd.options.display.max_columns = 0" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": { + "scrolled": false + }, + "outputs": [ + { + "data": { + "text/html": [ + "<div>\n", + "<style scoped>\n", + " .dataframe tbody tr th:only-of-type {\n", + " vertical-align: middle;\n", + " }\n", + "\n", + " .dataframe tbody tr th {\n", + " vertical-align: top;\n", + " }\n", + "\n", + " .dataframe thead th {\n", + " text-align: right;\n", + " }\n", + "</style>\n", + "<table border=\"1\" class=\"dataframe\">\n", + " <thead>\n", + " <tr style=\"text-align: right;\">\n", + " <th></th>\n", + " <th>runNumber</th>\n", + " <th>eventNumber</th>\n", + " <th>channelNumber</th>\n", + " <th>mcWeight</th>\n", + " <th>scaleFactor_PILEUP</th>\n", + " <th>scaleFactor_ELE</th>\n", + " <th>scaleFactor_MUON</th>\n", + " <th>scaleFactor_PHOTON</th>\n", + " <th>scaleFactor_TAU</th>\n", + " <th>scaleFactor_BTAG</th>\n", + " <th>scaleFactor_LepTRIGGER</th>\n", + " <th>scaleFactor_PhotonTRIGGER</th>\n", + " <th>trigE</th>\n", + " <th>trigM</th>\n", + " <th>trigP</th>\n", + " <th>lep_n</th>\n", + " <th>lep_truthMatched</th>\n", + " <th>lep_trigMatched</th>\n", + " <th>lep_pt</th>\n", + " <th>lep_eta</th>\n", + " <th>lep_phi</th>\n", + " <th>lep_E</th>\n", + " <th>lep_z0</th>\n", + " <th>lep_charge</th>\n", + " <th>lep_type</th>\n", + " <th>lep_isTightID</th>\n", + " <th>lep_ptcone30</th>\n", + " <th>lep_etcone20</th>\n", + " <th>lep_trackd0pvunbiased</th>\n", + " <th>lep_tracksigd0pvunbiased</th>\n", + " <th>met_et</th>\n", + " <th>met_phi</th>\n", + " <th>jet_n</th>\n", + " <th>jet_pt</th>\n", + " <th>jet_eta</th>\n", + " <th>jet_phi</th>\n", + " <th>jet_E</th>\n", + " <th>jet_jvt</th>\n", + " <th>jet_trueflav</th>\n", + " <th>jet_truthMatched</th>\n", + " <th>...</th>\n", + " <th>photon_n</th>\n", + " <th>photon_truthMatched</th>\n", + " <th>photon_trigMatched</th>\n", + " <th>photon_pt</th>\n", + " <th>photon_eta</th>\n", + " <th>photon_phi</th>\n", + " <th>photon_E</th>\n", + " <th>photon_isTightID</th>\n", + " <th>photon_ptcone30</th>\n", + " <th>photon_etcone20</th>\n", + " <th>photon_convType</th>\n", + " <th>tau_n</th>\n", + " <th>tau_pt</th>\n", + " <th>tau_eta</th>\n", + " <th>tau_phi</th>\n", + " <th>tau_E</th>\n", + " <th>tau_isTightID</th>\n", + " <th>tau_truthMatched</th>\n", + " <th>tau_trigMatched</th>\n", + " <th>tau_nTracks</th>\n", + " <th>tau_BDTid</th>\n", + " <th>ditau_m</th>\n", + " <th>lep_pt_syst</th>\n", + " <th>met_et_syst</th>\n", + " <th>jet_pt_syst</th>\n", + " <th>photon_pt_syst</th>\n", + " <th>tau_pt_syst</th>\n", + " <th>XSection</th>\n", + " <th>SumWeights</th>\n", + " <th>largeRjet_n</th>\n", + " <th>largeRjet_pt</th>\n", + " <th>largeRjet_eta</th>\n", + " <th>largeRjet_phi</th>\n", + " <th>largeRjet_E</th>\n", + " <th>largeRjet_m</th>\n", + " <th>largeRjet_truthMatched</th>\n", + " <th>largeRjet_D2</th>\n", + " <th>largeRjet_tau32</th>\n", + " <th>largeRjet_pt_syst</th>\n", + " <th>tau_charge</th>\n", + " </tr>\n", + " </thead>\n", + " <tbody>\n", + " <tr>\n", + " <th>0</th>\n", + " <td>284500</td>\n", + " <td>15440</td>\n", + " <td>361106</td>\n", + " <td>1940.819946</td>\n", + " <td>1.264594</td>\n", + " <td>0.941115</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999384</td>\n", + " <td>1.01508</td>\n", + " <td>1.0</td>\n", + " <td>True</td>\n", + " <td>False</td>\n", + " <td>False</td>\n", + " <td>2</td>\n", + " <td>[True, True]</td>\n", + " <td>[True, True]</td>\n", + " <td>[169903.48, 142998.28]</td>\n", + " <td>[-1.8025961, -1.29893]</td>\n", + " <td>[3.0199232, -2.9465997]</td>\n", + " <td>[529270.06, 281577.66]</td>\n", + " <td>[-0.006664276, -0.0015563965]</td>\n", + " <td>[1, -1]</td>\n", + " <td>[11, 11]</td>\n", + " <td>[True, True]</td>\n", + " <td>[0.0, 5793.038]</td>\n", + " <td>[-664.72797, 1438.3691]</td>\n", + " <td>[0.028476434, 0.05088436]</td>\n", + " <td>[2.7847657, 1.8824033]</td>\n", + " <td>22774.525391</td>\n", + " <td>0.666735</td>\n", + " <td>1</td>\n", + " <td>[304260.75]</td>\n", + " <td>[0.2824949]</td>\n", + " <td>[-0.060711365]</td>\n", + " <td>[320657.12]</td>\n", + " <td>[0.9991386]</td>\n", + " <td>[0]</td>\n", + " <td>[False]</td>\n", + " <td>...</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0.0</td>\n", + " <td>[2525.5752, 1415.3539]</td>\n", + " <td>4430.583984</td>\n", + " <td>[7431.259]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>2</td>\n", + " <td>[336389.22, 326111.62]</td>\n", + " <td>[-1.6073604, 0.34273493]</td>\n", + " <td>[-3.1196432, -0.05282484]</td>\n", + " <td>[878163.0, 352508.03]</td>\n", + " <td>[95666.125, 70169.38]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[0.06860469, 1.7374916]</td>\n", + " <td>[0.5651315, 0.6696122]</td>\n", + " <td>[17716.031, 4227.1562]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>1</th>\n", + " <td>284500</td>\n", + " <td>81085</td>\n", + " <td>361106</td>\n", + " <td>1941.910034</td>\n", + " <td>1.320948</td>\n", + " <td>0.992704</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.997777</td>\n", + " <td>1.00936</td>\n", + " <td>1.0</td>\n", + " <td>True</td>\n", + " <td>False</td>\n", + " <td>False</td>\n", + " <td>2</td>\n", + " <td>[True, True]</td>\n", + " <td>[True, True]</td>\n", + " <td>[211820.23, 78281.89]</td>\n", + " <td>[0.047025315, -0.12980008]</td>\n", + " <td>[-3.0646539, 2.5374086]</td>\n", + " <td>[212054.48, 78942.266]</td>\n", + " <td>[-0.0404315, -0.0126543045]</td>\n", + " <td>[-1, 1]</td>\n", + " <td>[11, 11]</td>\n", + " <td>[True, True]</td>\n", + " <td>[0.0, 1394.1083]</td>\n", + " <td>[1965.3975, -560.0946]</td>\n", + " <td>[0.02091093, 0.0027772957]</td>\n", + " <td>[2.136357, 0.23608878]</td>\n", + " <td>12231.850586</td>\n", + " <td>0.864738</td>\n", + " <td>4</td>\n", + " <td>[357777.12, 162996.27, 122650.02, 87809.67]</td>\n", + " <td>[-0.38215926, 0.7312151, 0.827546, -1.0698929]</td>\n", + " <td>[0.41023782, -2.3782916, -1.9060545, 0.34156215]</td>\n", + " <td>[386681.56, 208748.97, 167330.88, 143282.44]</td>\n", + " <td>[0.991028, 0.9987799, 0.9975304, 0.9996967]</td>\n", + " <td>[0, 0, 0, 0]</td>\n", + " <td>[False, False, False, False]</td>\n", + " <td>...</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0.0</td>\n", + " <td>[2346.7507, 302.21893]</td>\n", + " <td>14149.702148</td>\n", + " <td>[5875.099, 6026.0947, 4836.232, 2209.1172]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>3</td>\n", + " <td>[424948.16, 295925.38, 288773.56]</td>\n", + " <td>[-0.50103956, 0.00353557, 0.7989666]</td>\n", + " <td>[0.39738846, 3.0375211, -2.1750822]</td>\n", + " <td>[494285.72, 314167.06, 393288.84]</td>\n", + " <td>[120339.15, 105489.414, 75616.88]</td>\n", + " <td>[0.0, 0.0, 0.0]</td>\n", + " <td>[1.4906838, 0.1404522, 0.62883687]</td>\n", + " <td>[0.64305073, 0.19809924, 0.69501674]</td>\n", + " <td>[22269.969, 14174.656, 24763.469]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>2</th>\n", + " <td>284500</td>\n", + " <td>89512</td>\n", + " <td>361106</td>\n", + " <td>1940.760010</td>\n", + " <td>1.251332</td>\n", + " <td>0.993904</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999102</td>\n", + " <td>1.00632</td>\n", + " <td>1.0</td>\n", + " <td>True</td>\n", + " <td>False</td>\n", + " <td>False</td>\n", + " <td>2</td>\n", + " <td>[True, True]</td>\n", + " <td>[True, False]</td>\n", + " <td>[234039.14, 8108.9414]</td>\n", + " <td>[0.5145523, -0.5389483]</td>\n", + " <td>[-2.36613, -0.6087628]</td>\n", + " <td>[265711.38, 9315.408]</td>\n", + " <td>[-0.072546005, 0.023675919]</td>\n", + " <td>[1, -1]</td>\n", + " <td>[11, 11]</td>\n", + " <td>[True, True]</td>\n", + " <td>[0.0, 1282.8975]</td>\n", + " <td>[-1899.3718, -219.85101]</td>\n", + " <td>[-0.01494353, 0.019966032]</td>\n", + " <td>[1.7670692, 0.88004786]</td>\n", + " <td>10936.896484</td>\n", + " <td>-2.576727</td>\n", + " <td>1</td>\n", + " <td>[239787.16]</td>\n", + " <td>[1.6001471]</td>\n", + " <td>[0.7848504]</td>\n", + " <td>[618527.5]</td>\n", + " <td>[0.9953245]</td>\n", + " <td>[0]</td>\n", + " <td>[False]</td>\n", + " <td>...</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0.0</td>\n", + " <td>[1448.3579, 86.43445]</td>\n", + " <td>7753.230469</td>\n", + " <td>[8643.487]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1</td>\n", + " <td>[255974.6]</td>\n", + " <td>[1.6731088]</td>\n", + " <td>[0.80495465]</td>\n", + " <td>[706647.8]</td>\n", + " <td>[29652.18]</td>\n", + " <td>[0.0]</td>\n", + " <td>[3.0340846]</td>\n", + " <td>[0.8060808]</td>\n", + " <td>[0.0]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>3</th>\n", + " <td>284500</td>\n", + " <td>92009</td>\n", + " <td>361106</td>\n", + " <td>1941.709961</td>\n", + " <td>1.320948</td>\n", + " <td>0.957973</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.000000</td>\n", + " <td>1.00582</td>\n", + " <td>1.0</td>\n", + " <td>True</td>\n", + " <td>False</td>\n", + " <td>False</td>\n", + " <td>2</td>\n", + " <td>[True, True]</td>\n", + " <td>[True, False]</td>\n", + " <td>[282535.56, 60798.965]</td>\n", + " <td>[-0.15060467, -0.712555]</td>\n", + " <td>[-2.322139, -1.9640917]</td>\n", + " <td>[285745.8, 76898.055]</td>\n", + " <td>[-0.007938385, 0.06525421]</td>\n", + " <td>[-1, 1]</td>\n", + " <td>[11, 11]</td>\n", + " <td>[True, False]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[475.2024, 2495.0132]</td>\n", + " <td>[0.0026975302, -0.07869638]</td>\n", + " <td>[0.4592046, 3.6803665]</td>\n", + " <td>11597.808594</td>\n", + " <td>0.950277</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>...</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0.0</td>\n", + " <td>[3379.703, 160.1774]</td>\n", + " <td>21881.703125</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1</td>\n", + " <td>[357909.62]</td>\n", + " <td>[-0.23516123]</td>\n", + " <td>[-2.2556815]</td>\n", + " <td>[377280.16]</td>\n", + " <td>[83818.125]</td>\n", + " <td>[1.0]</td>\n", + " <td>[0.13190588]</td>\n", + " <td>[0.18124953]</td>\n", + " <td>[1438.0938]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>4</th>\n", + " <td>284500</td>\n", + " <td>102088</td>\n", + " <td>361106</td>\n", + " <td>1941.239990</td>\n", + " <td>1.320948</td>\n", + " <td>0.948454</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.997393</td>\n", + " <td>1.00632</td>\n", + " <td>1.0</td>\n", + " <td>True</td>\n", + " <td>False</td>\n", + " <td>False</td>\n", + " <td>2</td>\n", + " <td>[True, True]</td>\n", + " <td>[True, True]</td>\n", + " <td>[193694.88, 135181.3]</td>\n", + " <td>[0.43973032, -0.09355457]</td>\n", + " <td>[-1.3102051, -1.2278727]</td>\n", + " <td>[212725.27, 135773.31]</td>\n", + " <td>[0.06219864, 0.013455391]</td>\n", + " <td>[-1, 1]</td>\n", + " <td>[11, 11]</td>\n", + " <td>[True, True]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[-260.75305, 333.5913]</td>\n", + " <td>[0.023009095, 0.02202719]</td>\n", + " <td>[3.0678267, 3.2605937]</td>\n", + " <td>33095.066406</td>\n", + " <td>-2.607571</td>\n", + " <td>3</td>\n", + " <td>[331520.7, 50229.3, 21196.365]</td>\n", + " <td>[-1.121277, 0.30034274, 1.7375907]</td>\n", + " <td>[1.6749763, -2.7881513, 1.2454554]</td>\n", + " <td>[563785.44, 52941.465, 62252.97]</td>\n", + " <td>[0.9958764, 0.9900437, 0.99060017]</td>\n", + " <td>[0, 0, 0]</td>\n", + " <td>[False, False, False]</td>\n", + " <td>...</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0.0</td>\n", + " <td>[1733.7571, 575.47437]</td>\n", + " <td>330.886719</td>\n", + " <td>[4881.326, 1914.4854, 2238.7397]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>2</td>\n", + " <td>[339715.34, 337192.75]</td>\n", + " <td>[0.22991252, -1.1238769]</td>\n", + " <td>[-1.2795317, 1.6654552]</td>\n", + " <td>[362021.66, 574320.1]</td>\n", + " <td>[97183.07, 30168.135]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[0.042992596, 2.0432703]</td>\n", + " <td>[0.075145245, 0.71639764]</td>\n", + " <td>[13062.9375, 5456.125]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " </tbody>\n", + "</table>\n", + "<p>5 rows × 81 columns</p>\n", + "</div>" + ], + "text/plain": [ + " runNumber eventNumber ... largeRjet_pt_syst tau_charge\n", + "0 284500 15440 ... [17716.031, 4227.1562] []\n", + "1 284500 81085 ... [22269.969, 14174.656, 24763.469] []\n", + "2 284500 89512 ... [0.0] []\n", + "3 284500 92009 ... [1438.0938] []\n", + "4 284500 102088 ... [13062.9375, 5456.125] []\n", + "\n", + "[5 rows x 81 columns]" + ] + }, + "execution_count": 8, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "df.head()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Un poco de datos estadÃsticos" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "scrolled": false + }, + "outputs": [ + { + "data": { + "text/html": [ + "<div>\n", + "<style scoped>\n", + " .dataframe tbody tr th:only-of-type {\n", + " vertical-align: middle;\n", + " }\n", + "\n", + " .dataframe tbody tr th {\n", + " vertical-align: top;\n", + " }\n", + "\n", + " .dataframe thead th {\n", + " text-align: right;\n", + " }\n", + "</style>\n", + "<table border=\"1\" class=\"dataframe\">\n", + " <thead>\n", + " <tr style=\"text-align: right;\">\n", + " <th></th>\n", + " <th>runNumber</th>\n", + " <th>eventNumber</th>\n", + " <th>channelNumber</th>\n", + " <th>mcWeight</th>\n", + " <th>scaleFactor_PILEUP</th>\n", + " <th>scaleFactor_ELE</th>\n", + " <th>scaleFactor_MUON</th>\n", + " <th>scaleFactor_PHOTON</th>\n", + " <th>scaleFactor_TAU</th>\n", + " <th>scaleFactor_BTAG</th>\n", + " <th>scaleFactor_LepTRIGGER</th>\n", + " <th>scaleFactor_PhotonTRIGGER</th>\n", + " <th>lep_n</th>\n", + " <th>met_et</th>\n", + " <th>met_phi</th>\n", + " <th>jet_n</th>\n", + " <th>photon_n</th>\n", + " <th>tau_n</th>\n", + " <th>ditau_m</th>\n", + " <th>met_et_syst</th>\n", + " <th>XSection</th>\n", + " <th>SumWeights</th>\n", + " <th>largeRjet_n</th>\n", + " </tr>\n", + " </thead>\n", + " <tbody>\n", + " <tr>\n", + " <th>count</th>\n", + " <td>53653.0</td>\n", + " <td>5.365300e+04</td>\n", + " <td>53653.0</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.0</td>\n", + " <td>53653.0</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.00000</td>\n", + " <td>53653.0</td>\n", + " <td>53653.000000</td>\n", + " <td>5.365300e+04</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.0</td>\n", + " <td>53653.0</td>\n", + " <td>53653.0</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.000000</td>\n", + " <td>5.365300e+04</td>\n", + " <td>53653.000000</td>\n", + " </tr>\n", + " <tr>\n", + " <th>mean</th>\n", + " <td>284500.0</td>\n", + " <td>3.972488e+07</td>\n", + " <td>361106.0</td>\n", + " <td>1906.595947</td>\n", + " <td>0.993552</td>\n", + " <td>0.961276</td>\n", + " <td>0.999988</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.000510</td>\n", + " <td>1.00101</td>\n", + " <td>1.0</td>\n", + " <td>1.727098</td>\n", + " <td>2.999400e+04</td>\n", + " <td>-0.004278</td>\n", + " <td>2.494436</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>7676.959961</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1.458912</td>\n", + " </tr>\n", + " <tr>\n", + " <th>std</th>\n", + " <td>0.0</td>\n", + " <td>2.308588e+07</td>\n", + " <td>0.0</td>\n", + " <td>365.987671</td>\n", + " <td>0.372941</td>\n", + " <td>0.026062</td>\n", + " <td>0.001164</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.029347</td>\n", + " <td>0.03176</td>\n", + " <td>0.0</td>\n", + " <td>0.455918</td>\n", + " <td>2.267055e+04</td>\n", + " <td>1.814127</td>\n", + " <td>1.376218</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>7742.734375</td>\n", + " <td>0.000000</td>\n", + " <td>0.000000e+00</td>\n", + " <td>0.511383</td>\n", + " </tr>\n", + " <tr>\n", + " <th>min</th>\n", + " <td>284500.0</td>\n", + " <td>1.010000e+03</td>\n", + " <td>361106.0</td>\n", + " <td>-1943.000000</td>\n", + " <td>0.000000</td>\n", + " <td>0.712289</td>\n", + " <td>0.825380</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.836130</td>\n", + " <td>0.00000</td>\n", + " <td>1.0</td>\n", + " <td>1.000000</td>\n", + " <td>8.181367e+01</td>\n", + " <td>-3.141528</td>\n", + " <td>0.000000</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.000000</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1.000000</td>\n", + " </tr>\n", + " <tr>\n", + " <th>25%</th>\n", + " <td>284500.0</td>\n", + " <td>1.964419e+07</td>\n", + " <td>361106.0</td>\n", + " <td>1941.060059</td>\n", + " <td>0.738662</td>\n", + " <td>0.947944</td>\n", + " <td>1.000000</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.997438</td>\n", + " <td>0.99866</td>\n", + " <td>1.0</td>\n", + " <td>1.000000</td>\n", + " <td>1.572494e+04</td>\n", + " <td>-1.579283</td>\n", + " <td>1.000000</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>2188.400391</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1.000000</td>\n", + " </tr>\n", + " <tr>\n", + " <th>50%</th>\n", + " <td>284500.0</td>\n", + " <td>3.940982e+07</td>\n", + " <td>361106.0</td>\n", + " <td>1941.400024</td>\n", + " <td>1.196105</td>\n", + " <td>0.963610</td>\n", + " <td>1.000000</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.998668</td>\n", + " <td>1.00306</td>\n", + " <td>1.0</td>\n", + " <td>2.000000</td>\n", + " <td>2.586413e+04</td>\n", + " <td>-0.010504</td>\n", + " <td>2.000000</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>5281.721680</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1.000000</td>\n", + " </tr>\n", + " <tr>\n", + " <th>75%</th>\n", + " <td>284500.0</td>\n", + " <td>5.977106e+07</td>\n", + " <td>361106.0</td>\n", + " <td>1941.739990</td>\n", + " <td>1.251332</td>\n", + " <td>0.978670</td>\n", + " <td>1.000000</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999264</td>\n", + " <td>1.00632</td>\n", + " <td>1.0</td>\n", + " <td>2.000000</td>\n", + " <td>3.917682e+04</td>\n", + " <td>1.572061</td>\n", + " <td>3.000000</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>10652.117188</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>2.000000</td>\n", + " </tr>\n", + " <tr>\n", + " <th>max</th>\n", + " <td>284500.0</td>\n", + " <td>7.999848e+07</td>\n", + " <td>361106.0</td>\n", + " <td>1943.439941</td>\n", + " <td>3.348151</td>\n", + " <td>1.112705</td>\n", + " <td>1.026367</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.752337</td>\n", + " <td>1.21491</td>\n", + " <td>1.0</td>\n", + " <td>4.000000</td>\n", + " <td>1.986488e+06</td>\n", + " <td>3.141495</td>\n", + " <td>11.000000</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>91205.000000</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>4.000000</td>\n", + " </tr>\n", + " </tbody>\n", + "</table>\n", + "</div>" + ], + "text/plain": [ + " runNumber eventNumber ... SumWeights largeRjet_n\n", + "count 53653.0 5.365300e+04 ... 5.365300e+04 53653.000000\n", + "mean 284500.0 3.972488e+07 ... 1.502776e+11 1.458912\n", + "std 0.0 2.308588e+07 ... 0.000000e+00 0.511383\n", + "min 284500.0 1.010000e+03 ... 1.502776e+11 1.000000\n", + "25% 284500.0 1.964419e+07 ... 1.502776e+11 1.000000\n", + "50% 284500.0 3.940982e+07 ... 1.502776e+11 1.000000\n", + "75% 284500.0 5.977106e+07 ... 1.502776e+11 2.000000\n", + "max 284500.0 7.999848e+07 ... 1.502776e+11 4.000000\n", + "\n", + "[8 rows x 23 columns]" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "df.describe()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Iniciamos el análisis" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": {}, + "outputs": [], + "source": [ + "canvas = ROOT.TCanvas(\"Canvas\",\"c\",800,600)" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "53653" + ] + }, + "execution_count": 11, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "tree = file2.Get(\"mini\")\n", + "tree.GetEntries()" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "metadata": {}, + "outputs": [], + "source": [ + "hist = ROOT.TH1F(\"variable\",\"Mass of the Z boson; mass [GeV]; events\",30,40,140)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## H→Zγ→γ ll*" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "> Prestar atención al Cut número 4, porque se debe al artÃculo de la búsqueda del decaimiento de Dalitz observado en ATLAS (http://cdsweb.cern.ch/record/2750822), es uno de las restricciones en la búsqueda de " + ] + }, + { + "cell_type": "code", + "execution_count": 13, + "metadata": {}, + "outputs": [], + "source": [ + "leadLepton = ROOT.TLorentzVector()\n", + "trailLepton = ROOT.TLorentzVector()\n", + "\n", + "for event in tree:\n", + " \n", + " # Cut #1: At least 2 leptons\n", + " if tree.lep_n >= 2:\n", + " \n", + " # Cut #2: Leptons with opposite charge\n", + " if (tree.lep_charge[0] != tree.lep_charge[1]):\n", + " \n", + " # Cut #3: Leptons of the same family (2 electrons or 2 muons)\n", + " if (tree.lep_type[0] == tree.lep_type[1]):\n", + " \n", + " # Cut #4: E < 30 GeV\n", + " if tree.lep_E[0] < 30000:\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " # Let's define one TLorentz vector for each, e.i. two vectors!\n", + " leadLepton.SetPtEtaPhiE(tree.lep_pt[0]/1000., tree.lep_eta[0], tree.lep_phi[0], tree.lep_E[0]/1000.)\n", + " trailLepton.SetPtEtaPhiE(tree.lep_pt[1]/1000., tree.lep_eta[1], tree.lep_phi[1], tree.lep_E[1]/1000.)\n", + " # Next line: addition of two TLorentz vectors above --> ask mass very easy (devide by 1000 to get value in GeV)\n", + " invmass = leadLepton + trailLepton\n", + "\n", + " hist.Fill(invmass.M())" + ] + }, + { + "cell_type": "code", + "execution_count": 14, + "metadata": {}, + "outputs": [], + "source": [ + "hist.Draw()\n", + "hist.SetFillColor(3)" + ] + }, + { + "cell_type": "code", + "execution_count": 15, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "\n", + "<div id=\"root_plot_1\"\n", + " style=\"width: 796px; height: 572px\">\n", + "</div>\n", + "\n", + "<script>\n", + " requirejs.config({\n", + " paths: {\n", + " 'JSRootCore' : 'https://root.cern.ch/js/notebook//scripts/JSRootCore',\n", + " }\n", + " });\n", + " require(['JSRootCore'],\n", + " function(Core) {\n", + " var obj = Core.JSONR_unref({\"_typename\":\"TCanvas\",\"fUniqueID\":0,\"fBits\":53673992,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fLeftMargin\":0.1,\"fRightMargin\":0.1,\"fBottomMargin\":0.1,\"fTopMargin\":0.1,\"fXfile\":2,\"fYfile\":2,\"fAfile\":1,\"fXstat\":0.99,\"fYstat\":0.99,\"fAstat\":2,\"fFrameFillColor\":0,\"fFrameLineColor\":1,\"fFrameFillStyle\":1001,\"fFrameLineStyle\":1,\"fFrameLineWidth\":1,\"fFrameBorderSize\":1,\"fFrameBorderMode\":0,\"fX1\":27.4999990686774,\"fY1\":-0.262500019557775,\"fX2\":152.500000931323,\"fY2\":2.36250001955778,\"fXtoAbsPixelk\":-175.119941459846,\"fXtoPixelk\":-175.119941459846,\"fXtoPixel\":6.3679999051094,\"fYtoAbsPixelk\":514.800046590614,\"fYtoPixelk\":514.800046590614,\"fYtoPixel\":-217.904758657728,\"fUtoAbsPixelk\":5e-5,\"fUtoPixelk\":5e-5,\"fUtoPixel\":796,\"fVtoAbsPixelk\":572.00005,\"fVtoPixelk\":572,\"fVtoPixel\":-572,\"fAbsPixeltoXk\":27.4999990686774,\"fPixeltoXk\":27.4999990686774,\"fPixeltoX\":0.157035178219404,\"fAbsPixeltoYk\":2.36250001955778,\"fPixeltoYk\":-0.262500019557775,\"fPixeltoY\":-0.00458916090754467,\"fXlowNDC\":0,\"fYlowNDC\":0,\"fXUpNDC\":0,\"fYUpNDC\":0,\"fWNDC\":1,\"fHNDC\":1,\"fAbsXlowNDC\":0,\"fAbsYlowNDC\":0,\"fAbsWNDC\":1,\"fAbsHNDC\":1,\"fUxmin\":40,\"fUymin\":0,\"fUxmax\":140,\"fUymax\":2.1,\"fTheta\":30,\"fPhi\":30,\"fAspectRatio\":0,\"fNumber\":0,\"fTickx\":0,\"fTicky\":0,\"fLogx\":0,\"fLogy\":0,\"fLogz\":0,\"fPadPaint\":0,\"fCrosshair\":0,\"fCrosshairPos\":0,\"fBorderSize\":2,\"fBorderMode\":0,\"fModified\":false,\"fGridx\":false,\"fGridy\":false,\"fAbsCoord\":false,\"fEditable\":true,\"fFixedAspectRatio\":false,\"fPrimitives\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TFrame\",\"fUniqueID\":0,\"fBits\":50331656,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fX1\":40,\"fY1\":0,\"fX2\":140,\"fY2\":2.1,\"fBorderSize\":1,\"fBorderMode\":0},{\"_typename\":\"TH1F\",\"fUniqueID\":0,\"fBits\":50331656,\"fName\":\"variable\",\"fTitle\":\"Mass of the Z boson\",\"fLineColor\":602,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":3,\"fFillStyle\":1001,\"fMarkerColor\":1,\"fMarkerStyle\":1,\"fMarkerSize\":1,\"fNcells\":32,\"fXaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"xaxis\",\"fTitle\":\" mass [GeV]\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":30,\"fXmin\":40,\"fXmax\":140,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fYaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"yaxis\",\"fTitle\":\" events\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":0,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fZaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"zaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fBarOffset\":0,\"fBarWidth\":1000,\"fEntries\":5,\"fTsumw\":5,\"fTsumw2\":5,\"fTsumwx\":460.606963185387,\"fTsumwx2\":42465.3851897774,\"fMaximum\":-1111,\"fMinimum\":-1111,\"fNormFactor\":0,\"fContour\":[],\"fSumw2\":[],\"fOption\":\"\",\"fFunctions\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TPaveStats\",\"fUniqueID\":0,\"fBits\":50331657,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fX1\":125.000002533197,\"fY1\":1.771875026403,\"fX2\":150.000003278256,\"fY2\":2.19187502327375,\"fX1NDC\":0.780000016093254,\"fY1NDC\":0.775000005960464,\"fX2NDC\":0.980000019073486,\"fY2NDC\":0.935000002384186,\"fBorderSize\":1,\"fInit\":1,\"fShadowColor\":1,\"fCornerRadius\":0,\"fOption\":\"brNDC\",\"fName\":\"stats\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":12,\"fTextColor\":1,\"fTextFont\":42,\"fLabel\":\"\",\"fLongest\":18,\"fMargin\":0.05,\"fLines\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"variable\",\"fTextAngle\":0,\"fTextSize\":0.0368,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.0368000008165836},{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Entries = 5 \",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.04},{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Mean = 92.12\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.04},{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Std Dev = 2.593\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.04}],\"opt\":[\"\",\"\",\"\",\"\"]},\"fOptFit\":0,\"fOptStat\":1111,\"fFitFormat\":\"5.4g\",\"fStatFormat\":\"6.4g\",\"fParent\":{\"$ref\":3}}],\"opt\":[\"brNDC\"]},\"fBufferSize\":0,\"fBuffer\":[],\"fBinStatErrOpt\":0,\"fStatOverflows\":2,\"fArray\":[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,2,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0]},{\"_typename\":\"TPaveText\",\"fUniqueID\":0,\"fBits\":50331657,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":0,\"fX1\":66.8435926197668,\"fY1\":2.20500002757646,\"fX2\":113.156407380233,\"fY2\":2.34937503187917,\"fX1NDC\":0.314748743718593,\"fY1NDC\":0.940000003948808,\"fX2NDC\":0.685251256281407,\"fY2NDC\":0.995000004768372,\"fBorderSize\":0,\"fInit\":1,\"fShadowColor\":1,\"fCornerRadius\":0,\"fOption\":\"blNDC\",\"fName\":\"title\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":22,\"fTextColor\":1,\"fTextFont\":42,\"fLabel\":\"\",\"fLongest\":19,\"fMargin\":0.05,\"fLines\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Mass of the Z boson\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.0467500016093254}],\"opt\":[\"\"]}}],\"opt\":[\"\",\"\",\"blNDC\"]},\"fExecs\":null,\"fName\":\"Canvas\",\"fTitle\":\"c\",\"fNumPaletteColor\":0,\"fNextPaletteColor\":0,\"fDISPLAY\":\"$DISPLAY\",\"fDoubleBuffer\":0,\"fRetained\":true,\"fXsizeUser\":0,\"fYsizeUser\":0,\"fXsizeReal\":20,\"fYsizeReal\":15,\"fWindowTopX\":0,\"fWindowTopY\":0,\"fWindowWidth\":0,\"fWindowHeight\":0,\"fCw\":796,\"fCh\":572,\"fCatt\":{\"_typename\":\"TAttCanvas\",\"fXBetween\":2,\"fYBetween\":2,\"fTitleFromTop\":1.2,\"fXdate\":0.2,\"fYdate\":0.3,\"fAdate\":1},\"kMoveOpaque\":true,\"kResizeOpaque\":true,\"fHighLightColor\":2,\"fBatch\":true,\"kShowEventStatus\":false,\"kAutoExec\":true,\"kMenuBar\":true});\n", + " Core.draw(\"root_plot_1\", obj, \"\");\n", + " }\n", + " );\n", + "</script>\n" + ], + "text/plain": [ + "<IPython.core.display.HTML object>" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "canvas.Draw()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Missing transverse energy de los leptones en estado final" + ] + }, + { + "cell_type": "code", + "execution_count": 16, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "Warning in <TCanvas::Constructor>: Deleting canvas with same name: Canvas\n" + ] + } + ], + "source": [ + "canvas = ROOT.TCanvas(\"Canvas\",\"cz\",800,600)" + ] + }, + { + "cell_type": "code", + "execution_count": 17, + "metadata": {}, + "outputs": [], + "source": [ + "hist = ROOT.TH1F(\"X\",\"Y [GeV] ; events\",30,105,160)" + ] + }, + { + "cell_type": "code", + "execution_count": 18, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "53653" + ] + }, + "execution_count": 18, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "tree = file2.Get(\"mini\")\n", + "tree.GetEntries()" + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "metadata": {}, + "outputs": [], + "source": [ + "c = ROOT.TCanvas(\"testCanvas\",\"a first way to plot a variable\",800,600)" + ] + }, + { + "cell_type": "code", + "execution_count": 27, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "Warning in <TROOT::Append>: Replacing existing TH1: h_sig (Potential memory leak).\n", + "Warning in <TROOT::Append>: Replacing existing TH1: h2_sig (Potential memory leak).\n" + ] + } + ], + "source": [ + "h_sig = ROOT.TH1F(\"h_sig\",\"Example plot: Missing transverse energy\",20,0,200)\n", + "h2_sig = ROOT.TH1F(\"h2_sig\",\"Example plot: Number of Jets\",10,0,10)" + ] + }, + { + "cell_type": "code", + "execution_count": 28, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "10000\n", + "20000\n", + "30000\n", + "40000\n", + "50000\n" + ] + } + ], + "source": [ + "m=0 \n", + "for event in tree:\n", + " m += 1\n", + " ## printing the evolution in number of events\n", + " if(m%10000==0):\n", + " print(m)\n", + " h_sig.Fill((tree.met_et)/1000.)" + ] + }, + { + "cell_type": "code", + "execution_count": 29, + "metadata": {}, + "outputs": [], + "source": [ + "scale_sig = h_sig.Integral()\n", + "h_sig.Scale(1/scale_sig)" + ] + }, + { + "cell_type": "code", + "execution_count": 30, + "metadata": {}, + "outputs": [], + "source": [ + "h_sig.SetFillStyle(3003)\n", + "h_sig.SetFillColor(2)\n", + "h_sig.SetLineColor(2)" + ] + }, + { + "cell_type": "code", + "execution_count": 31, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "<ROOT.TLegendEntry object (\"TLegendEntry\") at 0x15c191c0>" + ] + }, + "execution_count": 31, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "legend=ROOT.TLegend(0.5,0.7,0.9,0.9)\n", + "legend.AddEntry(h_sig,\"Signal (Z #rightarrow ll*)\",\"l\")" + ] + }, + { + "cell_type": "code", + "execution_count": 32, + "metadata": {}, + "outputs": [], + "source": [ + "h_sig.SetStats(0)" + ] + }, + { + "cell_type": "code", + "execution_count": 33, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "\n", + "<div id=\"root_plot_3\"\n", + " style=\"width: 796px; height: 572px\">\n", + "</div>\n", + "\n", + "<script>\n", + " requirejs.config({\n", + " paths: {\n", + " 'JSRootCore' : 'https://root.cern.ch/js/notebook//scripts/JSRootCore',\n", + " }\n", + " });\n", + " require(['JSRootCore'],\n", + " function(Core) {\n", + " var obj = Core.JSONR_unref({\"_typename\":\"TCanvas\",\"fUniqueID\":0,\"fBits\":53673992,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fLeftMargin\":0.1,\"fRightMargin\":0.1,\"fBottomMargin\":0.1,\"fTopMargin\":0.1,\"fXfile\":2,\"fYfile\":2,\"fAfile\":1,\"fXstat\":0.99,\"fYstat\":0.99,\"fAstat\":2,\"fFrameFillColor\":0,\"fFrameLineColor\":1,\"fFrameFillStyle\":1001,\"fFrameLineStyle\":1,\"fFrameLineWidth\":1,\"fFrameBorderSize\":1,\"fFrameBorderMode\":0,\"fX1\":-25.0000018626453,\"fY1\":-0.0313701751731466,\"fX2\":225.000001862645,\"fY2\":0.282331557860238,\"fXtoAbsPixelk\":79.6000547445299,\"fXtoPixelk\":79.6000547445299,\"fXtoPixel\":3.1839999525547,\"fYtoAbsPixelk\":514.800046590614,\"fYtoPixelk\":514.800046590614,\"fYtoPixel\":-1823.38807780551,\"fUtoAbsPixelk\":5e-5,\"fUtoPixelk\":5e-5,\"fUtoPixel\":796,\"fVtoAbsPixelk\":572.00005,\"fVtoPixelk\":572,\"fVtoPixel\":-572,\"fAbsPixeltoXk\":-25.0000018626453,\"fPixeltoXk\":-25.0000018626453,\"fPixeltoX\":0.314070356438807,\"fAbsPixeltoYk\":0.282331557860238,\"fPixeltoYk\":-0.0313701751731466,\"fPixeltoY\":-5.48429603205219e-4,\"fXlowNDC\":0,\"fYlowNDC\":0,\"fXUpNDC\":0,\"fYUpNDC\":0,\"fWNDC\":1,\"fHNDC\":1,\"fAbsXlowNDC\":0,\"fAbsYlowNDC\":0,\"fAbsWNDC\":1,\"fAbsHNDC\":1,\"fUxmin\":0,\"fUymin\":0,\"fUxmax\":200,\"fUymax\":0.250961382687092,\"fTheta\":30,\"fPhi\":30,\"fAspectRatio\":0,\"fNumber\":0,\"fTickx\":0,\"fTicky\":0,\"fLogx\":0,\"fLogy\":0,\"fLogz\":0,\"fPadPaint\":0,\"fCrosshair\":0,\"fCrosshairPos\":0,\"fBorderSize\":2,\"fBorderMode\":0,\"fModified\":false,\"fGridx\":false,\"fGridy\":false,\"fAbsCoord\":false,\"fEditable\":true,\"fFixedAspectRatio\":false,\"fPrimitives\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TFrame\",\"fUniqueID\":0,\"fBits\":50331656,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fX1\":0,\"fY1\":0,\"fX2\":200,\"fY2\":0.250961382687092,\"fBorderSize\":1,\"fBorderMode\":0},{\"_typename\":\"TH1F\",\"fUniqueID\":0,\"fBits\":50332168,\"fName\":\"h_sig\",\"fTitle\":\"Example plot: Missing transverse energy\",\"fLineColor\":2,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":2,\"fFillStyle\":3003,\"fMarkerColor\":1,\"fMarkerStyle\":1,\"fMarkerSize\":1,\"fNcells\":22,\"fXaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"xaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":20,\"fXmin\":0,\"fXmax\":200,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fYaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"yaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":0,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fZaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"zaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fBarOffset\":0,\"fBarWidth\":1000,\"fEntries\":53653,\"fTsumw\":53621,\"fTsumw2\":53621,\"fTsumwx\":1598245.92758354,\"fTsumwx2\":68461690.7118449,\"fMaximum\":-1111,\"fMinimum\":-1111,\"fNormFactor\":0,\"fContour\":[],\"fSumw2\":[0,2.17340522377219e-6,4.45741100141852e-6,4.39028550802949e-6,3.1802876246076e-6,1.92890148360386e-6,1.1386987842263e-6,6.21867265179176e-7,3.33888464525732e-7,1.81551852585867e-7,9.77319359705528e-8,5.63436783887173e-8,3.23454450009303e-8,1.87812261295724e-8,1.2173016935834e-8,8.3472116131433e-9,5.9126082259765e-9,4.52140629045262e-9,3.82580532269068e-9,1.73900241940485e-9,1.39120193552388e-9,1.11296154841911e-8],\"fOption\":\"\",\"fFunctions\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[],\"opt\":[]},\"fBufferSize\":0,\"fBuffer\":[],\"fBinStatErrOpt\":0,\"fStatOverflows\":2,\"fArray\":[0,0.1165402,0.2390108,0.2354115,0.1705302,0.1034296,0.06105817,0.03334514,0.01790343,0.009734992,0.005240484,0.003021204,0.001734395,0.001007068,6.527294e-4,4.475858e-4,3.1704e-4,2.424423e-4,2.051435e-4,9.324705e-5,7.459764e-5,5.967811e-4]},{\"_typename\":\"TLegend\",\"fUniqueID\":0,\"fBits\":50331656,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fX1\":100,\"fY1\":0.188221037950223,\"fX2\":200.000001490116,\"fY2\":0.2509613845569,\"fX1NDC\":0.5,\"fY1NDC\":0.7,\"fX2NDC\":0.9,\"fY2NDC\":0.9,\"fBorderSize\":1,\"fInit\":1,\"fShadowColor\":1,\"fCornerRadius\":0,\"fOption\":\"brNDC\",\"fName\":\"TPave\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":12,\"fTextColor\":1,\"fTextFont\":42,\"fPrimitives\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TLegendEntry\",\"fUniqueID\":0,\"fBits\":50331656,\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":42,\"fLineColor\":2,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":0,\"fMarkerColor\":1,\"fMarkerStyle\":21,\"fMarkerSize\":1,\"fObject\":{\"$ref\":3},\"fLabel\":\"Signal (Z #rightarrow ll*)\",\"fOption\":\"l\"}],\"opt\":[\"l\"]},\"fEntrySeparation\":0.1,\"fMargin\":0.25,\"fNColumns\":1,\"fColumnSeparation\":0},{\"_typename\":\"TPaveText\",\"fUniqueID\":0,\"fBits\":50331657,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":0,\"fX1\":12.4999986961483,\"fY1\":0.261677700499376,\"fX2\":187.500001303852,\"fY2\":0.280763050690918,\"fX1NDC\":0.15,\"fY1NDC\":0.934160843929211,\"fX2NDC\":0.85,\"fY2NDC\":0.995000004768372,\"fBorderSize\":0,\"fInit\":1,\"fShadowColor\":1,\"fCornerRadius\":0,\"fOption\":\"blNDC\",\"fName\":\"title\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":22,\"fTextColor\":1,\"fTextFont\":42,\"fLabel\":\"\",\"fLongest\":39,\"fMargin\":0.05,\"fLines\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Example plot: Missing transverse energy\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.0517132878303528}],\"opt\":[\"\"]}}],\"opt\":[\"\",\"hist\",\"\",\"blNDC\"]},\"fExecs\":null,\"fName\":\"testCanvas\",\"fTitle\":\"a first way to plot a variable\",\"fNumPaletteColor\":0,\"fNextPaletteColor\":0,\"fDISPLAY\":\"$DISPLAY\",\"fDoubleBuffer\":0,\"fRetained\":true,\"fXsizeUser\":0,\"fYsizeUser\":0,\"fXsizeReal\":20,\"fYsizeReal\":15,\"fWindowTopX\":0,\"fWindowTopY\":0,\"fWindowWidth\":0,\"fWindowHeight\":0,\"fCw\":796,\"fCh\":572,\"fCatt\":{\"_typename\":\"TAttCanvas\",\"fXBetween\":2,\"fYBetween\":2,\"fTitleFromTop\":1.2,\"fXdate\":0.2,\"fYdate\":0.3,\"fAdate\":1},\"kMoveOpaque\":true,\"kResizeOpaque\":true,\"fHighLightColor\":2,\"fBatch\":true,\"kShowEventStatus\":false,\"kAutoExec\":true,\"kMenuBar\":true});\n", + " Core.draw(\"root_plot_3\", obj, \"\");\n", + " }\n", + " );\n", + "</script>\n" + ], + "text/plain": [ + "<IPython.core.display.HTML object>" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "h_sig.Draw(\"hist\")\n", + "legend.Draw()\n", + "c.Draw()" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3", + "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.6.9" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +} diff --git a/codigo/Data_Analysis.ipynb b/codigo/Data_Analysis.ipynb new file mode 100644 index 0000000000000000000000000000000000000000..852ecbee34ae9cefd02e0cfa03a5b388168e26e7 --- /dev/null +++ b/codigo/Data_Analysis.ipynb @@ -0,0 +1,1782 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Z -> ll (decaimiento del boson Z a dos leptones) \n", + "> Proveniente de https://atlas-opendata.web.cern.ch/atlas-opendata/samples/2020/2lep/MC/mc_361106.Zee.1largeRjet1lep.root" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Importamos los paquetes necesarios" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Welcome to JupyROOT 6.18/04\n" + ] + } + ], + "source": [ + "import ROOT\n", + "from ROOT import TMath\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Usamos pandas para leer los archivos .root" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [], + "source": [ + "import pandas as pd\n", + "from root_pandas import read_root" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [], + "source": [ + "%jsroot on" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Z boson invariant mass (reconstruction)" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": {}, + "outputs": [], + "source": [ + "file2 = ROOT.TFile.Open(\"/home/student/ejercicios-clase-08-datos/data-used/mc_361106.Zee.1largeRjet1lep.root\")\n", + "# file2 = ROOT.TFile.Open(\"/https://atlas-opendata.web.cern.ch/atlas-opendata/samples/2020/2lep/MC/mc_361106.Zee.1largeRjet1lep.root\")\n", + "# desde la web" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Visualizamos el archivo aprovechando que no es tan grande, podemos visualizarlo como si fuera un archivo csv" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [], + "source": [ + "df = read_root('/home/student/ejercicios-clase-08-datos/data-used/mc_361106.Zee.1largeRjet1lep.root')" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "scrolled": true + }, + "outputs": [ + { + "data": { + "text/html": [ + "<div>\n", + "<style scoped>\n", + " .dataframe tbody tr th:only-of-type {\n", + " vertical-align: middle;\n", + " }\n", + "\n", + " .dataframe tbody tr th {\n", + " vertical-align: top;\n", + " }\n", + "\n", + " .dataframe thead th {\n", + " text-align: right;\n", + " }\n", + "</style>\n", + "<table border=\"1\" class=\"dataframe\">\n", + " <thead>\n", + " <tr style=\"text-align: right;\">\n", + " <th></th>\n", + " <th>runNumber</th>\n", + " <th>eventNumber</th>\n", + " <th>channelNumber</th>\n", + " <th>mcWeight</th>\n", + " <th>scaleFactor_PILEUP</th>\n", + " <th>scaleFactor_ELE</th>\n", + " <th>scaleFactor_MUON</th>\n", + " <th>scaleFactor_PHOTON</th>\n", + " <th>scaleFactor_TAU</th>\n", + " <th>scaleFactor_BTAG</th>\n", + " <th>...</th>\n", + " <th>largeRjet_pt</th>\n", + " <th>largeRjet_eta</th>\n", + " <th>largeRjet_phi</th>\n", + " <th>largeRjet_E</th>\n", + " <th>largeRjet_m</th>\n", + " <th>largeRjet_truthMatched</th>\n", + " <th>largeRjet_D2</th>\n", + " <th>largeRjet_tau32</th>\n", + " <th>largeRjet_pt_syst</th>\n", + " <th>tau_charge</th>\n", + " </tr>\n", + " </thead>\n", + " <tbody>\n", + " <tr>\n", + " <th>0</th>\n", + " <td>284500</td>\n", + " <td>15440</td>\n", + " <td>361106</td>\n", + " <td>1940.819946</td>\n", + " <td>1.264594</td>\n", + " <td>0.941115</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999384</td>\n", + " <td>...</td>\n", + " <td>[336389.22, 326111.62]</td>\n", + " <td>[-1.6073604, 0.34273493]</td>\n", + " <td>[-3.1196432, -0.05282484]</td>\n", + " <td>[878163.0, 352508.03]</td>\n", + " <td>[95666.125, 70169.38]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[0.06860469, 1.7374916]</td>\n", + " <td>[0.5651315, 0.6696122]</td>\n", + " <td>[17716.031, 4227.1562]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>1</th>\n", + " <td>284500</td>\n", + " <td>81085</td>\n", + " <td>361106</td>\n", + " <td>1941.910034</td>\n", + " <td>1.320948</td>\n", + " <td>0.992704</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.997777</td>\n", + " <td>...</td>\n", + " <td>[424948.16, 295925.38, 288773.56]</td>\n", + " <td>[-0.50103956, 0.00353557, 0.7989666]</td>\n", + " <td>[0.39738846, 3.0375211, -2.1750822]</td>\n", + " <td>[494285.72, 314167.06, 393288.84]</td>\n", + " <td>[120339.15, 105489.414, 75616.88]</td>\n", + " <td>[0.0, 0.0, 0.0]</td>\n", + " <td>[1.4906838, 0.1404522, 0.62883687]</td>\n", + " <td>[0.64305073, 0.19809924, 0.69501674]</td>\n", + " <td>[22269.969, 14174.656, 24763.469]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>2</th>\n", + " <td>284500</td>\n", + " <td>89512</td>\n", + " <td>361106</td>\n", + " <td>1940.760010</td>\n", + " <td>1.251332</td>\n", + " <td>0.993904</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999102</td>\n", + " <td>...</td>\n", + " <td>[255974.6]</td>\n", + " <td>[1.6731088]</td>\n", + " <td>[0.80495465]</td>\n", + " <td>[706647.8]</td>\n", + " <td>[29652.18]</td>\n", + " <td>[0.0]</td>\n", + " <td>[3.0340846]</td>\n", + " <td>[0.8060808]</td>\n", + " <td>[0.0]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>3</th>\n", + " <td>284500</td>\n", + " <td>92009</td>\n", + " <td>361106</td>\n", + " <td>1941.709961</td>\n", + " <td>1.320948</td>\n", + " <td>0.957973</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.000000</td>\n", + " <td>...</td>\n", + " <td>[357909.62]</td>\n", + " <td>[-0.23516123]</td>\n", + " <td>[-2.2556815]</td>\n", + " <td>[377280.16]</td>\n", + " <td>[83818.125]</td>\n", + " <td>[1.0]</td>\n", + " <td>[0.13190588]</td>\n", + " <td>[0.18124953]</td>\n", + " <td>[1438.0938]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>4</th>\n", + " <td>284500</td>\n", + " <td>102088</td>\n", + " <td>361106</td>\n", + " <td>1941.239990</td>\n", + " <td>1.320948</td>\n", + " <td>0.948454</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.997393</td>\n", + " <td>...</td>\n", + " <td>[339715.34, 337192.75]</td>\n", + " <td>[0.22991252, -1.1238769]</td>\n", + " <td>[-1.2795317, 1.6654552]</td>\n", + " <td>[362021.66, 574320.1]</td>\n", + " <td>[97183.07, 30168.135]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[0.042992596, 2.0432703]</td>\n", + " <td>[0.075145245, 0.71639764]</td>\n", + " <td>[13062.9375, 5456.125]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>...</th>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " <td>...</td>\n", + " </tr>\n", + " <tr>\n", + " <th>53648</th>\n", + " <td>284500</td>\n", + " <td>5740044</td>\n", + " <td>361106</td>\n", + " <td>1941.760010</td>\n", + " <td>0.934886</td>\n", + " <td>0.986970</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999438</td>\n", + " <td>...</td>\n", + " <td>[272509.03, 264178.5]</td>\n", + " <td>[-0.11495946, -0.88012326]</td>\n", + " <td>[-2.6213317, 0.5289878]</td>\n", + " <td>[274337.56, 373833.25]</td>\n", + " <td>[3766.6348, 20426.285]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[0.0, 2.0461884]</td>\n", + " <td>[-1.0, 0.7430232]</td>\n", + " <td>[17955.438, 0.0]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>53649</th>\n", + " <td>284500</td>\n", + " <td>5752321</td>\n", + " <td>361106</td>\n", + " <td>1941.680054</td>\n", + " <td>0.596562</td>\n", + " <td>0.889910</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999002</td>\n", + " <td>...</td>\n", + " <td>[272032.44]</td>\n", + " <td>[1.6677606]</td>\n", + " <td>[-2.4962265]</td>\n", + " <td>[746732.94]</td>\n", + " <td>[14772.887]</td>\n", + " <td>[0.0]</td>\n", + " <td>[3.902746]</td>\n", + " <td>[0.58085346]</td>\n", + " <td>[21161.61]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>53650</th>\n", + " <td>284500</td>\n", + " <td>5757742</td>\n", + " <td>361106</td>\n", + " <td>1941.000000</td>\n", + " <td>1.251332</td>\n", + " <td>0.984580</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.998199</td>\n", + " <td>...</td>\n", + " <td>[262801.75]</td>\n", + " <td>[-0.25520322]</td>\n", + " <td>[-1.0126089]</td>\n", + " <td>[280588.1]</td>\n", + " <td>[71192.08]</td>\n", + " <td>[0.0]</td>\n", + " <td>[0.036484536]</td>\n", + " <td>[0.39032233]</td>\n", + " <td>[7288.8906]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>53651</th>\n", + " <td>284500</td>\n", + " <td>5745755</td>\n", + " <td>361106</td>\n", + " <td>1940.719971</td>\n", + " <td>1.196105</td>\n", + " <td>0.964722</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999469</td>\n", + " <td>...</td>\n", + " <td>[266810.53, 260362.53]</td>\n", + " <td>[-0.47211164, 0.4879255]</td>\n", + " <td>[-1.3692597, 1.7747238]</td>\n", + " <td>[298000.4, 312144.3]</td>\n", + " <td>[23127.982, 110384.94]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[1.6352673, 0.0088554425]</td>\n", + " <td>[0.65562195, -1.0]</td>\n", + " <td>[1491.4062, 8128.9688]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>53652</th>\n", + " <td>284500</td>\n", + " <td>5760130</td>\n", + " <td>361106</td>\n", + " <td>1940.819946</td>\n", + " <td>0.506784</td>\n", + " <td>1.006096</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.000000</td>\n", + " <td>...</td>\n", + " <td>[326995.22, 289713.94]</td>\n", + " <td>[0.7875121, 0.041406196]</td>\n", + " <td>[-1.3435338, 1.6638559]</td>\n", + " <td>[433746.4, 289980.38]</td>\n", + " <td>[1893.2017, 3234.1577]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[77.1253, 34.656742]</td>\n", + " <td>[-1.0, -1.0]</td>\n", + " <td>[788.78125, 122.65625]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " </tbody>\n", + "</table>\n", + "<p>53653 rows × 81 columns</p>\n", + "</div>" + ], + "text/plain": [ + " runNumber eventNumber channelNumber mcWeight scaleFactor_PILEUP \\\n", + "0 284500 15440 361106 1940.819946 1.264594 \n", + "1 284500 81085 361106 1941.910034 1.320948 \n", + "2 284500 89512 361106 1940.760010 1.251332 \n", + "3 284500 92009 361106 1941.709961 1.320948 \n", + "4 284500 102088 361106 1941.239990 1.320948 \n", + "... ... ... ... ... ... \n", + "53648 284500 5740044 361106 1941.760010 0.934886 \n", + "53649 284500 5752321 361106 1941.680054 0.596562 \n", + "53650 284500 5757742 361106 1941.000000 1.251332 \n", + "53651 284500 5745755 361106 1940.719971 1.196105 \n", + "53652 284500 5760130 361106 1940.819946 0.506784 \n", + "\n", + " scaleFactor_ELE scaleFactor_MUON scaleFactor_PHOTON scaleFactor_TAU \\\n", + "0 0.941115 1.0 1.0 1.0 \n", + "1 0.992704 1.0 1.0 1.0 \n", + "2 0.993904 1.0 1.0 1.0 \n", + "3 0.957973 1.0 1.0 1.0 \n", + "4 0.948454 1.0 1.0 1.0 \n", + "... ... ... ... ... \n", + "53648 0.986970 1.0 1.0 1.0 \n", + "53649 0.889910 1.0 1.0 1.0 \n", + "53650 0.984580 1.0 1.0 1.0 \n", + "53651 0.964722 1.0 1.0 1.0 \n", + "53652 1.006096 1.0 1.0 1.0 \n", + "\n", + " scaleFactor_BTAG ... largeRjet_pt \\\n", + "0 0.999384 ... [336389.22, 326111.62] \n", + "1 0.997777 ... [424948.16, 295925.38, 288773.56] \n", + "2 0.999102 ... [255974.6] \n", + "3 1.000000 ... [357909.62] \n", + "4 0.997393 ... [339715.34, 337192.75] \n", + "... ... ... ... \n", + "53648 0.999438 ... [272509.03, 264178.5] \n", + "53649 0.999002 ... [272032.44] \n", + "53650 0.998199 ... [262801.75] \n", + "53651 0.999469 ... [266810.53, 260362.53] \n", + "53652 1.000000 ... [326995.22, 289713.94] \n", + "\n", + " largeRjet_eta \\\n", + "0 [-1.6073604, 0.34273493] \n", + "1 [-0.50103956, 0.00353557, 0.7989666] \n", + "2 [1.6731088] \n", + "3 [-0.23516123] \n", + "4 [0.22991252, -1.1238769] \n", + "... ... \n", + "53648 [-0.11495946, -0.88012326] \n", + "53649 [1.6677606] \n", + "53650 [-0.25520322] \n", + "53651 [-0.47211164, 0.4879255] \n", + "53652 [0.7875121, 0.041406196] \n", + "\n", + " largeRjet_phi largeRjet_E \\\n", + "0 [-3.1196432, -0.05282484] [878163.0, 352508.03] \n", + "1 [0.39738846, 3.0375211, -2.1750822] [494285.72, 314167.06, 393288.84] \n", + "2 [0.80495465] [706647.8] \n", + "3 [-2.2556815] [377280.16] \n", + "4 [-1.2795317, 1.6654552] [362021.66, 574320.1] \n", + "... ... ... \n", + "53648 [-2.6213317, 0.5289878] [274337.56, 373833.25] \n", + "53649 [-2.4962265] [746732.94] \n", + "53650 [-1.0126089] [280588.1] \n", + "53651 [-1.3692597, 1.7747238] [298000.4, 312144.3] \n", + "53652 [-1.3435338, 1.6638559] [433746.4, 289980.38] \n", + "\n", + " largeRjet_m largeRjet_truthMatched \\\n", + "0 [95666.125, 70169.38] [0.0, 0.0] \n", + "1 [120339.15, 105489.414, 75616.88] [0.0, 0.0, 0.0] \n", + "2 [29652.18] [0.0] \n", + "3 [83818.125] [1.0] \n", + "4 [97183.07, 30168.135] [0.0, 0.0] \n", + "... ... ... \n", + "53648 [3766.6348, 20426.285] [0.0, 0.0] \n", + "53649 [14772.887] [0.0] \n", + "53650 [71192.08] [0.0] \n", + "53651 [23127.982, 110384.94] [0.0, 0.0] \n", + "53652 [1893.2017, 3234.1577] [0.0, 0.0] \n", + "\n", + " largeRjet_D2 \\\n", + "0 [0.06860469, 1.7374916] \n", + "1 [1.4906838, 0.1404522, 0.62883687] \n", + "2 [3.0340846] \n", + "3 [0.13190588] \n", + "4 [0.042992596, 2.0432703] \n", + "... ... \n", + "53648 [0.0, 2.0461884] \n", + "53649 [3.902746] \n", + "53650 [0.036484536] \n", + "53651 [1.6352673, 0.0088554425] \n", + "53652 [77.1253, 34.656742] \n", + "\n", + " largeRjet_tau32 \\\n", + "0 [0.5651315, 0.6696122] \n", + "1 [0.64305073, 0.19809924, 0.69501674] \n", + "2 [0.8060808] \n", + "3 [0.18124953] \n", + "4 [0.075145245, 0.71639764] \n", + "... ... \n", + "53648 [-1.0, 0.7430232] \n", + "53649 [0.58085346] \n", + "53650 [0.39032233] \n", + "53651 [0.65562195, -1.0] \n", + "53652 [-1.0, -1.0] \n", + "\n", + " largeRjet_pt_syst tau_charge \n", + "0 [17716.031, 4227.1562] [] \n", + "1 [22269.969, 14174.656, 24763.469] [] \n", + "2 [0.0] [] \n", + "3 [1438.0938] [] \n", + "4 [13062.9375, 5456.125] [] \n", + "... ... ... \n", + "53648 [17955.438, 0.0] [] \n", + "53649 [21161.61] [] \n", + "53650 [7288.8906] [] \n", + "53651 [1491.4062, 8128.9688] [] \n", + "53652 [788.78125, 122.65625] [] \n", + "\n", + "[53653 rows x 81 columns]" + ] + }, + "execution_count": 6, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "df" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Vemos que tiene 81 columnas, para poder verlas todas con la finalidad de ver que variables podemos usar, usamos la siguiente lÃnea." + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": {}, + "outputs": [], + "source": [ + "pd.options.display.max_columns = 0" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": { + "scrolled": false + }, + "outputs": [ + { + "data": { + "text/html": [ + "<div>\n", + "<style scoped>\n", + " .dataframe tbody tr th:only-of-type {\n", + " vertical-align: middle;\n", + " }\n", + "\n", + " .dataframe tbody tr th {\n", + " vertical-align: top;\n", + " }\n", + "\n", + " .dataframe thead th {\n", + " text-align: right;\n", + " }\n", + "</style>\n", + "<table border=\"1\" class=\"dataframe\">\n", + " <thead>\n", + " <tr style=\"text-align: right;\">\n", + " <th></th>\n", + " <th>runNumber</th>\n", + " <th>eventNumber</th>\n", + " <th>channelNumber</th>\n", + " <th>mcWeight</th>\n", + " <th>scaleFactor_PILEUP</th>\n", + " <th>scaleFactor_ELE</th>\n", + " <th>scaleFactor_MUON</th>\n", + " <th>scaleFactor_PHOTON</th>\n", + " <th>scaleFactor_TAU</th>\n", + " <th>scaleFactor_BTAG</th>\n", + " <th>scaleFactor_LepTRIGGER</th>\n", + " <th>scaleFactor_PhotonTRIGGER</th>\n", + " <th>trigE</th>\n", + " <th>trigM</th>\n", + " <th>trigP</th>\n", + " <th>lep_n</th>\n", + " <th>lep_truthMatched</th>\n", + " <th>lep_trigMatched</th>\n", + " <th>lep_pt</th>\n", + " <th>lep_eta</th>\n", + " <th>lep_phi</th>\n", + " <th>lep_E</th>\n", + " <th>lep_z0</th>\n", + " <th>lep_charge</th>\n", + " <th>lep_type</th>\n", + " <th>lep_isTightID</th>\n", + " <th>lep_ptcone30</th>\n", + " <th>lep_etcone20</th>\n", + " <th>lep_trackd0pvunbiased</th>\n", + " <th>lep_tracksigd0pvunbiased</th>\n", + " <th>met_et</th>\n", + " <th>met_phi</th>\n", + " <th>jet_n</th>\n", + " <th>jet_pt</th>\n", + " <th>jet_eta</th>\n", + " <th>jet_phi</th>\n", + " <th>jet_E</th>\n", + " <th>jet_jvt</th>\n", + " <th>jet_trueflav</th>\n", + " <th>jet_truthMatched</th>\n", + " <th>...</th>\n", + " <th>photon_n</th>\n", + " <th>photon_truthMatched</th>\n", + " <th>photon_trigMatched</th>\n", + " <th>photon_pt</th>\n", + " <th>photon_eta</th>\n", + " <th>photon_phi</th>\n", + " <th>photon_E</th>\n", + " <th>photon_isTightID</th>\n", + " <th>photon_ptcone30</th>\n", + " <th>photon_etcone20</th>\n", + " <th>photon_convType</th>\n", + " <th>tau_n</th>\n", + " <th>tau_pt</th>\n", + " <th>tau_eta</th>\n", + " <th>tau_phi</th>\n", + " <th>tau_E</th>\n", + " <th>tau_isTightID</th>\n", + " <th>tau_truthMatched</th>\n", + " <th>tau_trigMatched</th>\n", + " <th>tau_nTracks</th>\n", + " <th>tau_BDTid</th>\n", + " <th>ditau_m</th>\n", + " <th>lep_pt_syst</th>\n", + " <th>met_et_syst</th>\n", + " <th>jet_pt_syst</th>\n", + " <th>photon_pt_syst</th>\n", + " <th>tau_pt_syst</th>\n", + " <th>XSection</th>\n", + " <th>SumWeights</th>\n", + " <th>largeRjet_n</th>\n", + " <th>largeRjet_pt</th>\n", + " <th>largeRjet_eta</th>\n", + " <th>largeRjet_phi</th>\n", + " <th>largeRjet_E</th>\n", + " <th>largeRjet_m</th>\n", + " <th>largeRjet_truthMatched</th>\n", + " <th>largeRjet_D2</th>\n", + " <th>largeRjet_tau32</th>\n", + " <th>largeRjet_pt_syst</th>\n", + " <th>tau_charge</th>\n", + " </tr>\n", + " </thead>\n", + " <tbody>\n", + " <tr>\n", + " <th>0</th>\n", + " <td>284500</td>\n", + " <td>15440</td>\n", + " <td>361106</td>\n", + " <td>1940.819946</td>\n", + " <td>1.264594</td>\n", + " <td>0.941115</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999384</td>\n", + " <td>1.01508</td>\n", + " <td>1.0</td>\n", + " <td>True</td>\n", + " <td>False</td>\n", + " <td>False</td>\n", + " <td>2</td>\n", + " <td>[True, True]</td>\n", + " <td>[True, True]</td>\n", + " <td>[169903.48, 142998.28]</td>\n", + " <td>[-1.8025961, -1.29893]</td>\n", + " <td>[3.0199232, -2.9465997]</td>\n", + " <td>[529270.06, 281577.66]</td>\n", + " <td>[-0.006664276, -0.0015563965]</td>\n", + " <td>[1, -1]</td>\n", + " <td>[11, 11]</td>\n", + " <td>[True, True]</td>\n", + " <td>[0.0, 5793.038]</td>\n", + " <td>[-664.72797, 1438.3691]</td>\n", + " <td>[0.028476434, 0.05088436]</td>\n", + " <td>[2.7847657, 1.8824033]</td>\n", + " <td>22774.525391</td>\n", + " <td>0.666735</td>\n", + " <td>1</td>\n", + " <td>[304260.75]</td>\n", + " <td>[0.2824949]</td>\n", + " <td>[-0.060711365]</td>\n", + " <td>[320657.12]</td>\n", + " <td>[0.9991386]</td>\n", + " <td>[0]</td>\n", + " <td>[False]</td>\n", + " <td>...</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0.0</td>\n", + " <td>[2525.5752, 1415.3539]</td>\n", + " <td>4430.583984</td>\n", + " <td>[7431.259]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>2</td>\n", + " <td>[336389.22, 326111.62]</td>\n", + " <td>[-1.6073604, 0.34273493]</td>\n", + " <td>[-3.1196432, -0.05282484]</td>\n", + " <td>[878163.0, 352508.03]</td>\n", + " <td>[95666.125, 70169.38]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[0.06860469, 1.7374916]</td>\n", + " <td>[0.5651315, 0.6696122]</td>\n", + " <td>[17716.031, 4227.1562]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>1</th>\n", + " <td>284500</td>\n", + " <td>81085</td>\n", + " <td>361106</td>\n", + " <td>1941.910034</td>\n", + " <td>1.320948</td>\n", + " <td>0.992704</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.997777</td>\n", + " <td>1.00936</td>\n", + " <td>1.0</td>\n", + " <td>True</td>\n", + " <td>False</td>\n", + " <td>False</td>\n", + " <td>2</td>\n", + " <td>[True, True]</td>\n", + " <td>[True, True]</td>\n", + " <td>[211820.23, 78281.89]</td>\n", + " <td>[0.047025315, -0.12980008]</td>\n", + " <td>[-3.0646539, 2.5374086]</td>\n", + " <td>[212054.48, 78942.266]</td>\n", + " <td>[-0.0404315, -0.0126543045]</td>\n", + " <td>[-1, 1]</td>\n", + " <td>[11, 11]</td>\n", + " <td>[True, True]</td>\n", + " <td>[0.0, 1394.1083]</td>\n", + " <td>[1965.3975, -560.0946]</td>\n", + " <td>[0.02091093, 0.0027772957]</td>\n", + " <td>[2.136357, 0.23608878]</td>\n", + " <td>12231.850586</td>\n", + " <td>0.864738</td>\n", + " <td>4</td>\n", + " <td>[357777.12, 162996.27, 122650.02, 87809.67]</td>\n", + " <td>[-0.38215926, 0.7312151, 0.827546, -1.0698929]</td>\n", + " <td>[0.41023782, -2.3782916, -1.9060545, 0.34156215]</td>\n", + " <td>[386681.56, 208748.97, 167330.88, 143282.44]</td>\n", + " <td>[0.991028, 0.9987799, 0.9975304, 0.9996967]</td>\n", + " <td>[0, 0, 0, 0]</td>\n", + " <td>[False, False, False, False]</td>\n", + " <td>...</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0.0</td>\n", + " <td>[2346.7507, 302.21893]</td>\n", + " <td>14149.702148</td>\n", + " <td>[5875.099, 6026.0947, 4836.232, 2209.1172]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>3</td>\n", + " <td>[424948.16, 295925.38, 288773.56]</td>\n", + " <td>[-0.50103956, 0.00353557, 0.7989666]</td>\n", + " <td>[0.39738846, 3.0375211, -2.1750822]</td>\n", + " <td>[494285.72, 314167.06, 393288.84]</td>\n", + " <td>[120339.15, 105489.414, 75616.88]</td>\n", + " <td>[0.0, 0.0, 0.0]</td>\n", + " <td>[1.4906838, 0.1404522, 0.62883687]</td>\n", + " <td>[0.64305073, 0.19809924, 0.69501674]</td>\n", + " <td>[22269.969, 14174.656, 24763.469]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>2</th>\n", + " <td>284500</td>\n", + " <td>89512</td>\n", + " <td>361106</td>\n", + " <td>1940.760010</td>\n", + " <td>1.251332</td>\n", + " <td>0.993904</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999102</td>\n", + " <td>1.00632</td>\n", + " <td>1.0</td>\n", + " <td>True</td>\n", + " <td>False</td>\n", + " <td>False</td>\n", + " <td>2</td>\n", + " <td>[True, True]</td>\n", + " <td>[True, False]</td>\n", + " <td>[234039.14, 8108.9414]</td>\n", + " <td>[0.5145523, -0.5389483]</td>\n", + " <td>[-2.36613, -0.6087628]</td>\n", + " <td>[265711.38, 9315.408]</td>\n", + " <td>[-0.072546005, 0.023675919]</td>\n", + " <td>[1, -1]</td>\n", + " <td>[11, 11]</td>\n", + " <td>[True, True]</td>\n", + " <td>[0.0, 1282.8975]</td>\n", + " <td>[-1899.3718, -219.85101]</td>\n", + " <td>[-0.01494353, 0.019966032]</td>\n", + " <td>[1.7670692, 0.88004786]</td>\n", + " <td>10936.896484</td>\n", + " <td>-2.576727</td>\n", + " <td>1</td>\n", + " <td>[239787.16]</td>\n", + " <td>[1.6001471]</td>\n", + " <td>[0.7848504]</td>\n", + " <td>[618527.5]</td>\n", + " <td>[0.9953245]</td>\n", + " <td>[0]</td>\n", + " <td>[False]</td>\n", + " <td>...</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0.0</td>\n", + " <td>[1448.3579, 86.43445]</td>\n", + " <td>7753.230469</td>\n", + " <td>[8643.487]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1</td>\n", + " <td>[255974.6]</td>\n", + " <td>[1.6731088]</td>\n", + " <td>[0.80495465]</td>\n", + " <td>[706647.8]</td>\n", + " <td>[29652.18]</td>\n", + " <td>[0.0]</td>\n", + " <td>[3.0340846]</td>\n", + " <td>[0.8060808]</td>\n", + " <td>[0.0]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>3</th>\n", + " <td>284500</td>\n", + " <td>92009</td>\n", + " <td>361106</td>\n", + " <td>1941.709961</td>\n", + " <td>1.320948</td>\n", + " <td>0.957973</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.000000</td>\n", + " <td>1.00582</td>\n", + " <td>1.0</td>\n", + " <td>True</td>\n", + " <td>False</td>\n", + " <td>False</td>\n", + " <td>2</td>\n", + " <td>[True, True]</td>\n", + " <td>[True, False]</td>\n", + " <td>[282535.56, 60798.965]</td>\n", + " <td>[-0.15060467, -0.712555]</td>\n", + " <td>[-2.322139, -1.9640917]</td>\n", + " <td>[285745.8, 76898.055]</td>\n", + " <td>[-0.007938385, 0.06525421]</td>\n", + " <td>[-1, 1]</td>\n", + " <td>[11, 11]</td>\n", + " <td>[True, False]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[475.2024, 2495.0132]</td>\n", + " <td>[0.0026975302, -0.07869638]</td>\n", + " <td>[0.4592046, 3.6803665]</td>\n", + " <td>11597.808594</td>\n", + " <td>0.950277</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>...</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0.0</td>\n", + " <td>[3379.703, 160.1774]</td>\n", + " <td>21881.703125</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1</td>\n", + " <td>[357909.62]</td>\n", + " <td>[-0.23516123]</td>\n", + " <td>[-2.2556815]</td>\n", + " <td>[377280.16]</td>\n", + " <td>[83818.125]</td>\n", + " <td>[1.0]</td>\n", + " <td>[0.13190588]</td>\n", + " <td>[0.18124953]</td>\n", + " <td>[1438.0938]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " <tr>\n", + " <th>4</th>\n", + " <td>284500</td>\n", + " <td>102088</td>\n", + " <td>361106</td>\n", + " <td>1941.239990</td>\n", + " <td>1.320948</td>\n", + " <td>0.948454</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.997393</td>\n", + " <td>1.00632</td>\n", + " <td>1.0</td>\n", + " <td>True</td>\n", + " <td>False</td>\n", + " <td>False</td>\n", + " <td>2</td>\n", + " <td>[True, True]</td>\n", + " <td>[True, True]</td>\n", + " <td>[193694.88, 135181.3]</td>\n", + " <td>[0.43973032, -0.09355457]</td>\n", + " <td>[-1.3102051, -1.2278727]</td>\n", + " <td>[212725.27, 135773.31]</td>\n", + " <td>[0.06219864, 0.013455391]</td>\n", + " <td>[-1, 1]</td>\n", + " <td>[11, 11]</td>\n", + " <td>[True, True]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[-260.75305, 333.5913]</td>\n", + " <td>[0.023009095, 0.02202719]</td>\n", + " <td>[3.0678267, 3.2605937]</td>\n", + " <td>33095.066406</td>\n", + " <td>-2.607571</td>\n", + " <td>3</td>\n", + " <td>[331520.7, 50229.3, 21196.365]</td>\n", + " <td>[-1.121277, 0.30034274, 1.7375907]</td>\n", + " <td>[1.6749763, -2.7881513, 1.2454554]</td>\n", + " <td>[563785.44, 52941.465, 62252.97]</td>\n", + " <td>[0.9958764, 0.9900437, 0.99060017]</td>\n", + " <td>[0, 0, 0]</td>\n", + " <td>[False, False, False]</td>\n", + " <td>...</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>0.0</td>\n", + " <td>[1733.7571, 575.47437]</td>\n", + " <td>330.886719</td>\n", + " <td>[4881.326, 1914.4854, 2238.7397]</td>\n", + " <td>[]</td>\n", + " <td>[]</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>2</td>\n", + " <td>[339715.34, 337192.75]</td>\n", + " <td>[0.22991252, -1.1238769]</td>\n", + " <td>[-1.2795317, 1.6654552]</td>\n", + " <td>[362021.66, 574320.1]</td>\n", + " <td>[97183.07, 30168.135]</td>\n", + " <td>[0.0, 0.0]</td>\n", + " <td>[0.042992596, 2.0432703]</td>\n", + " <td>[0.075145245, 0.71639764]</td>\n", + " <td>[13062.9375, 5456.125]</td>\n", + " <td>[]</td>\n", + " </tr>\n", + " </tbody>\n", + "</table>\n", + "<p>5 rows × 81 columns</p>\n", + "</div>" + ], + "text/plain": [ + " runNumber eventNumber ... largeRjet_pt_syst tau_charge\n", + "0 284500 15440 ... [17716.031, 4227.1562] []\n", + "1 284500 81085 ... [22269.969, 14174.656, 24763.469] []\n", + "2 284500 89512 ... [0.0] []\n", + "3 284500 92009 ... [1438.0938] []\n", + "4 284500 102088 ... [13062.9375, 5456.125] []\n", + "\n", + "[5 rows x 81 columns]" + ] + }, + "execution_count": 8, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "df.head()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Un poco de datos estadÃsticos" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "scrolled": false + }, + "outputs": [ + { + "data": { + "text/html": [ + "<div>\n", + "<style scoped>\n", + " .dataframe tbody tr th:only-of-type {\n", + " vertical-align: middle;\n", + " }\n", + "\n", + " .dataframe tbody tr th {\n", + " vertical-align: top;\n", + " }\n", + "\n", + " .dataframe thead th {\n", + " text-align: right;\n", + " }\n", + "</style>\n", + "<table border=\"1\" class=\"dataframe\">\n", + " <thead>\n", + " <tr style=\"text-align: right;\">\n", + " <th></th>\n", + " <th>runNumber</th>\n", + " <th>eventNumber</th>\n", + " <th>channelNumber</th>\n", + " <th>mcWeight</th>\n", + " <th>scaleFactor_PILEUP</th>\n", + " <th>scaleFactor_ELE</th>\n", + " <th>scaleFactor_MUON</th>\n", + " <th>scaleFactor_PHOTON</th>\n", + " <th>scaleFactor_TAU</th>\n", + " <th>scaleFactor_BTAG</th>\n", + " <th>scaleFactor_LepTRIGGER</th>\n", + " <th>scaleFactor_PhotonTRIGGER</th>\n", + " <th>lep_n</th>\n", + " <th>met_et</th>\n", + " <th>met_phi</th>\n", + " <th>jet_n</th>\n", + " <th>photon_n</th>\n", + " <th>tau_n</th>\n", + " <th>ditau_m</th>\n", + " <th>met_et_syst</th>\n", + " <th>XSection</th>\n", + " <th>SumWeights</th>\n", + " <th>largeRjet_n</th>\n", + " </tr>\n", + " </thead>\n", + " <tbody>\n", + " <tr>\n", + " <th>count</th>\n", + " <td>53653.0</td>\n", + " <td>5.365300e+04</td>\n", + " <td>53653.0</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.0</td>\n", + " <td>53653.0</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.00000</td>\n", + " <td>53653.0</td>\n", + " <td>53653.000000</td>\n", + " <td>5.365300e+04</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.0</td>\n", + " <td>53653.0</td>\n", + " <td>53653.0</td>\n", + " <td>53653.000000</td>\n", + " <td>53653.000000</td>\n", + " <td>5.365300e+04</td>\n", + " <td>53653.000000</td>\n", + " </tr>\n", + " <tr>\n", + " <th>mean</th>\n", + " <td>284500.0</td>\n", + " <td>3.972488e+07</td>\n", + " <td>361106.0</td>\n", + " <td>1906.595947</td>\n", + " <td>0.993552</td>\n", + " <td>0.961276</td>\n", + " <td>0.999988</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.000510</td>\n", + " <td>1.00101</td>\n", + " <td>1.0</td>\n", + " <td>1.727098</td>\n", + " <td>2.999400e+04</td>\n", + " <td>-0.004278</td>\n", + " <td>2.494436</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>7676.959961</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1.458912</td>\n", + " </tr>\n", + " <tr>\n", + " <th>std</th>\n", + " <td>0.0</td>\n", + " <td>2.308588e+07</td>\n", + " <td>0.0</td>\n", + " <td>365.987671</td>\n", + " <td>0.372941</td>\n", + " <td>0.026062</td>\n", + " <td>0.001164</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.029347</td>\n", + " <td>0.03176</td>\n", + " <td>0.0</td>\n", + " <td>0.455918</td>\n", + " <td>2.267055e+04</td>\n", + " <td>1.814127</td>\n", + " <td>1.376218</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>7742.734375</td>\n", + " <td>0.000000</td>\n", + " <td>0.000000e+00</td>\n", + " <td>0.511383</td>\n", + " </tr>\n", + " <tr>\n", + " <th>min</th>\n", + " <td>284500.0</td>\n", + " <td>1.010000e+03</td>\n", + " <td>361106.0</td>\n", + " <td>-1943.000000</td>\n", + " <td>0.000000</td>\n", + " <td>0.712289</td>\n", + " <td>0.825380</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.836130</td>\n", + " <td>0.00000</td>\n", + " <td>1.0</td>\n", + " <td>1.000000</td>\n", + " <td>8.181367e+01</td>\n", + " <td>-3.141528</td>\n", + " <td>0.000000</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.000000</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1.000000</td>\n", + " </tr>\n", + " <tr>\n", + " <th>25%</th>\n", + " <td>284500.0</td>\n", + " <td>1.964419e+07</td>\n", + " <td>361106.0</td>\n", + " <td>1941.060059</td>\n", + " <td>0.738662</td>\n", + " <td>0.947944</td>\n", + " <td>1.000000</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.997438</td>\n", + " <td>0.99866</td>\n", + " <td>1.0</td>\n", + " <td>1.000000</td>\n", + " <td>1.572494e+04</td>\n", + " <td>-1.579283</td>\n", + " <td>1.000000</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>2188.400391</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1.000000</td>\n", + " </tr>\n", + " <tr>\n", + " <th>50%</th>\n", + " <td>284500.0</td>\n", + " <td>3.940982e+07</td>\n", + " <td>361106.0</td>\n", + " <td>1941.400024</td>\n", + " <td>1.196105</td>\n", + " <td>0.963610</td>\n", + " <td>1.000000</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.998668</td>\n", + " <td>1.00306</td>\n", + " <td>1.0</td>\n", + " <td>2.000000</td>\n", + " <td>2.586413e+04</td>\n", + " <td>-0.010504</td>\n", + " <td>2.000000</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>5281.721680</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>1.000000</td>\n", + " </tr>\n", + " <tr>\n", + " <th>75%</th>\n", + " <td>284500.0</td>\n", + " <td>5.977106e+07</td>\n", + " <td>361106.0</td>\n", + " <td>1941.739990</td>\n", + " <td>1.251332</td>\n", + " <td>0.978670</td>\n", + " <td>1.000000</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>0.999264</td>\n", + " <td>1.00632</td>\n", + " <td>1.0</td>\n", + " <td>2.000000</td>\n", + " <td>3.917682e+04</td>\n", + " <td>1.572061</td>\n", + " <td>3.000000</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>10652.117188</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>2.000000</td>\n", + " </tr>\n", + " <tr>\n", + " <th>max</th>\n", + " <td>284500.0</td>\n", + " <td>7.999848e+07</td>\n", + " <td>361106.0</td>\n", + " <td>1943.439941</td>\n", + " <td>3.348151</td>\n", + " <td>1.112705</td>\n", + " <td>1.026367</td>\n", + " <td>1.0</td>\n", + " <td>1.0</td>\n", + " <td>1.752337</td>\n", + " <td>1.21491</td>\n", + " <td>1.0</td>\n", + " <td>4.000000</td>\n", + " <td>1.986488e+06</td>\n", + " <td>3.141495</td>\n", + " <td>11.000000</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>0.0</td>\n", + " <td>91205.000000</td>\n", + " <td>1950.529541</td>\n", + " <td>1.502776e+11</td>\n", + " <td>4.000000</td>\n", + " </tr>\n", + " </tbody>\n", + "</table>\n", + "</div>" + ], + "text/plain": [ + " runNumber eventNumber ... SumWeights largeRjet_n\n", + "count 53653.0 5.365300e+04 ... 5.365300e+04 53653.000000\n", + "mean 284500.0 3.972488e+07 ... 1.502776e+11 1.458912\n", + "std 0.0 2.308588e+07 ... 0.000000e+00 0.511383\n", + "min 284500.0 1.010000e+03 ... 1.502776e+11 1.000000\n", + "25% 284500.0 1.964419e+07 ... 1.502776e+11 1.000000\n", + "50% 284500.0 3.940982e+07 ... 1.502776e+11 1.000000\n", + "75% 284500.0 5.977106e+07 ... 1.502776e+11 2.000000\n", + "max 284500.0 7.999848e+07 ... 1.502776e+11 4.000000\n", + "\n", + "[8 rows x 23 columns]" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "df.describe()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Iniciamos el análisis" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": {}, + "outputs": [], + "source": [ + "canvas = ROOT.TCanvas(\"Canvas\",\"c\",800,600)" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "53653" + ] + }, + "execution_count": 11, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "tree = file2.Get(\"mini\")\n", + "tree.GetEntries()" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "metadata": {}, + "outputs": [], + "source": [ + "hist = ROOT.TH1F(\"variable\",\"Mass of the Z boson; mass [GeV]; events\",30,40,140)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## H→Zγ→γ ll*" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "> Prestar atención al Cut número 4, porque se debe al artÃculo de la búsqueda del decaimiento de Dalitz observado en ATLAS (http://cdsweb.cern.ch/record/2750822), es uno de las restricciones en la búsqueda de " + ] + }, + { + "cell_type": "code", + "execution_count": 13, + "metadata": {}, + "outputs": [], + "source": [ + "leadLepton = ROOT.TLorentzVector()\n", + "trailLepton = ROOT.TLorentzVector()\n", + "\n", + "for event in tree:\n", + " \n", + " # Cut #1: At least 2 leptons\n", + " if tree.lep_n >= 2:\n", + " \n", + " # Cut #2: Leptons with opposite charge\n", + " if (tree.lep_charge[0] != tree.lep_charge[1]):\n", + " \n", + " # Cut #3: Leptons of the same family (2 electrons or 2 muons)\n", + " if (tree.lep_type[0] == tree.lep_type[1]):\n", + " \n", + " # Cut #4: E < 30 GeV\n", + " if tree.lep_E[0] < 30000:\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " # Let's define one TLorentz vector for each, e.i. two vectors!\n", + " leadLepton.SetPtEtaPhiE(tree.lep_pt[0]/1000., tree.lep_eta[0], tree.lep_phi[0], tree.lep_E[0]/1000.)\n", + " trailLepton.SetPtEtaPhiE(tree.lep_pt[1]/1000., tree.lep_eta[1], tree.lep_phi[1], tree.lep_E[1]/1000.)\n", + " # Next line: addition of two TLorentz vectors above --> ask mass very easy (devide by 1000 to get value in GeV)\n", + " invmass = leadLepton + trailLepton\n", + "\n", + " hist.Fill(invmass.M())" + ] + }, + { + "cell_type": "code", + "execution_count": 14, + "metadata": {}, + "outputs": [], + "source": [ + "hist.Draw()\n", + "hist.SetFillColor(3)" + ] + }, + { + "cell_type": "code", + "execution_count": 15, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "\n", + "<div id=\"root_plot_1\"\n", + " style=\"width: 796px; height: 572px\">\n", + "</div>\n", + "\n", + "<script>\n", + " requirejs.config({\n", + " paths: {\n", + " 'JSRootCore' : 'https://root.cern.ch/js/notebook//scripts/JSRootCore',\n", + " }\n", + " });\n", + " require(['JSRootCore'],\n", + " function(Core) {\n", + " var obj = Core.JSONR_unref({\"_typename\":\"TCanvas\",\"fUniqueID\":0,\"fBits\":53673992,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fLeftMargin\":0.1,\"fRightMargin\":0.1,\"fBottomMargin\":0.1,\"fTopMargin\":0.1,\"fXfile\":2,\"fYfile\":2,\"fAfile\":1,\"fXstat\":0.99,\"fYstat\":0.99,\"fAstat\":2,\"fFrameFillColor\":0,\"fFrameLineColor\":1,\"fFrameFillStyle\":1001,\"fFrameLineStyle\":1,\"fFrameLineWidth\":1,\"fFrameBorderSize\":1,\"fFrameBorderMode\":0,\"fX1\":27.4999990686774,\"fY1\":-0.262500019557775,\"fX2\":152.500000931323,\"fY2\":2.36250001955778,\"fXtoAbsPixelk\":-175.119941459846,\"fXtoPixelk\":-175.119941459846,\"fXtoPixel\":6.3679999051094,\"fYtoAbsPixelk\":514.800046590614,\"fYtoPixelk\":514.800046590614,\"fYtoPixel\":-217.904758657728,\"fUtoAbsPixelk\":5e-5,\"fUtoPixelk\":5e-5,\"fUtoPixel\":796,\"fVtoAbsPixelk\":572.00005,\"fVtoPixelk\":572,\"fVtoPixel\":-572,\"fAbsPixeltoXk\":27.4999990686774,\"fPixeltoXk\":27.4999990686774,\"fPixeltoX\":0.157035178219404,\"fAbsPixeltoYk\":2.36250001955778,\"fPixeltoYk\":-0.262500019557775,\"fPixeltoY\":-0.00458916090754467,\"fXlowNDC\":0,\"fYlowNDC\":0,\"fXUpNDC\":0,\"fYUpNDC\":0,\"fWNDC\":1,\"fHNDC\":1,\"fAbsXlowNDC\":0,\"fAbsYlowNDC\":0,\"fAbsWNDC\":1,\"fAbsHNDC\":1,\"fUxmin\":40,\"fUymin\":0,\"fUxmax\":140,\"fUymax\":2.1,\"fTheta\":30,\"fPhi\":30,\"fAspectRatio\":0,\"fNumber\":0,\"fTickx\":0,\"fTicky\":0,\"fLogx\":0,\"fLogy\":0,\"fLogz\":0,\"fPadPaint\":0,\"fCrosshair\":0,\"fCrosshairPos\":0,\"fBorderSize\":2,\"fBorderMode\":0,\"fModified\":false,\"fGridx\":false,\"fGridy\":false,\"fAbsCoord\":false,\"fEditable\":true,\"fFixedAspectRatio\":false,\"fPrimitives\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TFrame\",\"fUniqueID\":0,\"fBits\":50331656,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fX1\":40,\"fY1\":0,\"fX2\":140,\"fY2\":2.1,\"fBorderSize\":1,\"fBorderMode\":0},{\"_typename\":\"TH1F\",\"fUniqueID\":0,\"fBits\":50331656,\"fName\":\"variable\",\"fTitle\":\"Mass of the Z boson\",\"fLineColor\":602,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":3,\"fFillStyle\":1001,\"fMarkerColor\":1,\"fMarkerStyle\":1,\"fMarkerSize\":1,\"fNcells\":32,\"fXaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"xaxis\",\"fTitle\":\" mass [GeV]\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":30,\"fXmin\":40,\"fXmax\":140,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fYaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"yaxis\",\"fTitle\":\" events\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":0,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fZaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"zaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fBarOffset\":0,\"fBarWidth\":1000,\"fEntries\":5,\"fTsumw\":5,\"fTsumw2\":5,\"fTsumwx\":460.606963185387,\"fTsumwx2\":42465.3851897774,\"fMaximum\":-1111,\"fMinimum\":-1111,\"fNormFactor\":0,\"fContour\":[],\"fSumw2\":[],\"fOption\":\"\",\"fFunctions\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TPaveStats\",\"fUniqueID\":0,\"fBits\":50331657,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fX1\":125.000002533197,\"fY1\":1.771875026403,\"fX2\":150.000003278256,\"fY2\":2.19187502327375,\"fX1NDC\":0.780000016093254,\"fY1NDC\":0.775000005960464,\"fX2NDC\":0.980000019073486,\"fY2NDC\":0.935000002384186,\"fBorderSize\":1,\"fInit\":1,\"fShadowColor\":1,\"fCornerRadius\":0,\"fOption\":\"brNDC\",\"fName\":\"stats\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":12,\"fTextColor\":1,\"fTextFont\":42,\"fLabel\":\"\",\"fLongest\":18,\"fMargin\":0.05,\"fLines\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"variable\",\"fTextAngle\":0,\"fTextSize\":0.0368,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.0368000008165836},{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Entries = 5 \",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.04},{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Mean = 92.12\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.04},{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Std Dev = 2.593\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.04}],\"opt\":[\"\",\"\",\"\",\"\"]},\"fOptFit\":0,\"fOptStat\":1111,\"fFitFormat\":\"5.4g\",\"fStatFormat\":\"6.4g\",\"fParent\":{\"$ref\":3}}],\"opt\":[\"brNDC\"]},\"fBufferSize\":0,\"fBuffer\":[],\"fBinStatErrOpt\":0,\"fStatOverflows\":2,\"fArray\":[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,2,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0]},{\"_typename\":\"TPaveText\",\"fUniqueID\":0,\"fBits\":50331657,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":0,\"fX1\":66.8435926197668,\"fY1\":2.20500002757646,\"fX2\":113.156407380233,\"fY2\":2.34937503187917,\"fX1NDC\":0.314748743718593,\"fY1NDC\":0.940000003948808,\"fX2NDC\":0.685251256281407,\"fY2NDC\":0.995000004768372,\"fBorderSize\":0,\"fInit\":1,\"fShadowColor\":1,\"fCornerRadius\":0,\"fOption\":\"blNDC\",\"fName\":\"title\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":22,\"fTextColor\":1,\"fTextFont\":42,\"fLabel\":\"\",\"fLongest\":19,\"fMargin\":0.05,\"fLines\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Mass of the Z boson\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.0467500016093254}],\"opt\":[\"\"]}}],\"opt\":[\"\",\"\",\"blNDC\"]},\"fExecs\":null,\"fName\":\"Canvas\",\"fTitle\":\"c\",\"fNumPaletteColor\":0,\"fNextPaletteColor\":0,\"fDISPLAY\":\"$DISPLAY\",\"fDoubleBuffer\":0,\"fRetained\":true,\"fXsizeUser\":0,\"fYsizeUser\":0,\"fXsizeReal\":20,\"fYsizeReal\":15,\"fWindowTopX\":0,\"fWindowTopY\":0,\"fWindowWidth\":0,\"fWindowHeight\":0,\"fCw\":796,\"fCh\":572,\"fCatt\":{\"_typename\":\"TAttCanvas\",\"fXBetween\":2,\"fYBetween\":2,\"fTitleFromTop\":1.2,\"fXdate\":0.2,\"fYdate\":0.3,\"fAdate\":1},\"kMoveOpaque\":true,\"kResizeOpaque\":true,\"fHighLightColor\":2,\"fBatch\":true,\"kShowEventStatus\":false,\"kAutoExec\":true,\"kMenuBar\":true});\n", + " Core.draw(\"root_plot_1\", obj, \"\");\n", + " }\n", + " );\n", + "</script>\n" + ], + "text/plain": [ + "<IPython.core.display.HTML object>" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "canvas.Draw()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Missing transverse energy de los leptones en estado final" + ] + }, + { + "cell_type": "code", + "execution_count": 16, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "Warning in <TCanvas::Constructor>: Deleting canvas with same name: Canvas\n" + ] + } + ], + "source": [ + "canvas = ROOT.TCanvas(\"Canvas\",\"cz\",800,600)" + ] + }, + { + "cell_type": "code", + "execution_count": 17, + "metadata": {}, + "outputs": [], + "source": [ + "hist = ROOT.TH1F(\"X\",\"Y [GeV] ; events\",30,105,160)" + ] + }, + { + "cell_type": "code", + "execution_count": 18, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "53653" + ] + }, + "execution_count": 18, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "tree = file2.Get(\"mini\")\n", + "tree.GetEntries()" + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "metadata": {}, + "outputs": [], + "source": [ + "c = ROOT.TCanvas(\"testCanvas\",\"a first way to plot a variable\",800,600)" + ] + }, + { + "cell_type": "code", + "execution_count": 34, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "Warning in <TROOT::Append>: Replacing existing TH1: h_sig (Potential memory leak).\n" + ] + } + ], + "source": [ + "h_sig = ROOT.TH1F(\"h_sig\",\"Missing transverse energy\",20,0,200)" + ] + }, + { + "cell_type": "code", + "execution_count": 35, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "10000\n", + "20000\n", + "30000\n", + "40000\n", + "50000\n" + ] + } + ], + "source": [ + "m=0 \n", + "for event in tree:\n", + " m += 1\n", + " ## printing the evolution in number of events\n", + " if(m%10000==0):\n", + " print(m)\n", + " h_sig.Fill((tree.met_et)/1000.)" + ] + }, + { + "cell_type": "code", + "execution_count": 36, + "metadata": {}, + "outputs": [], + "source": [ + "scale_sig = h_sig.Integral()\n", + "h_sig.Scale(1/scale_sig)" + ] + }, + { + "cell_type": "code", + "execution_count": 37, + "metadata": {}, + "outputs": [], + "source": [ + "h_sig.SetFillStyle(3003)\n", + "h_sig.SetFillColor(2)\n", + "h_sig.SetLineColor(2)" + ] + }, + { + "cell_type": "code", + "execution_count": 38, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "<ROOT.TLegendEntry object (\"TLegendEntry\") at 0x15c191c0>" + ] + }, + "execution_count": 38, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "legend=ROOT.TLegend(0.5,0.7,0.9,0.9)\n", + "legend.AddEntry(h_sig,\"Signal (Z #rightarrow ll*)\",\"l\")" + ] + }, + { + "cell_type": "code", + "execution_count": 39, + "metadata": {}, + "outputs": [], + "source": [ + "h_sig.SetStats(0)" + ] + }, + { + "cell_type": "code", + "execution_count": 40, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "\n", + "<div id=\"root_plot_4\"\n", + " style=\"width: 796px; height: 572px\">\n", + "</div>\n", + "\n", + "<script>\n", + " requirejs.config({\n", + " paths: {\n", + " 'JSRootCore' : 'https://root.cern.ch/js/notebook//scripts/JSRootCore',\n", + " }\n", + " });\n", + " require(['JSRootCore'],\n", + " function(Core) {\n", + " var obj = Core.JSONR_unref({\"_typename\":\"TCanvas\",\"fUniqueID\":0,\"fBits\":53673992,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fLeftMargin\":0.1,\"fRightMargin\":0.1,\"fBottomMargin\":0.1,\"fTopMargin\":0.1,\"fXfile\":2,\"fYfile\":2,\"fAfile\":1,\"fXstat\":0.99,\"fYstat\":0.99,\"fAstat\":2,\"fFrameFillColor\":0,\"fFrameLineColor\":1,\"fFrameFillStyle\":1001,\"fFrameLineStyle\":1,\"fFrameLineWidth\":1,\"fFrameBorderSize\":1,\"fFrameBorderMode\":0,\"fX1\":-25.0000018626453,\"fY1\":-0.0313701751731466,\"fX2\":225.000001862645,\"fY2\":0.282331557860238,\"fXtoAbsPixelk\":79.6000547445299,\"fXtoPixelk\":79.6000547445299,\"fXtoPixel\":3.1839999525547,\"fYtoAbsPixelk\":514.800046590614,\"fYtoPixelk\":514.800046590614,\"fYtoPixel\":-1823.38807780551,\"fUtoAbsPixelk\":5e-5,\"fUtoPixelk\":5e-5,\"fUtoPixel\":796,\"fVtoAbsPixelk\":572.00005,\"fVtoPixelk\":572,\"fVtoPixel\":-572,\"fAbsPixeltoXk\":-25.0000018626453,\"fPixeltoXk\":-25.0000018626453,\"fPixeltoX\":0.314070356438807,\"fAbsPixeltoYk\":0.282331557860238,\"fPixeltoYk\":-0.0313701751731466,\"fPixeltoY\":-5.48429603205219e-4,\"fXlowNDC\":0,\"fYlowNDC\":0,\"fXUpNDC\":0,\"fYUpNDC\":0,\"fWNDC\":1,\"fHNDC\":1,\"fAbsXlowNDC\":0,\"fAbsYlowNDC\":0,\"fAbsWNDC\":1,\"fAbsHNDC\":1,\"fUxmin\":0,\"fUymin\":0,\"fUxmax\":200,\"fUymax\":0.250961382687092,\"fTheta\":30,\"fPhi\":30,\"fAspectRatio\":0,\"fNumber\":0,\"fTickx\":0,\"fTicky\":0,\"fLogx\":0,\"fLogy\":0,\"fLogz\":0,\"fPadPaint\":0,\"fCrosshair\":0,\"fCrosshairPos\":0,\"fBorderSize\":2,\"fBorderMode\":0,\"fModified\":false,\"fGridx\":false,\"fGridy\":false,\"fAbsCoord\":false,\"fEditable\":true,\"fFixedAspectRatio\":false,\"fPrimitives\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TFrame\",\"fUniqueID\":0,\"fBits\":50331656,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fX1\":0,\"fY1\":0,\"fX2\":200,\"fY2\":0.250961382687092,\"fBorderSize\":1,\"fBorderMode\":0},{\"_typename\":\"TH1F\",\"fUniqueID\":0,\"fBits\":50332168,\"fName\":\"h_sig\",\"fTitle\":\"Missing transverse energy\",\"fLineColor\":2,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":2,\"fFillStyle\":3003,\"fMarkerColor\":1,\"fMarkerStyle\":1,\"fMarkerSize\":1,\"fNcells\":22,\"fXaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"xaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":20,\"fXmin\":0,\"fXmax\":200,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fYaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"yaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":0,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fZaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"zaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fBarOffset\":0,\"fBarWidth\":1000,\"fEntries\":53653,\"fTsumw\":53621,\"fTsumw2\":53621,\"fTsumwx\":1598245.92758354,\"fTsumwx2\":68461690.7118449,\"fMaximum\":-1111,\"fMinimum\":-1111,\"fNormFactor\":0,\"fContour\":[],\"fSumw2\":[0,2.17340522377219e-6,4.45741100141852e-6,4.39028550802949e-6,3.1802876246076e-6,1.92890148360386e-6,1.1386987842263e-6,6.21867265179176e-7,3.33888464525732e-7,1.81551852585867e-7,9.77319359705528e-8,5.63436783887173e-8,3.23454450009303e-8,1.87812261295724e-8,1.2173016935834e-8,8.3472116131433e-9,5.9126082259765e-9,4.52140629045262e-9,3.82580532269068e-9,1.73900241940485e-9,1.39120193552388e-9,1.11296154841911e-8],\"fOption\":\"\",\"fFunctions\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[],\"opt\":[]},\"fBufferSize\":0,\"fBuffer\":[],\"fBinStatErrOpt\":0,\"fStatOverflows\":2,\"fArray\":[0,0.1165402,0.2390108,0.2354115,0.1705302,0.1034296,0.06105817,0.03334514,0.01790343,0.009734992,0.005240484,0.003021204,0.001734395,0.001007068,6.527294e-4,4.475858e-4,3.1704e-4,2.424423e-4,2.051435e-4,9.324705e-5,7.459764e-5,5.967811e-4]},{\"_typename\":\"TLegend\",\"fUniqueID\":0,\"fBits\":50331656,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fX1\":100,\"fY1\":0.188221037950223,\"fX2\":200.000001490116,\"fY2\":0.2509613845569,\"fX1NDC\":0.5,\"fY1NDC\":0.7,\"fX2NDC\":0.9,\"fY2NDC\":0.9,\"fBorderSize\":1,\"fInit\":1,\"fShadowColor\":1,\"fCornerRadius\":0,\"fOption\":\"brNDC\",\"fName\":\"TPave\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":12,\"fTextColor\":1,\"fTextFont\":42,\"fPrimitives\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TLegendEntry\",\"fUniqueID\":0,\"fBits\":50331656,\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":42,\"fLineColor\":2,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":0,\"fMarkerColor\":1,\"fMarkerStyle\":21,\"fMarkerSize\":1,\"fObject\":{\"$ref\":3},\"fLabel\":\"Signal (Z #rightarrow ll*)\",\"fOption\":\"l\"}],\"opt\":[\"l\"]},\"fEntrySeparation\":0.1,\"fMargin\":0.25,\"fNColumns\":1,\"fColumnSeparation\":0},{\"_typename\":\"TPaveText\",\"fUniqueID\":0,\"fBits\":50331657,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":0,\"fX1\":40.0251247344454,\"fY1\":0.261677700499376,\"fX2\":159.974875265555,\"fY2\":0.280763050690918,\"fX1NDC\":0.260100502512563,\"fY1NDC\":0.934160843929211,\"fX2NDC\":0.739899497487437,\"fY2NDC\":0.995000004768372,\"fBorderSize\":0,\"fInit\":1,\"fShadowColor\":1,\"fCornerRadius\":0,\"fOption\":\"blNDC\",\"fName\":\"title\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":22,\"fTextColor\":1,\"fTextFont\":42,\"fLabel\":\"\",\"fLongest\":25,\"fMargin\":0.05,\"fLines\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Missing transverse energy\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.0517132878303528}],\"opt\":[\"\"]}}],\"opt\":[\"\",\"hist\",\"\",\"blNDC\"]},\"fExecs\":null,\"fName\":\"testCanvas\",\"fTitle\":\"a first way to plot a variable\",\"fNumPaletteColor\":0,\"fNextPaletteColor\":0,\"fDISPLAY\":\"$DISPLAY\",\"fDoubleBuffer\":0,\"fRetained\":true,\"fXsizeUser\":0,\"fYsizeUser\":0,\"fXsizeReal\":20,\"fYsizeReal\":15,\"fWindowTopX\":0,\"fWindowTopY\":0,\"fWindowWidth\":0,\"fWindowHeight\":0,\"fCw\":796,\"fCh\":572,\"fCatt\":{\"_typename\":\"TAttCanvas\",\"fXBetween\":2,\"fYBetween\":2,\"fTitleFromTop\":1.2,\"fXdate\":0.2,\"fYdate\":0.3,\"fAdate\":1},\"kMoveOpaque\":true,\"kResizeOpaque\":true,\"fHighLightColor\":2,\"fBatch\":true,\"kShowEventStatus\":false,\"kAutoExec\":true,\"kMenuBar\":true});\n", + " Core.draw(\"root_plot_4\", obj, \"\");\n", + " }\n", + " );\n", + "</script>\n" + ], + "text/plain": [ + "<IPython.core.display.HTML object>" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "h_sig.Draw(\"hist\")\n", + "legend.Draw()\n", + "c.Draw()" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3", + "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.6.9" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +} diff --git a/data-used/mc_361106.Zee.1largeRjet1lep.root b/data-used/mc_361106.Zee.1largeRjet1lep.root new file mode 100644 index 0000000000000000000000000000000000000000..9a07662b61beb8760bbde9562dc5346c0d5e66b8 Binary files /dev/null and b/data-used/mc_361106.Zee.1largeRjet1lep.root differ diff --git a/reporte/reporte.pdf b/reporte/reporte.pdf new file mode 100644 index 0000000000000000000000000000000000000000..fddc5c15398e392cd50b5bdd2e896be2adce6f37 Binary files /dev/null and b/reporte/reporte.pdf differ