PartPropSvc.py

Go to the documentation of this file.

#!/usr/bin/env python3

"""Demo-file for interactive work with new Particle Property Service"""
 
__author__ = "Vanya BELYAEV Ivan.Belyaev@nikhef.nl"
import GaudiPartProp.decorators
 
# =============================================================================
import GaudiPartProp.PartPropAlg
 
# The following is misunderstood by flake8 - the import is needed as it
# has necessary side effects
import GaudiPartProp.Service  # noqa: F401
 
# =============================================================================

from GaudiPartProp.Nodes import (
    CC,
    Baryon,
    Bottom,
    EllMinus,
    EllPlus,
    FiveHalf,
    Hadron,
    Lepton,
    Meson,
    Nu,
    Tensor,
    ThreeHalf,
)
from GaudiPython.Bindings import AppMgr
 
gaudi = AppMgr()
gaudi.initialize()
 
 

def test():
    """Simple function for the test"""
    pps = gaudi.gaudiPartProp()
 
    for l in range(1, 6):
        # get
        lst1 = pps.get(lambda x: l == x.pid().lSpin())
        print(f" {l} == 2*lSpin+1 ")
        print(lst1)
 
    for s in range(0, 2):
        # get
        lst1 = pps.get(lambda x: s == x.pid().sSpin())
        print(f" {s} == 2*sSpin+1 ")
        print(lst1)
 
    for j in range(0, 7):
        # get
        lst1 = pps.get(lambda x: j == x.pid().jSpin())
        print(f" {j} == 2*jSpi+1 ")
        print(lst1)
 
    decays = [
        "B0   -> K- pi+",
        "B+   -> J/psi(1S) ( K*0 -> K+ pi- ) ",
        "D_s+ -> pi+ pi+ pi- ",
        "B_s0 -> ( J/psi(1S) -> mu+ mu- ) ( phi(1020) -> K+ K- ) ",
    ]
 
    for decay in decays:
        print(f"Decay&CC:  {decay} : {pps.cc(decay)} ")
 
    Item = pps.Item
    Items = pps.Items
    Decay = pps.Decay
    Decays = pps.Decays
 
    items = Items()
    items.push_back(Item("K+"))
    items.push_back(Item("K-"))
 
    decay = Decay(Item("B0"), items)
 
    print(decay.validate(pps.svc()))
 
    decays = Decays()
    decays.push_back(decay)
    decays.push_back(decay)
    decays.push_back(decay)
 
    print(decay)
    print(decays)
 
 
def test2():
    pps = gaudi.gaudiPartProp()
 
    nodes = [
        Lepton & ~Meson,
        EllPlus | Nu,
        Lepton & ~Nu,
        EllMinus | CC("B0"),
        Hadron & Tensor,
        Hadron & ThreeHalf,
        Hadron & FiveHalf,
        EllPlus | (Bottom & Baryon),
    ]
 
    for node in nodes:
        node.validate(pps.svc())
 
        lst1 = pps.get(node)
        print(f' selected by node: "{node}" ')
        print(lst1)
 
 
if "__main__" == __name__:
    test()
    test2()