def histo(variable, hname): h = None mtBinEdges = fakeu.mtBinEdges() ptBinEdges = fakeu.ptBinEdges() etaBinEdges = fakeu.etaBinEdges() if variable == 'mt0': h = r.TH1F(hname, ';m_{T}(tag,MET) [GeV]; entries/bin', len(mtBinEdges) - 1, mtBinEdges) elif variable == 'mt1': h = r.TH1F(hname, ';m_{T}(probe,MET) [GeV]; entries/bin', len(mtBinEdges) - 1, mtBinEdges) elif variable == 'pt0': h = r.TH1F(hname, ';p_{T,l0} [GeV]; entries/bin', len(ptBinEdges) - 1, ptBinEdges) elif variable == 'pt1': h = r.TH1F(hname, ';p_{T,l1} [GeV]; entries/bin', len(ptBinEdges) - 1, ptBinEdges) elif variable == 'eta1': h = r.TH1F(hname, ';#eta_{l1}; entries/bin', len(etaBinEdges) - 1, etaBinEdges) elif v == 'pt1_eta1': h = r.TH2F(hname, ';p_{T} [GeV]; #eta', len(ptBinEdges) - 1, ptBinEdges, len(etaBinEdges) - 1, etaBinEdges) else: print "unknown variable %s" % v h.SetDirectory(0) h.Sumw2() return h
def shiftWithinRange(pt, eta, mt, epsilon=1.0e-3): ptBins, etaBins, mtBins = fakeu.ptBinEdges(), fakeu.etaBinEdges(), fakeu.mtBinEdges() minPt, maxPt = min(ptBins), max(ptBins) minEta, maxEta = min(etaBins), max(etaBins) minMt, maxMt = min(mtBins), max(mtBins) pt = minPt*(1.0+epsilon) if pt<minPt else maxPt*(1.0-epsilon) if pt > maxPt else pt eta = minEta*(1.0+epsilon) if eta<minEta else maxEta*(1.0-epsilon) if eta > maxEta else eta mt = minMt*(1.0+epsilon) if mt<minMt else maxMt*(1.0-epsilon) if mt > maxMt else mt return pt, eta, mt
def histo(variable, hname): h = None mtBinEdges = fakeu.mtBinEdges() ptBinEdges = fakeu.ptBinEdges() etaBinEdges = fakeu.etaBinEdges() if variable=='mt0' : h = r.TH1F(hname, ';m_{T}(tag,MET) [GeV]; entries/bin', len(mtBinEdges)-1, mtBinEdges) elif variable=='mt1' : h = r.TH1F(hname, ';m_{T}(probe,MET) [GeV]; entries/bin', len(mtBinEdges)-1, mtBinEdges) elif variable=='pt0' : h = r.TH1F(hname, ';p_{T,l0} [GeV]; entries/bin', len(ptBinEdges)-1, ptBinEdges) elif variable=='pt1' : h = r.TH1F(hname, ';p_{T,l1} [GeV]; entries/bin', len(ptBinEdges)-1, ptBinEdges) elif variable=='eta1' : h = r.TH1F(hname, ';#eta_{l1}; entries/bin', len(etaBinEdges)-1, etaBinEdges) else : print "unknown variable %s"%v h.SetDirectory(0) h.Sumw2() return h
def histo(variable, hname): h = None mtBinEdges = fakeu.mtBinEdges() ptBinEdges = fakeu.ptBinEdges() etaBinEdges = fakeu.etaBinEdges() mdrBinEdges = fakeu.mdeltarBinEdges() if variable=='pt' : h = r.TH1F(hname, ';p_{T,l} [GeV]; entries/bin', len(ptBinEdges)-1, ptBinEdges) elif variable=='eta' : h = r.TH1F(hname, ';|#eta_{l}|; entries/bin', len(etaBinEdges)-1, etaBinEdges) elif variable=='mt' : h = r.TH1F(hname, ';m_{T}(l,MET) [GeV]; entries/bin', len(mtBinEdges)-1, mtBinEdges) elif variable=='mdeltar': h = r.TH1F(hname, ';m_{T}(l,MET) [GeV]; entries/bin', len(mdrBinEdges)-1, mdrBinEdges) elif variable=='pt_eta' : h = r.TH2F(hname, ';p_{T,l} [GeV]; #eta_{l};', len(ptBinEdges)-1, ptBinEdges, len(etaBinEdges)-1, etaBinEdges) else : print "unknown variable %s"%v h.SetDirectory(0) h.Sumw2() return h