Commit 85c5f6d0 authored by Benjamin Fischer's avatar Benjamin Fischer
Browse files

Merge branch 'miniaod' into dask

parents 10877401 c3b26424
......@@ -416,8 +416,7 @@ dataset_tt_sl = od.Dataset(
200,
campaign=campaign,
keys=[
"/TTToSemiLeptonic_TuneCP5_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v1/MINIAODSIM",
"/TTToSemiLeptonic_TuneCP5_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14_ext1-v1/MINIAODSIM",
"/TTToSemiLeptonic_TuneCP5_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v1/MINIAODSIM"
],
)
......@@ -426,8 +425,7 @@ dataset_tt_dl = od.Dataset(
201,
campaign=campaign,
keys=[
"/TTTo2L2Nu_TuneCP5_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v2/MINIAODSIM",
"/TTTo2L2Nu_TuneCP5_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-fixECALGT_LowPU_94X_mc2017_realistic_v10For2017H_v2_ext1-v1/MINIAODSIM",
"/TTTo2L2Nu_TuneCP5_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v1/MINIAODSIM",
"/TTTo2L2Nu_TuneCP5_PSweights_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v2/MINIAODSIM",
],
)
......@@ -448,7 +446,16 @@ dataset_dy_lep_10To50 = od.Dataset(
300,
campaign=campaign,
keys=[
"/DYJetsToLL_M-10to50_TuneCP5_13TeV-madgraphMLM-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v2/MINIAODSIM",
"/DYJetsToLL_M-10to50_TuneCP5_13TeV-madgraphMLM-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
dataset_dy_lep_50ToInf = od.Dataset(
"dy_lep_50ToInf",
299,
campaign=campaign,
keys=[
"/DYJetsToLL_M-50_TuneCP5_13TeV-amcatnloFXFX-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
"/DYJetsToLL_M-50_TuneCP5_13TeV-amcatnloFXFX-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14_ext1-v1/MINIAODSIM",
],
)
# c
......@@ -457,7 +464,7 @@ dataset_dy_lep_0j = od.Dataset(
301,
campaign=campaign,
keys=[
"/DYJetsToLL_0J_TuneCP5_13TeV-amcatnloFXFX-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v*/MINIAODSIM",
"/DYJetsToLL_0J_TuneCP5_13TeV-amcatnloFXFX-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
# c
......@@ -466,7 +473,7 @@ dataset_dy_lep_1j = od.Dataset(
302,
campaign=campaign,
keys=[
"/DYJetsToLL_1J_TuneCP5_13TeV-amcatnloFXFX-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v*/MINIAODSIM",
"/DYJetsToLL_1J_TuneCP5_13TeV-amcatnloFXFX-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
# c
......@@ -475,7 +482,7 @@ dataset_dy_lep_2j = od.Dataset(
303,
campaign=campaign,
keys=[
"/DYJetsToLL_2J_TuneCP5_13TeV-amcatnloFXFX-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v*/MINIAODSIM",
"/DYJetsToLL_2J_TuneCP5_13TeV-amcatnloFXFX-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
......@@ -652,99 +659,37 @@ dataset_dyjetstoll_m_50_incl.set_aux("kfactor", 1.23)
# Di-Boson VV
# Non-resonant WW
dataset_ww_sl = od.Dataset(
"ww_sl",
dataset_ww = od.Dataset(
"ww",
400,
campaign=campaign,
n_files=122,
keys=[
"/WWToLNuQQ_NNPDF31_TuneCP5_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
dataset_ww_dl = od.Dataset(
"ww_dl",
401,
campaign=campaign,
n_files=78,
keys=[
"/WWTo2L2Nu_NNPDF31_TuneCP5Down_PSweights_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
"/WW_TuneCP5_13TeV-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
dataset_ww_fh = od.Dataset(
"ww_fh",
4002,
campaign=campaign,
n_files=26,
keys=[
"/WWTo4Q_NNPDF31_TuneCP5_PSweights_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14_ext1-v1/MINIAODSIM",
],
)
# ZZ
dataset_zz_2l_2q = od.Dataset(
"zz_2l_2q",
402,
campaign=campaign,
n_files=409,
keys=[
"/ZZTo2L2Q_13TeV_amcatnloFXFX_madspin_pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
dataset_zz_2l_2nu = od.Dataset(
"zz_2l_2nu",
403,
dataset_zz = od.Dataset(
"zz",
401,
campaign=campaign,
n_files=112,
keys=[
"/ZZTo2L2Nu_13TeV_powheg_pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
"/ZZ_TuneCP5_13TeV-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
dataset_zz_4l = od.Dataset(
"zz_4l",
404,
campaign=campaign,
n_files=362 + 1128 + 1623,
keys=[
"/ZZTo4L_13TeV_powheg_pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v2/MINIAODSIM",
"/ZZTo4L_13TeV_powheg_pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14_ext1-v1/MINIAODSIM",
"/ZZTo4L_13TeV_powheg_pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14_ext2-v1/MINIAODSIM",
],
)
# WZ
dataset_wz_2l_2q = od.Dataset(
"wz_2l_2q",
405,
campaign=campaign,
n_files=406,
keys=[
"/WZTo2L2Q_13TeV_amcatnloFXFX_madspin_pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
dataset_wz_1l_3nu = od.Dataset(
"wz_1l_3nu",
406,
dataset_wz = od.Dataset(
"wz",
402,
campaign=campaign,
n_files=136,
keys=[
"/WZTo1L3Nu_13TeV_amcatnloFXFX_madspin_pythia8_v2/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
"/WZ_TuneCP5_13TeV-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
dataset_wz_3l_1nu = od.Dataset(
"wz_3l_1nu",
407,
campaign=campaign,
n_files=131,
keys=[
"/WZTo3LNu_TuneCP5_13TeV-amcatnloFXFX-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
# Single top
# t-channel
......@@ -754,7 +699,7 @@ dataset_st_t_channel_top = od.Dataset(
500,
campaign=campaign,
keys=[
"/ST_t-channel_top_4f_inclusiveDecays_TuneCP5_13TeV-powhegV2-madspin-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v*/MINIAODSIM",
"/ST_t-channel_top_4f_inclusiveDecays_TuneCP5_13TeV-powhegV2-madspin-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
# c
......@@ -774,7 +719,7 @@ dataset_st_s_channel = od.Dataset(
502,
campaign=campaign,
keys=[
"/ST_s-channel_4f_leptonDecays_TuneCP5_PSweights_13TeV-amcatnlo-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v1/MINIAODSIM",
"/ST_s-channel_4f_leptonDecays_TuneCP5_13TeV-amcatnlo-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v2/MINIAODSIM",
],
)
......@@ -785,7 +730,7 @@ dataset_st_tW_top = od.Dataset(
503,
campaign=campaign,
keys=[
"/ST_tW_top_5f_NoFullyHadronicDecays_TuneCP5_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v*/MINIAODSIM",
"/ST_tW_top_5f_inclusiveDecays_TuneCP5_PSweights_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
......@@ -794,90 +739,21 @@ dataset_st_tW_antitop = od.Dataset(
504,
campaign=campaign,
keys=[
"/ST_tW_antitop_5f_NoFullyHadronicDecays_TuneCP5_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v2/MINIAODSIM",
"/ST_tW_antitop_5f_inclusiveDecays_TuneCP5_PSweights_13TeV-powheg-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v2/MINIAODSIM",
],
)
# W+jets
# 100to200 exists only v2
# 2500to exists only v3
# rest is v1
# c
dataset_wjets_100to200 = od.Dataset(
"wjets_100to200",
dataset_wjets = od.Dataset(
"wjets",
600,
campaign=campaign,
keys=[
"/WJetsToLNu_HT-100To200_TuneCP5_13TeV-madgraphMLM-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v*/MINIAODSIM",
"/WJetsToLNu_TuneCP5_13TeV-madgraphMLM-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v2/MINIAODSIM",
"/WJetsToLNu_TuneCP5_13TeV-madgraphMLM-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14_ext1-v2/MINIAODSIM",
],
)
dataset_wjets_100to200.set_aux("kfactor", 1.21)
# c
dataset_wjets_200to400 = od.Dataset(
"wjets_200to400",
601,
campaign=campaign,
keys=[
"/WJetsToLNu_HT-200To400_TuneCP5_13TeV-madgraphMLM-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v*/MINIAODSIM",
],
)
dataset_wjets_200to400.set_aux("kfactor", 1.21)
# c
dataset_wjets_400to600 = od.Dataset(
"wjets_400to600",
602,
campaign=campaign,
keys=[
"/WJetsToLNu_HT-400To600_TuneCP5_13TeV-madgraphMLM-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v*/MINIAODSIM",
],
)
dataset_wjets_400to600.set_aux("kfactor", 1.21)
# c
dataset_wjets_600to800 = od.Dataset(
"wjets_600to800",
603,
campaign=campaign,
keys=[
"/WJetsToLNu_HT-600To800_TuneCP5_13TeV-madgraphMLM-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v*/MINIAODSIM",
],
)
dataset_wjets_600to800.set_aux("kfactor", 1.21)
# c
dataset_wjets_800to1200 = od.Dataset(
"wjets_800to1200",
604,
campaign=campaign,
keys=[
"/WJetsToLNu_HT-800To1200_TuneCP5_13TeV-madgraphMLM-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v*/MINIAODSIM",
],
)
dataset_wjets_800to1200.set_aux("kfactor", 1.21)
# c
dataset_wjets_1200to2500 = od.Dataset(
"wjets_1200to2500",
605,
campaign=campaign,
keys=[
"/WJetsToLNu_HT-1200To2500_TuneCP5_13TeV-madgraphMLM-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v*/MINIAODSIM",
],
)
dataset_wjets_1200to2500.set_aux("kfactor", 1.21)
# c
dataset_wjets_2500toInf = od.Dataset(
"wjets_2500toInf",
606,
campaign=campaign,
keys=[
"/WJetsToLNu_HT-2500ToInf_TuneCP5_13TeV-madgraphMLM-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v*/MINIAODSIM",
],
)
dataset_wjets_2500toInf.set_aux("kfactor", 1.21)
dataset_wbjets_100to200 = od.Dataset(
"wbjets_100to200",
......@@ -912,12 +788,12 @@ dataset_ttWjets_2q = od.Dataset(
],
)
# c
dataset_ttWjets_1l_2nu = od.Dataset(
"ttWjets_1l_2nu",
dataset_ttWjets_1l_1nu = od.Dataset(
"ttWjets_1l_1nu",
701,
campaign=campaign,
keys=[
"/TTWJetsToLNu_TuneCP5_13TeV-amcatnloFXFX-madspin-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v2/MINIAODSIM",
"/TTWJetsToLNu_TuneCP5_PSweights_13TeV-amcatnloFXFX-madspin-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_new_pmx_94X_mc2017_realistic_v14-v1/MINIAODSIM",
],
)
......@@ -927,13 +803,12 @@ dataset_ttZ_2q = od.Dataset(
campaign=campaign,
keys=[
"/TTZToQQ_TuneCP5_13TeV-amcatnlo-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v2/MINIAODSIM",
"/TTZToQQ_TuneCP5_13TeV-amcatnlo-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14_ext1-v1/MINIAODSIM",
],
)
# c
dataset_ttZ_2l_2nu = od.Dataset(
"ttZ_2l_2nu",
15886960,
703,
campaign=campaign,
keys=[
"/TTZToLLNuNu_M-10_TuneCP5_13TeV-amcatnlo-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v2/MINIAODSIM",
......@@ -1211,31 +1086,19 @@ dataset_dy_lep_1j.add_process(process_dy_lep_1j)
dataset_dy_lep_2j.add_process(process_dy_lep_2j)
dataset_dy_bjets_100to200.add_process(process_dy_bjets_100to200)
dataset_dy_bjets_200toInf.add_process(process_dy_bjets_200toInf)
dataset_ww_sl.add_process(process_ww_sl)
dataset_ww_dl.add_process(process_ww_dl)
dataset_ww_dl.add_process(process_ww_fh)
dataset_zz_2l_2q.add_process(process_zz_2l_2q)
dataset_zz_2l_2nu.add_process(process_zz_2l_2nu)
dataset_zz_4l.add_process(process_zz_4l)
dataset_wz_2l_2q.add_process(process_wz_2l_2q)
dataset_wz_1l_3nu.add_process(process_wz_1l_3nu)
dataset_wz_3l_1nu.add_process(process_wz_3l_1nu)
dataset_ww.add_process(process_ww)
dataset_zz.add_process(process_zz)
dataset_wz.add_process(process_wz)
dataset_st_t_channel_top.add_process(process_st_t_channel_top)
dataset_st_t_channel_antitop.add_process(process_st_t_channel_antitop)
dataset_st_s_channel.add_process(process_st_s_channel)
dataset_st_tW_top.add_process(process_st_tW_top)
dataset_st_tW_antitop.add_process(process_st_tW_antitop)
dataset_wjets_100to200.add_process(process_wjets_100to200)
dataset_wjets_200to400.add_process(process_wjets_200to400)
dataset_wjets_400to600.add_process(process_wjets_400to600)
dataset_wjets_600to800.add_process(process_wjets_600to800)
dataset_wjets_800to1200.add_process(process_wjets_800to1200)
dataset_wjets_1200to2500.add_process(process_wjets_1200to2500)
dataset_wjets_2500toInf.add_process(process_wjets_2500toInf)
dataset_wjets.add_process(process_wjets)
dataset_wbjets_100to200.add_process(process_wbjets_100to200)
dataset_wbjets_200toInf.add_process(process_wbjets_200toInf)
dataset_ttWjets_2q.add_process(process_ttWjets_2q)
dataset_ttWjets_1l_2nu.add_process(process_ttWjets_1l_2nu)
dataset_ttWjets_1l_1nu.add_process(process_ttWjets_1l_1nu)
dataset_ttZ_2q.add_process(process_ttZ_2q)
dataset_ttZ_2l_2nu.add_process(process_ttZ_2l_2nu)
dataset_qcd_ht100to200.add_process(process_qcd_ht100to200)
......@@ -1281,7 +1144,6 @@ dataset_dyjetstoll_m_50_ht_800to1200.add_process(process_dyjetstoll_m_50_ht_800t
dataset_dyjetstoll_m_50_ht_1200to2500.add_process(process_dyjetstoll_m_50_ht_1200to2500)
dataset_dyjetstoll_m_50_ht_2500toInf.add_process(process_dyjetstoll_m_50_ht_2500toInf)
dataset_dyjetstoll_m_50_incl.add_process(process_dyjetstoll_m_50_incl)
dataset_wjetstolnu_incl.add_process(process_wjetstolnu_incl)
dataset_dyjetstoll_bgenfilter_zpt_100to200.add_process(
process_dyjetstoll_bgenfilter_zpt_100to200
)
......
This diff is collapsed.
......@@ -17,3 +17,4 @@ BR_W_LEP = 1 - BR_W_HAD
BR_WW_SL = 2 * BR_W_HAD.mul(BR_W_LEP, rho=-1, inplace=False)
BR_WW_DL = BR_W_LEP ** 2
BR_WW_FH = BR_W_HAD ** 2
BR_Z_CLEP = sn.Number(0.033658, {"br_z_clep": 0.000023}) * N_LEPS
......@@ -12,7 +12,16 @@ Values taken from:
import order as od
import scinum as sn
from config.constants import N_LEPS, BR_WW_DL, BR_WW_SL, BR_H_BB, BR_H_WW, BR_WW_FH
from config.constants import (
N_LEPS,
BR_WW_DL,
BR_WW_SL,
BR_H_BB,
BR_H_WW,
BR_WW_FH,
BR_W_LEP,
BR_Z_CLEP,
)
process_data = od.Process("data", 0, is_data=True, label="data",)
process_data_dl = process_data.add_process("data_dl", 1, is_data=True, label="data",)
......@@ -131,6 +140,7 @@ process_tt = od.Process(
process_tt_sl = process_tt.add_process(
"tt_sl",
2100,
color=(255, 235, 59),
label=r"$t\bar{t}$ + Jets, SL",
xsecs={13: process_tt.get_xsec(13) * BR_WW_SL,},
)
......@@ -154,36 +164,28 @@ process_dy_lep = process_dy.add_process(
"dy_lep", 3100, label=r"Drell-Yan, $Z \rightarrow ll$",
)
# process_dy_lep_50ToInf = process_dy_lep.add_process(
# "dy_lep_50ToInf", 3160,
# label=r"Drell-Yan, $m_{ll} \gt 50$",
# xsecs={
# 13: sn.Number(1921.8, {
# "integration": 0.6,
# "pdf": 33.2,
# }) * N_LEPS,
# },
# )
# process_dy_lep_5To50 = process_dy_lep.add_process(
# "dy_lep_5To50", 3110,
# "dy_lep_5To50",
# 3101,
# label=r"Drell-Yan, $5 \leq m_{ll} \leq 50$",
# xsecs={
# 13: sn.Number(71310, {
# "scale": 70,
# }),
# },
# xsecs={13: sn.Number(71310, {"scale": 70})},
# )
#
# process_dy_lep_10To50 = process_dy_lep.add_process(
# "dy_lep_10To50", 3120,
# "dy_lep_10To50",
# 3102,
# label=r"Drell-Yan, $10 \leq m_{ll} \leq 50$",
# xsecs={
# 13: sn.Number(18610),
# },
# xsecs={13: sn.Number(18610)},
# )
#
# process_dy_lep_50ToInf = process_dy_lep.add_process(
# "dy_lep_50ToInf",
# 222,
# label=r"Drell-Yan, $m_{ll} > 50$",
# xsecs={13: sn.Number(1921.8, {"integration": 0.6, "pdf": 33.2,}) * N_LEPS,},
# )
process_dy_lep_0j = process_dy_lep.add_process(
"dy_lep_0j", 3130, label=r"Drell-Yan, $0 Jets$", xsecs={13: sn.Number(4758.9),},
)
......@@ -193,8 +195,13 @@ process_dy_lep_1j = process_dy_lep.add_process(
)
process_dy_lep_2j = process_dy_lep.add_process(
"dy_lep_2j", 3150, label=r"Drell-Yan, $2 Jets$", xsecs={13: sn.Number(337.1),},
"dy_lep_2j",
3150,
color=(76, 175, 80),
label=r"Drell-Yan, $2 Jets$",
xsecs={13: sn.Number(337.1),},
)
# Drell Yan b enriched
process_dyb = od.Process("dyb", 3001, color=(76, 175, 80), label="Drell-Yan bjets",)
# only b enriched
......@@ -229,7 +236,7 @@ process_dyjetstoll_bgenfilter_zpt_200toInf = process_dy_benriched.add_process(
)
# DY without b renriched samples
process_dynob = process_dyb.add_process(
"dynob", 3300, color=(76, 175, 80), label="Drell-Yan no benriched",
"dynob", 3300, color=(76, 175, 80), label="Drell-Yan",
)
process_dyjetstoll_m_4to50_ht_100to200 = process_dynob.add_process(
"DYJetsToLL_M_4to50_HT_100to200",
......@@ -316,45 +323,71 @@ process_ww = process_vv.add_process(
xsecs={13: sn.Number(118.7, {"scale": ("rel", 2.5, 2.2), "pdf": ("rel", 0.046),}),},
)
process_ww_sl = process_ww.add_process(
"ww_sl", 4110, label="WW (SL)", xsecs={13: process_ww.get_xsec(13) * BR_WW_SL,},
)
process_ww_dl = process_ww.add_process(
"ww_dl", 4120, label="WW (DL)", xsecs={13: process_ww.get_xsec(13) * BR_WW_DL,},
)
process_ww_fh = process_ww.add_process(
"ww_fh", 4130, label="WW (FH)", xsecs={13: process_ww.get_xsec(13) * BR_WW_FH,},
)
process_zz = process_vv.add_process("zz", 4200, label="ZZ",)
process_zz_2l_2q = process_zz.add_process(
"zz_2l_2q", 4210, label="ZZ (2l2q)", xsecs={13: sn.Number(3.22),},
)
# process_ww_sl = process_ww.add_process(
# "ww_sl",
# 4110,
# color=(33, 150, 243),
# label="WW (SL)",
# xsecs={13: process_ww.get_xsec(13) * BR_WW_SL,},
# )
#
# process_ww_dl = process_ww.add_process(
# "ww_dl", 4120, label="WW (DL)", xsecs={13: process_ww.get_xsec(13) * BR_WW_DL,},
# )
#
# process_ww_fh = process_ww.add_process(
# "ww_fh", 4130, label="WW (FH)", xsecs={13: process_ww.get_xsec(13) * BR_WW_FH,},
# )
process_zz_2l_2nu = process_zz.add_process(
"zz_2l_2nu", 4220, label="ZZ (2l2nu)", xsecs={13: sn.Number(0.564),},
process_zz = process_vv.add_process(
"zz",
430,
label="ZZ",
xsecs={
13: sn.Number(0.0719, {"scale": 0.0023, "pdf": 0.0027})
* ((N_LEPS / BR_Z_CLEP) ** 2 / 2.0),
},
)
process_zz_4l = process_zz.add_process(
"zz_4l", 4230, label="ZZ (4l)", xsecs={13: sn.Number(1.256),},
)
process_wz = process_vv.add_process("wz", 4300, label="WZ",)
# process_zz_2l_2q = process_zz.add_process(
# "zz_2l_2q", 4210, label="ZZ (2l2q)", xsecs={13: sn.Number(3.22),},
# )
#
# process_zz_2l_2nu = process_zz.add_process(
# "zz_2l_2nu", 4220, label="ZZ (2l2nu)", xsecs={13: sn.Number(0.564),},
# )
#
# process_zz_4l = process_zz.add_process(
# "zz_4l", 4230, label="ZZ (4l)", xsecs={13: sn.Number(1.256),},
# )
process_wz_2l_2q = process_wz.add_process(
"wz_2l_2q", 4310, label="WZ (2l2q)", xsecs={13: sn.Number(5.595),},
process_wz = process_vv.add_process(
"wz",
420,
label="WZ",
xsecs={
13: sn.Number.add(
sn.Number(0.106, {"scale": 0.0036, "pdf": 0.0050}),
sn.Number(0.0663, {"scale": 0.0029, "pdf": 0.0032}),
rho=1,
)
* (N_LEPS / BR_W_LEP * N_LEPS / BR_Z_CLEP),
},
)
process_wz_1l_3nu = process_wz.add_process(
"wz_1l_3nu", 4320, label="WZ (1l3nu)", xsecs={13: sn.Number(3.033),},
)
process_wz_3l_1nu = process_wz.add_process(
"wz_3l_1nu", 4330, label="WZ (3l1nu)", xsecs={13: sn.Number(4.42965),},
)
# process_wz_2l_2q = process_wz.add_process(
# "wz_2l_2q", 4310, label="WZ (2l2q)", xsecs={13: sn.Number(5.595),},
# )
#
# process_wz_1l_3nu = process_wz.add_process(
# "wz_1l_3nu", 4320, label="WZ (1l3nu)", xsecs={13: sn.Number(3.033),},
# )
#
# process_wz_3l_1nu = process_wz.add_process(
# "wz_3l_1nu", 4330, label="WZ (3l1nu)", xsecs={13: sn.Number(4.42965),},
# )
process_st = od.Process("st", 5000, color=(255, 152, 0), label="Single top",)
......@@ -377,11 +410,11 @@ process_st_s_channel = process_st.add_process(
)
process_st_tW_top = process_st.add_process(
"st_tW_top", 5040, label="tW", xsecs={13: sn.Number(19.5545),},
"st_tW_top", 5040, label="tW", xsecs={13: sn.Number(35.85),},
)
process_st_tW_antitop = process_st.add_process(
"st_tW_antitop", 5050, label=r"$\overline{t}W$", xsecs={13: sn.Number(19.5545),},
"st_tW_antitop", 5050, label=r"$\overline{t}W$", xsecs={13: sn.Number(35.85),},
)
# use wbjets for wjets +b enriched +bGenFilter
......@@ -420,49 +453,59 @@ process_wjetstolnu_bgenfilter_wpt_200toInf = process_wjets_benriched.add_process
)
# only inclusive samples
process_wjets = process_wbjets.add_process(
"wjets", 6000, color=(233, 30, 99), label="W+jets",
)
process_wjets_100to200 = process_wjets.add_process(
"wjets_100to200", 6010, label="W+jets (HT 100-200)", xsecs={13: sn.Number(1395.0),},
)
process_wjets_200to400 = process_wjets.add_process(
"wjets_200to400", 6020, label="W+jets (HT 200-400)", xsecs={13: sn.Number(407.9),},
)
process_wjets_400to600 = process_wjets.add_process(
"wjets_400to600", 6030, label="W+jets (HT 400-600)", xsecs={13: sn.Number(57.48),},
)
# process_wjets = process_wbjets.add_process(
# "wjets", 6000, color=(199, 21, 133), label="W+jets",
# )
process_wjets_600to800 = process_wjets.add_process(
"wjets_600to800", 6040, label="W+jets (HT 600-800)", xsecs={13: sn.Number(12.87),},
)
process_wjets_800to1200 = process_wjets.add_process(
"wjets_800to1200",
6050,
label="W+jets (HT 800-1200)",
xsecs={13: sn.Number(5.366),},
process_wjets = od.Process(
"wjets",
6000,
label=r"W + Jets, $W \rightarrow l \nu$",
color=(199, 21, 133),
xsecs={13: sn.Number(20508.9, {"scale": (165.7, 88.2), "pdf": 770.9,}) * N_LEPS},
)