add functions to collect particle yield in particleReader.py