-
diagGaussian_*figure.png : .png figure showing the (x,y) space before (top left), after (top right), the structure and after a 0.51-m drift
-
diagGaussian_*_init.npz : structured file with initial macroparticle distribution
-
diagGaussian_*_final.npz : structured file with macroparticle distribution at the exit of the structure
-
diagGaussian_*_scrn.npz : structured file with macroparticle distribution at the exit of the structure
import numpy as np
fileiN=np.load('diagGaussian_yoffum_00.0_tiltdeg_00.0_20210904_142002_scrn.npz')
# these are the particle positions (xs, ys, zs)
xs=fileiN['xs']
ys=fileiN['ys']
zs=fileiN['zs']
# these are the particle velocities (uxs, uys, uzs) following WARP convension u_i=c gamma*beta_i for i in[x,y,z]
uxs=fileiN['uxs']
uys=fileiN['uys']
uzs=fileiN['uzs']
An example is provided in the readexple file which can be run with the argument:
python readexple.py diagGaussian_yoffum_00.0_tiltdeg_02.0_20210904_143344require
- diagGaussian_noUpperSlab: Gaussian beam above a single slab
- diagGaussian: Gaussian beam injected in the structure
For all cases the beam was located at (x,y)=(0,0) [yoffum=0] and only the tilt was varied.
- Ls = 15e-2 # length of the structure
- W = 50e-3 # along x
- a = 0.914e-3/2. # half vacuum gap
- b = a+5e-3 # half conductor gap
beam_Np = 250000 x0 = 0.0 # initial beam position y0 = yoff # initial beam position z0 = z0Beam # initial beam position beam_weight= beam_Q / (beam_Np*echarge)
- beam_emit_x = 200e-6 # in m.rad
- beam_emit_y = 2e-6 # in m.rad
- beam_xrms = 2.0e-3 m
- beam_yrms = 0.10e-3 m
- beam_zrms = 640e-6 m
- beam_qbunch = 2e-9 C
_yoffum_00.0_tiltdeg_00.0_YYYYMMDD_HHMMSSfigure.png