Pythonic access to high energy particle data tables and PDG ID codes.
Pythonic access to high energy particle data tables and ID codes.
This is a rethink in Python of the commonly used HepPDT/HepPID packages for accessing particle data tables and PDG scheme particle ID codes.
Particle data table functions (cf. HepPDT)
HepPDT has no native Python interface, and it is not easy to naturally wrap its C++ interface for Python, so PyPDT is written from scratch in Python, with a lighter class count than the original: it’s just a convenience interface to a small ASCII data file, so who needs lots of classes?
>>> tbl = pypdt.ParticleDataTable('particle.tbl') >>> # or tbl = pypdt.ParticleDataTable() # to read from AFS default db >>> # or tbl = pypdt.PDT() # if you prefer minimal keystrokes to clarity :) >>> print tbl # print a summary of muon information >>> tbl.mass # access the top mass in GeV >>> tbl.width # access the Z width in GeV >>> tbl.lifetime # access the K0S lifetime in ps >>> tbl.ctau # access the K0S lifetime in mm >>> tbl.mean_disp(10.) # get the mean flight distance of a 10 GeV K0S in mm >>> for p in tbl: print p # print summaries for all known particles >>> # print info for all particles in asc lifetime order: >>> for t, p in sorted((p.ctau, p) for p in tbl): print p
The world is your bivalve of choice.
For maximum convenience, you can even use PyPDT from the command line:
$ python -m pypdt 13 310 mu^-: ID=13, m=1.06e-01 GeV, 3*q=-3, width=9.99e-19 GeV, tau=6.59e+05 ps, ctau=1.97e+05 mm K_S^0: ID=310, m=4.98e-01 GeV, 3*q=0, width=2.46e-14 GeV, tau=2.68e+01 ps, ctau=8.02e+00 mm
(If /afs/cern.ch is not present, add -f path/to/my/particle.tbl to the above.)
Particle ID functions (cf. HepPID)
PyPDT also contains Python versions of the HepPID library of functions for interpreting particle ID codes in the PDG scheme. The code of these functions has been directly translated to Python from HepPID’s ParticleIdMethods – please notify the author of this package in case of inaccuracies or evolution of the PID scheme.
>>> pypdt.isSUSY(2010011) True
- Add PID functions to the ParticleData class
- Use doctest