OPTION, VERSION = 20000;

OPTION, ECHO= FALSE;
OPTION, PSDUMPFREQ= 10;
OPTION, PSDUMPFRAME= REFERENCE;

//////////////////////////////////////////////////////////////////////////////
// input file for single particle tracking                  //
// the initial conditions are obtained from FIXPO GORBIT mode.              //
//////////////////////////////////////////////////////////////////////////////

TITLE, STRING= "Decy-13";

REAL Edes = 0.000040;
REAL pr0  = 0.00000;
REAL r0   = 22.3646;

REAL phi01 = 8.6;
REAL phi02 = phi01+146.0;

REAL volt1st = 0.040;

REAL r1   = 10.00;
REAL r2   = 500.00;

REAL f1 = 77.64;

REAL gamma=(Edes+PMASS)/PMASS;
REAL beta=sqrt(1-(1/gamma^2)); 
REAL gambet=gamma*beta; 
REAL P0 = gamma*beta*PMASS;

Ring: CYCLOTRON, TYPE="BANDRF", CYHARMON=4, PHIINIT=64.0, PRINIT=pr0, RINIT=-r0 , SYMMETRY=4.0, RFFREQ=f1, FMAPFN="sfield.dat",RFPHI=Pi, RFMAPFN="efield_CYC.h5part", ESCALE=1.0;

Rf0: RFCAVITY, VOLT=volt1st, FMAPFN="rffield1.dat", TYPE="SINGLEGAP", 
FREQ=f1, RMIN = r1, RMAX = r2, ANGLE=26.75,  PDIS = 0.0, GAPWIDTH = 0.0, PHI0=phi01; 

Rf1: RFCAVITY, VOLT=volt1st, FMAPFN="rffield1.dat", TYPE="SINGLEGAP", 
FREQ=f1, RMIN = r1, RMAX = r2, ANGLE=63.25, PDIS = 0.0, GAPWIDTH = 0.0, PHI0=phi02; 

Rf2: RFCAVITY, VOLT=volt1st, FMAPFN="rffield1.dat", TYPE="SINGLEGAP", 
FREQ=f1, RMIN = r1, RMAX = r2, ANGLE=206.75, PDIS = 0.0, GAPWIDTH = 0.0, PHI0=phi01; 

Rf3: RFCAVITY, VOLT=volt1st, FMAPFN="rffield1.dat", TYPE="SINGLEGAP", 
FREQ=f1, RMIN = r1, RMAX = r2, ANGLE=243.25, PDIS = 0.0, GAPWIDTH = 0.0, PHI0=phi02; 

L1:   LINE = (Ring);

Dist1:DISTRIBUTION, TYPE=FROMFILE,FNAME="dist2.dat"; 

Fs1:FIELDSOLVER, FSTYPE=NONE, MX=64, MY=64, MT=64, PARFFTX=true, PARFFTY=true, PARFFTT=false,
         BCFFTX=open, BCFFTY=open, BCFFTT=open;

Beam1: BEAM, PARTICLE=HMINUS, pc=P0, NPART=1, BCURRENT=1.0E-6, CHARGE=-1.0, BFREQ= f1;

SELECT, LINE=L1;

TRACK, LINE= L1, BEAM= Beam1, MAXSTEPS= 720*150, STEPSPERTURN= 720;
 RUN, METHOD = "CYCLOTRON-T", BEAM= Beam1, FIELDSOLVER= Fs1, DISTRIBUTION= Dist1;
ENDTRACK;
STOP;