def makeCombineStyleSeparationPlot(ifile, leg, ofile, rebin, xlow, xhigh):
print 'Plotting combine style separartion'
setTDRStyle()
extractSignificanceStats(options.unblind, leg,
options.outDir + '/' + ofile, ifile, rebin, xlow,
xhigh)
if not options.isBatch: raw_input('Looks ok?')
def makeCombineStyleSeparationPlot(ifile,leg,ofile,rebin,xlow,xhigh):
print 'Plotting combine style separartion'
setTDRStyle()
expCLs = r.Double(0)
obsCLs = r.Double(0)
extractSignificanceStats(expCLs,obsCLs,options.unblind,leg,options.outDir+'/'+ofile,ifile,rebin,xlow,xhigh,options.sqrtS)
if not options.isBatch: raw_input('Looks ok?')
return [expCLs,obsCLs]
range[1.00] = [-8.,8.]
elif options.sqrtS=='comb':
range[0.00] = [-15.,15.]
range[0.25] = [-12.,12.]
range[0.50] = [-10.,10.]
range[0.75] = [-10.,10.]
range[1.00] = [-12.,12.]
else:
sys.exit('Invalid sqrtS'+options.sqrtS)
for fqq in [0.0,0.25,0.5,0.75,1.0]:
fname = options.dir+'/testStat_fqq%4.2f.root'%fqq
if not os.path.exists(fname): fname = options.dir+'/qmu_qqbar%4.2f.root'%fqq
legname = '2^{+}_{m}(%3.0f%% q#bar{q})'%(fqq*100.)
altname = '2pm%4.2f'%fqq
expCLs = r.Double(0)
expSigSM = r.Double(0)
expSigALT = r.Double(0)
obsCLs = r.Double(0)
obsSigSM = r.Double(0)
obsSigALT = r.Double(0)
outname = options.dir+'/'+altname+'_blind'
if options.unblind: outname = outname.replace('blind','unblind')
extractSignificanceStats(expCLs,obsCLs,expSigSM,expSigALT,obsSigSM,obsSigALT,options.unblind,legname,outname,fname,options.rebin,range[fqq][0],range[fqq][1],options.sqrtS)
cls[fqq] = [expCLs,obsCLs]
sigma[fqq] = [expSigSM,expSigALT,obsSigSM,obsSigALT]
print 'fqq - 1-expCLs - expP(SM<ALT) - expP(ALT>SM) - 1-obsCLs - obsP(SM<Obs) - obsP(ALT>Obs)'
for fqq in [0.0,0.25,0.5,0.75,1.00]:
print '%4.2f - %5.3f - %4.2f - %4.2f - %5.3f - %4.2f - %4.2f'%(fqq,1-cls[fqq][0],r.RooStats.PValueToSignificance(sigma[fqq][0]),r.RooStats.PValueToSignificance(sigma[fqq][1]),1-cls[fqq][1],r.RooStats.PValueToSignificance(sigma[fqq][2]),r.RooStats.PValueToSignificance(sigma[fqq][3]))
def makeCombineStyleSeparationPlot(ifile,leg,ofile,xlow,xhigh):
print 'Plotting combine style separartion'
setTDRStyle()
extractSignificanceStats(options.unblind,leg,options.outDir+'/'+ofile,ifile,xlow,xhigh)
if not options.isBatch: raw_input('Looks ok?')