Here, we show a more complex example which may help after reading the rest of this manual. It is for a - collider. First define two electron beams and two lasers. The electron beams are converted into photons by non-linear Compton scattering with lasers. Then, the back-scattered photons and the spent electrons all come to the collision point.
! NLC Gamma-Gamma Collider ! Takes 2000 sec on HP workstation at KEK HEADER 'NLC G-G; G G;'; SET photon=1, electron=2, positron=3, mm=1D-3, micron=1D-6, nm=1D-9, mu0=4*Pi*1D-7, ee=250D9, gamma=ee/Emass, an=0.65D10, sigz=0.1*mm, betax=0.5*mm, betay=0.5*mm, emitx=5D-6/gamma, emity=8D-8/gamma, sigx=Sqrt(emitx*betax), sigy=Sqrt(emity*betay), off=sigy, laserwl=1.053*micron, pulseE=1.0, lambar=laserwl/(2*Pi), omegal=Hbarc/lambar, rlx=0.1*mm, rly=0.1*mm, sigt=0.23*mm, powerd=pulseE*Cvel/[Pi*lambar*sigt*Sqrt(2*Pi*rlx*rly)], xisq=powerd*mu0*Cvel*(lambar/Emass)^2, xi=Sqrt(xisq), eta=omegal*ee/Emass^2, lambda=4*eta, angle=0.0, dcp=5*mm ; ! dcp=CP-to-IP distance SET MsgLevel=1; PRINT 90*180*an^2/(4*Pi*sigx*sigy)/1D4, FORMAT=(' Lgeom=',1PD11.4); BEAM RIGHT, KIND=electron, NP=5000, AN=an, E0=ee, TXYS=(0,0,off/2,0), BETA=(betax,betay), EMIT=(emitx,emity), SIGT=sigz, SPIN=(0,0,-1); BEAM LEFT, KIND=electron, NP=5000, AN=an, E0=ee, TXYS=(0,0,-off/2,0), BETA=(betax,betay), EMIT=(emitx,emity), SIGT=sigz, SPIN=(0,0,1); LASER LEFT, WAVEL=laserwl, POWERD=powerd, TXYS=(-dcp,0,off/2,-dcp), E3=(0,-Sin(angle),-Cos(angle)), E1=(1,0,0), RAYLEIGH=(rlx,rly), SIGT=sigt, STOKES=(0,1,0) ; LASER RIGHT, WAVEL=laserwl, POWERD=powerd, TXYS=(-dcp,0,-off/2,dcp), E3=(0,-Sin(angle),Cos(angle)), E1=(1,0,0), RAYLEIGH=(rlx,rly), SIGT=sigt, STOKES=(0,1,0) ; LASERQED COMPTON, NPH=5, XIMAX=1.1*xi, LAMBDAMAX=1.1*lambda ; LASERQED BREITW, NPH=5, XIMAX=1.1*xi, ETAMAX=1.1*eta ; SET MsgLevel=0; FLAG OFF ECHO; SET Smesh=sigz/3; SET it=0; PRINT CPUTIME; PUSH Time=(-3*sigt-2.5*sigz-dcp,3*sigt+2.5*sigz-dcp,1000); IF Mod(it,100)=0; PRINT it, FORMAT=(F6.0,'-th time step'); PRINT STAT, SHORT; ENDIF; SET it=it+1; ENDPUSH; PRINT CPUTIME; DRIFT S=0; PRINT STAT; PLOT HIST, RIGHT, KIND=electron, H=En/1D9, HSCALE=(0,250.5,50), TITLE='Right-Going Electron Energy Spectrum after CP;', HTITLE='E0e1 (GeV); X X ;'; PLOT HIST, LEFT, KIND=electron, H=En/1D9, HSCALE=(0,250.5,50), TITLE='Left-Going Electron Energy Spectrum after CP;', HTITLE='E0e1 (GeV); X X ;'; PLOT HIST, KIND=photon, H=En/1D9, HSCALE=(0,250.5,50), TITLE='All Photon Energy Spectrum after CP;', HTITLE='E0G1 (GeV); XGX ;' ; PLOT HIST, RIGHT, KIND=photon, GEN=2, H=En/1D9, HSCALE=(0,250.5,50), PLOT HIST, RIGHT, KIND=photon, GEN=2, H=En/1D9, HSCALE=(0,250.5,50), TITLE='Right-Going Primary Photon Energy Spectrum after CP;', HTITLE='E0G1 (GeV); XGX ;' ; PLOT HIST, LEFT, KIND=photon, GEN=2, H=En/1D9, HSCALE=(0,250.5,50), TITLE='Left-Going Primary Photon Energy Spectrum after CP;', HTITLE='E0G1 (GeV); XGX ;' ; PLOT SCAT, RIGHT, KIND=photon, GEN=2, H=X/nm, V=Y/nm, HSCALE=(-5*sigx/nm,5*sigx/nm), VSCALE=(-10*sigy/nm,10*sigy/nm), TITLE='Right-Going Primary Photon IP Spot;', HTITLE='X (nm);', VTITLE='Y (nm);' ; PLOT SCAT, KIND=photon, RIGHT, GEN=2, H=En/1D9, V=Sqrt[X^2+Y^2]/nm, HSCALE=(0,250.5), VSCALE=(0,100), TITLE='Primary Photon Energy vs. R at IP;', HTITLE='E0G1 (GeV); XGX ;', VTITLE='r (nm);' ; PLOT SCAT, KIND=photon, RIGHT, GEN=2, H=Xi2, V=Sqrt[X^2+Y^2]/nm, HSCALE=(-1,1), VSCALE=(0,100), TITLE='Primary Photon Helicity vs. R at IP;', HTITLE='X021; X X;', VTITLE='r (nm);' ; PLOT SCAT, KIND=photon, RIGHT, GEN=2, H=En/1D9, V=Xi2, HSCALE=(0,250.5), VSCALE=(-1,1), TITLE='Primary Photon Energy vs. Helicity after CP;', HTITLE='E0G1 (GeV); XGX ;', VTITLE='X021; X X;',; SET MsgLevel=1; CLEAR LASER; CLEAR LASERQED; LUMINOSITY KIND=(electron,electron), W=(0,501D9,50), WX=(8*sigx/2,16*sigx),WY=(16*sigy,64*sigy), FREP=90*180 ; LUMINOSITY KIND=(photon,photon), W=(0,501D9,50); LUMINOSITY KIND=(photon,electron), W=(0,501D9,50); LUMINOSITY KIND=(photon,positron), W=(0,501D9,50); LUMINOSITY KIND=(electron,positron), W=(0,501D9,50); BBFIELD NX=32, NY=32, WX=8*sigx, R=sigx/sigy/8, WXMAX=16*sigx; CFQED BEAMSTRAHLUNG, POL; CFQED PAIR, POL; FLAG OFF SPIN; SET MsgLevel=0; SET it=0; PRINT CPUTIME; PUSH Time=(-2.5*sigz,2.5*sigz,300) ; IF Mod(it,10)=0; PRINT it, FORMAT=(F6.0,'-th time step'); PRINT STAT, SHORT; ENDIF; IF it=150; PLOT SCAT, KIND=electron, RIGHT, H=S/micron, V=Y/nm, HSCALE=(-2.5*sigz/micron,2.5*sigz/micron), VSCALE=(-50*sigy/nm,50*sigy/nm), TITLE='Electron at IP;', HTITLE='S (Mm); G ;', VTITLE='Y (nm);' ; PLOT BBFIELD, S=0; ENDIF; SET it=it+1; ENDPUSH; PRINT CPUTIME; DRIFT S=0; PLOT LUMINOSITY; PLOT SCAT, KIND=electron, H=En/1D9, V=Sqrt[(Px^2+Py^2)/Ps^2]/mm, HSCALE=(0,250.5), VSCALE=(0,10), TITLE='Electron Energy vs. Angle after IP;', HTITLE='E0e1 (GeV); X X ;', VTITLE='Q0e1 (mrad);GX X ;' ; STOP;