notebook.community

Edit and run



In [1]:

    
import hepmcanalysis
import ROOT
import rootnotes



In [2]:

    
h = ROOT.TH1F('finalstatept','finalstatept',50,0,500)



In [3]:

    
for e in hepmcanalysis.readhepmc('input.hepmc'):
    for p in e.fsParticles():
        h.Fill(p.momentum().perp())



In [4]:

    
c = rootnotes.default_canvas()



In [5]:

    
h.SetFillColor(38)
h.Draw()



In [6]:

    
c









    Out[6]:



In [6]: