From 0af018f0d0425a9ad22c7decb4406012de225a7d Mon Sep 17 00:00:00 2001 From: "Ryan@SOLARIS_testStation" Date: Mon, 10 Mar 2025 16:17:49 -0400 Subject: [PATCH] added convertXsec.py for DWBA.Xsec.txt --- Cleopatra/IAEANuclearData.py | 2 +- PyGUIQt6/convertXsec.py | 25 +++++++++++++++++++++++++ dwuck4/extractData.py | 1 - 3 files changed, 26 insertions(+), 2 deletions(-) create mode 100755 PyGUIQt6/convertXsec.py diff --git a/Cleopatra/IAEANuclearData.py b/Cleopatra/IAEANuclearData.py index 31d27f1..de94253 100644 --- a/Cleopatra/IAEANuclearData.py +++ b/Cleopatra/IAEANuclearData.py @@ -196,7 +196,7 @@ class IsotopeClass: [A, Z] = self.GetAZ(ASym) print("
========================= ", ASym) print("
A : %d, Z : %d, N : %d" % (A, Z, A-Z)) - print("
Jpi : ", self.GetJpi(ASym)) + print("
Jpi : ", self.GetJpi_ASym(ASym)) print("
half-live : %.2f sec" % (self.GetHalfLife(ASym))) print("
Mass : %9.2f MeV" % (self.GetMassFromSym(ASym) )) print("
Binding : %9.2f MeV/A" % (self.GetBindingPerA(ASym))) diff --git a/PyGUIQt6/convertXsec.py b/PyGUIQt6/convertXsec.py new file mode 100755 index 0000000..1fe563b --- /dev/null +++ b/PyGUIQt6/convertXsec.py @@ -0,0 +1,25 @@ +#!/usr/bin/env python3 + +import sys + +filename = sys.argv[1] + +if len(sys.argv) < 2: + print("Error: Not enough arguments provided.") + print("Usage: ./{sys.argv[0]} filename") + sys.exit(1) + + +import numpy as np + + +# Load data while skipping the first two lines +data = np.loadtxt(filename, skiprows=2) + +# Extract columns +angles, values = data[:, 0], data[:, 1] + +for a, b in zip(angles, values): + print(f"{{{a:5.1f}, {b:7.3f}}},", end="") + +print("\n") \ No newline at end of file diff --git a/dwuck4/extractData.py b/dwuck4/extractData.py index 55a01ba..8094b05 100755 --- a/dwuck4/extractData.py +++ b/dwuck4/extractData.py @@ -11,7 +11,6 @@ if len(sys.argv) < 2: ##################################################### import numpy as np import re -import numpy as np import matplotlib.pyplot as plt Lmax = None