c.add_category('H2Mu_gg', True)
c.add_dir_to_category('H2Mu_gg', path+'/H2Mu_ggH/', 0.009618, 'H2Mu_gg')
c.add_category('H2Mu_VBF', True)
c.add_dir_to_category('H2Mu_VBF', path+'/H2Mu_VBF/', 0.0008208, 'H2Mu_VBF')
c.add_category('ZJets_MG', False)
c.add_dir_to_category('ZJets_MG', path+'/ZJets_MG/', 5765.4, 'ZJets_MG')
c.add_category('tt_ll_AMC', False)
c.add_dir_to_category('tt_ll_AMC', path+'/tt_ll_AMC/', 85.656*0.9, 'tt_ll_AMC')
c.set_tree_path(treePath)
c.add_variable("muPairs.pt", 1)
c.add_variable("muPairs.eta", 1)
c.add_variable("muPairs.dEta", 1)
c.add_variable("muPairs.dPhi", 1)
c.add_variable("met.pt", 1)
# c.set_year("2016-2orMoreJets")
# var_set = "V2-2orMoreJets"
var_set = "V2"
c.set_year("2016")
if var_set == "V1":
c.add_variable("nJetsCent", 1)
c.add_variable("nJetsFwd", 1)
c.add_category(c.vbf_label, True)
c.add_dir_to_category(VBF_2017_powheg.name, VBF_2017_powheg.path,
VBF_2017_powheg.xSec, c.vbf_label)
c.add_category(c.dy_label, False)
c.add_dir_to_category(ZJets_aMC_2017_hiStat.name, ZJets_aMC_2017_hiStat.path,
ZJets_aMC_2017_hiStat.xSec, c.dy_label)
c.add_category(c.tt_label, False)
c.add_dir_to_category(tt_ll_POW_2017.name, tt_ll_POW_2017.path,
tt_ll_POW_2017.xSec, c.tt_label)
##########################################################
### ------ Raffaele's variables ------ ###
c.add_variable("muPairs.pt", 1)
c.add_variable("muPairs.eta", 1)
c.add_variable("muPairs.dEta", 1)
c.add_variable("muPairs.dPhi", 1)
c.add_variable("met.pt", 1)
c.add_variable("mu1_pt_Roch_over_mass", 1)
c.add_variable("mu2_pt_Roch_over_mass", 1)
c.add_variable('muons.eta', 2)
c.add_variable("min_dR_mu_jet", 1)
c.add_variable("nJets", 1)
c.add_variable("nBMed", 1)
c.add_variable("zeppenfeld", 1)
c.add_variable("jets.pt", 2)
c.add_variable("jetPairs.mass", 1)
# ['ZJets_MG_HT_1200_2500', "/DYJetsToLL_M-50_HT-1200to2500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/ZJets_MG_HT_1200_2500/180913_192208/0000/tuple_*root", 0.96*0.186222 ],
# ['ZJets_MG_HT_2500_inf', "/DYJetsToLL_M-50_HT-2500toInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/ZJets_MG_HT_2500_inf/180913_192230/0000/tuple_*root", 0.96*0.004385 ],
# ['ZJets_hiM', "/DYJetsToLL_M-100to200_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/ZJets_hiM/180913_191802/0000/tuple_*root", 5765.4*1.235 ]
]
for s in signal:
c.add_signal(s[0], mc_path+s[1], s[2])
for b in background:
c.add_background(b[0], bkg_path+b[1], b[2])
c.set_tree_path(treePath)
c.add_variable("muPairs.pt", 1) #second argument is the number of objects considered
c.add_variable("muPairs.eta", 1)
c.add_variable("muPairs.phi", 1)
c.add_variable("muPairs.dEta", 1)
c.add_variable("muPairs.dPhi", 1)
c.add_variable("muons.eta", 2)
c.add_variable("met.pt", 1)
c.add_variable("nJets", 1)
c.add_variable("nJetsCent", 1)
# c.add_variable("nJetsFwd", 1)
c.add_variable("nBMed", 1)
# c.add_variable("jets.eta", 2)
c.add_variable("jetPairs.dEta", 2)
c.add_variable("jetPairs.mass", 2)
# c.add_variable("muons.pt", 2) #two muons
['ttZ', 246.881972], ['tt_ll_AMC', 97340.321489],
['tt_ll_MG_1', 98109.617801], ['tt_ll_MG_2', 97747.879187]]
for s in signal:
c.add_signal(s[0], inFileDir + s[0] + ".root", s[1])
for b in background:
c.add_background(b[0], inFileDir + b[0] + ".root", b[1])
# c.add_signal('H2Mu_gg', inFileDir+"H2Mu_gg.root", 0.006343) # label, path, weight
# c.add_signal('H2Mu_VBF', inFileDir+"H2Mu_VBF.root", 0.000495)
# c.add_background('dy', inFileDir+"ZJets_AMC.root", 29.853717)
c.set_tree_path(treePath)
c.add_variable("muPairs.pt",
1) #second argument is the number of objects considered
c.add_variable("muPairs.eta", 1)
c.add_variable("muPairs.dEta", 1) #for example, one muon pair
c.add_variable("muPairs.dPhi", 1)
c.add_variable("met.pt", 1)
c.add_variable("nJetsCent", 1)
c.add_variable("nJetsFwd", 1)
c.add_variable("nBMed", 1)
# c.add_variable("muons.pt", 2) #two muons
# c.add_variable("muons.eta", 2)
# c.add_variable("muons.phi", 2)
c.add_package("TMVA")
# c.add_transf("N,G,P")
# c.add_package("Keras")
c.add_mc_to_evaluate(
'tt_ll_POW',
'/mnt/hadoop/store/user/dkondrat/TTTo2L2Nu_TuneCP5_PSweights_13TeV-powheg-pythia8/tt_ll_POW/190322_195819/0000/*root',
85.656)
c.add_mc_to_evaluate(
'H2Mu_gg',
'/mnt/hadoop/store/user/dkondrat/GluGluHToMuMu_M-125_TuneCP5_PSweights_13TeV_powheg_pythia8/H2Mu_gg_powheg/190323_193526/0000/*root',
0.009618)
c.add_mc_to_evaluate(
'H2Mu_VBF',
'/mnt/hadoop/store/user/dkondrat/VBFHToMuMu_M-125_TuneCP5_PSweights_13TeV_powheg_pythia8/H2Mu_VBF_powheg/190322_195711/0000/*root',
0.0008208)
c.set_tree_path(treePath)
c.add_variable("muPairs.pt", 1)
c.add_variable("muPairs.eta", 1)
c.add_variable("muPairs.dEta", 1)
c.add_variable("muPairs.dPhi", 1)
c.add_variable("met.pt", 1)
c.add_variable("nJetsCent", 1)
c.add_variable("nJetsFwd", 1)
c.add_variable("nBMed", 1)
c.add_variable("jets.eta", 2)
c.add_variable("jetPairs.dEta", 1)
c.add_variable("jetPairs.mass", 1)
c.set_year("2017")
c.weigh_by_event(True)
c.add_package("TMVA")
c.add_data_to_evaluate(SingleMu2017F.name, SingleMu2017F.path)
c.add_mc_dir_to_evaluate(ZJets_aMC_2017_hiStat.name,
ZJets_aMC_2017_hiStat.path,
ZJets_aMC_2017_hiStat.xSec)
# c.add_mc_dir_to_evaluate(ZJets_aMC_2017_hiStat.name, ZJets_aMC_2017_hiStat.path, ZJets_aMC_2017_hiStat.xSec)
c.add_mc_dir_to_evaluate(tt_ll_POW_2017.name, tt_ll_POW_2017.path,
tt_ll_POW_2017.xSec)
c.add_mc_dir_to_evaluate(ggH_2017_powheg.name, ggH_2017_powheg.path,
ggH_2017_powheg.xSec)
c.add_mc_dir_to_evaluate(VBF_2017_powheg.name, VBF_2017_powheg.path,
VBF_2017_powheg.xSec)
c.set_tree_path(treePath)
c.add_variable("muPairs.pt", 1)
c.add_variable("muPairs.eta", 1)
c.add_variable("muPairs.cosThetaCS", 1)
c.add_variable("muPairs.phiCS", 1)
c.add_variable("jets.pt", 2)
c.add_variable("jets.eta", 1)
c.add_variable("jetPairs.dEta", 1)
c.add_variable("jetPairs.dPhi", 1)
c.add_variable("jetPairs.mass", 1)
c.add_variable("met.pt", 1)
c.add_variable("zeppenfeld", 1)
c.add_variable("nBMed", 1)
c.add_variable("nJets", 1)
c.add_variable("min_dR_mu_jet", 1)
c.add_variable("mu1_pt_Roch_over_mass", 1)
c.add_variable("mu2_pt_Roch_over_mass", 1)
c.add_category(c.vbf_label, True)
c.add_dir_to_category(VBF_2017_powheg.name, VBF_2017_powheg.path,
VBF_2017_powheg.xSec, c.vbf_label)
c.add_category(c.dy_label, False)
c.add_dir_to_category(ZJets_aMC_2017_hiStat.name, ZJets_aMC_2017_hiStat.path,
ZJets_aMC_2017_hiStat.xSec, c.dy_label)
c.add_category(c.tt_label, False)
c.add_dir_to_category(tt_ll_POW_2017.name, tt_ll_POW_2017.path,
tt_ll_POW_2017.xSec, c.tt_label)
##########################################################
### ------ Raffaele's variables ------ ###
c.add_variable("muPairs.pt", 1)
c.add_variable("muPairs.eta", 1)
c.add_variable("muPairs.cosThetaCS", 1)
c.add_variable("muPairs.phiCS", 1)
c.add_variable("met.pt", 1)
# c.add_variable("mu1_pt_Roch_over_mass", 1)
# c.add_variable("mu2_pt_Roch_over_mass", 1)
# c.add_variable('muons.eta', 2)
c.add_variable("min_dR_mu_jet", 1)
c.add_variable("nJets", 1)
c.add_variable("nBMed", 1)
c.add_variable("zeppenfeld", 1)
c.add_variable("jets.pt", 2)
c.add_variable("jetPairs.mass", 1)
# ['ZJets_MG_HT_1200_2500', "/DYJetsToLL_M-50_HT-1200to2500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/ZJets_MG_HT_1200_2500/180913_192208/0000/", 0.96*0.186222 ],
# ['ZJets_MG_HT_2500_inf', "/DYJetsToLL_M-50_HT-2500toInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/ZJets_MG_HT_2500_inf/180913_192230/0000/", 0.96*0.004385 ],
# ['ZJets_hiM', "/DYJetsToLL_M-100to200_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/ZJets_hiM/180913_191802/0000/", 5765.4*1.235 ]
]
for s in signal:
c.add_signal_dir(s[0], mc_path+s[1], s[2])
for b in background:
c.add_background_dir(b[0], bkg_path+b[1], b[2])
c.set_tree_path(treePath)
c.add_variable("muPairs.pt", 1)
c.add_variable("muPairs.eta", 1)
c.add_variable("muPairs.dEta", 1)
c.add_variable("muPairs.dPhi", 1)
c.add_variable("met.pt", 1)
c.add_variable("nJetsCent", 1)
c.add_variable("nJetsFwd", 1)
c.add_variable("nBMed", 1)
c.add_variable("jets.eta", 2)
c.add_variable("jetPairs.dEta", 1)
c.add_variable("jetPairs.mass", 1)
c.add_variable("min_dR_mu_jet" , 1)
c.add_variable("max_dR_mu_jet" , 1)
c.add_variable("min_dR_mumu_jet", 1)
c.add_variable("max_dR_mumu_jet", 1)
c.add_variable("zepenfeld", 1)
c.add_category(c.dy_label, False)
c.add_dir_to_category(DY_ucsd_2016.name, DY_ucsd_2016.path, DY_ucsd_2016.xSec, c.dy_label)
c.add_dir_to_category(DY_ucsd_2017.name, DY_ucsd_2017.path, DY_ucsd_2017.xSec, c.dy_label)
c.add_dir_to_category(DY_ucsd_2018.name, DY_ucsd_2018.path, DY_ucsd_2018.xSec, c.dy_label)
c.add_category(c.tt_label, False)
c.add_dir_to_category(tt_ucsd_2016.name, tt_ucsd_2016.path, tt_ucsd_2016.xSec, c.tt_label)
c.add_dir_to_category(tt_ucsd_2017.name, tt_ucsd_2017.path, tt_ucsd_2017.xSec, c.tt_label)
c.add_dir_to_category(tt_ucsd_2018.name, tt_ucsd_2018.path, tt_ucsd_2018.xSec, c.tt_label)
##########################################################
### ------ Raffaele's variables ------ ###
c.add_variable("hmmpt", 1)
c.add_variable("hmmrap", 1)
c.add_variable("hmmthetacs", 1)
c.add_variable("hmmphics", 1)
c.add_variable("met", 1)
c.add_variable("m1ptOverMass", 1)
c.add_variable("m2ptOverMass", 1)
c.add_variable('m1eta', 1)
c.add_variable('m2eta', 1)
c.add_variable("drmj", 1)
c.add_variable("njets", 1)
c.add_variable("nbjets", 1)
c.add_variable("zepen", 1)
c.add_variable("j1pt", 1)