# set up some useful paths of locations on the MicroDST flavTagPath = locationRoot + "/" + mainLocation + "/FlavourTags" lhcbApp = LHCbApp() lhcbApp.DDDBtag = 'default' lhcbApp.CondDBtag = 'default' appMgr = AppMgr(outputlevel=4) appMgr.config(files=['$GAUDIPOOLDBROOT/options/GaudiPoolDbRoot.opts']) appMgr.initialize() appMgr.ExtSvc += ['LHCb::ParticlePropertySvc'] appMgr.HistogramPersistency = "ROOT" evtSvc = appMgr.evtSvc() toolSvc = appMgr.toolsvc() evtSel = appMgr.evtSel() nextEvent = Functors.NextEvent(appMgr) pp = Functors.PartPropSvc(appMgr) ppSvc = pp particleNameFunc = Functors.ParticleName(ppSvc) # open a DST or MicroDST evtSel.open(microDSTFile) histoPath = "MicroDST/Histos/" omegaPlots = {} # flavour tagging plot: plot mistag for each category. for cat in range(0, 6): omegaPlots[cat] = book(histoPath + "FlavourTags/MisTagCat" + str(cat), "Mis-tag fraction category " + str(cat), 100, 0.,
#pattern = '/lhcb/LHCb/Collision12/FMDST/00021391/0000/00021391_00000%03d_1.fmdst' from PRConfig import TestFileDB from GaudiConf.IOHelper import IOHelper IOHelper().inputFiles( TestFileDB.test_file_db["2012_CaloFemtoDST"].filenames) #evtSel.open ( 'KaliPi0_2k+11.fmdst_test') from GaudiPython.Bindings import AppMgr gaudi = AppMgr() ## === OR with test fmDSTs evtSel = gaudi.evtSel() ## run over events gaudi.run(5000) from KaliCalo.FitUtils import fitPi0, getPi0Params, s2b import GaudiPython.GaudiAlgs pi0 = gaudi.algorithm('KaliPi0') histos = pi0.Histos() if histos: keys = histos.keys() keys.sort() for h in keys: histo = histos[h] if hasattr(histo, 'dump'): print histo.dump(60, 25, True)