signals = ["BulkGWW"]
for sig in signals:
for p in purities:
cat = '_'.join(['JJ', sig, p, '13TeV'])
card = DataCardMaker('', p, '13TeV', 35900, 'JJ', cat)
cmd = cmd + " " + cat + '=datacard_' + cat + '.txt '
#SIGNAL
card.addMVVSignalParametricShape2("%s_MVV" % sig, "MJJ",
"JJ_%s_MVV.json" % sig,
{'CMS_scale_j': 1},
{'CMS_res_j': 1.0})
card.addMJJSignalParametricShapeNOEXP(
"Wqq1", "MJ1", "JJ_%s_MJl1_" % sig + p + ".json",
{'CMS_scale_prunedj': 1}, {'CMS_res_prunedj': 1.0})
card.addMJJSignalParametricShapeNOEXP(
"Wqq2", "MJ2", "JJ_%s_MJl1_" % sig + p + ".json",
{'CMS_scale_prunedj': 1}, {'CMS_res_prunedj': 1.0})
card.addParametricYieldWithUncertainty(
"%s" % sig, 0, "JJ_%s_" % sig + p + "_yield.json", 1,
'CMS_tau21_PtDependence', 'log(MH/600)', 0.041)
card.product3D("%s" % sig, "Wqq1", "Wqq2", "%s_MVV" % sig)
#Vjets
if p == 'HPHP':
from JJ_VJets_HPHP import JJ_VJets__MVV, JJ_VJets__Res_l1, JJ_VJets__ratio_l1, JJ_VJets__Res_l2, JJ_VJets__ratio_l2
if p == 'HPLP':
from JJ_VJets_HPLP import JJ_VJets__MVV, JJ_VJets__Res_l1, JJ_VJets__ratio_l1, JJ_VJets__Res_l2, JJ_VJets__ratio_l2
if p == 'LPLP':
{'CMS_scale_prunedj': 1},
{'CMS_res_prunedj': 1.0})
card.addMJJSignalParametricShape("Zqq", "MJ",
"JJ_XqZ_MJJ_" + p + ".json",
{'CMS_scale_prunedj': 1},
{'CMS_res_prunedj': 1.0})
card.addParametricYieldWithUncertainty(
"XqW", 0, "JJ_XqW_" + p + "_yield.json", 1,
'CMS_tau21_PtDependence', '((0.054/0.041)*(-log(MH/600)))', 0.041)
card.addParametricYieldWithUncertainty(
"XqZ", 0, "JJ_XqZ_" + p + "_yield.json", 1,
'CMS_tau21_PtDependence', '((0.054/0.041)*(-log(MH/600)))', 0.041)
else:
card.addMJJSignalParametricShapeNOEXP("Wqq", "MJ",
"JJ_XqW_MJJ_" + p + ".json",
{'CMS_scale_prunedj': 1},
{'CMS_res_prunedj': 1.0}) #
card.addMJJSignalParametricShapeNOEXP("Zqq", "MJ",
"JJ_XqZ_MJJ_" + p + ".json",
{'CMS_scale_prunedj': 1},
{'CMS_res_prunedj': 1.0}) #
card.addParametricYieldWithUncertainty("XqW", 0,
"JJ_XqW_" + p + "_yield.json",
1, 'CMS_tau21_PtDependence',
'log(MH/600)', 0.041)
card.addParametricYieldWithUncertainty("XqZ", 0,
"JJ_XqZ_" + p + "_yield.json",
1, 'CMS_tau21_PtDependence',
'log(MH/600)', 0.041)
#WW signal-MVV
card.addMVVSignalParametricShape(
"XWW_MVV", "MLNuJ", "LNuJJ_XWW_MVV_" + lepton + ".json", {
'CMS_scale_j': 1,
'CMS_scale_MET': 1.0
}, {
'CMS_res_j': 1.0,
'CMS_res_MET': 1.0
})
if purity == 'LP':
card.addMJJSignalParametricShape(
"Wqq", "MJ", "LNuJJ_XWW_MJJ_" + purity + ".json",
{'CMS_scale_prunedj': 1}, {'CMS_res_prunedj': 1.0})
else:
card.addMJJSignalParametricShapeNOEXP(
"Wqq", "MJ", "LNuJJ_XWW_MJJ_" + purity + ".json",
{'CMS_scale_prunedj': 1}, {'CMS_res_prunedj': 1.0}) #
card.product("XWW", "Wqq", "XWW_MVV")
#WZ signal-MVV
# card.addMVVSignalParametricShape("XWZ_MVV","MLNuJ","LNuJJ_XWZ_MVV_"+lepton+".json",{'CMS_scale_j':1,'CMS_scale_MET':1.0},{'CMS_res_j':1.0,'CMS_res_MET':1.0})
# if purity=='LP':
# card.addMJJSignalParametricShape("Zqq","MJ","LNuJJ_XWZ_MJJ_"+purity+".json",{'CMS_scale_prunedj':1},{'CMS_res_prunedj':1.0})
# else:
# card.addMJJSignalParametricShapeNOEXP("Zqq","MJ","LNuJJ_XWZ_MJJ_"+purity+".json",{'CMS_scale_prunedj':1},{'CMS_res_prunedj':1.0})
# card.product("XWZ","Zqq","XWZ_MVV")
#WH signal-MVV
# card.addMVVSignalParametricShape("XWH_MVV","MLNuJ","LNuJJ_XWH_MVV_"+lepton+".json",{'CMS_scale_j':1,'CMS_scale_MET':1.0},{'CMS_res_j':1.0,'CMS_res_MET':1.0})
# if purity=='LP':
'CMS_scale_MET_' + YEAR: 1.0,
'CMS_scale_' + LYtag: 1.0
}, {
'CMS_res_j_' + YEAR: 1.0,
'CMS_res_MET_' + YEAR: 1.0
})
if purity == 'LP':
card.addMJJSignalParametricShape(
sig + "_MJJ", varMJJ,
inputDir + "LNuJJ_" + sig + "_MJJ_" + PCEtag + ".json",
{'CMS_scale_prunedj_WPeak_' + PYtag: '0.05'},
{'CMS_res_prunedj_WPeak_' + PYtag: '0.25'})
else:
card.addMJJSignalParametricShapeNOEXP(
sig + "_MJJ", varMJJ,
inputDir + "LNuJJ_" + sig + "_MJJ_" + PCEtag + ".json",
{'CMS_scale_prunedj_WPeak_' + PYtag: '0.05'},
{'CMS_res_prunedj_WPeak_' + PYtag: '0.25'})
card.product(sig, sig + "_MJJ", sig + "_MVV")
if purity == 'HP':
card.addParametricYieldWithUncertainty(
sig, 0, inputDir + "LNuJJ_" + sig + "_" + LPCEtag +
"_yield.json", 1, 'CMS_Vtag_PtDependence_' + YEAR,
'(4.95e-5*(MH-650))', 1.)
else:
card.addParametricYieldWithUncertainty(
sig, 0, inputDir + "LNuJJ_" + sig + "_" + LPCEtag +
"_yield.json", 1, 'CMS_Vtag_PtDependence_' + YEAR,
'(-3.54e-5*(MH-650))', 1.)