StudyXClean/rootScripts/MakeStudy.py

import ROOT
import ConfigParser, re
import sys
sys.path.append( '/afs/naf.desy.de/user/n/nowaf/CMSSW_3_3_6/src/LeptonSel/LeptonSel/rootScripts' )
import definitions as Def
from math import sqrt

def main(argv):
    #import sys

    basepath = ''

    try:
        options, arguments = getopt.getopt(argv, "hp:", ["help","path="])
        
    except getopt.GetoptError:
        usage()
        sys.exit(2)

    for opt, arg in options:
        if opt in ("-h", "--help"):
            usage()
            sys.exit(2)
        elif opt in ("-p", "--path"):
            basepath = arg

    if basepath == '':
        print "Error! Need path to .cfg file to continue."
        sys.exit(2)

    return basepath

def usage():

    print "Script to plot all XC study plots for the di lepton analysis"
    print "Options and arguments:"
    print "-h, --help\t\t: print this help message and exit"
    print "-p, --path\t\t: relative path to the .cfg file"

    pass


def readTheConfig(name):
    """
    read in the config file
    """
    config = ConfigParser.ConfigParser()
    config.read( name )

    return config

def getStudiedCut( name ):
    """
    gets the studied cut out of the section name
    """
    studied = ""
    if re.search( "EJ", name ) or re.search( "PJ", name ):
        studied = "cut on E_{shared}/E_{Jet}"
        if re.search( "dR", name ):
            studied = "max #delta R (non iso elec, jet)"
    elif re.search( "MJ", name ):
        studied = "#delta R (muon, jet)"
    
    return studied

def MakeDivision( h_1, h_2, option ):
    """
    divides two histograms, computes errors
    and returns the result
    """
    #h_1.Sumw2()
    #h_2.Sumw2()
    hist = h_2.Clone()
    hist.Divide( h_1, h_2, 1., 1., option )

    return hist

def Fill( hist, option, weight  ):
    """
    h.SetBinContent sets the content right, but has only
    one wheighted entry then. This gives problems in dividing the
    histograms. This fill function takes care of that
    """
    for i in range( 0, weight ):
        hist.Fill( float( option ) )

    return hist
    

if __name__ == "__main__":

    import getopt

    basepath = main(sys.argv[1:])


    Def.SetGlobalStyles( "emr" )

    saves = {}

    config = readTheConfig( basepath )

    for section in config.sections():

        print "section: ", section

        studied_cut = getStudiedCut( section )

        nbin = 9
        min = 0.05
        max = 0.95
        if( re.search( "EJ_dR", section ) ):
            min = 0.15
            max = 1.05
            pass
       
        h_mu_m    = ROOT.TH1F( "h_mu_m", "h_mu_m",nbin,min,max )
        h_mu_gen  = ROOT.TH1F( "h_mu_gen", "h_mu_gen",nbin,min,max )
        h_mu_nm   = ROOT.TH1F( "h_mu_nm", "h_mu_nm",nbin,min,max )
        h_mu_all  = ROOT.TH1F( "h_mu_all", "h_mu_all",nbin,min,max )
        h_el_m    = ROOT.TH1F( "h_el_m", "h_el_m",nbin,min,max )
        h_el_gen  = ROOT.TH1F( "h_el_gen", "h_el_gen",nbin,min,max )
        h_el_nm   = ROOT.TH1F( "h_el_nm", "h_el_nm",nbin,min,max )
        h_el_all  = ROOT.TH1F( "h_el_all", "h_el_all",nbin,min,max )
        h_ph_m    = ROOT.TH1F( "h_ph_m", "h_ph_m",nbin,min,max )
        h_ph_gen  = ROOT.TH1F( "h_ph_gen", "h_ph_gen",nbin,min,max )
        h_ph_nm   = ROOT.TH1F( "h_ph_nm", "h_ph_nm",nbin,min,max )
        h_ph_all  = ROOT.TH1F( "h_ph_all", "h_ph_all",nbin,min,max )
        h_je_m    = ROOT.TH1F( "h_je_m", "h_je_m",nbin,min,max )
        h_je_gen  = ROOT.TH1F( "h_je_gen", "h_je_gen",nbin,min,max )
        h_je_nm   = ROOT.TH1F( "h_je_nm", "h_je_nm",nbin,min,max )
        h_je_all  = ROOT.TH1F( "h_je_all", "h_je_all",nbin,min,max )

        h_mu_m.Sumw2()
        h_mu_gen.Sumw2()
        h_mu_nm.Sumw2()
        h_mu_all.Sumw2()
        h_el_m.Sumw2()
        h_el_gen.Sumw2()
        h_el_nm.Sumw2()
        h_el_all.Sumw2()
        h_ph_m.Sumw2()
        h_ph_gen.Sumw2()
        h_ph_nm.Sumw2()
        h_ph_all.Sumw2()
        h_je_m.Sumw2()
        h_je_gen.Sumw2()
        h_je_nm.Sumw2()
        h_je_all.Sumw2()  
        

        #print type( h_mu_eff)
        
        for option in config.options( section ):

            print "option: ", option

            if re.search( "afs", config.get( section, option ) ) or re.search( "current", config.get( section, option ) ):
                fin = Def.ReadInRootFile( config.get( section, option ) + "/studyXC.root" )
            else:
                fin = Def.ReadInRootFile( "./../test/" + config.get( section, option ) + "/studyXC.root" )

               
            Tree_ = fin.Get( 'StudyXC' )
            mu_all = Def.defineArray( 'i', 1 )
            mu_gen = Def.defineArray( 'i', 1 )
            mu_m   = Def.defineArray( 'i', 1 )
            mu_nm  = Def.defineArray( 'i', 1 )
            el_all = Def.defineArray( 'i', 1 )
            el_gen = Def.defineArray( 'i', 1 )
            el_m   = Def.defineArray( 'i', 1 )
            el_nm  = Def.defineArray( 'i', 1 )
            ph_all = Def.defineArray( 'i', 1 )
            ph_gen = Def.defineArray( 'i', 1 )
            ph_m   = Def.defineArray( 'i', 1 )
            ph_nm  = Def.defineArray( 'i', 1 )
            je_all = Def.defineArray( 'i', 1 )
            je_gen = Def.defineArray( 'i', 1 )
            je_m   = Def.defineArray( 'i', 1 )
            je_nm  = Def.defineArray( 'i', 1 )
            
            Tree_.SetBranchAddress( "MuonsAll", mu_all )
            Tree_.SetBranchAddress( "MuonsGen", mu_gen )
            Tree_.SetBranchAddress( "MuonsNotMatched", mu_nm )
            Tree_.SetBranchAddress( "MuonsMatched", mu_m )
            Tree_.SetBranchAddress( "ElectronsAll", el_all )
            Tree_.SetBranchAddress( "ElectronsGen", el_gen )
            Tree_.SetBranchAddress( "ElectronsNotMatched", el_nm )
            Tree_.SetBranchAddress( "ElectronsMatched", el_m )
            Tree_.SetBranchAddress( "PhotonsAll", ph_all )
            Tree_.SetBranchAddress( "PhotonsGen", ph_gen )
            Tree_.SetBranchAddress( "PhotonsNotMatched", ph_nm )
            Tree_.SetBranchAddress( "PhotonsMatched", ph_m )
            Tree_.SetBranchAddress( "JetsAll", je_all )
            Tree_.SetBranchAddress( "JetsGen", je_gen )
            Tree_.SetBranchAddress( "JetsNotMatched", je_nm )
            Tree_.SetBranchAddress( "JetsMatched", je_m )
            
            n_mu_all = 0
            n_mu_gen = 0
            n_mu_m   = 0
            n_mu_nm  = 0
            n_el_all = 0
            n_el_gen = 0
            n_el_m   = 0
            n_el_nm  = 0
            n_ph_all = 0
            n_ph_gen = 0
            n_ph_m   = 0
            n_ph_nm  = 0
            n_je_all = 0
            n_je_gen = 0
            n_je_m   = 0
            n_je_nm  = 0

            n_bytes = 0
            for ent in range( 0, Tree_.GetEntries() ):
                
                n_bytes += Tree_.GetEntry( ent )
                n_mu_all += mu_all[0]
                n_mu_gen += mu_gen[0]
                n_mu_m   += mu_m[0]
                n_mu_nm  += mu_nm[0]
                n_el_all += el_all[0]
                n_el_gen += el_gen[0]
                n_el_m   += el_m[0]
                n_el_nm  += el_nm[0]
                n_ph_all += ph_all[0]
                n_ph_gen += ph_gen[0]
                n_ph_m   += ph_m[0]
                n_ph_nm  += ph_nm[0]
                n_je_all += je_all[0]
                n_je_gen += je_gen[0]
                n_je_m   += je_m[0]
                n_je_nm  += je_nm[0]
                pass


            Fill( h_mu_m, option, n_mu_m ) 
            Fill( h_mu_gen, option, n_mu_gen )
            Fill( h_mu_nm, option, n_mu_nm )
            Fill( h_mu_all, option, n_mu_all )
            Fill( h_el_m, option, n_el_m ) 
            Fill( h_el_gen, option, n_el_gen )
            Fill( h_el_nm, option, n_el_nm )
            Fill( h_el_all, option, n_el_all )
            Fill( h_ph_m, option, n_ph_m ) 
            Fill( h_ph_gen, option, n_ph_gen )
            Fill( h_ph_nm, option, n_ph_nm )
            Fill( h_ph_all, option, n_ph_all )
            Fill( h_je_m, option, n_je_m ) 
            Fill( h_je_gen, option, n_je_gen )
            Fill( h_je_nm, option, n_je_nm )
            Fill( h_je_all, option, n_je_all )

            print "n_mu_all: ", n_mu_all
            print "n_mu_m  : ", n_mu_m
            print "n_mu_nm : ", n_mu_nm
            print "n_mu_gen: ", n_mu_gen
           
            pass

        h_mu_eff = MakeDivision( h_mu_m, h_mu_gen, "B" )
        h_mu_cont = MakeDivision( h_mu_nm, h_mu_all, "B" )
        h_el_eff = MakeDivision( h_el_m, h_el_gen, "B" )
        h_el_cont = MakeDivision( h_el_nm, h_el_all, "B" )
        h_ph_eff = MakeDivision( h_ph_m, h_ph_gen, "B" )
        h_ph_cont = MakeDivision( h_ph_nm, h_ph_all, "B" )
        h_je_eff = MakeDivision( h_je_m, h_je_gen, "B" )
        h_je_cont = MakeDivision( h_je_nm, h_je_all, "B" )

        #for i in range( 0, 9 ):
        #    print "error in bin ", i, ": ",  h_mu_eff.GetBinError( i )


        ### wrong errors here
        h_mu_ratio = MakeDivision( h_mu_eff, h_mu_cont, "" )
        h_el_ratio = MakeDivision( h_el_eff, h_el_cont, "" )
        h_ph_ratio = MakeDivision( h_ph_eff, h_ph_cont, "" )
        h_je_ratio = MakeDivision( h_je_eff, h_je_cont, "" )     
        

        ## draw the hists for each section
        c_mu = ROOT.TCanvas( "c_mu", section + "_Muon" ,1,71,300,600 )
        c_mu.Divide(1,2)
        c_mu.cd(1)
        Def.DrawSmall( c_mu, \
                       [h_mu_eff, h_mu_cont], [8,50], ["eff","cont"], \
                       studied_cut, "", "", \
                       False, 0, False, "ur", section + "muons", saves)

        c_mu.cd(2)
        Def.DrawSmall( c_mu, \
                       [h_mu_ratio], [4], ["ratio"], \
                       studied_cut, "ratio eff / sqrt(cont)", "", \
                       False, 0, False, "ur", section + "muons_ratio", saves)
        
        c_el = ROOT.TCanvas( "c_el", section + "_Elecs" ,1,71,300,600 )
        c_el.Divide(1,2)
        c_el.cd(1)
        Def.DrawSmall( c_el, \
                       [h_el_eff, h_el_cont], [8,50], ["eff","cont"], \
                       studied_cut, "", "", \
                       False, 0, False, "ur", section + "elecs", saves)

        c_el.cd(2)
        Def.DrawSmall( c_el, \
                       [h_el_ratio], [4], ["ratio"], \
                       studied_cut, "ratio eff / sqrt(cont)", "", \
                       False, 0, False, "ur", section + "elecs_ratio", saves)
        
        c_ph = ROOT.TCanvas( "c_ph", section + "_Photons" ,1,71,300,600 )
        c_ph.Divide(1,2)
        c_ph.cd(1)
        Def.DrawSmall( c_ph, \
                       [h_ph_eff, h_ph_cont], [8,50], ["eff","cont"], \
                       studied_cut, "", "", \
                       False, 0, False, "ur", section + "phots", saves)

        c_ph.cd(2)
        Def.DrawSmall( c_ph, \
                       [h_ph_ratio], [4], ["ratio"], \
                       studied_cut, "ratio eff / sqrt(cont)", "", \
                       False, 0, False, "ur", section + "phots_ratio", saves)
        
        c_je = ROOT.TCanvas( "c_je", section + "_Jets" ,1,71,300,600 )
        c_je.Divide(1,2)
        c_je.cd(1)
        Def.DrawSmall( c_je, \
                       [h_je_eff, h_je_cont], [8,50], ["eff","cont"], \
                       studied_cut, "", "", \
                       False, 0, False, "ur", section + "jets", saves)

        c_je.cd(2)
        Def.DrawSmall( c_je, \
                       [h_je_ratio], [4], ["ratio"], \
                       studied_cut, "ratio eff / sqrt(cont)", "", \
                       False, 0, False, "ur", section + "jets_ratio", saves)


        pass
    
    
    Def.DontQuit()
Revision:	1.2
Committed:	Thu Feb 25 06:27:10 2010 UTC (15 years, 2 months ago) by nowak
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	HEAD
Changes since 1.1:	+206 -81 lines
Log Message:	more stuff
#	User	Rev	Content
1	nowak	1.1	import ROOT
2			import ConfigParser, re
3			import sys
4			sys.path.append( '/afs/naf.desy.de/user/n/nowaf/CMSSW_3_3_6/src/LeptonSel/LeptonSel/rootScripts' )
5			import definitions as Def
6	nowak	1.2	from math import sqrt
7
8			def main(argv):
9			#import sys
10
11			basepath = ''
12
13			try:
14			options, arguments = getopt.getopt(argv, "hp:", ["help","path="])
15
16			except getopt.GetoptError:
17			usage()
18			sys.exit(2)
19
20			for opt, arg in options:
21			if opt in ("-h", "--help"):
22			usage()
23			sys.exit(2)
24			elif opt in ("-p", "--path"):
25			basepath = arg
26
27			if basepath == '':
28			print "Error! Need path to .cfg file to continue."
29			sys.exit(2)
30
31			return basepath
32
33			def usage():
34
35			print "Script to plot all XC study plots for the di lepton analysis"
36			print "Options and arguments:"
37			print "-h, --help\t\t: print this help message and exit"
38			print "-p, --path\t\t: relative path to the .cfg file"
39
40			pass
41
42
43	nowak	1.1
44			def readTheConfig(name):
45			"""
46			read in the config file
47			"""
48			config = ConfigParser.ConfigParser()
49			config.read( name )
50
51			return config
52
53			def getStudiedCut( name ):
54			"""
55			gets the studied cut out of the section name
56			"""
57			studied = ""
58			if re.search( "EJ", name ) or re.search( "PJ", name ):
59			studied = "cut on E_{shared}/E_{Jet}"
60	nowak	1.2	if re.search( "dR", name ):
61			studied = "max #delta R (non iso elec, jet)"
62	nowak	1.1	elif re.search( "MJ", name ):
63			studied = "#delta R (muon, jet)"
64
65			return studied
66
67	nowak	1.2	def MakeDivision( h_1, h_2, option ):
68			"""
69			divides two histograms, computes errors
70			and returns the result
71			"""
72			#h_1.Sumw2()
73			#h_2.Sumw2()
74			hist = h_2.Clone()
75			hist.Divide( h_1, h_2, 1., 1., option )
76
77			return hist
78
79			def Fill( hist, option, weight ):
80			"""
81			h.SetBinContent sets the content right, but has only
82			one wheighted entry then. This gives problems in dividing the
83			histograms. This fill function takes care of that
84			"""
85			for i in range( 0, weight ):
86			hist.Fill( float( option ) )
87
88			return hist
89
90
91	nowak	1.1	if __name__ == "__main__":
92
93	nowak	1.2	import getopt
94
95			basepath = main(sys.argv[1:])
96
97
98	nowak	1.1	Def.SetGlobalStyles( "emr" )
99
100			saves = {}
101
102	nowak	1.2	config = readTheConfig( basepath )
103	nowak	1.1
104			for section in config.sections():
105
106	nowak	1.2	print "section: ", section
107
108	nowak	1.1	studied_cut = getStudiedCut( section )
109	nowak	1.2
110			nbin = 9
111			min = 0.05
112			max = 0.95
113			if( re.search( "EJ_dR", section ) ):
114			min = 0.15
115			max = 1.05
116			pass
117
118			h_mu_m = ROOT.TH1F( "h_mu_m", "h_mu_m",nbin,min,max )
119			h_mu_gen = ROOT.TH1F( "h_mu_gen", "h_mu_gen",nbin,min,max )
120			h_mu_nm = ROOT.TH1F( "h_mu_nm", "h_mu_nm",nbin,min,max )
121			h_mu_all = ROOT.TH1F( "h_mu_all", "h_mu_all",nbin,min,max )
122			h_el_m = ROOT.TH1F( "h_el_m", "h_el_m",nbin,min,max )
123			h_el_gen = ROOT.TH1F( "h_el_gen", "h_el_gen",nbin,min,max )
124			h_el_nm = ROOT.TH1F( "h_el_nm", "h_el_nm",nbin,min,max )
125			h_el_all = ROOT.TH1F( "h_el_all", "h_el_all",nbin,min,max )
126			h_ph_m = ROOT.TH1F( "h_ph_m", "h_ph_m",nbin,min,max )
127			h_ph_gen = ROOT.TH1F( "h_ph_gen", "h_ph_gen",nbin,min,max )
128			h_ph_nm = ROOT.TH1F( "h_ph_nm", "h_ph_nm",nbin,min,max )
129			h_ph_all = ROOT.TH1F( "h_ph_all", "h_ph_all",nbin,min,max )
130			h_je_m = ROOT.TH1F( "h_je_m", "h_je_m",nbin,min,max )
131			h_je_gen = ROOT.TH1F( "h_je_gen", "h_je_gen",nbin,min,max )
132			h_je_nm = ROOT.TH1F( "h_je_nm", "h_je_nm",nbin,min,max )
133			h_je_all = ROOT.TH1F( "h_je_all", "h_je_all",nbin,min,max )
134
135			h_mu_m.Sumw2()
136			h_mu_gen.Sumw2()
137			h_mu_nm.Sumw2()
138			h_mu_all.Sumw2()
139			h_el_m.Sumw2()
140			h_el_gen.Sumw2()
141			h_el_nm.Sumw2()
142			h_el_all.Sumw2()
143			h_ph_m.Sumw2()
144			h_ph_gen.Sumw2()
145			h_ph_nm.Sumw2()
146			h_ph_all.Sumw2()
147			h_je_m.Sumw2()
148			h_je_gen.Sumw2()
149			h_je_nm.Sumw2()
150			h_je_all.Sumw2()
151	nowak	1.1
152	nowak	1.2
153			#print type( h_mu_eff)
154	nowak	1.1
155			for option in config.options( section ):
156
157	nowak	1.2	print "option: ", option
158
159			if re.search( "afs", config.get( section, option ) ) or re.search( "current", config.get( section, option ) ):
160			fin = Def.ReadInRootFile( config.get( section, option ) + "/studyXC.root" )
161			else:
162			fin = Def.ReadInRootFile( "./../test/" + config.get( section, option ) + "/studyXC.root" )
163	nowak	1.1
164	nowak	1.2
165	nowak	1.1	Tree_ = fin.Get( 'StudyXC' )
166			mu_all = Def.defineArray( 'i', 1 )
167			mu_gen = Def.defineArray( 'i', 1 )
168			mu_m = Def.defineArray( 'i', 1 )
169			mu_nm = Def.defineArray( 'i', 1 )
170			el_all = Def.defineArray( 'i', 1 )
171			el_gen = Def.defineArray( 'i', 1 )
172			el_m = Def.defineArray( 'i', 1 )
173			el_nm = Def.defineArray( 'i', 1 )
174			ph_all = Def.defineArray( 'i', 1 )
175			ph_gen = Def.defineArray( 'i', 1 )
176			ph_m = Def.defineArray( 'i', 1 )
177			ph_nm = Def.defineArray( 'i', 1 )
178			je_all = Def.defineArray( 'i', 1 )
179			je_gen = Def.defineArray( 'i', 1 )
180			je_m = Def.defineArray( 'i', 1 )
181			je_nm = Def.defineArray( 'i', 1 )
182
183			Tree_.SetBranchAddress( "MuonsAll", mu_all )
184			Tree_.SetBranchAddress( "MuonsGen", mu_gen )
185			Tree_.SetBranchAddress( "MuonsNotMatched", mu_nm )
186			Tree_.SetBranchAddress( "MuonsMatched", mu_m )
187			Tree_.SetBranchAddress( "ElectronsAll", el_all )
188			Tree_.SetBranchAddress( "ElectronsGen", el_gen )
189			Tree_.SetBranchAddress( "ElectronsNotMatched", el_nm )
190			Tree_.SetBranchAddress( "ElectronsMatched", el_m )
191			Tree_.SetBranchAddress( "PhotonsAll", ph_all )
192			Tree_.SetBranchAddress( "PhotonsGen", ph_gen )
193			Tree_.SetBranchAddress( "PhotonsNotMatched", ph_nm )
194			Tree_.SetBranchAddress( "PhotonsMatched", ph_m )
195			Tree_.SetBranchAddress( "JetsAll", je_all )
196			Tree_.SetBranchAddress( "JetsGen", je_gen )
197			Tree_.SetBranchAddress( "JetsNotMatched", je_nm )
198			Tree_.SetBranchAddress( "JetsMatched", je_m )
199
200			n_mu_all = 0
201			n_mu_gen = 0
202			n_mu_m = 0
203			n_mu_nm = 0
204			n_el_all = 0
205			n_el_gen = 0
206			n_el_m = 0
207			n_el_nm = 0
208			n_ph_all = 0
209			n_ph_gen = 0
210			n_ph_m = 0
211			n_ph_nm = 0
212			n_je_all = 0
213			n_je_gen = 0
214			n_je_m = 0
215			n_je_nm = 0
216
217			n_bytes = 0
218			for ent in range( 0, Tree_.GetEntries() ):
219
220			n_bytes += Tree_.GetEntry( ent )
221			n_mu_all += mu_all[0]
222			n_mu_gen += mu_gen[0]
223			n_mu_m += mu_m[0]
224			n_mu_nm += mu_nm[0]
225			n_el_all += el_all[0]
226			n_el_gen += el_gen[0]
227			n_el_m += el_m[0]
228			n_el_nm += el_nm[0]
229			n_ph_all += ph_all[0]
230			n_ph_gen += ph_gen[0]
231			n_ph_m += ph_m[0]
232			n_ph_nm += ph_nm[0]
233			n_je_all += je_all[0]
234			n_je_gen += je_gen[0]
235			n_je_m += je_m[0]
236			n_je_nm += je_nm[0]
237			pass
238
239
240	nowak	1.2	Fill( h_mu_m, option, n_mu_m )
241			Fill( h_mu_gen, option, n_mu_gen )
242			Fill( h_mu_nm, option, n_mu_nm )
243			Fill( h_mu_all, option, n_mu_all )
244			Fill( h_el_m, option, n_el_m )
245			Fill( h_el_gen, option, n_el_gen )
246			Fill( h_el_nm, option, n_el_nm )
247			Fill( h_el_all, option, n_el_all )
248			Fill( h_ph_m, option, n_ph_m )
249			Fill( h_ph_gen, option, n_ph_gen )
250			Fill( h_ph_nm, option, n_ph_nm )
251			Fill( h_ph_all, option, n_ph_all )
252			Fill( h_je_m, option, n_je_m )
253			Fill( h_je_gen, option, n_je_gen )
254			Fill( h_je_nm, option, n_je_nm )
255			Fill( h_je_all, option, n_je_all )
256
257			print "n_mu_all: ", n_mu_all
258			print "n_mu_m : ", n_mu_m
259			print "n_mu_nm : ", n_mu_nm
260			print "n_mu_gen: ", n_mu_gen
261
262	nowak	1.1	pass
263
264	nowak	1.2	h_mu_eff = MakeDivision( h_mu_m, h_mu_gen, "B" )
265			h_mu_cont = MakeDivision( h_mu_nm, h_mu_all, "B" )
266			h_el_eff = MakeDivision( h_el_m, h_el_gen, "B" )
267			h_el_cont = MakeDivision( h_el_nm, h_el_all, "B" )
268			h_ph_eff = MakeDivision( h_ph_m, h_ph_gen, "B" )
269			h_ph_cont = MakeDivision( h_ph_nm, h_ph_all, "B" )
270			h_je_eff = MakeDivision( h_je_m, h_je_gen, "B" )
271			h_je_cont = MakeDivision( h_je_nm, h_je_all, "B" )
272
273			#for i in range( 0, 9 ):
274			# print "error in bin ", i, ": ", h_mu_eff.GetBinError( i )
275
276
277			### wrong errors here
278			h_mu_ratio = MakeDivision( h_mu_eff, h_mu_cont, "" )
279			h_el_ratio = MakeDivision( h_el_eff, h_el_cont, "" )
280			h_ph_ratio = MakeDivision( h_ph_eff, h_ph_cont, "" )
281			h_je_ratio = MakeDivision( h_je_eff, h_je_cont, "" )
282
283
284	nowak	1.1	## draw the hists for each section
285	nowak	1.2	c_mu = ROOT.TCanvas( "c_mu", section + "_Muon" ,1,71,300,600 )
286			c_mu.Divide(1,2)
287			c_mu.cd(1)
288	nowak	1.1	Def.DrawSmall( c_mu, \
289			[h_mu_eff, h_mu_cont], [8,50], ["eff","cont"], \
290			studied_cut, "", "", \
291	nowak	1.2	False, 0, False, "ur", section + "muons", saves)
292
293			c_mu.cd(2)
294			Def.DrawSmall( c_mu, \
295			[h_mu_ratio], [4], ["ratio"], \
296			studied_cut, "ratio eff / sqrt(cont)", "", \
297			False, 0, False, "ur", section + "muons_ratio", saves)
298	nowak	1.1
299	nowak	1.2	c_el = ROOT.TCanvas( "c_el", section + "_Elecs" ,1,71,300,600 )
300			c_el.Divide(1,2)
301			c_el.cd(1)
302	nowak	1.1	Def.DrawSmall( c_el, \
303			[h_el_eff, h_el_cont], [8,50], ["eff","cont"], \
304			studied_cut, "", "", \
305	nowak	1.2	False, 0, False, "ur", section + "elecs", saves)
306
307			c_el.cd(2)
308			Def.DrawSmall( c_el, \
309			[h_el_ratio], [4], ["ratio"], \
310			studied_cut, "ratio eff / sqrt(cont)", "", \
311			False, 0, False, "ur", section + "elecs_ratio", saves)
312	nowak	1.1
313	nowak	1.2	c_ph = ROOT.TCanvas( "c_ph", section + "_Photons" ,1,71,300,600 )
314			c_ph.Divide(1,2)
315			c_ph.cd(1)
316	nowak	1.1	Def.DrawSmall( c_ph, \
317			[h_ph_eff, h_ph_cont], [8,50], ["eff","cont"], \
318			studied_cut, "", "", \
319	nowak	1.2	False, 0, False, "ur", section + "phots", saves)
320
321			c_ph.cd(2)
322			Def.DrawSmall( c_ph, \
323			[h_ph_ratio], [4], ["ratio"], \
324			studied_cut, "ratio eff / sqrt(cont)", "", \
325			False, 0, False, "ur", section + "phots_ratio", saves)
326	nowak	1.1
327	nowak	1.2	c_je = ROOT.TCanvas( "c_je", section + "_Jets" ,1,71,300,600 )
328			c_je.Divide(1,2)
329			c_je.cd(1)
330	nowak	1.1	Def.DrawSmall( c_je, \
331			[h_je_eff, h_je_cont], [8,50], ["eff","cont"], \
332			studied_cut, "", "", \
333	nowak	1.2	False, 0, False, "ur", section + "jets", saves)
334
335			c_je.cd(2)
336			Def.DrawSmall( c_je, \
337			[h_je_ratio], [4], ["ratio"], \
338			studied_cut, "ratio eff / sqrt(cont)", "", \
339			False, 0, False, "ur", section + "jets_ratio", saves)
340
341
342	nowak	1.1	pass
343
344
345
346
347
348			Def.DontQuit()