Example #1
0
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?')
Example #2
0
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]
Example #3
0
import JPsi.MuMu.common.energyScaleChains as esChains
from ROOT import *

#ebselection = 'r9 > 0.94'
#eeselection = 'r9 > 0.95'

#ebselection = 'r9 < 0.94'
#eeselection = 'r9 < 0.95'

ebselection = '1'
eeselection = '1'

gROOT.LoadMacro("tdrstyle.C")
from ROOT import setTDRStyle

setTDRStyle()
#gStyle.SetPadRightMargin(0.05)
gStyle.SetOptFit(1111)
gStyle.SetOptTitle(1)
gStyle.SetOptStat(1111)
#gStyle.SetStatFontSize(0.05)
gStyle.SetStatW(0.2)
gStyle.SetStatH(0.2)

canvases = []
report = []

oplus = lambda x, y: math.sqrt(x * x + y * y)


def customizeStats(stats):
Example #4
0
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?')
Example #5
0
def prob(fin,name,var="pt",profile=True,rebin=4):

    prb = fin.Get("vtxprob_%s_cat0_tot" % var)

    if rebin > 0:
        prb.RebinX(rebin)
    
    if profile:
        return prb.ProfileX(name)

    return prb.Clone(name)
    
gROOT.LoadMacro("../Macros/rootglobestyle.C")
from ROOT import setTDRStyle, gStyle
setTDRStyle()
gStyle.SetOptTitle(0)

## file="histograms_vtxeff.root"
file1=argv[1]
## file2=argv[2]

fin1=TFile.Open(file1)
## fin2=TFile.Open(file2)


fout = TFile.Open("vertex_efficiency.root","recreate")
fout.cd()
rescaleEff(fin1).Write("efficiency")
fout.Close()