rtang/PtolemyGUI
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

small update

Browse Source
This commit is contained in:
Ryan Tang 2025-02-25 19:04:28 -05:00
parent b1433f3d98
commit 6a125d4a64
2 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
Raphael/dwba_zr.py
View File

@ -438,7 +438,7 @@ class DWBA_ZR:
[dummy, index1, indexL2, index2] = self.ConvertLJ2RadialIndex(L1, J1, L2, J2) [dummy, index1, indexL2, index2] = self.ConvertLJ2RadialIndex(L1, J1, L2, J2)
ri = self.radialInt[int(L1)][index1][indexL2][index2] 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 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): for m in np.arange(-self.j + mb - ma, self.j + mb -ma + 1, 1):
haha = self.Beta(m, ma, mb) haha = self.Beta(m, ma, mb)
xsec += np.abs(haha)**2 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 return xsec * self.xsecScalingfactor * 10 # factor 10 for fm^2 = 10 mb

2
dwuck4/inFileCreatorDW.py
View File

@ -65,7 +65,7 @@ BindingEnergy = reactionData.BindingEnergy
#=================== outfile name #=================== outfile name
fileOutName = str(sym_A) + str(A_A) + "_" + str(nu_a) + str(nu_b) + "_" \ 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 #=================== find the maximum L for partial wave