evtpdl.py

#!/usr/bin/env python3
 

 
# this is a test executable, not a module so we don't need doxygen warnings
# @cond SUPPRESS_DOXYGEN
 
"""
Check reading and writing of evt.pdl files as well as errors when creating new particles.
 
This script tests that
 - the ROOT.TDatabasePDG instance is in fact an EvtGenDatabasePDG
 - reading a pdl file raises a warning if it's not the original one
 - creating a particle which already exists is forbidden
 - creating particles with whitespace in the name is forbidden
 - writing an evt.pdl results in the same output as the original pdl, except for
   whitespace and additionally defined particles
"""
 
 
from tempfile import NamedTemporaryFile
import basf2
import ROOT
import os
import sys
import subprocess
import pdg
 
db = ROOT.TDatabasePDG.Instance()
 
retcode = 0
 
# Test particle data.
if not pdg._get_instance().testParticleData():
    retcode = 1
 
# create particle which already exists, should print error
db.AddParticle("e-", "electron", 0, False, 0, 0, "duplicate", 11, 0, 0, 0, 0, 0, 0)
# create particle with whitespace in the name, should print error
db.AddParticle("invalid name", "invalid name", 0, False, 0, 0, "custom", 2345, 0, 0, 0, 0, 0, 0)
# and also create a particle with the python function: should give error due to
# duplicate pdg code
pdg.add_particle("test2", 1, 5.28, 0.1, 0, 0.5)
pdg.add_particle("foo\tbar", 10001, 0, 0, 0, 0)
 
# default evt.pdl filename
default_evtpdl = basf2.find_file(os.path.join("data", "framework", "particledb", "evt.pdl"))
# create a temporary one and compare
with NamedTemporaryFile() as tempfile:
    # write a evt.pdl from the EventGenDatabasePDG
    db.WriteEvtGenTable(tempfile.name)
    # and make sure that it's the same as the original one except for
    # whitespace, empty lines and lines starting with a "*"
    with open(default_evtpdl) as default:
        old_lines = [x.split() for x in default if x[0] != '*' and not x.startswith("end")]
    with open(tempfile.name) as changed:
        new_lines = [x.split() for x in changed if x[0] != '*' and not x.startswith("end")]
 
    if len(old_lines) != len(new_lines):
        basf2.B2ERROR("Different number of lines when writing pdl file")
        retcode = 1
    else:
        # field types for comparison
        types = [str] * 4 + [int, float, float, float, int, int, float, int]
        # check all lines for equality but don't care about formatting
        for old, new in zip(old_lines, new_lines):
            # make sure the lines have the same number of fields
            if len(old) != len(new):
                basf2.B2ERROR("Number of field different between old and new:",
                              old=old, new=new)
                retcode = 1
            # and that they have enough fields for comparison
            elif len(new) < len(types):
                basf2.B2ERROR("Not enough fields", old=old, new=new)
                retcode = 1
            # and then make sure they're the same compariny the correct types
            else:
                if not all(t(o) == t(n) for t, o, n in zip(types, old, new)):
                    basf2.B2ERROR("Particles differ", old=old, new=new)
                    retcode = 1
 
    with NamedTemporaryFile() as tempfile2:
        # create a new particle
        db.AddParticle("testparticle", "testparticle", 1, False, 2, 3, "custom", 10001, 4, 5, 6, 7, 8, 9)
 
        # also test creating particle with the python functions
        pdg.add_particle("test2", 10002, 5.28, 0.1, -1, 0.5)
        # write the table again
        db.WriteEvtGenTable(tempfile2.name)
        # and make sure that it has the additional particle
        subprocess.call(["diff", "--old-line-format=- %L", "--new-line-format=+ %L",
                         "--unchanged-line-format=", tempfile.name, tempfile2.name])
 
db.ReadEvtGenTable("/no/such/file/(I hope)")
# and also test the python function
pdg.load("/no/such/file/(again)")
 
print("now we should have zero entries: ", db.ParticleList().GetEntries())
print("so reread default file...")
db.ReadEvtGenTable()
print("number of entries is now", db.ParticleList().GetEntries() and "positive" or "zero")
sys.exit(retcode)
 
# @endcond