ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/VHbb/python/HistoMaker.py
(Generate patch)

Comparing UserCode/VHbb/python/HistoMaker.py (file contents):
Revision 1.7 by peller, Thu Oct 11 13:40:18 2012 UTC vs.
Revision 1.9 by nmohr, Thu Oct 18 13:30:32 2012 UTC

# Line 18 | Line 18 | class HistoMaker:
18          self.region = region
19          self.lumi=0.
20  
21 <    def getScale(self,job,subsample=-1):
21 >    def getScale(self,job,subsample=-1,MC_rescale_factor=1):
22          anaTag=self.config.get('Analysis','tag')
23          input = TFile.Open(self.path+'/'+job.getpath())
24          CountWithPU = input.Get("CountWithPU")
# Line 32 | Line 32 | class HistoMaker:
32              sf=float(job.sf)
33          theScale = 1.
34          if anaTag == '7TeV':
35 <            theScale = float(self.lumi)*xsec*sf/(0.46502*CountWithPU.GetBinContent(1)+0.53498*CountWithPU2011B.GetBinContent(1))*self.rescale/float(job.split)
35 >            theScale = float(self.lumi)*xsec*sf/(0.46502*CountWithPU.GetBinContent(1)+0.53498*CountWithPU2011B.GetBinContent(1))*MC_rescale_factor/float(job.split)
36          elif anaTag == '8TeV':
37 <            theScale = float(self.lumi)*xsec*sf/(CountWithPU.GetBinContent(1))*self.rescale/float(job.split)
37 >            theScale = float(self.lumi)*xsec*sf/(CountWithPU.GetBinContent(1))*MC_rescale_factor/float(job.split)
38          return theScale
39  
40  
# Line 46 | Line 46 | class HistoMaker:
46  
47          #get the conversion rate in case of BDT plots
48          TrainFlag = eval(self.config.get('Analysis','TrainFlag'))
49 –        if TrainFlag:
50 –            MC_rescale_factor=2.
51 –            print 'I RESCALE BY 2.0'
52 –        else: MC_rescale_factor = 1.
53 –
49          BDT_add_cut='EventForTraining == 0'
50  
51  
# Line 103 | Line 98 | class HistoMaker:
98              if job.type != 'DATA':
99                  if CuttedTree.GetEntries():
100                      
101 <                    if 'BDT' in treeVar: drawoption = '(%s)*(%s)'%(weightF,BDT_add_cut)
101 >                    if 'RTight' in treeVar or 'RMed' in treeVar: drawoption = '(%s)*(%s)'%(weightF,BDT_add_cut)
102                      else: drawoption = '%s'%(weightF)
103                      output.cd()
104                      CuttedTree.Draw('%s>>%s(%s,%s,%s)' %(treeVar,name,nBins,xMin,xMax), drawoption, "goff,e")
# Line 114 | Line 109 | class HistoMaker:
109                  if options[11] == 'blind':
110                      output.cd()
111                      if treeVar == 'H.mass':
112 <                        CuttedTree.Draw('%s>>%s(%s,%s,%s)' %(treeVar,name,nBins,xMin,xMax),treeVar+'<80. || '+treeVar + '>150.' , "goff,e")
112 >                        CuttedTree.Draw('%s>>%s(%s,%s,%s)' %(treeVar,name,nBins,xMin,xMax),treeVar+'<90. || '+treeVar + '>150.' , "goff,e")
113                      else:
114                          CuttedTree.Draw('%s>>%s(%s,%s,%s)' %(treeVar,name,nBins,xMin,xMax),treeVar+'<0', "goff,e")
115  
# Line 129 | Line 124 | class HistoMaker:
124                  hTree = ROOT.TH1F('%s'%name,'%s'%name,nBins,xMin,xMax)
125                  hTree.Sumw2()
126              if job.type != 'DATA':
127 <                if 'BDT' in treeVar: ScaleFactor = self.getScale(job,subsample,MC_rescale_factor)
127 >                if 'RTight' in treeVar or 'RMed' in treeVar:
128 >                    if TrainFlag:
129 >                        MC_rescale_factor=2.
130 >                        print 'I RESCALE BY 2.0'
131 >                    else: MC_rescale_factor = 1.
132 >                    ScaleFactor = self.getScale(job,subsample,MC_rescale_factor)
133                  else: ScaleFactor = self.getScale(job,subsample)
134                  if ScaleFactor != 0:
135                      hTree.Scale(ScaleFactor)

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines