diff --git a/data/xsecs.json b/data/xsecs.json index 5e6e9f5c..137a4a20 100644 --- a/data/xsecs.json +++ b/data/xsecs.json @@ -22,19 +22,19 @@ "QCD_Pt_3200toInf": 0.0001353, "ST_s-channel_4f_leptonDecays": 1.188915, - "ST_t-channel_antitop_4f_InclusiveDecays": 67.93, - "ST_t-channel_top_4f_InclusiveDecays": 113.4, + "ST_t-channel_antitop_4f_InclusiveDecays": 80.0, + "ST_t-channel_top_4f_InclusiveDecays": 134.2, "ST_t-channel_antitop_5f_InclusiveDecays": 71.74, "ST_t-channel_top_5f_InclusiveDecays": 119.7, - "ST_tW_antitop_5f_inclusiveDecays": 32.51, + "ST_tW_antitop_5f_inclusiveDecays": 39.65, "ST_tW_antitop_5f_NoFullyHadronicDecays": 10.890849999999999, - "ST_tW_top_5f_inclusiveDecays": 32.45, + "ST_tW_top_5f_inclusiveDecays": 39.65, "ST_tW_top_5f_NoFullyHadronicDecays": 10.87075, - "TTTo2L2Nu": 77.10979749999998, - "TTToHadronic": 303.8527975, - "TTToSemiLeptonic": 306.137405, - "TTToSemiLeptonic_ext1": 306.137405, + "TTTo2L2Nu": "670.3*1.24088*0.333*0.333", + "TTToHadronic": "670.3*1.24088*0.667*0.667", + "TTToSemiLeptonic": "670.3*1.24088*2*(0.667*0.333)", + "TTToSemiLeptonic_ext1": "670.3*1.24088*2*(0.667*0.333)", "TT_MTT1000ToInf": 76.61, "TT_MTT700To1000": 20.58, @@ -57,9 +57,9 @@ "ZJetsToQQ_HT-600to800": 25.41, "ZJetsToQQ_HT-800toInf": 12.91, - "WW": 75.83, - "WZ": 27.56, - "ZZ": 12.14, + "WW": 118.7, + "WZ": "28.55 + 18.19", + "ZZ": 16.91, "GluGluToHHTobbVV_node_cHHH0": "0.069725 * 5.824e-01 * (0.2154 * 0.676 ** 2 + 0.02643 * 0.692 ** 2) * 2", "GluGluToHHTobbVV_node_cHHH1": "31.05e-3 * 5.824e-01 * (0.2154 * 0.676 ** 2 + 0.02643 * 0.692 ** 2) * 2", @@ -94,5 +94,4 @@ "ZH_HToBB_ZToQQ": 0.326, "ggZH_HToBB_ZToQQ": 0.025, "ttHToBB": 0.2912 - } diff --git a/paper/latex_tables.ipynb b/paper/latex_tables.ipynb index 69e80d4b..a1b6b8f8 100644 --- a/paper/latex_tables.ipynb +++ b/paper/latex_tables.ipynb @@ -111,7 +111,7 @@ }, { "cell_type": "code", - "execution_count": 3, + "execution_count": 24, "metadata": {}, "outputs": [], "source": [ @@ -119,12 +119,46 @@ " if isinstance(xsec, str):\n", " xsec = eval(xsec)\n", "\n", + " xsec = str(np.round(xsec, 3))\n", + "\n", + " # remove trailing zeros\n", + " for i in range(len(xsec)):\n", + " if xsec[-1] == \"0\":\n", + " xsec = xsec[:-1]\n", + " else:\n", + " break\n", + "\n", + " # remove trailing \".\"\n", + " if xsec[-1] == \".\":\n", + " xsec = xsec[:-1]\n", + "\n", " return xsec" ] }, { "cell_type": "code", - "execution_count": 4, + "execution_count": 12, + "metadata": {}, + "outputs": [], + "source": [ + "# process or sample and corresponding citation in paper\n", + "refs = {\n", + " \"WW\": \"Gehrmann:2014fva\",\n", + " \"WZ\": \"Campbell:2011bn\",\n", + " \"ZZ\": \"Cascioli:2014yka\",\n", + " \"SingleTop\": \"SingleTopNNLORef\",\n", + " \"TTbar\": \"TtbarNNLO\",\n", + " \"HH\": \"LHCHiggsHH\",\n", + " # \"HWW (inclusive)\": \"YR4\",\n", + " # \"Hbb (inclusive)\": \"YR4\",\n", + " # \"QCD\": \"xsdb\",\n", + " # \"V+Jets\": \"xsdb\"\n", + "}" + ] + }, + { + "cell_type": "code", + "execution_count": 25, "metadata": {}, "outputs": [], "source": [ @@ -137,12 +171,17 @@ "lines = []\n", "\n", "for sample in np.sort(list(samples.keys())):\n", + " sample_text = sample\n", + " if sample in refs:\n", + " sample_text += rf\"~\\cite{{{refs[sample]}}}\"\n", + "\n", " for key, dataset in samples[sample][\"datasets\"].items():\n", " dname = dataset.split(\"/\")[1].replace(\"_\", \"\\_\")\n", - " lines.append(rf\" & {dname} & {process_xsec(xsecs[key]):.2f} \\\\\" + \"\\n\")\n", + " citation = \"\" if not key in refs else rf\"~\\cite{{{refs[key]}}}\"\n", + " lines.append(rf\" & {dname} & {process_xsec(xsecs[key])}{citation} \\\\\" + \"\\n\")\n", "\n", " lend = len(samples[sample][\"datasets\"])\n", - " lines[-lend] = rf\"\\multirow{{{lend}}}{{*}}{{{sample}}}\" + lines[-lend]\n", + " lines[-lend] = rf\"\\multirow{{{lend}}}{{*}}{{{sample_text}}}\" + lines[-lend]\n", " lines.append(r\"\\hline\" + \"\\n\")\n", "\n", "# remove trailing \"\\\\\" and \"\\hline\"\n", @@ -155,7 +194,7 @@ }, { "cell_type": "code", - "execution_count": 12, + "execution_count": 11, "metadata": {}, "outputs": [], "source": [ @@ -169,13 +208,24 @@ "\n", "for sample in np.sort(list(samples.keys())):\n", " for key, dataset in samples[sample][\"datasets\"].items():\n", - " key = key.split(\"_pn4q\")[0]\n", - " cl = key.split(\"cHHH\")[1].split(\"_\")[0].replace(\"p\", \".\")\n", - " sm = \"SM\" if cl == \"1\" else \"BSM\"\n", - " pname = rf\"ggHH {sm} ($\\kappa_{{\\lambda}} = {cl}$)\"\n", + " if key.startswith(\"GluGluToHHTo\"):\n", + " key = key.split(\"_pn4q\")[0]\n", + " cl = key.split(\"cHHH\")[1].split(\"_\")[0].replace(\"p\", \".\")\n", + " sm = \"SM\" if cl == \"1\" else \"BSM\"\n", + " pname = rf\"ggHH {sm} ($\\kappa_{{\\lambda}} = {cl}$)\"\n", + " else:\n", + " cv = key.split(\"CV_\")[1].split(\"_C2V\")[0].replace(\"_\", \".\")\n", + " c2v = key.split(\"C2V_\")[1].split(\"_C3\")[0].replace(\"_\", \".\")\n", + " c3 = key.split(\"C3_\")[1].replace(\"_\", \".\")\n", + " sm = \"SM\" if cl == \"1\" and c2v == \"1\" and c3 == \"1\" else \"BSM\"\n", + " pname = rf\"VBFHH {sm} ($\\kappa_\\mathrm{{V}} = {cv}$, $\\kappa_\\mathrm{{2V}} = {c2v}$, $\\kappa_{{\\lambda}} = {c3}$)\"\n", + "\n", " dname = dataset.split(\"/\")[1].replace(\"_\", \"\\_\")\n", " lines.append(rf\"{pname} & {dname} & {process_xsec(xsecs[key]) * 1000:.2f} \\\\\" + \"\\n\")\n", "\n", + "# remove trailing \"\\\\\"\n", + "lines[-1] = lines[-1][:-4]\n", + "\n", "with open(\"tables/datasets_mcsig_nonres.tex\", \"w\") as f:\n", " f.writelines(lines)" ] diff --git a/paper/samples/MC_bg.yaml b/paper/samples/MC_bg.yaml index 9fd26869..da708c99 100644 --- a/paper/samples/MC_bg.yaml +++ b/paper/samples/MC_bg.yaml @@ -23,9 +23,11 @@ samples: V+Jets: datasets: + WJetsToQQ_HT-200to400: /WJetsToQQ_HT-200to400_TuneCP5_13TeV-madgraphMLM-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM WJetsToQQ_HT-400to600: /WJetsToQQ_HT-400to600_TuneCP5_13TeV-madgraphMLM-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM WJetsToQQ_HT-600to800: /WJetsToQQ_HT-600to800_TuneCP5_13TeV-madgraphMLM-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM WJetsToQQ_HT-800toInf: /WJetsToQQ_HT-800toInf_TuneCP5_13TeV-madgraphMLM-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM + ZJetsToQQ_HT-200to400: /ZJetsToQQ_HT-200to600_TuneCP5_13TeV-madgraphMLM-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM ZJetsToQQ_HT-400to600: /ZJetsToQQ_HT-400to600_TuneCP5_13TeV-madgraphMLM-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM ZJetsToQQ_HT-600to800: /ZJetsToQQ_HT-600to800_TuneCP5_13TeV-madgraphMLM-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM ZJetsToQQ_HT-800toInf: /ZJetsToQQ_HT-800toInf_TuneCP5_13TeV-madgraphMLM-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM @@ -70,17 +72,4 @@ samples: HH: datasets: - VBF_HHTobbVV_CV_1_C2V_1_C3_1: /VBF_HHTobbVV_CV_1_C2V_1_C3_1_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM - GluGluToHHTo4B_node_cHHH1_preUL: /GluGluToHHTo4B_node_cHHH1_TuneCP5_PSWeights_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM - - # GluGluToHHTobbVV_node_cHHH1_pn4q: /GluGluToHHTobbVV_node_cHHH1_TuneCP5_13TeV-powheg-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v1/MINIAODSIM - # GluGluToHHTobbVV_node_cHHH2p45_pn4q: /GluGluToHHTobbVV_node_cHHH2p45_TuneCP5_13TeV-powheg-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v1/MINIAODSIM - # GluGluToHHTobbVV_node_cHHH0_pn4q: /GluGluToHHTobbVV_node_cHHH0_TuneCP5_13TeV-powheg-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v1/MINIAODSIM - # GluGluToHHTobbVV_node_cHHH5_pn4q: /GluGluToHHTobbVV_node_cHHH5_TuneCP5_13TeV-powheg-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v1/MINIAODSIM - # - # VBF_HHTobbVV_CV_1_C2V_2_C3_1: /VBF_HHTobbVV_CV_1_C2V_2_C3_1_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM - # VBF_HHTobbVV_CV_1_C2V_1_C3_2: /VBF_HHTobbVV_CV_1_C2V_1_C3_2_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM - # VBF_HHTobbVV_CV_1_C2V_1_C3_0: /VBF_HHTobbVV_CV_1_C2V_1_C3_0_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM - # VBF_HHTobbVV_CV_1_C2V_0_C3_1: /VBF_HHTobbVV_CV_1_C2V_0_C3_1_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM - # VBF_HHTobbVV_CV_1_5_C2V_1_C3_1: /VBF_HHTobbVV_CV_1_5_C2V_1_C3_1_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM - # VBF_HHTobbVV_CV_0_5_C2V_1_C3_1: /VBF_HHTobbVV_CV_0_5_C2V_1_C3_1_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM \ No newline at end of file + GluGluToHHTo4B_node_cHHH1_preUL: /GluGluToHHTo4B_node_cHHH1_TuneCP5_PSWeights_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM \ No newline at end of file diff --git a/paper/samples/MC_sig.yaml b/paper/samples/MC_sig.yaml index 3208669a..22724bc2 100644 --- a/paper/samples/MC_sig.yaml +++ b/paper/samples/MC_sig.yaml @@ -4,4 +4,10 @@ samples: GluGluToHHTobbVV_node_cHHH1_pn4q: /GluGluToHHTobbVV_node_cHHH1_TuneCP5_13TeV-powheg-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v1/MINIAODSIM GluGluToHHTobbVV_node_cHHH0_pn4q: /GluGluToHHTobbVV_node_cHHH0_TuneCP5_13TeV-powheg-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v1/MINIAODSIM GluGluToHHTobbVV_node_cHHH2p45_pn4q: /GluGluToHHTobbVV_node_cHHH2p45_TuneCP5_13TeV-powheg-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v1/MINIAODSIM - GluGluToHHTobbVV_node_cHHH5_pn4q: /GluGluToHHTobbVV_node_cHHH5_TuneCP5_13TeV-powheg-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v1/MINIAODSIM \ No newline at end of file + GluGluToHHTobbVV_node_cHHH5_pn4q: /GluGluToHHTobbVV_node_cHHH5_TuneCP5_13TeV-powheg-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v1/MINIAODSIM + VBF_HHTobbVV_CV_1_C2V_2_C3_1: /VBF_HHTobbVV_CV_1_C2V_2_C3_1_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM + VBF_HHTobbVV_CV_1_C2V_1_C3_2: /VBF_HHTobbVV_CV_1_C2V_1_C3_2_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM + VBF_HHTobbVV_CV_1_C2V_1_C3_0: /VBF_HHTobbVV_CV_1_C2V_1_C3_0_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM + VBF_HHTobbVV_CV_1_C2V_0_C3_1: /VBF_HHTobbVV_CV_1_C2V_0_C3_1_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM + VBF_HHTobbVV_CV_1_5_C2V_1_C3_1: /VBF_HHTobbVV_CV_1_5_C2V_1_C3_1_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM + VBF_HHTobbVV_CV_0_5_C2V_1_C3_1: /VBF_HHTobbVV_CV_0_5_C2V_1_C3_1_TuneCP5_13TeV-madgraph-pythia8/RunIISummer20UL17MiniAODv2-106X_mc2017_realistic_v9-v2/MINIAODSIM \ No newline at end of file diff --git a/paper/tables/datasets_mcbg.tex b/paper/tables/datasets_mcbg.tex index 64b80e09..4e61642c 100644 --- a/paper/tables/datasets_mcbg.tex +++ b/paper/tables/datasets_mcbg.tex @@ -1,45 +1,46 @@ -\multirow{3}{*}{Diboson} & WW\_TuneCP5\_13TeV-pythia8 & 75.83 \\ - & WZ\_TuneCP5\_13TeV-pythia8 & 27.56 \\ - & ZZ\_TuneCP5\_13TeV-pythia8 & 12.14 \\ +\multirow{3}{*}{Diboson} & WW\_TuneCP5\_13TeV-pythia8 & 118.7~\cite{Gehrmann:2014fva} \\ + & WZ\_TuneCP5\_13TeV-pythia8 & 46.74~\cite{Campbell:2011bn} \\ + & ZZ\_TuneCP5\_13TeV-pythia8 & 16.91~\cite{Cascioli:2014yka} \\ \hline -\multirow{2}{*}{HH} & VBF\_HHTobbVV\_CV\_1\_C2V\_1\_C3\_1\_TuneCP5\_13TeV-madgraph-pythia8 & 0.00 \\ - & GluGluToHHTo4B\_node\_cHHH1\_TuneCP5\_PSWeights\_13TeV-powheg-pythia8 & 0.01 \\ +\multirow{1}{*}{HH~\cite{LHCHiggsHH}} & GluGluToHHTo4B\_node\_cHHH1\_TuneCP5\_PSWeights\_13TeV-powheg-pythia8 & 0.011 \\ \hline -\multirow{5}{*}{HWW (inclusive)} & GluGluHToWW\_Pt-200ToInf\_M-125\_TuneCP5\_MINLO\_13TeV-powheg-pythia8 & 0.10 \\ - & HWminusJ\_HToWW\_M-125\_TuneCP5\_13TeV-powheg-jhugen727-pythia8 & 0.12 \\ - & HWplusJ\_HToWW\_M-125\_TuneCP5\_13TeV-powheg-jhugen727-pythia8 & 0.19 \\ - & HZJ\_HToWW\_M-125\_TuneCP5\_13TeV-powheg-jhugen727-pythia8 & 0.20 \\ +\multirow{5}{*}{HWW (inclusive)} & GluGluHToWW\_Pt-200ToInf\_M-125\_TuneCP5\_MINLO\_13TeV-powheg-pythia8 & 0.101 \\ + & HWminusJ\_HToWW\_M-125\_TuneCP5\_13TeV-powheg-jhugen727-pythia8 & 0.116 \\ + & HWplusJ\_HToWW\_M-125\_TuneCP5\_13TeV-powheg-jhugen727-pythia8 & 0.186 \\ + & HZJ\_HToWW\_M-125\_TuneCP5\_13TeV-powheg-jhugen727-pythia8 & 0.205 \\ & ttHToNonbb\_M125\_TuneCP5\_13TeV-powheg-pythia8 & 0.21 \\ \hline -\multirow{7}{*}{Hbb (inclusive)} & GluGluHToBB\_Pt-200ToInf\_M-125\_TuneCP5\_MINLO\_13TeV-powheg-pythia8 & 0.27 \\ +\multirow{7}{*}{Hbb (inclusive)} & GluGluHToBB\_Pt-200ToInf\_M-125\_TuneCP5\_MINLO\_13TeV-powheg-pythia8 & 0.274 \\ & VBFHToBB\_M-125\_TuneCP5\_13TeV-powheg-pythia8 & 2.25 \\ - & ttHTobb\_M125\_TuneCP5\_13TeV-powheg-pythia8 & 0.29 \\ - & WplusH\_HToBB\_WToQQ\_M-125\_TuneCP5\_13TeV-powheg-pythia8 & 0.34 \\ - & WminusH\_HToBB\_WToQQ\_M-125\_TuneCP5\_13TeV-powheg-pythia8 & 0.21 \\ - & ZH\_HToBB\_ZToQQ\_M-125\_TuneCP5\_13TeV-powheg-pythia8 & 0.33 \\ - & ggZH\_HToBB\_ZToQQ\_M-125\_TuneCP5\_13TeV-powheg-pythia8 & 0.03 \\ + & ttHTobb\_M125\_TuneCP5\_13TeV-powheg-pythia8 & 0.291 \\ + & WplusH\_HToBB\_WToQQ\_M-125\_TuneCP5\_13TeV-powheg-pythia8 & 0.342 \\ + & WminusH\_HToBB\_WToQQ\_M-125\_TuneCP5\_13TeV-powheg-pythia8 & 0.213 \\ + & ZH\_HToBB\_ZToQQ\_M-125\_TuneCP5\_13TeV-powheg-pythia8 & 0.326 \\ + & ggZH\_HToBB\_ZToQQ\_M-125\_TuneCP5\_13TeV-powheg-pythia8 & 0.025 \\ \hline -\multirow{9}{*}{QCD} & QCD\_HT50to100\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 248600000.00 \\ - & QCD\_HT100to200\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 27990000.00 \\ - & QCD\_HT200to300\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 1712000.00 \\ - & QCD\_HT300to500\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 347700.00 \\ - & QCD\_HT500to700\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 30330.00 \\ - & QCD\_HT700to1000\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 6412.00 \\ - & QCD\_HT1000to1500\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 1118.00 \\ - & QCD\_HT1500to2000\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 108.50 \\ +\multirow{9}{*}{QCD} & QCD\_HT50to100\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 248600000 \\ + & QCD\_HT100to200\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 27990000 \\ + & QCD\_HT200to300\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 1712 \\ + & QCD\_HT300to500\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 3477 \\ + & QCD\_HT500to700\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 30330 \\ + & QCD\_HT700to1000\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 6412 \\ + & QCD\_HT1000to1500\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 1118 \\ + & QCD\_HT1500to2000\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 108.5 \\ & QCD\_HT2000toInf\_TuneCP5\_PSWeights\_13TeV-madgraph-pythia8 & 21.94 \\ \hline -\multirow{4}{*}{SingleTop} & ST\_t-channel\_top\_4f\_InclusiveDecays\_TuneCP5\_13TeV-powheg-madspin-pythia8 & 113.40 \\ - & ST\_t-channel\_antitop\_4f\_InclusiveDecays\_TuneCP5\_13TeV-powheg-madspin-pythia8 & 67.93 \\ - & ST\_tW\_antitop\_5f\_inclusiveDecays\_TuneCP5\_13TeV-powheg-pythia8 & 32.51 \\ - & ST\_tW\_top\_5f\_inclusiveDecays\_TuneCP5\_13TeV-powheg-pythia8 & 32.45 \\ +\multirow{4}{*}{SingleTop~\cite{SingleTopNNLORef}} & ST\_t-channel\_top\_4f\_InclusiveDecays\_TuneCP5\_13TeV-powheg-madspin-pythia8 & 134.2 \\ + & ST\_t-channel\_antitop\_4f\_InclusiveDecays\_TuneCP5\_13TeV-powheg-madspin-pythia8 & 80 \\ + & ST\_tW\_antitop\_5f\_inclusiveDecays\_TuneCP5\_13TeV-powheg-pythia8 & 39.65 \\ + & ST\_tW\_top\_5f\_inclusiveDecays\_TuneCP5\_13TeV-powheg-pythia8 & 39.65 \\ \hline -\multirow{2}{*}{TTbar} & TTToHadronic\_TuneCP5\_13TeV-powheg-pythia8 & 303.85 \\ - & TTToSemiLeptonic\_TuneCP5\_13TeV-powheg-pythia8 & 306.14 \\ +\multirow{2}{*}{TTbar~\cite{TtbarNNLO}} & TTToHadronic\_TuneCP5\_13TeV-powheg-pythia8 & 370.042 \\ + & TTToSemiLeptonic\_TuneCP5\_13TeV-powheg-pythia8 & 369.487 \\ \hline -\multirow{6}{*}{V+Jets} & WJetsToQQ\_HT-400to600\_TuneCP5\_13TeV-madgraphMLM-pythia8 & 277.00 \\ +\multirow{8}{*}{V+Jets} & WJetsToQQ\_HT-200to400\_TuneCP5\_13TeV-madgraphMLM-pythia8 & 2549 \\ + & WJetsToQQ\_HT-400to600\_TuneCP5\_13TeV-madgraphMLM-pythia8 & 277 \\ & WJetsToQQ\_HT-600to800\_TuneCP5\_13TeV-madgraphMLM-pythia8 & 59.06 \\ & WJetsToQQ\_HT-800toInf\_TuneCP5\_13TeV-madgraphMLM-pythia8 & 28.75 \\ - & ZJetsToQQ\_HT-400to600\_TuneCP5\_13TeV-madgraphMLM-pythia8 & 114.50 \\ + & ZJetsToQQ\_HT-200to600\_TuneCP5\_13TeV-madgraphMLM-pythia8 & 1012 \\ + & ZJetsToQQ\_HT-400to600\_TuneCP5\_13TeV-madgraphMLM-pythia8 & 114.5 \\ & ZJetsToQQ\_HT-600to800\_TuneCP5\_13TeV-madgraphMLM-pythia8 & 25.41 \\ & ZJetsToQQ\_HT-800toInf\_TuneCP5\_13TeV-madgraphMLM-pythia8 & 12.91 \ No newline at end of file diff --git a/paper/tables/datasets_mcsig_nonres.tex b/paper/tables/datasets_mcsig_nonres.tex index 8ec3de0d..d0d64e17 100644 --- a/paper/tables/datasets_mcsig_nonres.tex +++ b/paper/tables/datasets_mcsig_nonres.tex @@ -2,3 +2,9 @@ ggHH BSM ($\kappa_{\lambda} = 0$) & GluGluToHHTobbVV\_node\_cHHH0\_TuneCP5\_13TeV-powheg-pythia8 & 9.02 \\ ggHH BSM ($\kappa_{\lambda} = 2.45$) & GluGluToHHTobbVV\_node\_cHHH2p45\_TuneCP5\_13TeV-powheg-pythia8 & 1.70 \\ ggHH BSM ($\kappa_{\lambda} = 5$) & GluGluToHHTobbVV\_node\_cHHH5\_TuneCP5\_13TeV-powheg-pythia8 & 11.80 \\ +VBFHH BSM ($\kappa_\mathrm{V} = 1$, $\kappa_\mathrm{2V} = 2$, $\kappa_{\lambda} = 1$) & VBF\_HHTobbVV\_CV\_1\_C2V\_2\_C3\_1\_TuneCP5\_13TeV-madgraph-pythia8 & 1.84 \\ +VBFHH BSM ($\kappa_\mathrm{V} = 1$, $\kappa_\mathrm{2V} = 1$, $\kappa_{\lambda} = 2$) & VBF\_HHTobbVV\_CV\_1\_C2V\_1\_C3\_2\_TuneCP5\_13TeV-madgraph-pythia8 & 0.18 \\ +VBFHH BSM ($\kappa_\mathrm{V} = 1$, $\kappa_\mathrm{2V} = 1$, $\kappa_{\lambda} = 0$) & VBF\_HHTobbVV\_CV\_1\_C2V\_1\_C3\_0\_TuneCP5\_13TeV-madgraph-pythia8 & 0.60 \\ +VBFHH BSM ($\kappa_\mathrm{V} = 1$, $\kappa_\mathrm{2V} = 0$, $\kappa_{\lambda} = 1$) & VBF\_HHTobbVV\_CV\_1\_C2V\_0\_C3\_1\_TuneCP5\_13TeV-madgraph-pythia8 & 3.50 \\ +VBFHH BSM ($\kappa_\mathrm{V} = 1.5$, $\kappa_\mathrm{2V} = 1$, $\kappa_{\lambda} = 1$) & VBF\_HHTobbVV\_CV\_1\_5\_C2V\_1\_C3\_1\_TuneCP5\_13TeV-madgraph-pythia8 & 8.54 \\ +VBFHH BSM ($\kappa_\mathrm{V} = 0.5$, $\kappa_\mathrm{2V} = 1$, $\kappa_{\lambda} = 1$) & VBF\_HHTobbVV\_CV\_0\_5\_C2V\_1\_C3\_1\_TuneCP5\_13TeV-madgraph-pythia8 & 1.40 \ No newline at end of file diff --git a/src/HHbbVV/postprocessing/PostProcess.ipynb b/src/HHbbVV/postprocessing/PostProcess.ipynb index 82e46657..eeca032e 100644 --- a/src/HHbbVV/postprocessing/PostProcess.ipynb +++ b/src/HHbbVV/postprocessing/PostProcess.ipynb @@ -132,8 +132,12 @@ "cutflow = pd.DataFrame(index=list(samples.keys()) + list(nonres_samples.keys()))\n", "\n", "# utils.remove_empty_parquets(samples_dir, year)\n", - "events_dict = utils.load_samples(signal_samples_dir, nonres_samples, year, filters)\n", - "events_dict |= utils.load_samples(samples_dir, samples, year, filters)\n", + "\n", + "# no HEM cleaning for now because it wasn't done for BDT samples\n", + "events_dict = utils.load_samples(\n", + " signal_samples_dir, nonres_samples, year, filters, hem_cleaning=False\n", + ")\n", + "events_dict |= utils.load_samples(samples_dir, samples, year, filters, hem_cleaning=False)\n", "\n", "utils.add_to_cutflow(events_dict, \"BDTPreselection\", \"weight\", cutflow)\n", "\n", @@ -330,6 +334,50 @@ ")" ] }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "events = events_dict[\"HHbbVV\"]\n", + "bb_mask = bb_masks[\"HHbbVV\"]\n", + "weight = events[\"finalWeight\"].values.squeeze()\n", + "weight_lp = weight * events[\"VV_lp_sf_nom\"].values.squeeze()\n", + "weight_lp_sys_up = weight * events[\"VV_lp_sf_sys_up\"].values.squeeze()\n", + "weight_lp_sys_down = weight * events[\"VV_lp_sf_sys_down\"].values.squeeze()\n", + "\n", + "plt.hist(\n", + " utils.get_feat(events, \"bbFatJetPt\", bb_mask),\n", + " np.linspace(250, 1200, 31),\n", + " weights=weight,\n", + " histtype=\"step\",\n", + " label=\"Pre-LP\",\n", + ")\n", + "plt.hist(\n", + " utils.get_feat(events, \"bbFatJetPt\", bb_mask),\n", + " np.linspace(250, 1200, 31),\n", + " weights=weight_lp,\n", + " histtype=\"step\",\n", + " label=\"Post-LP\",\n", + ")\n", + "plt.title(\"2018 HHbbVV\")\n", + "plt.xlabel(r\"$p_T^{VV}$ (GeV)\")\n", + "plt.ylabel(\"Events\")\n", + "# plt.hist(utils.get_feat(events, \"VVFatJetPt\", bb_mask), np.linspace(250, 2000, 31), weights=weight_lp_sys_up, histtype=\"step\", label=\"Post-LP Sys Up\")\n", + "# plt.hist(utils.get_feat(events, \"VVFatJetPt\", bb_mask), np.linspace(250, 2000, 31), weights=weight_lp_sys_down, histtype=\"step\", label=\"Post-LP Sys Down\")\n", + "plt.legend()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "list(events_dict[\"HHbbVV\"].columns)" + ] + }, { "attachments": {}, "cell_type": "markdown", diff --git a/src/HHbbVV/triton/export_and_check.ipynb b/src/HHbbVV/triton/export_and_check.ipynb index c0f9ed16..5f433864 100644 --- a/src/HHbbVV/triton/export_and_check.ipynb +++ b/src/HHbbVV/triton/export_and_check.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "metadata": {}, "outputs": [], "source": [