gwm17/spspy
1
0
Fork 0
mirror of https://github.com/gwm17/spspy.git synced 2024-11-22 18:18:52 -05:00
Code Issues Projects Releases Wiki Activity

Remove spurious debug prints in SPSTarget

Browse Source
This commit is contained in:
Gordon McCann 2023-01-25 09:18:44 -05:00
parent 1da81cb7ad
commit 4f622097cd
1 changed files with 0 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
spspy/SPSTarget.py
View File

@ -114,7 +114,6 @@ class SPSTarget:
if angle == pi*0.5: if angle == pi*0.5:
return e_initial return e_initial
print(f"z{zp} a{ap}")
projectile = catima.Projectile(ap, zp) projectile = catima.Projectile(ap, zp)
e_current = e_initial/ap e_current = e_initial/ap
for (idx, layer) in enumerate(self.layer_details): for (idx, layer) in enumerate(self.layer_details):
@ -123,12 +122,10 @@ class SPSTarget:
if idx == rxn_layer: if idx == rxn_layer:
material.thickness(self.layer_details[idx].thickness * self.UG2G / (2.0 * abs(cos(angle)))) material.thickness(self.layer_details[idx].thickness * self.UG2G / (2.0 * abs(cos(angle))))
e_current -= get_energyloss(projectile, material) e_current -= get_energyloss(projectile, material)
print("e_current: ", e_current*ap)
return e_initial - e_current*ap return e_initial - e_current*ap
else: else:
material.thickness(self.layer_details[idx].thickness * self.UG2G / abs(cos(angle))) material.thickness(self.layer_details[idx].thickness * self.UG2G / abs(cos(angle)))
e_current -= get_energyloss(projectile, material) e_current -= get_energyloss(projectile, material)
print("e_current2: ", e_current*ap)
return e_initial - e_current*ap return e_initial - e_current*ap