# GaudiPython configuration. gaudi = GaudiPython.AppMgr() tes = gaudi.evtsvc() # Run. import sys, ROOT from math import floor evtmax = -1 #1000#-1 #try: evtmax = int(sys.argv[1]) #except: evtmax = float('inf') evtnum = 0 ntuple = Ntuple('output.root', tes, gaudi.toolsvc(), gaudi.detSvc(), recJB.Output, recSVs.outputLocation(), recMus.algorithm().Output, addBackTags=True) while evtmax < 0 or evtnum < evtmax: gaudi.run(1) if not bool(tes['/Event']): break evtnum += 1 ntuple.clear() # Fill event info. try: ntuple.ntuple['evt_pvr_n'][0] = len(tes['Rec/Vertex/Primary']) except: continue try:
from Ntuple import Ntuple # GaudiPython configuration. gaudi = GaudiPython.AppMgr() tes = gaudi.evtsvc() # Run. import sys, ROOT from math import floor evtmax = -1 #' #try: evtmax = int(sys.argv[1]) #except: evtmax = float('inf') evtnum = 0 ntuple = Ntuple('output.root', tes, gaudi.toolsvc(), gaudi.detSvc(), recJB.Output, recSVs.outputLocation(), recMus.algorithm().Output) #gaudi.run(39000)#TODO while evtmax < 0 or evtnum < evtmax: gaudi.run(1) if not bool(tes['/Event']): break evtnum += 1 ntuple.clear() # Fill event info. try: ntuple.ntuple['evt_pvr_n'][0] = len(tes['Rec/Vertex/Primary']) except: continue try: ntuple.addEventInfo(); except: continue fill = False;
# GaudiPython configuration. gaudi = GaudiPython.AppMgr() tes = gaudi.evtsvc() # Run. import sys, ROOT from math import floor #evtmax = 10000 #TODO evtmax = -1 #TODO #try: evtmax = int(sys.argv[1]) #except: evtmax = float('inf') evtnum = 0 ntuple = Ntuple('output.root', tes, gaudi.toolsvc(), gaudi.detSvc(), primaries='/Event/' + stream + '/Rec/Vertex/Primary') while evtmax < 0 or evtnum < evtmax: gaudi.run(1) if not bool(tes['/Event']): break evtnum += 1 ntuple.clear() # Fill event info. try: ntuple.ntuple['evt_pvr_n'][0] = len(tes['/Event/' + stream + '/Rec/Vertex/Primary']) except: continue
from Ntuple import Ntuple # GaudiPython configuration. gaudi = GaudiPython.AppMgr() tes = gaudi.evtsvc() # Run. import sys, ROOT from math import floor evtmax = -1 #try: evtmax = int(sys.argv[1]) #except: evtmax = float('inf') evtnum = 0 ntuple = Ntuple('output.root', tes, gaudi.toolsvc(), gaudi.detSvc(), recJB.Output, recSVs.outputLocation(), recMus.algorithm().Output) while evtmax < 0 or evtnum < evtmax: gaudi.run(1) if not bool(tes['/Event']): break evtnum += 1 ntuple.clear() # Fill event info. try: ntuple.ntuple['evt_pvr_n'][0] = len(tes['Rec/Vertex/Primary']) except: continue try: ntuple.ntuple['evt_trk_n'][0] = len( tes['Phys/StdAllNoPIDsPions/Particles'])
LHCB = gbl.LHCb from Ntuple import Ntuple # GaudiPython configuration. gaudi = GaudiPython.AppMgr() tes = gaudi.evtsvc() # Run. import sys, ROOT from math import floor evtmax = -1 #try: evtmax = int(sys.argv[1]) #except: evtmax = float('inf') evtnum = 0 ntuple = Ntuple('output.root', tes, gaudi.toolsvc(), gaudi.detSvc(), recJB.Output, recSVs.outputLocation(), recMus.algorithm().Output) while evtmax < 0 or evtnum < evtmax: gaudi.run(1) if not bool(tes['/Event']): break evtnum += 1 ntuple.clear() # Fill event info. try: ntuple.ntuple['evt_pvr_n'][0] = len(tes['Rec/Vertex/Primary']) except: continue try: ntuple.addEventInfo(); except: continue # Fill generator level info. fill = False; gens = tes['MC/Particles']
#!/usr/bin/env python import ROOT from Ntuple import Ntuple ROOT.gROOT.SetBatch(True) n = Ntuple("DataTree") n.add("data.root") entries = n.GetEntries() print "Entries: ", entries gauss_h = ROOT.TH1F("gauss", "Gaussian distribution", 100, -5., 5.) uniform_h = ROOT.TH1I("uniform", "Uniform distribution", 101, 0, 101) for entry in range(entries): n.GetEntry(entry) gauss_h.Fill(n.val("gauss_0")) uniform_h.Fill(n.val("uniform_0")) canvas = ROOT.TCanvas("c", "c", 400, 400) gauss_h.Draw() canvas.SaveAs("gauss.ntuple_py.pdf") canvas.Clear() uniform_h.Draw() canvas.SaveAs("uniform.ntuple_py.pdf")
LHCB = gbl.LHCb from Ntuple import Ntuple # GaudiPython configuration. gaudi = GaudiPython.AppMgr() tes = gaudi.evtsvc() # Run. import sys, ROOT from math import floor evtmax = -1 #TODO #try: evtmax = int(sys.argv[1]) #except: evtmax = float('inf') evtnum = 0 ntuple = Ntuple('output.root', tes, gaudi.toolsvc(), gaudi.detSvc()) while evtmax < 0 or evtnum < evtmax: gaudi.run(1) if not bool(tes['/Event']): break evtnum += 1 ntuple.clear() # Fill event info. try: ntuple.ntuple['evt_pvr_n'][0] = len(tes['Rec/Vertex/Primary']) except: continue # Fill generator level info. fill = False