#ePlotter.addHistogram('HppHmm1.0TeV', sigMap['HppHmm1.0TeV'],style={'name':'#Phi^{++}#Phi^{--} (1 TeV)', 'linecolor':sigcolors[8], 'fillcolor':sigcolors[8]})
plots = {
'dxy': {
'legendpos': 34,
'numcol': 3,
'invert': True
},
'dz': {
'legendpos': 34,
'numcol': 3,
'invert': True
},
'mvaTrig': {
'legendpos': 34,
'numcol': 3,
'invert': False
},
}
for plot, kwargs in plots.iteritems():
sig = {
'HppHmm500GeVBarrel': '{0}_barrel'.format(plot),
'HppHmm500GeVEndcap': '{0}_endcap'.format(plot)
}
bg = {
'HppHmm500GeVBarrel': '{0}_barrel_fake'.format(plot),
'HppHmm500GeVEndcap': '{0}_endcap_fake'.format(plot)
}
ePlotter.plotROC(sig, bg, 'roc_{0}'.format(plot), **kwargs)
}
for mass in masses:
name = 'HppHm{0}GeV'.format(mass)
sigOrder += [name]
bgOrder += ['BG']
hpp3lPlotter.addHistogram(name,sigMap[name],signal=True,style={'linecolor': sigColors[mass]})
workingPoints['mllMinusMZ'][name] = {'Z Veto 5 GeV': 5, 'Z Veto 5 GeV': 10, 'Z Veto 5 GeV': 50}
workingPoints['hppMass'][name] = {'0.2*m_{#Phi}': 0.2*mass, '0.3*m_{#Phi}': 0.3*mass, '0.4*m_{#Phi}': 0.4*mass, '0.5*m_{#Phi}': 0.5*mass, '0.9*m_{#Phi}': 0.9*mass,}
workingPoints['hppMt'][name] = {'0.2*m_{#Phi}': 0.2*mass, '0.3*m_{#Phi}': 0.3*mass, '0.4*m_{#Phi}': 0.4*mass, '0.5*m_{#Phi}': 0.5*mass, '0.9*m_{#Phi}': 0.9*mass,}
workingPoints['st'][name] = {'0.2*m_{#Phi}': 0.2*mass, '0.4*m_{#Phi}': 0.4*mass, '0.6*m_{#Phi}': 0.6*mass, '0.8*m_{#Phi}': 0.8*mass, '1.0*m_{#Phi}': 1.0*mass,}
for plot in roc_cust:
kwargs = deepcopy(plots[plot])
kwargs.update(roc_cust[plot])
for higgsChan in ['ee','em','et','mm','mt','tt']:
plotnames = []
bgnames = []
for gen in genRecoMap:
if len(gen)==4: continue
for reco in genRecoMap[gen]:
if gen[:2]==higgsChan:
plotnames += ['default/{0}/gen_{1}/{2}'.format(reco,gen,plot)]
bgnames += ['default/{0}/{1}'.format(reco,plot)]
if plotnames:
savename = 'signal/{0}/{1}_genMatched_roc'.format(higgsChan,plot)
wp = workingPoints[plot] if plot in workingPoints else {}
hpp3lPlotter.plotROC(plotnames,bgnames,savename,sigOrder=sigOrder,bgOrder=bgOrder,workingPoints=wp,**kwargs)
ROOT.TColor.GetColor('#FFB2B2'),
ROOT.TColor.GetColor('#FFCCCC'),
]
#ePlotter.addHistogram('HppHmm500GeV',sigMap['HppHmm500GeV'])
ePlotter.addHistogram('HppHmm500GeVBarrel',sigMap['HppHmm500GeV'],style={'name':'#Phi^{++}#Phi^{--} (Barrel)','linecolor':sigcolors[3],'fillcolor':sigcolors[3]})
ePlotter.addHistogram('HppHmm500GeVEndcap',sigMap['HppHmm500GeV'],style={'name':'#Phi^{++}#Phi^{--} (Endcap)','linecolor':sigcolors[6],'fillcolor':sigcolors[6]})
#
#ePlotter.addHistogram('HppHmm200GeV', sigMap['HppHmm200GeV'],style={'name':'#Phi^{++}#Phi^{--} (200 GeV)', 'linecolor':sigcolors[0], 'fillcolor':sigcolors[0]})
#ePlotter.addHistogram('HppHmm300GeV', sigMap['HppHmm300GeV'],style={'name':'#Phi^{++}#Phi^{--} (300 GeV)', 'linecolor':sigcolors[1], 'fillcolor':sigcolors[1]})
#ePlotter.addHistogram('HppHmm400GeV', sigMap['HppHmm400GeV'],style={'name':'#Phi^{++}#Phi^{--} (400 GeV)', 'linecolor':sigcolors[2], 'fillcolor':sigcolors[2]})
#ePlotter.addHistogram('HppHmm500GeV', sigMap['HppHmm500GeV'],style={'name':'#Phi^{++}#Phi^{--} (500 GeV)', 'linecolor':sigcolors[3], 'fillcolor':sigcolors[3]})
#ePlotter.addHistogram('HppHmm600GeV', sigMap['HppHmm600GeV'],style={'name':'#Phi^{++}#Phi^{--} (600 GeV)', 'linecolor':sigcolors[4], 'fillcolor':sigcolors[4]})
#ePlotter.addHistogram('HppHmm700GeV', sigMap['HppHmm700GeV'],style={'name':'#Phi^{++}#Phi^{--} (700 GeV)', 'linecolor':sigcolors[5], 'fillcolor':sigcolors[5]})
#ePlotter.addHistogram('HppHmm800GeV', sigMap['HppHmm800GeV'],style={'name':'#Phi^{++}#Phi^{--} (800 GeV)', 'linecolor':sigcolors[6], 'fillcolor':sigcolors[6]})
#ePlotter.addHistogram('HppHmm900GeV', sigMap['HppHmm900GeV'],style={'name':'#Phi^{++}#Phi^{--} (900 GeV)', 'linecolor':sigcolors[7], 'fillcolor':sigcolors[7]})
#ePlotter.addHistogram('HppHmm1.0TeV', sigMap['HppHmm1.0TeV'],style={'name':'#Phi^{++}#Phi^{--} (1 TeV)', 'linecolor':sigcolors[8], 'fillcolor':sigcolors[8]})
plots = {
'dxy' : {'legendpos':34,'numcol':3,'invert':True},
'dz' : {'legendpos':34,'numcol':3,'invert':True},
'mvaTrig' : {'legendpos':34,'numcol':3,'invert':False},
}
for plot,kwargs in plots.iteritems():
sig = {'HppHmm500GeVBarrel':'{0}_barrel'.format(plot), 'HppHmm500GeVEndcap':'{0}_endcap'.format(plot)}
bg = {'HppHmm500GeVBarrel':'{0}_barrel_fake'.format(plot), 'HppHmm500GeVEndcap':'{0}_endcap_fake'.format(plot)}
ePlotter.plotROC(sig, bg, 'roc_{0}'.format(plot), **kwargs)
'st' : {},
}
for mass in masses:
name = 'HppHmm{0}GeV'.format(mass)
sigOrder += [name]
bgOrder += ['BG']
hpp4lPlotter.addHistogram(name,sigMap[name],signal=True,style={'linecolor': sigColors[mass]})
workingPoints['mllMinusMZ'][name] = {'Z Veto 5 GeV': 5, 'Z Veto 5 GeV': 10, 'Z Veto 5 GeV': 50}
workingPoints['hppMass'][name] = {'0.2*m_{#Phi}': 0.2*mass, '0.3*m_{#Phi}': 0.3*mass, '0.4*m_{#Phi}': 0.4*mass, '0.5*m_{#Phi}': 0.5*mass, '0.9*m_{#Phi}': 0.9*mass,}
workingPoints['hppMt'][name] = {'0.2*m_{#Phi}': 0.2*mass, '0.3*m_{#Phi}': 0.3*mass, '0.4*m_{#Phi}': 0.4*mass, '0.5*m_{#Phi}': 0.5*mass, '0.9*m_{#Phi}': 0.9*mass,}
workingPoints['st'][name] = {'0.2*m_{#Phi}': 0.2*mass, '0.4*m_{#Phi}': 0.4*mass, '0.6*m_{#Phi}': 0.6*mass, '0.8*m_{#Phi}': 0.8*mass, '1.0*m_{#Phi}': 1.0*mass,}
for plot in roc_cust:
kwargs = deepcopy(plots[plot])
kwargs.update(roc_cust[plot])
for higgsChan in ['ee','em','et','mm','mt','tt']:
plotnames = []
bgnames = []
for gen in genRecoMap:
for reco in genRecoMap[gen]:
if gen[:2]==higgsChan:
plotnames += ['default/{0}/gen_{1}/{2}'.format(reco,gen,plot)]
bgnames += ['default/{0}/{1}'.format(reco,plot)]
if gen[2:]==higgsChan and 'hpp' in plot:
plotnames += ['default/{0}/gen_{1}/{2}'.format(reco,gen,plot.replace('hpp','hmm'))]
if plotnames:
savename = 'signal/{0}/{1}_genMatched_roc'.format(higgsChan,plot)
wp = workingPoints[plot] if plot in workingPoints else {}
hpp4lPlotter.plotROC(plotnames,bgnames,savename,sigOrder=sigOrder,bgOrder=bgOrder,workingPoints=wp,**kwargs)