Fix phasor_evol

Correct time decay !
Now ok versus for loop (decay then beam loading)

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