From 6a125d4a64fc28c223eb77246fefde7a1f207d22 Mon Sep 17 00:00:00 2001 From: "Ryan@SOLARIS_testStation" Date: Tue, 25 Feb 2025 19:04:28 -0500 Subject: [PATCH] small update --- Raphael/dwba_zr.py | 3 ++- dwuck4/inFileCreatorDW.py | 2 +- 2 files changed, 3 insertions(+), 2 deletions(-) diff --git a/Raphael/dwba_zr.py b/Raphael/dwba_zr.py index 1c21146..8650f51 100755 --- a/Raphael/dwba_zr.py +++ b/Raphael/dwba_zr.py @@ -438,7 +438,7 @@ class DWBA_ZR: [dummy, index1, indexL2, index2] = self.ConvertLJ2RadialIndex(L1, J1, L2, J2) ri = self.radialInt[int(L1)][index1][indexL2][index2] - # print(f"{L1:2d}, {J1:4.1f}({index1:d}), {L2:2d}({indexL2:d}), {J2:4.1f}({index2:d}), {gg:10.6f}, {ri *self.ffactor :.10f}, {lp:10.6f}") + # print(f"{L1:2d}, {J1:4.1f}({index1:d}), {L2:2d}({indexL2:d}), {J2:4.1f}({index2:d}), {gg:10.6f}, {ri :.10f}, {lp:10.6f}") result += gg * lp * ri @@ -459,6 +459,7 @@ class DWBA_ZR: for m in np.arange(-self.j + mb - ma, self.j + mb -ma + 1, 1): haha = self.Beta(m, ma, mb) xsec += np.abs(haha)**2 + # print(f"{ma:4.1f}, {mb:4.1f}, {m:4.0f}, {haha:.6f}, {np.abs(haha)**2:.6e}, {xsec:.6e}") return xsec * self.xsecScalingfactor * 10 # factor 10 for fm^2 = 10 mb diff --git a/dwuck4/inFileCreatorDW.py b/dwuck4/inFileCreatorDW.py index c8eda80..69635e6 100755 --- a/dwuck4/inFileCreatorDW.py +++ b/dwuck4/inFileCreatorDW.py @@ -65,7 +65,7 @@ BindingEnergy = reactionData.BindingEnergy #=================== outfile name fileOutName = str(sym_A) + str(A_A) + "_" + str(nu_a) + str(nu_b) + "_" \ - + str(node) + str(l_sym) + str(int(2*j)) + "_" + str(Ex) + "_" + str(ELab) + ".in" + + str(node) + str(l_sym) + str(int(2*j)) + "_" + str(Ex) + "_" + str(ELab) + "_" + JB_pi +".in" #=================== find the maximum L for partial wave