'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':
# card.addMJJSignalParametricShape("Hqq","MJ","LNuJJ_XWH_MJJ_"+purity+".json",{'CMS_scale_prunedj':1},{'CMS_res_prunedj':1.0})
# else:
# card.addMJJSignalParametricShapeNOEXP("Hqq","MJ","LNuJJ_XWH_MJJ_"+purity+".json",{'CMS_scale_prunedj':1},{'CMS_res_prunedj':1.0})
{'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)
card.product("XqW", "Wqq", "XqW_MVV")
card.product("XqZ", "Zqq", "XqZ_MVV")
#QCD
rootFile = "JJ_nonRes_2D_" + p + ".root"
qcdTag = "_".join([p])
card.addHistoShapeFromFile("nonRes", ["MJJ", "MJ"], rootFile, "histo", [
'PTX:CMS_VV_JJ_nonRes_PTX_' + qcdTag, 'OPTX:CMS_VV_JJ_nonRes_OPTX_' +
qcdTag, 'OPTY:CMS_VV_JJ_nonRes_OPTY_' + qcdTag,
'PTY:CMS_VV_JJ_nonRes_PTY_' + qcdTag
], False, 0)
card.addFixedYieldFromFile("nonRes", 1, "JJ_" + p + ".root", "nonRes")
#DATA
card.importBinnedData("JJ_" + p + ".root", "data", ["MJJ", "MJ"])
card=DataCardMaker(lepton,purity,'13TeV',7650,category)
cat='_'.join([category,lepton,purity,'13TeV'])
cmd=cmd+" "+cat+'=datacard_'+cat+'.txt '
#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})
# card.addMJJSignalParametricShape("Wqq","MJ","LNuJJ_XWW_MJJ_"+purity+".json",{'CMS_scale_prunedj_0':1,'CMS_scale_prunedj_1':'MH'},{'CMS_res_prunedj_0':1.0,'CMS_res_prunedj_1':'MH'})
# card.product("XWW","Wqq","XWW_MVV")
# card.addParametricYield("XWW",0,"LNuJJ_XWW_"+lepton+"_"+purity+"_"+category+"_yield.json")
# card.addParametricYieldWithCrossSection("XWW",0,"LNuJJ_XWW_"+lepton+"_"+purity+"_"+category+"_yield.json",'sigma_hvt.json','sigma0','BRWW')
#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})
card.addMJJSignalParametricShape("Zqq","MJ","LNuJJ_XWZ_MJJ_"+purity+".json",{'CMS_scale_prunedj_0':1,'CMS_scale_prunedj_1':'MH'},{'CMS_res_prunedj_0':1.0,'CMS_res_prunedj_1':'MH'})
card.product("XWZ","Zqq","XWZ_MVV")
card.addParametricYield("XWZ",0,"LNuJJ_XWZ_"+lepton+"_"+purity+"_"+category+"_yield.json")
#W+jets
rootFile="LNuJJ_MVVHist_Wjets_"+lepton+"_"+purity+".root"
card.addHistoShapeFromFile("Wjets",["MLNuJ","MJ"],rootFile,"histo",["slopeSyst_Wjets_"+lepton,"widthSyst_Wjets_"+lepton,"meanSyst0_Wjets_"+lepton,"meanSyst1_Wjets_"+lepton,"slopeSystMJJ_Wjets_"+purity,"widthSystMJJ_Wjets_"+purity,"meanSystMJJ_Wjets_"+purity],False,2)
card.addFixedYieldFromFile("Wjets",1,"LNuJJ_"+lepton+"_"+purity+"_"+category+".root","Wjets")
#TOP RES
card.addMJJSignalParametricShapeCB("WqqTop","MJ","LNuJJ_MJJ_topRes_"+purity+".json",{'CMS_scale_prunedj_0':1,'CMS_scale_prunedj_1':'MLNuJ'},{'CMS_res_prunedj_0':1.0,'CMS_res_prunedj_1':'MLNuJ'},"MLNuJ")
rootFile="LNuJJ_MVVHist_topRes_"+lepton+"_"+category+".root"
card.addHistoShapeFromFile("topRes_MVV",["MLNuJ"],rootFile,"histo",["slopeSyst_topRes_"+lepton+"_"+category,"widthSyst_topRes_"+lepton+"_"+category,"meanSyst_topRes_"+lepton+"_"+category],False,2)
'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.)
## Non-resonant bkgd
rootFile = inputDir + "LNuJJ_nonRes" + sfx_rgn + "_2D_" + LPCEtag + ".root"
card.addHistoShapeFromFile(