"""PDF set writer."""
import io
import pathlib
import re
import numpy as np
import yaml
[docs]
def list_to_str(ls, fmt="%.6e"):
"""Convert a list of numbers to a string.
Parameters
----------
ls : list(float)
list
fmt : str
format string
Returns
-------
str :
final string
"""
return " ".join([fmt % x for x in ls])
[docs]
def array_to_str(ar):
"""Convert an array of numbers to a string.
Parameters
----------
ar : list(list(float))
array
Returns
-------
str :
final string
"""
table = ""
for line in ar:
table += f"{line[0]:.8e} " + list_to_str(line[1:], fmt="%.8e") + "\n"
return table
[docs]
def dump_blocks(name, member, blocks, pdf_type=None):
"""Write LHAPDF data file.
Parameters
----------
name : str or os.PathLike
target name or path
member : int
PDF member
blocks : list(dict)
pdf blocks of data
inherit : str
str to be copied in the head of member files
"""
path_name = pathlib.Path(name)
target = path_name / f"{path_name.stem}_{member:04d}.dat"
target.parent.mkdir(exist_ok=True)
with open(target, "w", encoding="utf-8") as o:
if pdf_type is None:
if member == 0:
o.write("PdfType: central\n")
else:
o.write("PdfType: replica\n")
else:
o.write(pdf_type)
o.write("Format: lhagrid1\n---\n")
for b in blocks:
o.write(list_to_str(b["xgrid"]) + "\n")
o.write(list_to_str(list(np.sqrt(b["mu2grid"]))) + "\n")
o.write(list_to_str(b["pids"], "%d") + "\n")
o.write(array_to_str(b["data"]))
o.write("---\n")
[docs]
def dump_info(name, info):
"""Write LHAPDF info file.
NOTE: Since LHAPDF info files are not truely yaml files,
we have to use a slightly more complicated function to
dump the info file.
Parameters
----------
name : str or os.Pathlike
target name or path
info : dict
info dictionary
"""
path_name = pathlib.Path(name)
target = path_name / f"{path_name.stem}.info"
target.parent.mkdir(exist_ok=True)
# write on string stream to capture output
stream = io.StringIO()
yaml.safe_dump(info, stream, default_flow_style=True, width=100000, line_break="\n")
cnt = stream.getvalue()
# now insert some newlines for each key
new_cnt = re.sub(r", ([A-Za-z_]+\d?):", r"\n\1:", cnt.strip()[1:-1])
with open(target, "w", encoding="utf-8") as o:
o.write(new_cnt)
[docs]
def dump_set(name, info, member_blocks, pdf_type_list=None):
"""Dump a whole set.
Parameters
----------
name : str
target name
info : dict
info dictionary
member_blocks : list(list(dict))
blocks for all members
pdf_type : list(str)
list of strings to be copied in the head of member files
"""
dump_info(name, info)
for mem, blocks in enumerate(member_blocks):
if not isinstance(pdf_type_list, list) or len(pdf_type_list) == 0:
dump_blocks(name, mem, blocks)
else:
dump_blocks(name, mem, blocks, pdf_type=pdf_type_list[mem])