(options,args) = parser.parse_args()
YEAR=options.year
inputDir='Inputs_'+YEAR+'/'
intlumi=35920+41530+59740
inCR = True
sfx_rgn = "CR" if inCR else ""
sfx_year = "_"+YEAR if YEAR!="Run2" else ""
for lepton in ['allL']: #['e','mu']:
for purity in ['HP','LP','NP']:
for category in ['allC']: #['bb','nobb']:
card=DataCardMaker(lepton,purity,YEAR,intlumi,category)
cat='_'.join([category,lepton,purity,YEAR])
cmd=cmd+" "+cat+'=datacard_'+cat+'.txt '
varMJJ = "MJ"
PCtag = "_".join([purity,category])
PYtag = "_".join([purity,YEAR])
LPCtag = "_".join([lepton,purity,category])
PCYtag = "_".join([purity,category,YEAR])
LPCYtag = "_".join([lepton,purity,category,YEAR])
rootFile=inputDir+"LNuJJ_res"+sfx_rgn+"_MJJ_"+LPCtag+".root"
import ROOT
ROOT.gSystem.Load("libHiggsAnalysisCombinedLimit")
from CMGTools.VVResonances.statistics.DataCardMaker import DataCardMaker
cmd = 'combineCards.py '
for purity in ['HP', 'LP']:
for category in ['WW', 'WZ', 'ZZ']:
card = DataCardMaker('dijet', purity, '13TeV', 2630, category)
cat = '_'.join([category, 'dijet', purity, '13TeV'])
cmd = cmd + " " + cat + '=datacard_' + cat + '.txt '
#WW signal-MVV
card.addMVVSignalParametricShape("XWW", "MVV", "JJ_XWW.json", {
'CMS_scale_j': 1,
'CMS_scale_MET': 1.0
}, {
'CMS_res_j': 1.0,
'CMS_res_MET': 1.0
})
card.addParametricYield(
"XWW", 0, "JJ_XWW_" + purity + "_" + category + "_yield.json")
#QCD function
card.addMVVBackgroundShapeQCD("QCD", "MVV", logTerm=False)
card.addFloatingYield("QCD", 1, 1000, mini=0, maxi=1e+9)
card.importBinnedData("JJ_" + purity + "_" + category + ".root",
"data", ["MVV"])
#####
#####SYSTEMATICS
import sys
import ROOT
ROOT.gSystem.Load("libHiggsAnalysisCombinedLimit")
from CMGTools.VVResonances.statistics.DataCardMaker import DataCardMaker
cmd = 'combineCards.py '
purities = ['HPHP', 'HPLP']
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)
import ROOT
ROOT.gSystem.Load("libHiggsAnalysisCombinedLimit")
from CMGTools.VVResonances.statistics.DataCardMaker import DataCardMaker
cmd='combineCards.py '
for category in ['nob','b']:
if category=='vbf':
pur=['NP']
else:
pur=['HP','LP']
for purity in pur:
for lepton in ['e','mu']:
card=DataCardMaker(lepton,purity,'13TeV',35900,category)
cat='_'.join([category,lepton,purity,'13TeV'])
cmd=cmd+" "+cat+'=datacard_'+cat+'.txt '
#WW signal-MVV
card.addMVVSignalParametricShape("XWW","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.addParametricYieldWithUncertainty("XWW",0,"LNuJJ_XWW_"+lepton+"_"+purity+"_"+category+"_yield.json",1,'CMS_tau21_PtDependence','log(MH/600)',0.041)
# card.addMVVSignalParametricShape("XWZ","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.addParametricYieldWithUncertainty("XWZ",0,"LNuJJ_XWZ_"+lepton+"_"+purity+"_"+category+"_yield.json",1,'CMS_tau21_PtDependence','log(MH/600)',0.041)
rootFile="LNuJJ_"+lepton+"_"+purity+"_"+category+".root"
qcdTag ="_".join([lepton,purity,category])
card.addHistoShapeFromFile("nonRes",["MLNuJ"],rootFile,"nonRes",["shape_nonRes:shape_nonRes"+qcdTag],False,0)
# card.addHistoShapeFromFile("nonRes",["MLNuJ"],rootFile,"nonRes",[],False,0)
card.addFixedYieldFromFile("nonRes",1,"LNuJJ_"+lepton+"_"+purity+"_"+category+".root","nonRes")
import ROOT
ROOT.gSystem.Load("libHiggsAnalysisCombinedLimit")
from CMGTools.VVResonances.statistics.DataCardMaker import DataCardMaker
cmd = 'combineCards.py '
for category in ['nob']:
if category == 'vbf':
pur = ['NP']
else:
pur = ['HP', 'LP']
for purity in pur:
for lepton in ['e', 'mu']:
card = DataCardMaker(lepton, purity, '13TeV', 35900, 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
})
if purity == 'LP':
card.addMJJSignalParametricShape(
"Wqq", "MJ", "LNuJJ_XWW_MJJ_" + purity + ".json",
{'CMS_scale_prunedj': 1}, {'CMS_res_prunedj': 1.0})
else:
card.addMJJSignalParametricShapeNOEXP(
import ROOT
ROOT.gSystem.Load("libHiggsAnalysisCombinedLimit")
from CMGTools.VVResonances.statistics.DataCardMaker import DataCardMaker
cmd='combineCards.py '
for lepton in ['e','mu']:
for purity in ['HP','LP']:
for category in ['nob','b']:
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
import sys
import ROOT
ROOT.gSystem.Load("libHiggsAnalysisCombinedLimit")
from CMGTools.VVResonances.statistics.DataCardMaker import DataCardMaker
cmd = 'combineCards.py '
purities = ['HP', 'LP']
for p in purities:
card = DataCardMaker('', p, '13TeV', 35900, 'JJ')
cat = '_'.join(['JJ', p, '13TeV'])
cmd = cmd + " " + cat + '=datacard_' + cat + '.txt '
#SIGNAL
card.addMVVSignalParametricShape2("XqW_MVV", "MJJ", "JJ_XqW_MVV.json", {
'CMS_scale_j': 1,
'CMS_scale_MET': 1.0
}, {
'CMS_res_j': 1.0,
'CMS_res_MET': 1.0
})
card.addMVVSignalParametricShape2("XqZ_MVV", "MJJ", "JJ_XqZ_MVV.json", {
'CMS_scale_j': 1,
'CMS_scale_MET': 1.0
}, {
'CMS_res_j': 1.0,
'CMS_res_MET': 1.0
})
import ROOT
ROOT.gSystem.Load("libHiggsAnalysisCombinedLimit")
from CMGTools.VVResonances.statistics.DataCardMaker import DataCardMaker
cmd='combineCards.py '
for lepton in ['e','mu']:
for purity in ['HP','LP']:
for category in ['nob','b']:
card=DataCardMaker(lepton,purity,'13TeV',2630,category)
cat='_'.join([category,lepton,purity,'13TeV'])
cmd=cmd+" "+cat+'=datacard_'+cat+'.txt '
#WW signal-MVV
card.addMVVSignalParametricShape("XWW_MVV","MLNuJ","LNuJJ_XWW_"+lepton+"_"+purity+".json",{'CMS_scale_j':1,'CMS_scale_MET':1.0},{'CMS_res_j':1.0,'CMS_res_MET':1.0})
#WWsignal M-JJ
card.addMJJSignalParametricShape("Wqq","MJ","LNuJJ_XWW_"+lepton+"_"+purity+".json",{'CMS_scale_prunedj_0':1,'CMS_scale_prunedj_1':'MH'},{'CMS_res_prunedj_0':1.0,'CMS_res_prunedj_1':'MH'})
card.conditionalProduct("XWW","Wqq","MLNuJ","XWW_MVV")
#parametric yield for signal
card.addParametricYield("XWW",0,"LNuJJ_XWW_"+lepton+"_"+purity+"_"+category+"_yield.json")
#W+jets - quark
card.addMJJParametricBackgroundShapeErfExp("Wjets_quark_MJ","MJ","LNuJJ_MJJ_Wjets_quark_"+lepton+"_"+purity+".json",{'wMJJSlopeQ_'+lepton+"_"+purity:1.0},{'wMJJMeanQ_'+lepton+"_"+purity:1.0},{'wMJJSigmaQ_'+lepton+"_"+purity:1.0})
card.addParametricMVVBKGShapeErfPow("Wjets_quark_MVV","MLNuJ","MJ","LNuJJ_MVV_Wjets_quark_"+lepton+"_"+purity+".json","",{'wSlope_'+lepton+"_"+purity:1.0},{'wMean0_'+lepton+"_"+purity:1.0,'wMean1_'+lepton+"_"+purity:'MJ'},{'wSigma_'+lepton+"_"+purity:1.0})
card.conditionalProduct("Wjets_quark","Wjets_quark_MVV","MJ","Wjets_quark_MJ")
card.addFixedYieldFromFile("Wjets_quark",1,"LNuJJ_"+lepton+"_"+purity+"_"+category+".root","Wjets_quark")
sfx_year = "_" + YEAR if YEAR != "Run2" else ""
sig = options.signalType
if sig not in [
'GbuToWW', 'RadToWW', 'ZprToWW', 'WprToWZ', 'WprToWH', 'VBFGbuToWW',
'VBFRadToWW', 'VBFZprToWW', 'VBFWprToWZ', 'VBFWprToWH'
]:
sys.exit('Error: unrecognized signal')
for lepton in ['e', 'mu']:
for purity in ['HP', 'LP']:
for category in ['bb', 'nobb', 'vbf']:
for dy in ['LDy', 'HDy']:
card = DataCardMaker(lepton, purity, YEAR, intlumi,
category + "_" + dy)
cat = '_'.join([lepton, purity, category, dy, YEAR])
cmd = cmd + " " + cat + '=datacard_' + cat + '.txt '
varMVV = "MLNuJ"
varMJJ = "MJ"
if options.differentBinning and category in ['bb', 'vbf']:
varMVV = "MLNuJ_coarse"
varMJJ = "MJ_coarse"
PCEtag = "_".join([purity, category, dy])
LPCEtag = "_".join([lepton, purity, category, dy])
LYtag = "_".join([lepton, YEAR])
PYtag = "_".join([purity, YEAR])
PCEYtag = "_".join([purity, category, dy, YEAR])
LPCEYtag = "_".join([lepton, purity, category, dy, YEAR])