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?')