title, 'test of rf-multipole';
option, debug, -echo;

bunchcharge := 2.e10;

twss_ebc1 : beta0, 
            betx=3.3999, bety=3.7645,
            alfx=-1.0879, alfy= 1.1907;

beam_ebc1 : beam, energy=5, npart=bunchcharge, 
            sige=1.5e-3, sigt=6e-03, 
            ex=8e-6 * emass/5, ey=20e-9 * emass/5;

!!!
!!! it creates a sequence with an RF-Multipole and multipole that cancel each other
!!!

k:=sqrt(pi);
rfm1 : rfmultipole, l=0, freq=12, knl:={0, sqrt(1), -sqrt(3)}, ksl:={0, sqrt(5), -sqrt(7), sqrt(11), -sqrt(13) }, pnl:={ 0, 0, 0,0,0,0,0,0,0 }, psl:={0,0,0,0,0,0,0,0};
m1 : multipole, l=0.0, knl:={0, -sqrt(1), sqrt(3)}, ksl:={0, -sqrt(5), sqrt(7), -sqrt(11), sqrt(13) };
q1 : quadrupole, l=1.0, k1=0.1;

test_line : line=( m1, rfm1, q1 );

!!! 
!!! TEST 1: prints some of the RF-Multipole attributes
!!!

use, sequence=test_line;

track,onepass,file=test1_track;
start,x=1,t=0;
observe,place=test_line$END;
run,turns=10;
endtrack;