"""Collection of utility functions for the :mod:`hepi` package."""
import json
from typing import List, Tuple
import numpy as np
import hashlib
from particle.converters.bimap import DirectionalMaps
from particle import PDGID
import lhapdf
import pandas as pd
import warnings
[docs]class DictData:
[docs] def __str__(self):
"""Returns attributes as dict as string"""
return str(self.__dict__)
[docs]def LD2DL(l: List, actual_dict=False) -> dict:
"""
Convert a list of objects into a dictionary of lists.
The values of each object are first converted to a `dict` through the `__dict__` attribute.
Args:
l (List) : list of objects.
actual_dict (bool) : objects are already dicts
Returns:
dict : dictionary of numpy arrays.
Examples:
>>> class Param:
... def __init__(self,a,b,c):
... self.a = a
... self.b = b
... self.c = c
>>> LD2DL([ Param(1,2,3), Param(4,5,6) , Param(7,8,9) ])
{'a': array([1, 4, 7]), 'b': array([2, 5, 8]), 'c': array([3, 6, 9])}
"""
# Check l[0] keys in all dictionaries.
for m in l:
md = m if actual_dict else m.__dict__
for k in l[0] if actual_dict else l[0].__dict__:
assert k in md
# switch them
return {
k: np.array([dic[k] if actual_dict else dic.__dict__[k] for dic in l])
for k in (l[0] if actual_dict else l[0].__dict__)
}
[docs]def LD2DF(ld: dict) -> pd.DataFrame:
"""
Convert a `dict` of `list`s to a `pandas.DataFrame`.
"""
return pd.DataFrame.from_dict(ld)
PDG2LaTeXNameMap, LaTeX2PDGNameMap = DirectionalMaps("PDGID",
"LaTexName",
converters=(PDGID, str))
PDG2Name2IDMap, PDGID2NameMap = DirectionalMaps("PDGName",
"PDGID",
converters=(str, PDGID))
[docs]def get_name(pid: int) -> str:
"""
Get the latex name of a particle.
Args:
pid (int) : PDG Monte Carlo identifier for the particle.
Returns:
str: Latex name.
Examples:
>>> get_name(21)
'g'
>>> get_name(1000022)
'\\\\tilde{\\\\chi}_{1}^{0}'
"""
global PDG2LaTeXNameMap
pdgid = PDG2LaTeXNameMap[pid]
return pdgid
[docs]def get_LR_partner(pid: int) -> Tuple[int, int]:
"""Transforms a PDG id to it's left-right partner.
Args:
pid (int) : PDG Monte Carlo identifier for the particle.
Returns:
tuple : First int is -1 for Left and 1 for Right. Second int is the PDG id.
Examples:
>>> get_LR_partner(1000002)
(-1, 2000002)
"""
n = PDGID2NameMap[pid]
if "L" in n:
n = n.replace("L", "R")
return -1, int(PDG2Name2IDMap[n])
if "R" in n:
n = n.replace("R", "L")
return 1, int(PDG2Name2IDMap[n])
if "1" in n:
n = n.replace("1", "2")
return -1, int(PDG2Name2IDMap[n])
if "2" in n:
n = n.replace("2", "1")
return 1, int(PDG2Name2IDMap[n])
return None
[docs]def namehash(n: any) -> str:
"""
Creates a sha256 hash from the objects string representation.
Args:
n (any) : object.
Returns:
str: sha256 of object.
Examples:
>>> p = {'a':1,'b':2}
>>> str(p)
"{'a': 1, 'b': 2}"
>>> namehash(str(p))
'3dffaea891e5dbadb390da33bad65f509dd667779330e2720df8165a253462b8'
>>> namehash(p)
'3dffaea891e5dbadb390da33bad65f509dd667779330e2720df8165a253462b8'
"""
m = hashlib.sha256()
m.update(str(n).encode('utf-8'))
return m.hexdigest()
[docs]def lhapdf_name_to_id(name: str) -> int:
"""
Converts a LHAPDF name to the sets id.
Args:
name (str) : LHAPDF set name.
Returns:
int: id of the LHAPDF set.
Examples:
>>> lhapdf_name_to_id("CT14lo")
13200
"""
if not name in lhapdf.availablePDFSets():
warnings.warn("PDF set '" + name + "' not installed!")
return 0
return lhapdf.getPDFSet(name).lhapdfID