コード例 #1
0
 def endJob(self):
     #Normalise histograms at the end to xsec/sumw
     if self.doNorm:
         for obj in self.objs:
             if isinstance(obj, ROOT.TH1) and obj!=self.h_count:
                 scale = 1./self.genEventSumw
                 if self.xsec > 0:
                     scale *= self.xsec
                 obj.Scale(scale)
         #store x-sec and sum-of-gen-weights in count histogram
         self.h_count.SetBinContent(0,self.genEventSumw) #underflow
         self.h_count.SetBinContent(self.h_count.GetNbinsX()+1,self.xsec) #overflow
     # event count -> efficiency
     nbins = self.h_count.GetNbinsX()
     n_init = self.h_count.GetBinContent(1)
     e2_init = pow(self.h_count.GetBinError(1),2)
     for ibin in range(1,nbins+1):
         n_ibin = self.h_count.GetBinContent(ibin)
         self.h_eff.SetBinContent(ibin,n_ibin/n_init)
         #binomial error
         err = 0
         if ibin!=1 and n_ibin!=0:
             e2_ibin = pow(self.h_count.GetBinError(ibin),2)
             err = math.sqrt(abs(((1.-2.*n_init/n_ibin)*e2_init+pow(n_init,2)*e2_ibin/pow(n_ibin,2))/pow(n_ibin,2)))
         self.h_eff.SetBinError(ibin,err)
     Module.endJob(self)
コード例 #2
0
 def endJob(self):
     print "\nendJob()"
     if self.h_Stat1:
         print "nendJob():: printCutFlowTable::"
         printCutFlowTable(self.h_Stat1)
         
     Module.endJob(self)
コード例 #3
0
    def endJob(self):
        Module.endJob(self)

        print 'Events passing:', self.count

        if nSelect == 2.5 or nSelect == 5:
            self.diphoton_file.Write()
            self.diphoton_file.Close()
コード例 #4
0
 def endJob(self):
     """Calculate the correlated and uncorrelated errors.
     """
     if not self.writePredDist:
         for cat in self.anacatvals : 
             for hist in [self.predJetP[cat], self.predJetMTTBAR[cat],self.predJetMTTBARMod[cat] ] : 
                 hist.SetCalculatedErrors()
     Module.endJob(self)
     pass
コード例 #5
0
 def endJob(self):
     Module.endJob(self)
     pass
コード例 #6
0
ファイル: signalStudy.py プロジェクト: jwill24/PPSAnalyzer
    def endJob(self):
        Module.endJob(self)

        print 'Signal Efficiency:', str(
            float(100) * float(self.passing_wp90) / float(self.total)) + r'%'
        '''
コード例 #7
0
 def endJob(self):
     """end Job"""
     Module.endJob(self)
コード例 #8
0
 def endJob(self):
     """end Job"""
     #self.evelurun.close()
     Module.endJob(self)
コード例 #9
0
 def endJob(self):
     Module.endJob(self)
     print "+ Module ended successfully,"  #, self.h_events.GetEntries(), "events analyzed"
     pass
コード例 #10
0
 def endJob(self):
     # print("nEvents: {0}\t nPositivEvents: {1}\t nNegativeEvents: {2}\t nZeroEvents: {3}".format(
     #     self.nEvents, self.nPositiveEvents, self.nNegativeEvents, self.nZeroEvents)
     # )
     Module.endJob(self)
コード例 #11
0
 def endJob(self):
     Module.endJob(self)
     print "Module ended successfully,", self.nEvent, "events analyzed"
     pass
コード例 #12
0
    def endJob(self):
        self.MakeRatio("gen_ele_efficiency", self.h_gen_ele_pt_match,
                       self.h_gen_ele_pt)

        Module.endJob(self)
コード例 #13
0
ファイル: SkimmerDiMuon.py プロジェクト: nurfikri89/PUjetID
 def endJob(self):
   Module.endJob(self)
   print "SkimmerDiMuon module ended successfully"
   pass