PDGHelpers Namespace Reference

Classes
class	ExtraParticle
class	PDGParser

Functions
	getPDGTABLE (table)
	getExtraParticleAcceptList (acceptlist)
	updateExtraParticleAcceptList (listName='G4particle_acceptlist_ExtraParticles.txt', pdgcodes=[])

Detailed Description

Parser for the PDGTABLE.MeV file

    Creates a compact object that holds configuration
    for the ExtraParticlesPhysicsTool.

    Gaudi cannot parse nested std::map objects, so the extra
    particles are converted into a

      std::map< std::string, std::vector< double > >

    object and ordering of the values is therefore important.

    Rows are:
        name
        mass [MeV]
        width [MeV]
        pdg
        lifetime [ns]

    Charge, spin, parity, and isospin3 are set to zero.

    AUTHOR: miha muskinja
    DATE: August 2019

Function Documentation

getExtraParticleAcceptList()

PDGHelpers.getExtraParticleAcceptList

(

acceptlist

)

Definition at line 47 of file PDGHelpers.py.

def getExtraParticleAcceptList(acceptlist):
    # Delete a local file if present
    if os.path.isfile(acceptlist):
        os.remove(acceptlist)
    #create blank file
    blank = open('G4particle_acceptlist_ExtraParticles.txt', 'x')
    blank.close()
    return True
 
 

getPDGTABLE()

PDGHelpers.getPDGTABLE

(

table

)

Definition at line 37 of file PDGHelpers.py.

def getPDGTABLE(table):
    # Delete a local file if present
    if os.path.isfile(table):
        os.remove(table)
    # Grab the file
    os.system('get_files -data %s' % table)
    return True
 
 
@lru_cache

updateExtraParticleAcceptList()

PDGHelpers.updateExtraParticleAcceptList	(		listName = 'G4particle_acceptlist_ExtraParticles.txt',
			pdgcodes = [] )

Definition at line 57 of file PDGHelpers.py.

def updateExtraParticleAcceptList(listName='G4particle_acceptlist_ExtraParticles.txt', pdgcodes=[]):
    if getExtraParticleAcceptList(listName):
        import shutil
        shutil.copy(listName, listName+'.org')
    existingpdgcodes = [int(x) for x in open(listName).readlines()]
    newpdgcodes = list(set(pdgcodes).difference(existingpdgcodes))
    # update the acceptlist for GenParticleSimAcceptList
    with open(listName, 'a') as writer:
        for pdg in newpdgcodes:
            writer.write('%s\n' % pdg)