Check.py

Go to the documentation of this file.

#!/usr/bin/env python3

"""
The trivial example to inspect the particle data table from the data base
"""
 
# =============================================================================
__author__ = "Vanya BELYAEV vanya@nikhef.nl"
__verison__ = " "
# =============================================================================
 
from Configurables import ApplicationMgr, Gaudi__ParticlePropertySvc
 

ApplicationMgr().ExtSvc += [Gaudi__ParticlePropertySvc()]
 
# The following is misunderstood by flake8 - the import is needed as it
# has necessary side effects
import GaudiPartProp.Service  # noqa: F401
from GaudiKernel.SystemOfUnits import micrometer
 

from GaudiPython.Bindings import AppMgr
 
# instantiate the actual (C++) application manager
gaudi = AppMgr()
 
# switch off all algorithms
gaudi.TopAlg = []
 
# initialize everything properly
gaudi.run(1)
 
# get the actual (C++) particle properties service
ppsvc = gaudi.gaudiPartProp()
 
# finally: play with the service
 
# get all a self-conjugated particles
print(" Self-charge conjugated particles ")
selfcc = ppsvc.get(lambda x: x.selfcc())
print(selfcc)
 
# set self-charge conjugated conjugated mesons:
print(" CC-neutral mesons")
mesonscc = ppsvc.get(lambda x: (x.selfcc() and x.isMeson()))
print(mesonscc)
 
# set self-charge conjugated conjugated mesons:
print(" CC-neutral beauty mesons")
mesonscc = ppsvc.get(lambda x: (x.isMeson() and x.hasBottom()))
print(mesonscc)
 
# get "stable" particles ( lifetime in excess of 1 micrometer
print(" STABLE particles (c*tau> 1 um)")
stable = ppsvc.get(lambda x: (x.ctau() > micrometer))
print(stable)