From 3ead6b4aab85c0bd2f55a834998d607e989c3346 Mon Sep 17 00:00:00 2001
From: Gamelin Alexis <gamelin@synchrotron-soleil.fr>
Date: Tue, 6 Apr 2021 10:56:28 +0200
Subject: [PATCH] Fix phasor_evol
Correct time decay !
Now ok versus for loop (decay then beam loading)
---
tracking/rf.py | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/tracking/rf.py b/tracking/rf.py
index a12098a..b956380 100644
--- a/tracking/rf.py
+++ b/tracking/rf.py
@@ -66,7 +66,7 @@ class CavityResonator():
m : int or float
Harmonic number of the cavity.
Rs : float
- Shunt impedance of the cavity in [Ohm].
+ Shunt impedance of the cavity in [Ohm], defined as 0.5*Vc*Vc/Pc.
Q : float
Quality factor of the cavity.
QL : float
@@ -226,7 +226,6 @@ class CavityResonator():
if index != self.bunch_index:
self.phasor_evol(profile, bin_length, charge_per_mp, ref_frame=ref_frame)
- print('here')
else:
# modify beam phasor
for i, center0 in enumerate(center):
@@ -334,7 +333,7 @@ class CavityResonator():
# Phasor evolution due to induced voltage by marco-particles
k = np.arange(0, n_bin)
var = np.exp( (-1/self.filling_time + 1j*delta) *
- (n_bin-k) * deltaT )
+ (n_bin-k) * bin_length )
sum_tot = np.sum(profile * var)
sum_val = -2 * sum_tot * charge_per_mp * self.loss_factor
self.beam_phasor += sum_val
--
GitLab