add-xsec-uncert-quadrature-N.py File Reference

Go to the source code of this file.

Namespaces
namespace	add-xsec-uncert-quadrature-N

Variables
str	add-xsec-uncert-quadrature-N.filename_pattern = "pwg-*-st2-stat.dat"
list	add-xsec-uncert-quadrature-N.N_values = [50, 100,200,400,800,1000,1100,1200,1250,1260,1270,1275,1278,1279,1280]
str	add-xsec-uncert-quadrature-N.plot_file_postfix = ""
int	add-xsec-uncert-quadrature-N.num_repeats = 30
	add-xsec-uncert-quadrature-N.file_list = glob.glob(filename_pattern)
list	add-xsec-uncert-quadrature-N.empty_files = [f for f in file_list if os.path.getsize(f) == 0]
	add-xsec-uncert-quadrature-N.num_empty_files = len(empty_files)
dict	add-xsec-uncert-quadrature-N.results = {N: [] for N in N_values}
	add-xsec-uncert-quadrature-N.selected_files = random.sample(file_list, min(N, len(file_list)))
list	add-xsec-uncert-quadrature-N.value_list = []
	add-xsec-uncert-quadrature-N.match = re.search(r"^ grand total total \‍(pos\.-\|neg\.\|\‍):\s+(\d+\.\d+|\d+)\s+\+-\s+(\d+\.\d+E[-+]\d+|\d+\.\d+)", line)
	add-xsec-uncert-quadrature-N.value = float(match.group(2))
	add-xsec-uncert-quadrature-N.quadrature_sum = np.sqrt(np.sum(np.array(value_list) ** 2))/len(value_list)
list	add-xsec-uncert-quadrature-N.upper_bounds = []
list	add-xsec-uncert-quadrature-N.lower_bounds = []
list	add-xsec-uncert-quadrature-N.upper_bound_value_list = value_list[-N:]
list	add-xsec-uncert-quadrature-N.lower_bound_value_list = value_list[:N]
list	add-xsec-uncert-quadrature-N.mean_values = []
list	add-xsec-uncert-quadrature-N.std_dev_values = []
	add-xsec-uncert-quadrature-N.figsize
	add-xsec-uncert-quadrature-N.label
	add-xsec-uncert-quadrature-N.yerr
	add-xsec-uncert-quadrature-N.fmt
	add-xsec-uncert-quadrature-N.capsize
	add-xsec-uncert-quadrature-N.alpha
	add-xsec-uncert-quadrature-N.color
	add-xsec-uncert-quadrature-N.format
	add-xsec-uncert-quadrature-N.dpi
	add-xsec-uncert-quadrature-N.bbox_inches
	add-xsec-uncert-quadrature-N.markersize