diff --git a/beta_spectrum/spectrum.py b/beta_spectrum/spectrum.py index f244f7c..62a99ea 100644 --- a/beta_spectrum/spectrum.py +++ b/beta_spectrum/spectrum.py @@ -15,6 +15,7 @@ from beta_spectrum.components.screening import ScreeningCorrection from beta_spectrum.components.exchange import ExchangeCorrection from beta_spectrum.components.radiative import RadiativeCorrection from beta_spectrum.utils import T_to_W +from beta_spectrum.constants import ME_MEV @dataclass @@ -250,11 +251,13 @@ class BetaSpectrum: W0 = float(T_to_W(config.endpoint_MeV)) if config.use_phase_space: + # Convert neutrino mass from MeV to natural units + m_nu_natural = config.neutrino_mass_MeV / ME_MEV components.append( PhaseSpace( W0=W0, transition_type=config.transition_type, - m_nu=config.neutrino_mass_MeV, + m_nu=m_nu_natural, logger=logger, ) ) diff --git a/beta_spectrum/visualize/__init__.py b/beta_spectrum/visualize/__init__.py index 52ff51a..7bbe793 100644 --- a/beta_spectrum/visualize/__init__.py +++ b/beta_spectrum/visualize/__init__.py @@ -628,11 +628,9 @@ class BetaSpectrumAnalyzer: total = self.total_spectrum(normalize=True) components = self.components + # Only include enabled corrections that are not basic components + # Basic components (phase_space and fermi) are always used and should not be listed enabled = [] - if self.config.use_phase_space: - enabled.append("phase_space") - if self.config.use_fermi: - enabled.append("fermi") if self.config.use_screening: enabled.append("screening") if self.config.use_finite_size: @@ -647,10 +645,19 @@ class BetaSpectrumAnalyzer: parent_symbol = self._element_symbol(self.config.Z_parent) daughter_symbol = self._element_symbol(self.config.Z_daughter) + # Get git commit hash + git_commit = get_git_short_hash() + + # Format source line with paceENSDF URL if applicable + if source_type == "paceENSDF": + source_line = "# source: paceENSDF (https://pypi.org/project/paceENSDF/)" + else: + source_line = f"# source: {source_type}" + header_lines = [ - f"# beta-spectrum v{__import__('beta_spectrum').__version__}", + f"# beta-spectrum v{__import__('beta_spectrum').__version__}(build.{git_commit})", f"# timestamp: {datetime.now(timezone.utc).isoformat()}", - f"# source: {source_type}", + source_line, f"# nuclide: {parent_symbol}{self.config.A_number} -> {daughter_symbol}{self.config.A_number} (Z={self.config.Z_parent}->{self.config.Z_daughter})", f"# basic_endpoint: {self.config.endpoint_MeV * 1000:.1f} keV", ] @@ -667,49 +674,82 @@ class BetaSpectrumAnalyzer: header_lines.append(f"# corrections: {', '.join(enabled)}") header_lines.append(f"# e_step: {self.config.e_step_MeV:.4f} MeV") + header_lines.append(f"# n_points: {len(self.energies_MeV)}") + + # Add normalization mode + header_lines.append("# branch normalization: intensity") + + # Separate branches into active and excluded lists + active_branches = [] + excluded_branches = [] - # Add transition types for all branches - transition_types = [] for i, branch in enumerate(self._unified_branches): - transition_types.append(f"branch_{i + 1}:{branch.transition_type}") - header_lines.append(f"# transition_types: {'; '.join(transition_types)}") + branch_endpoint_keV = branch.endpoint_MeV * 1000 + neutrino_mass_keV = self.config.neutrino_mass_MeV * 1000 - # Branch info — unified view - for i, branch in enumerate(self._unified_branches): - header_lines.append( - f"# branch_{i + 1}: E₀={branch.endpoint_MeV * 1000:.1f} keV, " - f"transition={branch.transition_type}, " - f"intensity={branch.intensity:.4f}" - ) + # Check if branch is excluded due to neutrino mass + if branch_endpoint_keV <= neutrino_mass_keV: + excluded_branches.append( + f" branch_{i + 1}: E₀={branch_endpoint_keV:.1f} keV, transition={branch.transition_type}, intensity={branch.intensity:.4f} - excluded: endpoint < m_nu" + ) + # Check if branch is excluded due to intensity cutoff + elif branch.intensity < self.config.intensity_cutoff: + excluded_branches.append( + f" branch_{i + 1}: E₀={branch_endpoint_keV:.1f} keV, transition={branch.transition_type}, intensity={branch.intensity:.4f} - excluded: intensity < cutoff" + ) + else: + active_branches.append( + f" branch_{i + 1}: E₀={branch_endpoint_keV:.1f} keV, transition={branch.transition_type}, intensity={branch.intensity:.4f}" + ) - header_lines.append(f"# git_commit: {get_git_short_hash()}") + # Add active branches list + if active_branches: + header_lines.append("# active branches:") + for branch_info in active_branches: + header_lines.append(f"# {branch_info}") + # Add excluded branches list if there are any + if excluded_branches: + header_lines.append("# excluded branches:") + for branch_info in excluded_branches: + header_lines.append(f"# {branch_info}") + + # Create data dictionary with concise column headers data: Dict[str, np.ndarray] = { "energy_MeV": self.energies_MeV, - "spectrum": total, + "total_spectrum": total, } if csv_verbosity == "full": - # Per-branch spectra and components + # Per-branch spectra and components with concise branch naming for i in range(self._n_branches): - data[f"branch_{i + 1}_spectrum"] = self._get_branch_spectrum(i) + data[f"branch {i + 1}"] = self._get_branch_spectrum(i) for comp_name in ["PhaseSpace", "Radiative"]: key = self._get_component_key(i, comp_name) if key in components: - data[f"branch_{i + 1}_{comp_name}"] = components[key] + data[f"branch {i + 1}_{comp_name}"] = components[key] # Universal components (no prefix) for name, values in components.items(): if not name.startswith("branch_"): data[name] = values elif csv_verbosity == "concise": - # Only total spectrum and normalized branch components + # Only total spectrum and normalized branch components with concise naming for i in range(self._n_branches): - data[f"branch_{i + 1}_spectrum_normalized"] = ( - self._get_branch_spectrum(i) / np.sum(self._get_branch_spectrum(i)) - if np.sum(self._get_branch_spectrum(i)) > 0 - else self._get_branch_spectrum(i) - ) + branch_spectrum = self._get_branch_spectrum(i) + # Check if branch is degenerate (neutrino mass >= branch endpoint) + branch_endpoint_keV = self._unified_branches[i].endpoint_MeV * 1000 + neutrino_mass_keV = self.config.neutrino_mass_MeV * 1000 + if branch_endpoint_keV <= neutrino_mass_keV: + # Branch is degenerate, spectrum should be all zeros + data[f"branch {i + 1}"] = np.zeros_like(branch_spectrum) + else: + # Branch is valid, normalize as usual + sum_spectrum = np.sum(branch_spectrum) + if sum_spectrum > 0: + data[f"branch {i + 1}"] = branch_spectrum / sum_spectrum + else: + data[f"branch {i + 1}"] = branch_spectrum df = pd.DataFrame(data) df.to_csv(filename, index=False, float_format="%.4e") diff --git a/create_comparison_plot.py b/create_comparison_plot.py deleted file mode 100644 index 88a2123..0000000 --- a/create_comparison_plot.py +++ /dev/null @@ -1,144 +0,0 @@ -#!/usr/bin/env python3 -""" -Script to create a comparison plot of total beta spectra for edge neutrino masses (1.5 keV and 2.5 keV) -""" - -import pandas as pd -import matplotlib.pyplot as plt -import numpy as np - -# Input files -mass_1_5_csv = "/home/unzhakov/workspace/massive_nu_pr/pr144_nu_mass_1_50keV.csv" -mass_2_45_csv = "/home/unzhakov/workspace/massive_nu_pr/pr144_nu_mass_2_45keV.csv" - -# Output file -output_png = ( - "/home/unzhakov/workspace/massive_nu_pr/pr144_comparison_1_5_vs_2_45_keV.png" -) - -print(f"Reading data from {mass_1_5_csv} and {mass_2_45_csv}...") - - -# Read CSV files, skipping metadata headers -# The metadata headers start with '#' -def read_csv_with_metadata(filepath): - with open(filepath, "r") as f: - lines = f.readlines() - - # Find the line where data starts (not starting with '#') - data_start_line = 0 - for i, line in enumerate(lines): - if not line.startswith("#"): - data_start_line = i - break - - # Read the data - df = pd.read_csv(filepath, skiprows=data_start_line) - return df - - -# Read the data -df_1_5 = read_csv_with_metadata(mass_1_5_csv) -df_2_45 = read_csv_with_metadata(mass_2_45_csv) - -print("Data loaded successfully.") - -# Extract energy and spectrum data -energies_1_5 = df_1_5["energy_MeV"].values -spectrum_1_5 = df_1_5["spectrum"].values - -energies_2_45 = df_2_45["energy_MeV"].values -spectrum_2_45 = df_2_45["spectrum"].values - -print( - f"Energy range: {energies_1_5[0] * 1000:.1f} eV to {energies_1_5[-1] * 1000:.1f} eV" -) - -# Create comparison plot -fig, ax = plt.subplots(figsize=(12, 8)) - -# Plot both spectra -ax.plot(energies_1_5 * 1000, spectrum_1_5, "b-", lw=2, label="m_ν = 1.5 keV") -ax.plot(energies_2_45 * 1000, spectrum_2_45, "r--", lw=2, label="m_ν = 2.45 keV") - -# Set labels and title -ax.set_xlabel(r"Electron kinetic energy $E$ [keV]", fontsize=12) -ax.set_ylabel("Normalized Spectrum", fontsize=12) -ax.set_title( - "Pr144 -> Nd144 Beta Decay: Comparison of Neutrino Mass Effects\n" - "(1.5 keV vs 2.45 keV)", - fontsize=14, -) - -# Set log scale for y-axis to see the endpoint region clearly -ax.set_yscale("log") - -# Add grid -ax.grid(True, alpha=0.3) - -# Add legend -ax.legend(fontsize=12) - -# Focus on the endpoint region (last part of the spectrum) -# Find the energy range near the endpoint -max_energy_keV = max(energies_1_5) * 1000 -# Show last 10 keV of the spectrum -start_energy_keV = max_energy_keV - 10.0 -ax.set_xlim(start_energy_keV, max_energy_keV) - -# Save the plot -plt.tight_layout() -plt.savefig(output_png, dpi=150, bbox_inches="tight") -print(f"Comparison plot saved to {output_png}") - -# Also create a linear scale plot for the endpoint region -fig_linear, ax_linear = plt.subplots(figsize=(12, 6)) - -# Plot both spectra in linear scale -ax_linear.plot(energies_1_5 * 1000, spectrum_1_5, "b-", lw=2, label="m_ν = 1.5 keV") -ax_linear.plot(energies_2_45 * 1000, spectrum_2_45, "r--", lw=2, label="m_ν = 2.45 keV") - -# Set labels and title -ax_linear.set_xlabel(r"Electron kinetic energy $E$ [keV]", fontsize=12) -ax_linear.set_ylabel("Normalized Spectrum", fontsize=12) -ax_linear.set_title( - "Pr144 -> Nd144 Beta Decay: Comparison of Neutrino Mass Effects (Endpoint Region)\n" - "(1.5 keV vs 2.45 keV) - Linear Scale", - fontsize=14, -) - -# Set linear scale for y-axis -ax_linear.set_yscale("linear") - -# Set y-axis limits to show the variation in the spectrum -# Find the min and max values in the endpoint region to set appropriate limits -min_spec_1_5 = np.min(spectrum_1_5[-100:]) # Last 100 points -max_spec_1_5 = np.max(spectrum_1_5[-100:]) -min_spec_2_45 = np.min(spectrum_2_45[-100:]) -max_spec_2_45 = np.max(spectrum_2_45[-100:]) - -overall_min = ( - min(min_spec_1_5, min_spec_2_45) * 0.9 - if min(min_spec_1_5, min_spec_2_45) > 0 - else 0 -) -overall_max = max(max_spec_1_5, max_spec_2_45) * 1.1 - -ax_linear.set_ylim(overall_min, overall_max) - -# Add grid -ax_linear.grid(True, alpha=0.3) - -# Add legend -ax_linear.legend(fontsize=12) - -# Focus on the endpoint region (last part of the spectrum) -ax_linear.set_xlim(start_energy_keV, max_energy_keV) - -# Save the linear plot -plt.tight_layout() -linear_output_png = ( - "/home/unzhakov/workspace/massive_nu_pr/pr144_comparison_1_5_vs_2_45_keV_linear.png" -) -plt.savefig(linear_output_png, dpi=150, bbox_inches="tight") -print(f"Linear comparison plot saved to {linear_output_png}") diff --git a/run_pr144_calculations.py b/run_pr144_calculations.py deleted file mode 100644 index 6abe9f1..0000000 --- a/run_pr144_calculations.py +++ /dev/null @@ -1,86 +0,0 @@ -#!/usr/bin/env python3 -""" -Script to run Pr144 -> Nd144 beta decay calculations for neutrino mass range 1.5-2.5 keV -""" - -import subprocess -import os - -# Create output directory -output_dir = "/home/unzhakov/workspace/massive_nu_pr" -os.makedirs(output_dir, exist_ok=True) - -# Neutrino mass range: 1.5 to 2.5 keV with step of 50 eV -# Convert to MeV: 1.5 keV = 0.0015 MeV, 2.5 keV = 0.0025 MeV, 50 eV = 0.00005 MeV -start_mass_meV = 1.5 # keV -end_mass_meV = 2.5 # keV -step_mass_meV = 0.05 # keV (50 eV) - -# Convert to MeV -start_mass_MeV = start_mass_meV / 1000.0 -end_mass_MeV = end_mass_meV / 1000.0 -step_mass_MeV = step_mass_meV / 1000.0 - -# Energy step: 100 eV = 0.0001 MeV -e_step_MeV = 0.0001 - -# Intensity cutoff: 0.001 -intensity_cutoff = 0.001 - -# Generate neutrino mass values -mass_values = [] -mass = start_mass_MeV -while mass <= end_mass_MeV: - mass_values.append(mass) - mass += step_mass_MeV - -print(f"Running calculations for {len(mass_values)} neutrino mass values...") -print( - f"Neutrino mass range: {start_mass_meV} - {end_mass_meV} keV with step {step_mass_meV} keV" -) -print(f"Energy step: {e_step_MeV * 1000} eV") -print(f"Intensity cutoff: {intensity_cutoff}") - -# Run calculations for each neutrino mass -for i, mass_MeV in enumerate(mass_values): - mass_keV = mass_MeV * 1000.0 - print( - f"\n[{i + 1}/{len(mass_values)}] Calculating for neutrino mass = {mass_keV:.2f} keV ({mass_MeV:.6f} MeV)" - ) - - # Generate filename - mass_str = f"{mass_keV:.2f}".replace(".", "_") - csv_filename = f"{output_dir}/pr144_nu_mass_{mass_str}keV.csv" - png_filename = f"{output_dir}/pr144_nu_mass_{mass_str}keV.png" - - # Command to run - cmd = [ - "python", - "-m", - "beta_spectrum.cli", - "--nuclide", - "Pr144", - "--neutrino-mass", - str(mass_MeV), - "--e-step", - str(e_step_MeV), - "--intensity-cutoff", - str(intensity_cutoff), - "--output", - csv_filename, - "--plot", - png_filename, - "-q", # quiet mode - ] - - print(f"Running: {' '.join(cmd)}") - result = subprocess.run(cmd, capture_output=True, text=True) - - if result.returncode != 0: - print(f"Error in calculation for mass {mass_keV:.2f} keV:") - print(f"stdout: {result.stdout}") - print(f"stderr: {result.stderr}") - else: - print(f"Successfully completed calculation for mass {mass_keV:.2f} keV") - -print("\nAll calculations completed!") diff --git a/test_neutrino_mass.csv b/test_neutrino_mass.csv deleted file mode 100644 index ffdb224..0000000 --- a/test_neutrino_mass.csv +++ /dev/null @@ -1,310 +0,0 @@ -# beta-spectrum v0.4.0 -# timestamp: 2026-07-01T13:25:57.578522+00:00 -# source: paceENSDF -# nuclide: Tc99 -> Ru99 (Z=43->44) -# endpoint: 297.5 keV -# transition: F2 -# corrections: phase_space, fermi, screening, finite_size, charge_dist, radiative, exchange -# e_step: 0.0010 MeV -# neutrino_mass_MeV: 0.010000 -# branch_1: E₀=297.5 keV, transition=F2, intensity=1.0000 -# branch_2: E₀=208.0 keV, transition=F2U, intensity=0.0000 -# git_commit: 406fd57 -energy_MeV,spectrum,branch_1_spectrum,branch_1_PhaseSpace,branch_1_Radiative,branch_2_spectrum,branch_2_PhaseSpace,branch_2_Radiative,Fermi,FiniteSizeL0,ChargeDistributionU,Screening,Exchange -1.0000e-03,5.0134e+00,7.2723e-03,7.1882e-03,1.0117e+00,2.8198e-04,2.7884e-04,1.0113e+00,6.2469e+01,1.0154e+00,1.0000e+00,9.9894e-01,1.0496e+00 -2.0000e-03,5.0002e+00,1.0289e-02,1.0171e-02,1.0116e+00,3.9175e-04,3.8744e-04,1.0111e+00,4.4215e+01,1.0154e+00,1.0000e+00,9.9760e-01,1.0468e+00 -3.0000e-03,4.9742e+00,1.2607e-02,1.2464e-02,1.0115e+00,4.7190e-04,4.6675e-04,1.0110e+00,3.6137e+01,1.0154e+00,1.0000e+00,9.9594e-01,1.0416e+00 -4.0000e-03,4.9540e+00,1.4564e-02,1.4399e-02,1.0114e+00,5.3660e-04,5.3077e-04,1.0110e+00,3.1325e+01,1.0154e+00,1.0000e+00,9.9395e-01,1.0380e+00 -5.0000e-03,4.9369e+00,1.6290e-02,1.6107e-02,1.0114e+00,5.9152e-04,5.8513e-04,1.0109e+00,2.8045e+01,1.0154e+00,1.0000e+00,9.9164e-01,1.0353e+00 -6.0000e-03,4.9210e+00,1.7853e-02,1.7654e-02,1.0113e+00,6.3970e-04,6.3283e-04,1.0108e+00,2.5626e+01,1.0154e+00,1.0000e+00,9.8901e-01,1.0333e+00 -7.0000e-03,4.9054e+00,1.9293e-02,1.9078e-02,1.0113e+00,6.8302e-04,6.7572e-04,1.0108e+00,2.3748e+01,1.0154e+00,1.0000e+00,9.8609e-01,1.0316e+00 -8.0000e-03,4.8897e+00,2.0634e-02,2.0405e-02,1.0112e+00,7.2274e-04,7.1506e-04,1.0107e+00,2.2235e+01,1.0154e+00,1.0000e+00,9.8290e-01,1.0302e+00 -9.0000e-03,4.8802e+00,2.1895e-02,2.1653e-02,1.0112e+00,7.5978e-04,7.5174e-04,1.0107e+00,2.0983e+01,1.0154e+00,1.0000e+00,9.8079e-01,1.0290e+00 -1.0000e-02,4.8762e+00,2.3089e-02,2.2836e-02,1.0111e+00,7.9484e-04,7.8647e-04,1.0107e+00,1.9925e+01,1.0154e+00,1.0000e+00,9.7968e-01,1.0279e+00 -1.1000e-02,4.8732e+00,2.4227e-02,2.3961e-02,1.0111e+00,8.2847e-04,8.1977e-04,1.0106e+00,1.9016e+01,1.0154e+00,1.0000e+00,9.7867e-01,1.0269e+00 -1.2000e-02,4.8709e+00,2.5315e-02,2.5038e-02,1.0110e+00,8.6110e-04,8.5209e-04,1.0106e+00,1.8224e+01,1.0154e+00,1.0000e+00,9.7776e-01,1.0260e+00 -1.3000e-02,4.8695e+00,2.6359e-02,2.6072e-02,1.0110e+00,8.9308e-04,8.8378e-04,1.0105e+00,1.7525e+01,1.0154e+00,1.0000e+00,9.7694e-01,1.0252e+00 -1.4000e-02,4.8688e+00,2.7365e-02,2.7068e-02,1.0110e+00,9.2474e-04,9.1514e-04,1.0105e+00,1.6904e+01,1.0153e+00,1.0000e+00,9.7623e-01,1.0244e+00 -1.5000e-02,4.8686e+00,2.8336e-02,2.8030e-02,1.0109e+00,9.5631e-04,9.4642e-04,1.0105e+00,1.6347e+01,1.0153e+00,1.0000e+00,9.7560e-01,1.0237e+00 -1.6000e-02,4.8691e+00,2.9276e-02,2.8961e-02,1.0109e+00,9.8802e-04,9.7784e-04,1.0104e+00,1.5843e+01,1.0153e+00,1.0000e+00,9.7505e-01,1.0230e+00 -1.7000e-02,4.8700e+00,3.0188e-02,2.9864e-02,1.0109e+00,1.0201e-03,1.0096e-03,1.0104e+00,1.5385e+01,1.0153e+00,1.0000e+00,9.7457e-01,1.0223e+00 -1.8000e-02,4.8714e+00,3.1074e-02,3.0741e-02,1.0108e+00,1.0526e-03,1.0418e-03,1.0103e+00,1.4966e+01,1.0153e+00,1.0000e+00,9.7416e-01,1.0217e+00 -1.9000e-02,4.8731e+00,3.1936e-02,3.1595e-02,1.0108e+00,1.0858e-03,1.0747e-03,1.0103e+00,1.4581e+01,1.0153e+00,1.0000e+00,9.7381e-01,1.0211e+00 -2.0000e-02,4.8752e+00,3.2776e-02,3.2427e-02,1.0108e+00,1.1197e-03,1.1083e-03,1.0103e+00,1.4226e+01,1.0153e+00,1.0000e+00,9.7350e-01,1.0205e+00 -2.1000e-02,4.8775e+00,3.3595e-02,3.3238e-02,1.0107e+00,1.1545e-03,1.1428e-03,1.0102e+00,1.3897e+01,1.0153e+00,1.0000e+00,9.7324e-01,1.0200e+00 -2.2000e-02,4.8800e+00,3.4396e-02,3.4031e-02,1.0107e+00,1.1903e-03,1.1782e-03,1.0102e+00,1.3591e+01,1.0153e+00,1.0000e+00,9.7301e-01,1.0194e+00 -2.3000e-02,4.8828e+00,3.5178e-02,3.4806e-02,1.0107e+00,1.2271e-03,1.2147e-03,1.0102e+00,1.3305e+01,1.0153e+00,1.0000e+00,9.7281e-01,1.0189e+00 -2.4000e-02,4.8857e+00,3.5944e-02,3.5565e-02,1.0107e+00,1.2650e-03,1.2523e-03,1.0102e+00,1.3039e+01,1.0153e+00,1.0000e+00,9.7263e-01,1.0184e+00 -2.5000e-02,4.8888e+00,3.6694e-02,3.6308e-02,1.0106e+00,1.3042e-03,1.2911e-03,1.0101e+00,1.2788e+01,1.0153e+00,1.0000e+00,9.7248e-01,1.0180e+00 -2.6000e-02,4.8920e+00,3.7428e-02,3.7036e-02,1.0106e+00,1.3446e-03,1.3311e-03,1.0101e+00,1.2553e+01,1.0153e+00,1.0000e+00,9.7235e-01,1.0175e+00 -2.7000e-02,4.8953e+00,3.8149e-02,3.7750e-02,1.0106e+00,1.3863e-03,1.3725e-03,1.0101e+00,1.2331e+01,1.0153e+00,1.0000e+00,9.7223e-01,1.0170e+00 -2.8000e-02,4.8986e+00,3.8856e-02,3.8450e-02,1.0106e+00,1.4293e-03,1.4151e-03,1.0100e+00,1.2121e+01,1.0153e+00,1.0000e+00,9.7213e-01,1.0166e+00 -2.9000e-02,4.9021e+00,3.9550e-02,3.9138e-02,1.0105e+00,1.4738e-03,1.4591e-03,1.0100e+00,1.1923e+01,1.0153e+00,1.0000e+00,9.7204e-01,1.0162e+00 -3.0000e-02,4.9055e+00,4.0232e-02,3.9813e-02,1.0105e+00,1.5196e-03,1.5046e-03,1.0100e+00,1.1735e+01,1.0153e+00,1.0000e+00,9.7196e-01,1.0158e+00 -3.1000e-02,4.9090e+00,4.0901e-02,4.0477e-02,1.0105e+00,1.5669e-03,1.5514e-03,1.0100e+00,1.1557e+01,1.0153e+00,1.0000e+00,9.7188e-01,1.0154e+00 -3.2000e-02,4.9125e+00,4.1560e-02,4.1130e-02,1.0105e+00,1.6156e-03,1.5997e-03,1.0099e+00,1.1387e+01,1.0153e+00,1.0000e+00,9.7182e-01,1.0150e+00 -3.3000e-02,4.9161e+00,4.2207e-02,4.1772e-02,1.0104e+00,1.6658e-03,1.6495e-03,1.0099e+00,1.1225e+01,1.0153e+00,1.0000e+00,9.7176e-01,1.0146e+00 -3.4000e-02,4.9196e+00,4.2844e-02,4.2403e-02,1.0104e+00,1.7175e-03,1.7007e-03,1.0099e+00,1.1071e+01,1.0153e+00,1.0000e+00,9.7170e-01,1.0143e+00 -3.5000e-02,4.9231e+00,4.3471e-02,4.3024e-02,1.0104e+00,1.7707e-03,1.7534e-03,1.0098e+00,1.0924e+01,1.0152e+00,1.0000e+00,9.7165e-01,1.0139e+00 -3.6000e-02,4.9266e+00,4.4088e-02,4.3636e-02,1.0104e+00,1.8254e-03,1.8076e-03,1.0098e+00,1.0783e+01,1.0152e+00,1.0000e+00,9.7161e-01,1.0136e+00 -3.7000e-02,4.9301e+00,4.4695e-02,4.4238e-02,1.0103e+00,1.8816e-03,1.8633e-03,1.0098e+00,1.0648e+01,1.0152e+00,1.0000e+00,9.7157e-01,1.0132e+00 -3.8000e-02,4.9336e+00,4.5293e-02,4.4831e-02,1.0103e+00,1.9392e-03,1.9205e-03,1.0098e+00,1.0518e+01,1.0152e+00,1.0000e+00,9.7153e-01,1.0129e+00 -3.9000e-02,4.9370e+00,4.5882e-02,4.5415e-02,1.0103e+00,1.9984e-03,1.9791e-03,1.0097e+00,1.0394e+01,1.0152e+00,1.0000e+00,9.7150e-01,1.0126e+00 -4.0000e-02,4.9404e+00,4.6462e-02,4.5990e-02,1.0103e+00,2.0590e-03,2.0391e-03,1.0097e+00,1.0275e+01,1.0152e+00,1.0000e+00,9.7147e-01,1.0123e+00 -4.1000e-02,4.9438e+00,4.7033e-02,4.6556e-02,1.0102e+00,2.1210e-03,2.1007e-03,1.0097e+00,1.0160e+01,1.0152e+00,1.0000e+00,9.7144e-01,1.0120e+00 -4.2000e-02,4.9471e+00,4.7597e-02,4.7115e-02,1.0102e+00,2.1845e-03,2.1636e-03,1.0097e+00,1.0050e+01,1.0152e+00,1.0000e+00,9.7142e-01,1.0117e+00 -4.3000e-02,4.9503e+00,4.8151e-02,4.7665e-02,1.0102e+00,2.2494e-03,2.2280e-03,1.0096e+00,9.9439e+00,1.0152e+00,1.0000e+00,9.7140e-01,1.0114e+00 -4.4000e-02,4.9534e+00,4.8698e-02,4.8208e-02,1.0102e+00,2.3158e-03,2.2937e-03,1.0096e+00,9.8415e+00,1.0152e+00,1.0000e+00,9.7138e-01,1.0111e+00 -4.5000e-02,4.9565e+00,4.9237e-02,4.8742e-02,1.0102e+00,2.3834e-03,2.3608e-03,1.0096e+00,9.7428e+00,1.0152e+00,1.0000e+00,9.7136e-01,1.0108e+00 -4.6000e-02,4.9596e+00,4.9768e-02,4.9269e-02,1.0101e+00,2.4525e-03,2.4292e-03,1.0096e+00,9.6474e+00,1.0152e+00,1.0000e+00,9.7135e-01,1.0106e+00 -4.7000e-02,4.9625e+00,5.0292e-02,4.9788e-02,1.0101e+00,2.5228e-03,2.4989e-03,1.0095e+00,9.5552e+00,1.0152e+00,1.0000e+00,9.7133e-01,1.0103e+00 -4.8000e-02,4.9653e+00,5.0808e-02,5.0300e-02,1.0101e+00,2.5944e-03,2.5700e-03,1.0095e+00,9.4661e+00,1.0152e+00,1.0000e+00,9.7132e-01,1.0100e+00 -4.9000e-02,4.9681e+00,5.1317e-02,5.0805e-02,1.0101e+00,2.6673e-03,2.6422e-03,1.0095e+00,9.3799e+00,1.0152e+00,1.0000e+00,9.7132e-01,1.0098e+00 -5.0000e-02,4.9708e+00,5.1818e-02,5.1302e-02,1.0101e+00,2.7414e-03,2.7156e-03,1.0095e+00,9.2965e+00,1.0152e+00,1.0000e+00,9.7131e-01,1.0096e+00 -5.1000e-02,4.9733e+00,5.2312e-02,5.1793e-02,1.0100e+00,2.8166e-03,2.7903e-03,1.0094e+00,9.2157e+00,1.0152e+00,1.0000e+00,9.7130e-01,1.0093e+00 -5.2000e-02,4.9758e+00,5.2799e-02,5.2276e-02,1.0100e+00,2.8930e-03,2.8660e-03,1.0094e+00,9.1374e+00,1.0152e+00,1.0000e+00,9.7130e-01,1.0091e+00 -5.3000e-02,4.9782e+00,5.3279e-02,5.2752e-02,1.0100e+00,2.9705e-03,2.9428e-03,1.0094e+00,9.0614e+00,1.0152e+00,1.0000e+00,9.7130e-01,1.0089e+00 -5.4000e-02,4.9804e+00,5.3752e-02,5.3222e-02,1.0100e+00,3.0490e-03,3.0207e-03,1.0094e+00,8.9877e+00,1.0151e+00,1.0000e+00,9.7130e-01,1.0087e+00 -5.5000e-02,4.9825e+00,5.4218e-02,5.3684e-02,1.0099e+00,3.1285e-03,3.0995e-03,1.0093e+00,8.9161e+00,1.0151e+00,1.0000e+00,9.7130e-01,1.0084e+00 -5.6000e-02,4.9845e+00,5.4678e-02,5.4140e-02,1.0099e+00,3.2090e-03,3.1793e-03,1.0093e+00,8.8466e+00,1.0151e+00,1.0000e+00,9.7130e-01,1.0082e+00 -5.7000e-02,4.9864e+00,5.5130e-02,5.4590e-02,1.0099e+00,3.2903e-03,3.2600e-03,1.0093e+00,8.7791e+00,1.0151e+00,1.0000e+00,9.7131e-01,1.0080e+00 -5.8000e-02,4.9882e+00,5.5576e-02,5.5032e-02,1.0099e+00,3.3726e-03,3.3416e-03,1.0093e+00,8.7134e+00,1.0151e+00,1.0000e+00,9.7132e-01,1.0078e+00 -5.9000e-02,4.9898e+00,5.6015e-02,5.5468e-02,1.0099e+00,3.4556e-03,3.4239e-03,1.0093e+00,8.6496e+00,1.0151e+00,1.0000e+00,9.7132e-01,1.0076e+00 -6.0000e-02,4.9912e+00,5.6448e-02,5.5897e-02,1.0098e+00,3.5394e-03,3.5070e-03,1.0092e+00,8.5874e+00,1.0151e+00,1.0000e+00,9.7133e-01,1.0074e+00 -6.1000e-02,4.9926e+00,5.6873e-02,5.6320e-02,1.0098e+00,3.6239e-03,3.5908e-03,1.0092e+00,8.5269e+00,1.0151e+00,1.0000e+00,9.7134e-01,1.0073e+00 -6.2000e-02,4.9938e+00,5.7293e-02,5.6737e-02,1.0098e+00,3.7090e-03,3.6753e-03,1.0092e+00,8.4680e+00,1.0151e+00,1.0000e+00,9.7135e-01,1.0071e+00 -6.3000e-02,4.9948e+00,5.7705e-02,5.7146e-02,1.0098e+00,3.7947e-03,3.7603e-03,1.0092e+00,8.4106e+00,1.0151e+00,1.0000e+00,9.7137e-01,1.0069e+00 -6.4000e-02,4.9957e+00,5.8111e-02,5.7550e-02,1.0098e+00,3.8810e-03,3.8458e-03,1.0091e+00,8.3546e+00,1.0151e+00,1.0000e+00,9.7138e-01,1.0067e+00 -6.5000e-02,4.9964e+00,5.8511e-02,5.7947e-02,1.0097e+00,3.9677e-03,3.9318e-03,1.0091e+00,8.3000e+00,1.0151e+00,1.0000e+00,9.7139e-01,1.0066e+00 -6.6000e-02,4.9969e+00,5.8904e-02,5.8337e-02,1.0097e+00,4.0548e-03,4.0183e-03,1.0091e+00,8.2468e+00,1.0151e+00,1.0000e+00,9.7141e-01,1.0064e+00 -6.7000e-02,4.9973e+00,5.9291e-02,5.8721e-02,1.0097e+00,4.1422e-03,4.1050e-03,1.0091e+00,8.1949e+00,1.0151e+00,1.0000e+00,9.7143e-01,1.0062e+00 -6.8000e-02,4.9976e+00,5.9671e-02,5.9099e-02,1.0097e+00,4.2299e-03,4.1921e-03,1.0090e+00,8.1442e+00,1.0151e+00,1.0000e+00,9.7145e-01,1.0061e+00 -6.9000e-02,4.9976e+00,6.0045e-02,5.9470e-02,1.0097e+00,4.3179e-03,4.2793e-03,1.0090e+00,8.0947e+00,1.0151e+00,1.0000e+00,9.7146e-01,1.0059e+00 -7.0000e-02,4.9975e+00,6.0412e-02,5.9835e-02,1.0096e+00,4.4060e-03,4.3668e-03,1.0090e+00,8.0463e+00,1.0151e+00,1.0000e+00,9.7148e-01,1.0058e+00 -7.1000e-02,4.9972e+00,6.0773e-02,6.0193e-02,1.0096e+00,4.4942e-03,4.4543e-03,1.0090e+00,7.9991e+00,1.0151e+00,1.0000e+00,9.7150e-01,1.0056e+00 -7.2000e-02,4.9967e+00,6.1127e-02,6.0546e-02,1.0096e+00,4.5824e-03,4.5418e-03,1.0089e+00,7.9529e+00,1.0150e+00,1.0000e+00,9.7153e-01,1.0055e+00 -7.3000e-02,4.9960e+00,6.1475e-02,6.0891e-02,1.0096e+00,4.6706e-03,4.6293e-03,1.0089e+00,7.9078e+00,1.0150e+00,1.0000e+00,9.7155e-01,1.0054e+00 -7.4000e-02,4.9951e+00,6.1816e-02,6.1231e-02,1.0096e+00,4.7587e-03,4.7168e-03,1.0089e+00,7.8637e+00,1.0150e+00,1.0000e+00,9.7157e-01,1.0052e+00 -7.5000e-02,4.9940e+00,6.2151e-02,6.1563e-02,1.0095e+00,4.8466e-03,4.8040e-03,1.0089e+00,7.8205e+00,1.0150e+00,1.0000e+00,9.7159e-01,1.0051e+00 -7.6000e-02,4.9927e+00,6.2479e-02,6.1890e-02,1.0095e+00,4.9342e-03,4.8910e-03,1.0088e+00,7.7782e+00,1.0150e+00,1.0000e+00,9.7162e-01,1.0050e+00 -7.7000e-02,4.9913e+00,6.2801e-02,6.2210e-02,1.0095e+00,5.0216e-03,4.9777e-03,1.0088e+00,7.7369e+00,1.0150e+00,1.0000e+00,9.7164e-01,1.0048e+00 -7.8000e-02,4.9896e+00,6.3117e-02,6.2524e-02,1.0095e+00,5.1086e-03,5.0640e-03,1.0088e+00,7.6964e+00,1.0150e+00,1.0000e+00,9.7167e-01,1.0047e+00 -7.9000e-02,4.9877e+00,6.3425e-02,6.2831e-02,1.0095e+00,5.1951e-03,5.1499e-03,1.0088e+00,7.6568e+00,1.0150e+00,1.0000e+00,9.7170e-01,1.0046e+00 -8.0000e-02,4.9856e+00,6.3728e-02,6.3132e-02,1.0094e+00,5.2811e-03,5.2353e-03,1.0087e+00,7.6180e+00,1.0150e+00,1.0000e+00,9.7172e-01,1.0045e+00 -8.1000e-02,4.9833e+00,6.4024e-02,6.3426e-02,1.0094e+00,5.3665e-03,5.3201e-03,1.0087e+00,7.5799e+00,1.0150e+00,1.0000e+00,9.7175e-01,1.0044e+00 -8.2000e-02,4.9807e+00,6.4313e-02,6.3714e-02,1.0094e+00,5.4512e-03,5.4042e-03,1.0087e+00,7.5427e+00,1.0150e+00,1.0000e+00,9.7178e-01,1.0043e+00 -8.3000e-02,4.9780e+00,6.4596e-02,6.3995e-02,1.0094e+00,5.5352e-03,5.4876e-03,1.0087e+00,7.5061e+00,1.0150e+00,1.0000e+00,9.7181e-01,1.0041e+00 -8.4000e-02,4.9750e+00,6.4872e-02,6.4270e-02,1.0094e+00,5.6184e-03,5.5703e-03,1.0086e+00,7.4703e+00,1.0150e+00,1.0000e+00,9.7184e-01,1.0040e+00 -8.5000e-02,4.9718e+00,6.5141e-02,6.4538e-02,1.0093e+00,5.7008e-03,5.6520e-03,1.0086e+00,7.4352e+00,1.0150e+00,1.0000e+00,9.7187e-01,1.0039e+00 -8.6000e-02,4.9683e+00,6.5404e-02,6.4800e-02,1.0093e+00,5.7822e-03,5.7329e-03,1.0086e+00,7.4007e+00,1.0150e+00,1.0000e+00,9.7190e-01,1.0038e+00 -8.7000e-02,4.9647e+00,6.5660e-02,6.5055e-02,1.0093e+00,5.8625e-03,5.8127e-03,1.0086e+00,7.3669e+00,1.0150e+00,1.0000e+00,9.7193e-01,1.0037e+00 -8.8000e-02,4.9607e+00,6.5910e-02,6.5304e-02,1.0093e+00,5.9418e-03,5.8914e-03,1.0085e+00,7.3338e+00,1.0150e+00,1.0000e+00,9.7196e-01,1.0036e+00 -8.9000e-02,4.9566e+00,6.6153e-02,6.5546e-02,1.0093e+00,6.0199e-03,5.9691e-03,1.0085e+00,7.3012e+00,1.0150e+00,1.0000e+00,9.7199e-01,1.0035e+00 -9.0000e-02,4.9522e+00,6.6389e-02,6.5781e-02,1.0092e+00,6.0968e-03,6.0455e-03,1.0085e+00,7.2693e+00,1.0149e+00,1.0000e+00,9.7203e-01,1.0034e+00 -9.1000e-02,4.9476e+00,6.6619e-02,6.6010e-02,1.0092e+00,6.1724e-03,6.1206e-03,1.0085e+00,7.2379e+00,1.0149e+00,1.0000e+00,9.7206e-01,1.0034e+00 -9.2000e-02,4.9427e+00,6.6841e-02,6.6232e-02,1.0092e+00,6.2466e-03,6.1943e-03,1.0084e+00,7.2071e+00,1.0149e+00,1.0000e+00,9.7209e-01,1.0033e+00 -9.3000e-02,4.9375e+00,6.7057e-02,6.6447e-02,1.0092e+00,6.3194e-03,6.2666e-03,1.0084e+00,7.1768e+00,1.0149e+00,1.0000e+00,9.7213e-01,1.0032e+00 -9.4000e-02,4.9322e+00,6.7266e-02,6.6656e-02,1.0092e+00,6.3907e-03,6.3375e-03,1.0084e+00,7.1471e+00,1.0149e+00,1.0000e+00,9.7216e-01,1.0031e+00 -9.5000e-02,4.9265e+00,6.7469e-02,6.6857e-02,1.0091e+00,6.4604e-03,6.4068e-03,1.0084e+00,7.1179e+00,1.0149e+00,1.0000e+00,9.7219e-01,1.0030e+00 -9.6000e-02,4.9206e+00,6.7664e-02,6.7052e-02,1.0091e+00,6.5284e-03,6.4744e-03,1.0083e+00,7.0892e+00,1.0149e+00,1.0000e+00,9.7223e-01,1.0029e+00 -9.7000e-02,4.9145e+00,6.7853e-02,6.7241e-02,1.0091e+00,6.5947e-03,6.5403e-03,1.0083e+00,7.0610e+00,1.0149e+00,1.0000e+00,9.7226e-01,1.0029e+00 -9.8000e-02,4.9080e+00,6.8034e-02,6.7422e-02,1.0091e+00,6.6593e-03,6.6045e-03,1.0083e+00,7.0332e+00,1.0149e+00,1.0000e+00,9.7230e-01,1.0028e+00 -9.9000e-02,4.9013e+00,6.8209e-02,6.7596e-02,1.0091e+00,6.7219e-03,6.6668e-03,1.0083e+00,7.0060e+00,1.0149e+00,1.0000e+00,9.7233e-01,1.0027e+00 -1.0000e-01,4.8944e+00,6.8377e-02,6.7764e-02,1.0090e+00,6.7827e-03,6.7272e-03,1.0082e+00,6.9792e+00,1.0149e+00,1.0000e+00,9.7237e-01,1.0026e+00 -1.0100e-01,4.8872e+00,6.8537e-02,6.7925e-02,1.0090e+00,6.8415e-03,6.7857e-03,1.0082e+00,6.9528e+00,1.0149e+00,1.0000e+00,9.7240e-01,1.0026e+00 -1.0200e-01,4.8797e+00,6.8691e-02,6.8079e-02,1.0090e+00,6.8982e-03,6.8421e-03,1.0082e+00,6.9269e+00,1.0149e+00,1.0000e+00,9.7244e-01,1.0025e+00 -1.0300e-01,4.8719e+00,6.8838e-02,6.8225e-02,1.0090e+00,6.9528e-03,6.8965e-03,1.0082e+00,6.9013e+00,1.0149e+00,1.0000e+00,9.7248e-01,1.0024e+00 -1.0400e-01,4.8638e+00,6.8978e-02,6.8365e-02,1.0090e+00,7.0052e-03,6.9486e-03,1.0081e+00,6.8762e+00,1.0149e+00,1.0000e+00,9.7251e-01,1.0024e+00 -1.0500e-01,4.8555e+00,6.9110e-02,6.8498e-02,1.0089e+00,7.0553e-03,6.9986e-03,1.0081e+00,6.8515e+00,1.0149e+00,1.0000e+00,9.7255e-01,1.0023e+00 -1.0600e-01,4.8469e+00,6.9236e-02,6.8624e-02,1.0089e+00,7.1032e-03,7.0462e-03,1.0081e+00,6.8272e+00,1.0149e+00,1.0000e+00,9.7259e-01,1.0022e+00 -1.0700e-01,4.8380e+00,6.9354e-02,6.8742e-02,1.0089e+00,7.1486e-03,7.0915e-03,1.0081e+00,6.8032e+00,1.0148e+00,1.0000e+00,9.7262e-01,1.0022e+00 -1.0800e-01,4.8289e+00,6.9465e-02,6.8854e-02,1.0089e+00,7.1917e-03,7.1344e-03,1.0080e+00,6.7797e+00,1.0148e+00,1.0000e+00,9.7266e-01,1.0021e+00 -1.0900e-01,4.8194e+00,6.9569e-02,6.8958e-02,1.0089e+00,7.2322e-03,7.1748e-03,1.0080e+00,6.7564e+00,1.0148e+00,1.0000e+00,9.7270e-01,1.0020e+00 -1.1000e-01,4.8097e+00,6.9666e-02,6.9056e-02,1.0088e+00,7.2703e-03,7.2127e-03,1.0080e+00,6.7336e+00,1.0148e+00,1.0000e+00,9.7274e-01,1.0020e+00 -1.1100e-01,4.7996e+00,6.9755e-02,6.9146e-02,1.0088e+00,7.3057e-03,7.2481e-03,1.0079e+00,6.7111e+00,1.0148e+00,1.0000e+00,9.7277e-01,1.0019e+00 -1.1200e-01,4.7893e+00,6.9838e-02,6.9229e-02,1.0088e+00,7.3384e-03,7.2807e-03,1.0079e+00,6.6889e+00,1.0148e+00,1.0000e+00,9.7281e-01,1.0019e+00 -1.1300e-01,4.7787e+00,6.9913e-02,6.9305e-02,1.0088e+00,7.3684e-03,7.3107e-03,1.0079e+00,6.6670e+00,1.0148e+00,1.0000e+00,9.7285e-01,1.0018e+00 -1.1400e-01,4.7678e+00,6.9980e-02,6.9373e-02,1.0088e+00,7.3957e-03,7.3380e-03,1.0079e+00,6.6455e+00,1.0148e+00,1.0000e+00,9.7289e-01,1.0018e+00 -1.1500e-01,4.7565e+00,7.0041e-02,6.9434e-02,1.0087e+00,7.4202e-03,7.3625e-03,1.0078e+00,6.6243e+00,1.0148e+00,1.0000e+00,9.7293e-01,1.0017e+00 -1.1600e-01,4.7450e+00,7.0093e-02,6.9488e-02,1.0087e+00,7.4417e-03,7.3841e-03,1.0078e+00,6.6034e+00,1.0148e+00,1.0000e+00,9.7297e-01,1.0017e+00 -1.1700e-01,4.7332e+00,7.0139e-02,6.9535e-02,1.0087e+00,7.4604e-03,7.4028e-03,1.0078e+00,6.5828e+00,1.0148e+00,1.0000e+00,9.7300e-01,1.0016e+00 -1.1800e-01,4.7211e+00,7.0177e-02,6.9574e-02,1.0087e+00,7.4761e-03,7.4186e-03,1.0078e+00,6.5624e+00,1.0148e+00,1.0000e+00,9.7304e-01,1.0016e+00 -1.1900e-01,4.7087e+00,7.0208e-02,6.9606e-02,1.0086e+00,7.4888e-03,7.4314e-03,1.0077e+00,6.5424e+00,1.0148e+00,1.0000e+00,9.7308e-01,1.0015e+00 -1.2000e-01,4.6960e+00,7.0231e-02,6.9631e-02,1.0086e+00,7.4985e-03,7.4412e-03,1.0077e+00,6.5226e+00,1.0148e+00,1.0000e+00,9.7312e-01,1.0015e+00 -1.2100e-01,4.6830e+00,7.0247e-02,6.9648e-02,1.0086e+00,7.5051e-03,7.4480e-03,1.0077e+00,6.5032e+00,1.0148e+00,1.0000e+00,9.7316e-01,1.0014e+00 -1.2200e-01,4.6697e+00,7.0255e-02,6.9657e-02,1.0086e+00,7.5085e-03,7.4516e-03,1.0076e+00,6.4840e+00,1.0148e+00,1.0000e+00,9.7320e-01,1.0014e+00 -1.2300e-01,4.6560e+00,7.0256e-02,6.9660e-02,1.0086e+00,7.5088e-03,7.4521e-03,1.0076e+00,6.4650e+00,1.0148e+00,1.0000e+00,9.7324e-01,1.0013e+00 -1.2400e-01,4.6421e+00,7.0249e-02,6.9654e-02,1.0085e+00,7.5059e-03,7.4495e-03,1.0076e+00,6.4463e+00,1.0147e+00,1.0000e+00,9.7328e-01,1.0013e+00 -1.2500e-01,4.6278e+00,7.0235e-02,6.9642e-02,1.0085e+00,7.4998e-03,7.4436e-03,1.0075e+00,6.4279e+00,1.0147e+00,1.0000e+00,9.7332e-01,1.0013e+00 -1.2600e-01,4.6133e+00,7.0213e-02,6.9622e-02,1.0085e+00,7.4904e-03,7.4345e-03,1.0075e+00,6.4097e+00,1.0147e+00,1.0000e+00,9.7335e-01,1.0012e+00 -1.2700e-01,4.5984e+00,7.0184e-02,6.9594e-02,1.0085e+00,7.4777e-03,7.4221e-03,1.0075e+00,6.3917e+00,1.0147e+00,1.0000e+00,9.7339e-01,1.0012e+00 -1.2800e-01,4.5832e+00,7.0147e-02,6.9559e-02,1.0085e+00,7.4616e-03,7.4064e-03,1.0075e+00,6.3740e+00,1.0147e+00,1.0000e+00,9.7343e-01,1.0012e+00 -1.2900e-01,4.5677e+00,7.0102e-02,6.9516e-02,1.0084e+00,7.4423e-03,7.3874e-03,1.0074e+00,6.3565e+00,1.0147e+00,1.0000e+00,9.7347e-01,1.0011e+00 -1.3000e-01,4.5519e+00,7.0050e-02,6.9466e-02,1.0084e+00,7.4196e-03,7.3651e-03,1.0074e+00,6.3393e+00,1.0147e+00,1.0000e+00,9.7351e-01,1.0011e+00 -1.3100e-01,4.5358e+00,6.9990e-02,6.9408e-02,1.0084e+00,7.3935e-03,7.3395e-03,1.0074e+00,6.3223e+00,1.0147e+00,1.0000e+00,9.7355e-01,1.0010e+00 -1.3200e-01,4.5194e+00,6.9922e-02,6.9342e-02,1.0084e+00,7.3640e-03,7.3104e-03,1.0073e+00,6.3054e+00,1.0147e+00,1.0000e+00,9.7359e-01,1.0010e+00 -1.3300e-01,4.5027e+00,6.9847e-02,6.9269e-02,1.0083e+00,7.3312e-03,7.2780e-03,1.0073e+00,6.2888e+00,1.0147e+00,1.0000e+00,9.7363e-01,1.0010e+00 -1.3400e-01,4.4856e+00,6.9764e-02,6.9188e-02,1.0083e+00,7.2949e-03,7.2422e-03,1.0073e+00,6.2724e+00,1.0147e+00,1.0000e+00,9.7367e-01,1.0009e+00 -1.3500e-01,4.4682e+00,6.9674e-02,6.9100e-02,1.0083e+00,7.2551e-03,7.2030e-03,1.0072e+00,6.2563e+00,1.0147e+00,1.0000e+00,9.7371e-01,1.0009e+00 -1.3600e-01,4.4505e+00,6.9575e-02,6.9004e-02,1.0083e+00,7.2120e-03,7.1604e-03,1.0072e+00,6.2403e+00,1.0147e+00,1.0000e+00,9.7375e-01,1.0009e+00 -1.3700e-01,4.4325e+00,6.9469e-02,6.8900e-02,1.0083e+00,7.1654e-03,7.1144e-03,1.0072e+00,6.2245e+00,1.0147e+00,1.0000e+00,9.7379e-01,1.0009e+00 -1.3800e-01,4.4142e+00,6.9356e-02,6.8789e-02,1.0082e+00,7.1154e-03,7.0649e-03,1.0071e+00,6.2089e+00,1.0147e+00,1.0000e+00,9.7383e-01,1.0008e+00 -1.3900e-01,4.3956e+00,6.9234e-02,6.8670e-02,1.0082e+00,7.0619e-03,7.0121e-03,1.0071e+00,6.1935e+00,1.0147e+00,1.0000e+00,9.7387e-01,1.0008e+00 -1.4000e-01,4.3766e+00,6.9105e-02,6.8544e-02,1.0082e+00,7.0051e-03,6.9559e-03,1.0071e+00,6.1783e+00,1.0146e+00,1.0000e+00,9.7391e-01,1.0008e+00 -1.4100e-01,4.3573e+00,6.8968e-02,6.8409e-02,1.0082e+00,6.9448e-03,6.8963e-03,1.0070e+00,6.1632e+00,1.0146e+00,1.0000e+00,9.7395e-01,1.0007e+00 -1.4200e-01,4.3378e+00,6.8824e-02,6.8268e-02,1.0081e+00,6.8812e-03,6.8333e-03,1.0070e+00,6.1484e+00,1.0146e+00,1.0000e+00,9.7399e-01,1.0007e+00 -1.4300e-01,4.3179e+00,6.8671e-02,6.8118e-02,1.0081e+00,6.8142e-03,6.7670e-03,1.0070e+00,6.1337e+00,1.0146e+00,1.0000e+00,9.7403e-01,1.0007e+00 -1.4400e-01,4.2976e+00,6.8511e-02,6.7961e-02,1.0081e+00,6.7438e-03,6.6973e-03,1.0069e+00,6.1192e+00,1.0146e+00,1.0000e+00,9.7406e-01,1.0007e+00 -1.4500e-01,4.2771e+00,6.8343e-02,6.7796e-02,1.0081e+00,6.6701e-03,6.6244e-03,1.0069e+00,6.1048e+00,1.0146e+00,1.0000e+00,9.7410e-01,1.0006e+00 -1.4600e-01,4.2563e+00,6.8168e-02,6.7623e-02,1.0081e+00,6.5931e-03,6.5481e-03,1.0069e+00,6.0907e+00,1.0146e+00,1.0000e+00,9.7414e-01,1.0006e+00 -1.4700e-01,4.2351e+00,6.7985e-02,6.7443e-02,1.0080e+00,6.5129e-03,6.4687e-03,1.0068e+00,6.0766e+00,1.0146e+00,1.0000e+00,9.7418e-01,1.0006e+00 -1.4800e-01,4.2136e+00,6.7794e-02,6.7255e-02,1.0080e+00,6.4294e-03,6.3860e-03,1.0068e+00,6.0628e+00,1.0146e+00,1.0000e+00,9.7422e-01,1.0006e+00 -1.4900e-01,4.1919e+00,6.7595e-02,6.7060e-02,1.0080e+00,6.3427e-03,6.3001e-03,1.0068e+00,6.0491e+00,1.0146e+00,1.0000e+00,9.7426e-01,1.0005e+00 -1.5000e-01,4.1698e+00,6.7389e-02,6.6857e-02,1.0080e+00,6.2529e-03,6.2112e-03,1.0067e+00,6.0356e+00,1.0146e+00,1.0000e+00,9.7430e-01,1.0005e+00 -1.5100e-01,4.1474e+00,6.7175e-02,6.6646e-02,1.0079e+00,6.1600e-03,6.1191e-03,1.0067e+00,6.0222e+00,1.0146e+00,1.0000e+00,9.7434e-01,1.0005e+00 -1.5200e-01,4.1246e+00,6.6953e-02,6.6428e-02,1.0079e+00,6.0641e-03,6.0240e-03,1.0066e+00,6.0089e+00,1.0146e+00,1.0000e+00,9.7438e-01,1.0005e+00 -1.5300e-01,4.1016e+00,6.6724e-02,6.6202e-02,1.0079e+00,5.9651e-03,5.9260e-03,1.0066e+00,5.9959e+00,1.0146e+00,1.0000e+00,9.7442e-01,1.0005e+00 -1.5400e-01,4.0783e+00,6.6487e-02,6.5968e-02,1.0079e+00,5.8633e-03,5.8250e-03,1.0066e+00,5.9829e+00,1.0146e+00,1.0000e+00,9.7446e-01,1.0004e+00 -1.5500e-01,4.0546e+00,6.6242e-02,6.5727e-02,1.0078e+00,5.7586e-03,5.7212e-03,1.0065e+00,5.9701e+00,1.0146e+00,1.0000e+00,9.7450e-01,1.0004e+00 -1.5600e-01,4.0307e+00,6.5990e-02,6.5478e-02,1.0078e+00,5.6511e-03,5.6147e-03,1.0065e+00,5.9574e+00,1.0145e+00,1.0000e+00,9.7454e-01,1.0004e+00 -1.5700e-01,4.0064e+00,6.5730e-02,6.5222e-02,1.0078e+00,5.5409e-03,5.5054e-03,1.0065e+00,5.9449e+00,1.0145e+00,1.0000e+00,9.7457e-01,1.0004e+00 -1.5800e-01,3.9819e+00,6.5463e-02,6.4958e-02,1.0078e+00,5.4281e-03,5.3935e-03,1.0064e+00,5.9325e+00,1.0145e+00,1.0000e+00,9.7461e-01,1.0004e+00 -1.5900e-01,3.9570e+00,6.5188e-02,6.4687e-02,1.0077e+00,5.3127e-03,5.2791e-03,1.0064e+00,5.9203e+00,1.0145e+00,1.0000e+00,9.7465e-01,1.0003e+00 -1.6000e-01,3.9319e+00,6.4905e-02,6.4408e-02,1.0077e+00,5.1949e-03,5.1623e-03,1.0063e+00,5.9081e+00,1.0145e+00,1.0000e+00,9.7469e-01,1.0003e+00 -1.6100e-01,3.9064e+00,6.4615e-02,6.4121e-02,1.0077e+00,5.0748e-03,5.0431e-03,1.0063e+00,5.8962e+00,1.0145e+00,1.0000e+00,9.7473e-01,1.0003e+00 -1.6200e-01,3.8806e+00,6.4317e-02,6.3828e-02,1.0077e+00,4.9523e-03,4.9216e-03,1.0062e+00,5.8843e+00,1.0145e+00,1.0000e+00,9.7477e-01,1.0003e+00 -1.6300e-01,3.8546e+00,6.4012e-02,6.3527e-02,1.0076e+00,4.8278e-03,4.7980e-03,1.0062e+00,5.8725e+00,1.0145e+00,1.0000e+00,9.7481e-01,1.0003e+00 -1.6400e-01,3.8283e+00,6.3700e-02,6.3218e-02,1.0076e+00,4.7011e-03,4.6724e-03,1.0062e+00,5.8609e+00,1.0145e+00,1.0000e+00,9.7485e-01,1.0003e+00 -1.6500e-01,3.8016e+00,6.3380e-02,6.2902e-02,1.0076e+00,4.5726e-03,4.5448e-03,1.0061e+00,5.8494e+00,1.0145e+00,1.0000e+00,9.7488e-01,1.0002e+00 -1.6600e-01,3.7747e+00,6.3053e-02,6.2579e-02,1.0076e+00,4.4423e-03,4.4155e-03,1.0061e+00,5.8380e+00,1.0145e+00,1.0000e+00,9.7492e-01,1.0002e+00 -1.6700e-01,3.7475e+00,6.2718e-02,6.2248e-02,1.0075e+00,4.3102e-03,4.2845e-03,1.0060e+00,5.8267e+00,1.0145e+00,1.0000e+00,9.7496e-01,1.0002e+00 -1.6800e-01,3.7200e+00,6.2376e-02,6.1910e-02,1.0075e+00,4.1767e-03,4.1519e-03,1.0060e+00,5.8156e+00,1.0145e+00,1.0000e+00,9.7500e-01,1.0002e+00 -1.6900e-01,3.6922e+00,6.2027e-02,6.1565e-02,1.0075e+00,4.0417e-03,4.0179e-03,1.0059e+00,5.8045e+00,1.0145e+00,1.0000e+00,9.7504e-01,1.0002e+00 -1.7000e-01,3.6642e+00,6.1671e-02,6.1213e-02,1.0075e+00,3.9055e-03,3.8827e-03,1.0059e+00,5.7936e+00,1.0145e+00,1.0000e+00,9.7508e-01,1.0002e+00 -1.7100e-01,3.6358e+00,6.1307e-02,6.0854e-02,1.0074e+00,3.7682e-03,3.7464e-03,1.0058e+00,5.7828e+00,1.0144e+00,1.0000e+00,9.7511e-01,1.0002e+00 -1.7200e-01,3.6072e+00,6.0936e-02,6.0488e-02,1.0074e+00,3.6299e-03,3.6091e-03,1.0058e+00,5.7721e+00,1.0144e+00,1.0000e+00,9.7515e-01,1.0002e+00 -1.7300e-01,3.5783e+00,6.0558e-02,6.0114e-02,1.0074e+00,3.4909e-03,3.4710e-03,1.0057e+00,5.7614e+00,1.0144e+00,1.0000e+00,9.7519e-01,1.0001e+00 -1.7400e-01,3.5492e+00,6.0174e-02,5.9734e-02,1.0074e+00,3.3512e-03,3.3323e-03,1.0057e+00,5.7509e+00,1.0144e+00,1.0000e+00,9.7523e-01,1.0001e+00 -1.7500e-01,3.5197e+00,5.9782e-02,5.9346e-02,1.0073e+00,3.2111e-03,3.1932e-03,1.0056e+00,5.7405e+00,1.0144e+00,1.0000e+00,9.7527e-01,1.0001e+00 -1.7600e-01,3.4901e+00,5.9383e-02,5.8952e-02,1.0073e+00,3.0707e-03,3.0538e-03,1.0056e+00,5.7302e+00,1.0144e+00,1.0000e+00,9.7530e-01,1.0001e+00 -1.7700e-01,3.4601e+00,5.8977e-02,5.8550e-02,1.0073e+00,2.9303e-03,2.9143e-03,1.0055e+00,5.7200e+00,1.0144e+00,1.0000e+00,9.7534e-01,1.0001e+00 -1.7800e-01,3.4299e+00,5.8565e-02,5.8142e-02,1.0073e+00,2.7900e-03,2.7749e-03,1.0054e+00,5.7099e+00,1.0144e+00,1.0000e+00,9.7538e-01,1.0001e+00 -1.7900e-01,3.3995e+00,5.8145e-02,5.7728e-02,1.0072e+00,2.6501e-03,2.6359e-03,1.0054e+00,5.6999e+00,1.0144e+00,1.0000e+00,9.7542e-01,1.0001e+00 -1.8000e-01,3.3688e+00,5.7719e-02,5.7306e-02,1.0072e+00,2.5107e-03,2.4974e-03,1.0053e+00,5.6900e+00,1.0144e+00,1.0000e+00,9.7545e-01,1.0001e+00 -1.8100e-01,3.3378e+00,5.7287e-02,5.6878e-02,1.0072e+00,2.3720e-03,2.3596e-03,1.0053e+00,5.6802e+00,1.0144e+00,1.0000e+00,9.7549e-01,1.0001e+00 -1.8200e-01,3.3066e+00,5.6847e-02,5.6443e-02,1.0072e+00,2.2344e-03,2.2228e-03,1.0052e+00,5.6705e+00,1.0144e+00,1.0000e+00,9.7553e-01,1.0000e+00 -1.8300e-01,3.2752e+00,5.6401e-02,5.6002e-02,1.0071e+00,2.0979e-03,2.0872e-03,1.0051e+00,5.6608e+00,1.0144e+00,1.0000e+00,9.7557e-01,1.0000e+00 -1.8400e-01,3.2435e+00,5.5949e-02,5.5555e-02,1.0071e+00,1.9629e-03,1.9530e-03,1.0051e+00,5.6513e+00,1.0144e+00,1.0000e+00,9.7560e-01,1.0000e+00 -1.8500e-01,3.2116e+00,5.5490e-02,5.5101e-02,1.0071e+00,1.8297e-03,1.8206e-03,1.0050e+00,5.6418e+00,1.0144e+00,1.0000e+00,9.7564e-01,1.0000e+00 -1.8600e-01,3.1795e+00,5.5025e-02,5.4641e-02,1.0070e+00,1.6983e-03,1.6900e-03,1.0049e+00,5.6324e+00,1.0143e+00,1.0000e+00,9.7568e-01,1.0000e+00 -1.8700e-01,3.1471e+00,5.4554e-02,5.4174e-02,1.0070e+00,1.5692e-03,1.5617e-03,1.0049e+00,5.6232e+00,1.0143e+00,1.0000e+00,9.7571e-01,1.0000e+00 -1.8800e-01,3.1145e+00,5.4077e-02,5.3701e-02,1.0070e+00,1.4426e-03,1.4358e-03,1.0048e+00,5.6139e+00,1.0143e+00,1.0000e+00,9.7575e-01,9.9999e-01 -1.8900e-01,3.0817e+00,5.3593e-02,5.3223e-02,1.0070e+00,1.3188e-03,1.3126e-03,1.0047e+00,5.6048e+00,1.0143e+00,1.0000e+00,9.7579e-01,9.9998e-01 -1.9000e-01,3.0487e+00,5.3104e-02,5.2738e-02,1.0069e+00,1.1980e-03,1.1925e-03,1.0046e+00,5.5958e+00,1.0143e+00,1.0000e+00,9.7582e-01,9.9998e-01 -1.9100e-01,3.0155e+00,5.2608e-02,5.2248e-02,1.0069e+00,1.0806e-03,1.0757e-03,1.0045e+00,5.5868e+00,1.0143e+00,1.0000e+00,9.7586e-01,9.9997e-01 -1.9200e-01,2.9821e+00,5.2107e-02,5.1751e-02,1.0069e+00,9.6677e-04,9.6250e-04,1.0044e+00,5.5780e+00,1.0143e+00,1.0000e+00,9.7590e-01,9.9996e-01 -1.9300e-01,2.9485e+00,5.1600e-02,5.1249e-02,1.0068e+00,8.5690e-04,8.5320e-04,1.0043e+00,5.5692e+00,1.0143e+00,1.0000e+00,9.7593e-01,9.9995e-01 -1.9400e-01,2.9147e+00,5.1087e-02,5.0741e-02,1.0068e+00,7.5129e-04,7.4811e-04,1.0042e+00,5.5605e+00,1.0143e+00,1.0000e+00,9.7597e-01,9.9995e-01 -1.9500e-01,2.8807e+00,5.0569e-02,5.0228e-02,1.0068e+00,6.5023e-04,6.4755e-04,1.0041e+00,5.5518e+00,1.0143e+00,1.0000e+00,9.7601e-01,9.9994e-01 -1.9600e-01,2.8466e+00,5.0045e-02,4.9709e-02,1.0068e+00,5.5403e-04,5.5182e-04,1.0040e+00,5.5432e+00,1.0143e+00,1.0000e+00,9.7604e-01,9.9993e-01 -1.9700e-01,2.8122e+00,4.9516e-02,4.9185e-02,1.0067e+00,4.6300e-04,4.6121e-04,1.0039e+00,5.5348e+00,1.0143e+00,1.0000e+00,9.7608e-01,9.9993e-01 -1.9800e-01,2.7777e+00,4.8982e-02,4.8656e-02,1.0067e+00,3.7740e-04,3.7599e-04,1.0038e+00,5.5263e+00,1.0143e+00,1.0000e+00,9.7611e-01,9.9992e-01 -1.9900e-01,2.7430e+00,4.8442e-02,4.8121e-02,1.0067e+00,2.9744e-04,2.9637e-04,1.0036e+00,5.5180e+00,1.0143e+00,1.0000e+00,9.7615e-01,9.9992e-01 -2.0000e-01,2.7082e+00,4.7897e-02,4.7582e-02,1.0066e+00,2.2316e-04,2.2239e-04,1.0035e+00,5.5097e+00,1.0143e+00,1.0000e+00,9.7619e-01,9.9991e-01 -2.0100e-01,2.6732e+00,4.7348e-02,4.7037e-02,1.0066e+00,1.5417e-04,1.5367e-04,1.0033e+00,5.5015e+00,1.0142e+00,1.0000e+00,9.7622e-01,9.9990e-01 -2.0200e-01,2.6381e+00,4.6793e-02,4.6488e-02,1.0066e+00,8.8228e-05,8.7957e-05,1.0031e+00,5.4934e+00,1.0142e+00,1.0000e+00,9.7626e-01,9.9990e-01 -2.0300e-01,2.6028e+00,4.6234e-02,4.5934e-02,1.0065e+00,0.0000e+00,0.0000e+00,1.0028e+00,5.4853e+00,1.0142e+00,1.0000e+00,9.7629e-01,9.9989e-01 -2.0400e-01,2.5674e+00,4.5670e-02,4.5375e-02,1.0065e+00,0.0000e+00,0.0000e+00,1.0026e+00,5.4773e+00,1.0142e+00,1.0000e+00,9.7633e-01,9.9989e-01 -2.0500e-01,2.5318e+00,4.5102e-02,4.4812e-02,1.0065e+00,0.0000e+00,0.0000e+00,1.0022e+00,5.4694e+00,1.0142e+00,1.0000e+00,9.7636e-01,9.9988e-01 -2.0600e-01,2.4961e+00,4.4530e-02,4.4244e-02,1.0064e+00,0.0000e+00,0.0000e+00,1.0017e+00,5.4615e+00,1.0142e+00,1.0000e+00,9.7640e-01,9.9988e-01 -2.0700e-01,2.4603e+00,4.3953e-02,4.3672e-02,1.0064e+00,0.0000e+00,0.0000e+00,1.0009e+00,5.4537e+00,1.0142e+00,1.0000e+00,9.7644e-01,9.9987e-01 -2.0800e-01,2.4244e+00,4.3372e-02,4.3097e-02,1.0064e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.4460e+00,1.0142e+00,1.0000e+00,9.7647e-01,9.9987e-01 -2.0900e-01,2.3884e+00,4.2787e-02,4.2517e-02,1.0064e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.4383e+00,1.0142e+00,1.0000e+00,9.7651e-01,9.9986e-01 -2.1000e-01,2.3523e+00,4.2198e-02,4.1933e-02,1.0063e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.4307e+00,1.0142e+00,1.0000e+00,9.7654e-01,9.9986e-01 -2.1100e-01,2.3161e+00,4.1606e-02,4.1346e-02,1.0063e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.4232e+00,1.0142e+00,1.0000e+00,9.7658e-01,9.9985e-01 -2.1200e-01,2.2798e+00,4.1009e-02,4.0754e-02,1.0063e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.4157e+00,1.0142e+00,1.0000e+00,9.7661e-01,9.9985e-01 -2.1300e-01,2.2435e+00,4.0410e-02,4.0160e-02,1.0062e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.4083e+00,1.0142e+00,1.0000e+00,9.7665e-01,9.9985e-01 -2.1400e-01,2.2071e+00,3.9807e-02,3.9562e-02,1.0062e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.4009e+00,1.0142e+00,1.0000e+00,9.7668e-01,9.9984e-01 -2.1500e-01,2.1706e+00,3.9201e-02,3.8961e-02,1.0062e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3936e+00,1.0142e+00,1.0000e+00,9.7671e-01,9.9984e-01 -2.1600e-01,2.1341e+00,3.8592e-02,3.8357e-02,1.0061e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3864e+00,1.0141e+00,1.0000e+00,9.7675e-01,9.9984e-01 -2.1700e-01,2.0975e+00,3.7980e-02,3.7751e-02,1.0061e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3792e+00,1.0141e+00,1.0000e+00,9.7678e-01,9.9983e-01 -2.1800e-01,2.0608e+00,3.7366e-02,3.7141e-02,1.0061e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3720e+00,1.0141e+00,1.0000e+00,9.7682e-01,9.9983e-01 -2.1900e-01,2.0242e+00,3.6749e-02,3.6529e-02,1.0060e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3650e+00,1.0141e+00,1.0000e+00,9.7685e-01,9.9982e-01 -2.2000e-01,1.9875e+00,3.6129e-02,3.5915e-02,1.0060e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3579e+00,1.0141e+00,1.0000e+00,9.7689e-01,9.9982e-01 -2.2100e-01,1.9508e+00,3.5508e-02,3.5298e-02,1.0059e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3510e+00,1.0141e+00,1.0000e+00,9.7692e-01,9.9982e-01 -2.2200e-01,1.9142e+00,3.4884e-02,3.4679e-02,1.0059e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3441e+00,1.0141e+00,1.0000e+00,9.7695e-01,9.9982e-01 -2.2300e-01,1.8775e+00,3.4259e-02,3.4059e-02,1.0059e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3372e+00,1.0141e+00,1.0000e+00,9.7699e-01,9.9981e-01 -2.2400e-01,1.8408e+00,3.3632e-02,3.3437e-02,1.0058e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3304e+00,1.0141e+00,1.0000e+00,9.7702e-01,9.9981e-01 -2.2500e-01,1.8042e+00,3.3003e-02,3.2813e-02,1.0058e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3236e+00,1.0141e+00,1.0000e+00,9.7706e-01,9.9981e-01 -2.2600e-01,1.7675e+00,3.2373e-02,3.2188e-02,1.0058e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3169e+00,1.0141e+00,1.0000e+00,9.7709e-01,9.9980e-01 -2.2700e-01,1.7310e+00,3.1742e-02,3.1562e-02,1.0057e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3103e+00,1.0141e+00,1.0000e+00,9.7712e-01,9.9980e-01 -2.2800e-01,1.6944e+00,3.1110e-02,3.0934e-02,1.0057e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.3037e+00,1.0141e+00,1.0000e+00,9.7716e-01,9.9980e-01 -2.2900e-01,1.6580e+00,3.0478e-02,3.0306e-02,1.0056e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2971e+00,1.0141e+00,1.0000e+00,9.7719e-01,9.9980e-01 -2.3000e-01,1.6216e+00,2.9844e-02,2.9678e-02,1.0056e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2906e+00,1.0141e+00,1.0000e+00,9.7722e-01,9.9979e-01 -2.3100e-01,1.5852e+00,2.9211e-02,2.9049e-02,1.0056e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2842e+00,1.0140e+00,1.0000e+00,9.7726e-01,9.9979e-01 -2.3200e-01,1.5490e+00,2.8577e-02,2.8420e-02,1.0055e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2778e+00,1.0140e+00,1.0000e+00,9.7729e-01,9.9979e-01 -2.3300e-01,1.5129e+00,2.7943e-02,2.7791e-02,1.0055e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2714e+00,1.0140e+00,1.0000e+00,9.7732e-01,9.9979e-01 -2.3400e-01,1.4768e+00,2.7310e-02,2.7162e-02,1.0055e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2651e+00,1.0140e+00,1.0000e+00,9.7735e-01,9.9979e-01 -2.3500e-01,1.4409e+00,2.6677e-02,2.6533e-02,1.0054e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2588e+00,1.0140e+00,1.0000e+00,9.7739e-01,9.9978e-01 -2.3600e-01,1.4052e+00,2.6045e-02,2.5906e-02,1.0054e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2526e+00,1.0140e+00,1.0000e+00,9.7742e-01,9.9978e-01 -2.3700e-01,1.3695e+00,2.5413e-02,2.5279e-02,1.0053e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2465e+00,1.0140e+00,1.0000e+00,9.7745e-01,9.9978e-01 -2.3800e-01,1.3340e+00,2.4783e-02,2.4653e-02,1.0053e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2403e+00,1.0140e+00,1.0000e+00,9.7749e-01,9.9978e-01 -2.3900e-01,1.2987e+00,2.4154e-02,2.4028e-02,1.0052e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2342e+00,1.0140e+00,1.0000e+00,9.7752e-01,9.9978e-01 -2.4000e-01,1.2635e+00,2.3527e-02,2.3405e-02,1.0052e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2282e+00,1.0140e+00,1.0000e+00,9.7755e-01,9.9978e-01 -2.4100e-01,1.2286e+00,2.2902e-02,2.2784e-02,1.0052e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2222e+00,1.0140e+00,1.0000e+00,9.7758e-01,9.9977e-01 -2.4200e-01,1.1938e+00,2.2278e-02,2.2165e-02,1.0051e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2162e+00,1.0140e+00,1.0000e+00,9.7762e-01,9.9977e-01 -2.4300e-01,1.1592e+00,2.1657e-02,2.1548e-02,1.0051e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2103e+00,1.0140e+00,1.0000e+00,9.7765e-01,9.9977e-01 -2.4400e-01,1.1249e+00,2.1039e-02,2.0934e-02,1.0050e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.2045e+00,1.0140e+00,1.0000e+00,9.7768e-01,9.9977e-01 -2.4500e-01,1.0908e+00,2.0423e-02,2.0322e-02,1.0050e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1986e+00,1.0139e+00,1.0000e+00,9.7771e-01,9.9977e-01 -2.4600e-01,1.0569e+00,1.9811e-02,1.9713e-02,1.0049e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1928e+00,1.0139e+00,1.0000e+00,9.7774e-01,9.9977e-01 -2.4700e-01,1.0233e+00,1.9201e-02,1.9108e-02,1.0049e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1871e+00,1.0139e+00,1.0000e+00,9.7778e-01,9.9977e-01 -2.4800e-01,9.8995e-01,1.8596e-02,1.8506e-02,1.0048e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1814e+00,1.0139e+00,1.0000e+00,9.7781e-01,9.9977e-01 -2.4900e-01,9.5688e-01,1.7994e-02,1.7908e-02,1.0048e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1757e+00,1.0139e+00,1.0000e+00,9.7784e-01,9.9976e-01 -2.5000e-01,9.2411e-01,1.7396e-02,1.7314e-02,1.0047e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1701e+00,1.0139e+00,1.0000e+00,9.7787e-01,9.9976e-01 -2.5100e-01,8.9165e-01,1.6803e-02,1.6724e-02,1.0047e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1645e+00,1.0139e+00,1.0000e+00,9.7790e-01,9.9976e-01 -2.5200e-01,8.5950e-01,1.6214e-02,1.6139e-02,1.0046e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1590e+00,1.0139e+00,1.0000e+00,9.7793e-01,9.9976e-01 -2.5300e-01,8.2769e-01,1.5630e-02,1.5559e-02,1.0046e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1534e+00,1.0139e+00,1.0000e+00,9.7796e-01,9.9976e-01 -2.5400e-01,7.9623e-01,1.5052e-02,1.4984e-02,1.0045e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1480e+00,1.0139e+00,1.0000e+00,9.7800e-01,9.9976e-01 -2.5500e-01,7.6512e-01,1.4479e-02,1.4414e-02,1.0045e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1425e+00,1.0139e+00,1.0000e+00,9.7803e-01,9.9976e-01 -2.5600e-01,7.3439e-01,1.3911e-02,1.3850e-02,1.0044e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1371e+00,1.0139e+00,1.0000e+00,9.7806e-01,9.9976e-01 -2.5700e-01,7.0405e-01,1.3350e-02,1.3292e-02,1.0044e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1318e+00,1.0139e+00,1.0000e+00,9.7809e-01,9.9976e-01 -2.5800e-01,6.7411e-01,1.2795e-02,1.2741e-02,1.0043e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1264e+00,1.0139e+00,1.0000e+00,9.7812e-01,9.9976e-01 -2.5900e-01,6.4459e-01,1.2247e-02,1.2196e-02,1.0042e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1211e+00,1.0139e+00,1.0000e+00,9.7815e-01,9.9976e-01 -2.6000e-01,6.1550e-01,1.1706e-02,1.1658e-02,1.0042e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1159e+00,1.0138e+00,1.0000e+00,9.7818e-01,9.9976e-01 -2.6100e-01,5.8686e-01,1.1173e-02,1.1127e-02,1.0041e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1106e+00,1.0138e+00,1.0000e+00,9.7821e-01,9.9976e-01 -2.6200e-01,5.5867e-01,1.0647e-02,1.0604e-02,1.0041e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1054e+00,1.0138e+00,1.0000e+00,9.7824e-01,9.9976e-01 -2.6300e-01,5.3096e-01,1.0129e-02,1.0088e-02,1.0040e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.1003e+00,1.0138e+00,1.0000e+00,9.7827e-01,9.9976e-01 -2.6400e-01,5.0375e-01,9.6190e-03,9.5813e-03,1.0039e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0952e+00,1.0138e+00,1.0000e+00,9.7830e-01,9.9975e-01 -2.6500e-01,4.7704e-01,9.1179e-03,9.0828e-03,1.0039e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0901e+00,1.0138e+00,1.0000e+00,9.7833e-01,9.9975e-01 -2.6600e-01,4.5085e-01,8.6257e-03,8.5931e-03,1.0038e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0850e+00,1.0138e+00,1.0000e+00,9.7836e-01,9.9975e-01 -2.6700e-01,4.2520e-01,8.1429e-03,8.1126e-03,1.0037e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0800e+00,1.0138e+00,1.0000e+00,9.7839e-01,9.9975e-01 -2.6800e-01,4.0011e-01,7.6697e-03,7.6417e-03,1.0037e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0750e+00,1.0138e+00,1.0000e+00,9.7842e-01,9.9975e-01 -2.6900e-01,3.7559e-01,7.2065e-03,7.1807e-03,1.0036e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0701e+00,1.0138e+00,1.0000e+00,9.7846e-01,9.9975e-01 -2.7000e-01,3.5165e-01,6.7536e-03,6.7300e-03,1.0035e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0651e+00,1.0138e+00,1.0000e+00,9.7848e-01,9.9975e-01 -2.7100e-01,3.2832e-01,6.3114e-03,6.2898e-03,1.0034e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0602e+00,1.0138e+00,1.0000e+00,9.7851e-01,9.9975e-01 -2.7200e-01,3.0560e-01,5.8803e-03,5.8606e-03,1.0034e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0554e+00,1.0138e+00,1.0000e+00,9.7854e-01,9.9975e-01 -2.7300e-01,2.8353e-01,5.4606e-03,5.4428e-03,1.0033e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0505e+00,1.0138e+00,1.0000e+00,9.7857e-01,9.9975e-01 -2.7400e-01,2.6211e-01,5.0527e-03,5.0367e-03,1.0032e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0457e+00,1.0137e+00,1.0000e+00,9.7860e-01,9.9975e-01 -2.7500e-01,2.4135e-01,4.6570e-03,4.6426e-03,1.0031e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0410e+00,1.0137e+00,1.0000e+00,9.7863e-01,9.9975e-01 -2.7600e-01,2.2129e-01,4.2738e-03,4.2610e-03,1.0030e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0362e+00,1.0137e+00,1.0000e+00,9.7866e-01,9.9975e-01 -2.7700e-01,2.0193e-01,3.9035e-03,3.8922e-03,1.0029e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0315e+00,1.0137e+00,1.0000e+00,9.7869e-01,9.9975e-01 -2.7800e-01,1.8330e-01,3.5465e-03,3.5366e-03,1.0028e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0268e+00,1.0137e+00,1.0000e+00,9.7872e-01,9.9975e-01 -2.7900e-01,1.6541e-01,3.2033e-03,3.1946e-03,1.0027e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0222e+00,1.0137e+00,1.0000e+00,9.7875e-01,9.9975e-01 -2.8000e-01,1.4828e-01,2.8741e-03,2.8666e-03,1.0026e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0175e+00,1.0137e+00,1.0000e+00,9.7878e-01,9.9975e-01 -2.8100e-01,1.3192e-01,2.5593e-03,2.5530e-03,1.0025e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0129e+00,1.0137e+00,1.0000e+00,9.7881e-01,9.9975e-01 -2.8200e-01,1.1636e-01,2.2595e-03,2.2541e-03,1.0024e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0084e+00,1.0137e+00,1.0000e+00,9.7884e-01,9.9975e-01 -2.8300e-01,1.0161e-01,1.9748e-03,1.9704e-03,1.0022e+00,0.0000e+00,0.0000e+00,1.0000e+00,5.0038e+00,1.0137e+00,1.0000e+00,9.7887e-01,9.9975e-01 -2.8400e-01,8.7693e-02,1.7058e-03,1.7022e-03,1.0021e+00,0.0000e+00,0.0000e+00,1.0000e+00,4.9993e+00,1.0137e+00,1.0000e+00,9.7890e-01,9.9975e-01 -2.8500e-01,7.4619e-02,1.4528e-03,1.4499e-03,1.0020e+00,0.0000e+00,0.0000e+00,1.0000e+00,4.9948e+00,1.0137e+00,1.0000e+00,9.7892e-01,9.9975e-01 -2.8600e-01,6.2405e-02,1.2160e-03,1.2138e-03,1.0018e+00,0.0000e+00,0.0000e+00,1.0000e+00,4.9904e+00,1.0137e+00,1.0000e+00,9.7895e-01,9.9975e-01 -2.8700e-01,5.1064e-02,9.9590e-04,9.9425e-04,1.0017e+00,0.0000e+00,0.0000e+00,1.0000e+00,4.9859e+00,1.0137e+00,1.0000e+00,9.7898e-01,9.9975e-01 -2.8800e-01,4.0602e-02,7.9253e-04,7.9136e-04,1.0015e+00,0.0000e+00,0.0000e+00,1.0000e+00,4.9815e+00,1.0136e+00,1.0000e+00,9.7901e-01,9.9975e-01 -2.8900e-01,3.1012e-02,6.0587e-04,6.0509e-04,1.0013e+00,0.0000e+00,0.0000e+00,1.0000e+00,4.9771e+00,1.0136e+00,1.0000e+00,9.7904e-01,9.9975e-01 -2.9000e-01,2.2259e-02,4.3523e-04,4.3477e-04,1.0011e+00,0.0000e+00,0.0000e+00,1.0000e+00,4.9728e+00,1.0136e+00,1.0000e+00,9.7907e-01,9.9975e-01 -2.9100e-01,1.4207e-02,2.7804e-04,2.7781e-04,1.0008e+00,0.0000e+00,0.0000e+00,1.0000e+00,4.9685e+00,1.0136e+00,1.0000e+00,9.7910e-01,9.9975e-01 -2.9200e-01,6.1264e-03,1.1999e-04,1.1993e-04,1.0005e+00,0.0000e+00,0.0000e+00,1.0000e+00,4.9642e+00,1.0136e+00,1.0000e+00,9.7912e-01,9.9975e-01 -2.9300e-01,0.0000e+00,0.0000e+00,0.0000e+00,1.0002e+00,0.0000e+00,0.0000e+00,1.0000e+00,4.9599e+00,1.0136e+00,1.0000e+00,9.7915e-01,9.9976e-01 -2.9400e-01,0.0000e+00,0.0000e+00,0.0000e+00,9.9978e-01,0.0000e+00,0.0000e+00,1.0000e+00,4.9556e+00,1.0136e+00,1.0000e+00,9.7918e-01,9.9976e-01 -2.9500e-01,0.0000e+00,0.0000e+00,0.0000e+00,9.9923e-01,0.0000e+00,0.0000e+00,1.0000e+00,4.9514e+00,1.0136e+00,1.0000e+00,9.7921e-01,9.9976e-01 -2.9600e-01,0.0000e+00,0.0000e+00,0.0000e+00,9.9841e-01,0.0000e+00,0.0000e+00,1.0000e+00,4.9472e+00,1.0136e+00,1.0000e+00,9.7924e-01,9.9976e-01 -2.9700e-01,0.0000e+00,0.0000e+00,0.0000e+00,9.9673e-01,0.0000e+00,0.0000e+00,1.0000e+00,4.9430e+00,1.0136e+00,1.0000e+00,9.7926e-01,9.9976e-01 diff --git a/tests/physics/test_multi_branch.py b/tests/physics/test_multi_branch.py index efc9eeb..9618648 100644 --- a/tests/physics/test_multi_branch.py +++ b/tests/physics/test_multi_branch.py @@ -464,11 +464,11 @@ class TestMultiBranchCSVExport: df = pd.read_csv(csv_path, comment="#") # Should have total spectrum column - assert "spectrum" in df.columns + assert "total_spectrum" in df.columns # Should have per-branch spectrum columns (intensity/endpoint/transition are in the header) - assert "branch_1_spectrum" in df.columns - assert "branch_2_spectrum" in df.columns + assert "branch 1" in df.columns + assert "branch 2" in df.columns # Should NOT have repeated per-branch metadata columns assert "branch_1_intensity" not in df.columns @@ -502,7 +502,7 @@ class TestMultiBranchCSVExport: df = pd.read_csv(csv_path, comment="#") - assert "spectrum" in df.columns + assert "total_spectrum" in df.columns assert "Fermi" in df.columns diff --git a/tests/quality/test_logging_utils.py b/tests/quality/test_logging_utils.py index 3482ec8..e84038e 100644 --- a/tests/quality/test_logging_utils.py +++ b/tests/quality/test_logging_utils.py @@ -154,7 +154,7 @@ class TestCSVMetadataHeader: lines = f.readlines() assert any("# beta-spectrum v" in line for line in lines) - assert any("0.3.0" in line for line in lines) + assert any("0.4.0" in line for line in lines) os.unlink(filename) def test_header_contains_timestamp(self): @@ -236,7 +236,9 @@ class TestCSVMetadataHeader: with open(filename, "r") as f: content = f.read() - assert "# git_commit:" in content + # Check for version with build hash instead of git_commit line + assert "# beta-spectrum v" in content + assert "(build." in content os.unlink(filename) def test_header_contains_corrections(self): @@ -264,8 +266,12 @@ class TestCSVMetadataHeader: content = f.read() assert "# corrections:" in content - assert "phase_space" in content - assert "fermi" in content + # phase_space and fermi are basic components and should not be listed + assert "screening" in content + assert "finite_size" in content + assert "charge_dist" in content + assert "radiative" in content + assert "exchange" in content os.unlink(filename)