OPTION, ECHO=FALSE;
OPTION, INFO=FALSE;
OPTION, PSDUMPFREQ=50000000;
OPTION, TELL=TRUE;

// OPAL Input forLA1line
TDeg = 1e-12;
Edes = 0.249489;
gamma = (Edes+PMASS)/PMASS;
beta = sqrt(1-(1/gamma^2));
gambet = gamma*beta;
P0 = gamma*beta*PMASS;
brho = (PMASS*1.0e9*gambet) / CLIGHT;
VALUE,{gamma,brho,Edes,beta,gambet};

rf = 50.6328e6;
NProt = 5e5;

DISTRIB1: DISTRIBUTION,
SIGMAX = 0.005, SIGMAPX = 0.0,
SIGMAY = 0.005, SIGMAPY = 0.0,
SIGMAZ = 1e-05, SIGMAPZ = 0.0,
CORRX = 0.503, CORRY = 0.858, CORRZ = 0,
R51 = 0, R52 = 0, R61 = 0, R62 = 0,
CUTOFFX = 5, CUTOFFPX = 5, CUTOFFY = 5,
CUTOFFPY = 5, CUTOFFLONG = 5, CUTOFFPZ = 5,
OFFSETX = 0.0, OFFSETY =0.0, OFFSETZ = 0.0,
DISTRIBUTION = GAUSS, INPUTMOUNITS = NONE; 

D15: DRIFT, L=0.02, ELEMEDGE=0.0;

MDegD0: MONITOR, OUTFN="MDegD0.h5", ELEMEDGE=0.01;

WEDGE_HLEN=0.0197293;
WEDGE_FLEN=2*WEDGE_HLEN;
DRIFT_LENGTH=0.04111-WEDGE_FLEN;
START = 0.02;

VALUE,{WEDGE_HLEN, WEDGE_FLEN, DRIFT_LENGTH};

//surface physics
DEGPHYS_Wedge :  SURFACEPHYSICS, TYPE="DEGRADER", MATERIAL="GraphiteR6710";

// DMAD
Wedge1: DEGRADER, L=WEDGE_HLEN, ZSIZE= WEDGE_HLEN, OUTFN="Wedge1.h5", SURFACEPHYSICS=DEGPHYS_Wedge, ELEMEDGE=START;
DegD1: DRIFT, L=DRIFT_LENGTH, ELEMEDGE=START+WEDGE_HLEN;
MDegD1: MONITOR, OUTFN="MDegD1.h5", ELEMEDGE=START+WEDGE_HLEN+DRIFT_LENGTH/2;

Wedge2: DEGRADER, L=WEDGE_FLEN, ZSIZE= 2*WEDGE_HLEN, OUTFN="Wedge2.h5", SURFACEPHYSICS=DEGPHYS_Wedge, ELEMEDGE=START+WEDGE_HLEN+DRIFT_LENGTH;
DegD2: DRIFT, L=DRIFT_LENGTH, ELEMEDGE=START+3*WEDGE_HLEN+DRIFT_LENGTH;
MDegD2: MONITOR, OUTFN="MDegD2.h5", ELEMEDGE=START+3*WEDGE_HLEN+DRIFT_LENGTH + DRIFT_LENGTH/2;

Wedge3: DEGRADER, L=WEDGE_FLEN, ZSIZE= 2*WEDGE_HLEN, OUTFN="Wedge3.h5", SURFACEPHYSICS=DEGPHYS_Wedge, ELEMEDGE=START+3*WEDGE_HLEN+2*DRIFT_LENGTH;
DegD3: DRIFT, L=DRIFT_LENGTH, ELEMEDGE=START+5*WEDGE_HLEN+2*DRIFT_LENGTH;
MDegD3: MONITOR, OUTFN="MDegD3.h5", ELEMEDGE=START+5*WEDGE_HLEN+2*DRIFT_LENGTH + DRIFT_LENGTH/2;

Wedge4: DEGRADER, L=WEDGE_FLEN, ZSIZE= 2*WEDGE_HLEN, OUTFN="Wedge4.h5", SURFACEPHYSICS=DEGPHYS_Wedge, ELEMEDGE=START+5*WEDGE_HLEN+3*DRIFT_LENGTH;
DegD4: DRIFT, L=DRIFT_LENGTH, ELEMEDGE=START+7*WEDGE_HLEN+3*DRIFT_LENGTH;
MDegD4: MONITOR, OUTFN="MDegD4.h5", ELEMEDGE=START+7*WEDGE_HLEN+3*DRIFT_LENGTH + DRIFT_LENGTH/2;

Wedge5: DEGRADER, L=WEDGE_FLEN, ZSIZE= 2*WEDGE_HLEN, OUTFN="Wedge5.h5", SURFACEPHYSICS=DEGPHYS_Wedge, ELEMEDGE=START+7*WEDGE_HLEN+4*DRIFT_LENGTH;
DegD5: DRIFT, L=DRIFT_LENGTH, ELEMEDGE=START+9*WEDGE_HLEN+4*DRIFT_LENGTH;
MDegD5: MONITOR, OUTFN="MDegD5.h5", ELEMEDGE=START+9*WEDGE_HLEN+4*DRIFT_LENGTH + DRIFT_LENGTH/2;

Wedge6: DEGRADER, L=WEDGE_HLEN, ZSIZE= WEDGE_HLEN, OUTFN="Wedge6.h5", SURFACEPHYSICS=DEGPHYS_Wedge, ELEMEDGE=START+9*WEDGE_HLEN+5*DRIFT_LENGTH;
END: DRIFT, L=4.0, ELEMEDGE=START+9*WEDGE_HLEN+5*DRIFT_LENGTH+WEDGE_HLEN;
MDegD6: MONITOR, OUTFN="MDegD6.h5", ELEMEDGE=DRIFT_LENGTH/2+START+9*WEDGE_HLEN+5*DRIFT_LENGTH+WEDGE_HLEN;

M_Moliere: MONITOR, OUTFN="M_Moliere.h5", ELEMEDGE=1.22555;


BEAMLINE_A1: LINE=(D15,MDegD0,Wedge1,Deg_D1,MDegD1,Wedge2,Deg_D2,MDegD2,Wedge3,Deg_D3,MDegD3,Wedge4,Deg_D4,MDegD4,Wedge5,Deg_D5,MDegD5,Wedge6,END,MDegD6,M_Moliere);

FS1:FIELDSOLVER, FSTYPE=NONE, MX=64, MY=64, MT=64, PARFFTX=true, PARFFTY=true, PARFFTT=true, BCFFTX=open, BCFFTY=open, BCFFTT=open, BBOXINCR=1, GREENSF=STANDARD;
BEAM_A1: BEAM, PARTICLE=PROTON, PC=P0, NPART=NProt, BCURRENT=NProt*1.6e-19*rf, BFREQ=rf,  CHARGE=1;

SELECT, LINE=BEAMLINE_A1;
TRACK, LINE=BEAMLINE_A1, BEAM=BEAM_A1, MAXSTEPS= 5e18, DT=TDeg, ZSTOP=4.22555;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_A1, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STOP;