Option, ECHO=FALSE;
Option, STATDUMPFREQ=30000000;
Option, PSDUMPFREQ=100000000;
Option, SPTDUMPFREQ = 1000000000;
Option, PSDUMPEACHTURN=true;
Option, VERSION=100600;
//Option, PSDUMPLOCALFRAME=false;

Title,string="Test test";

REAL turns = 40;
REAL steps = 1800;

REAL Edes= 2.95681E-01/1000.0; //GeV
REAL phi0 = -1.19999E+01;    //deg  par loc
REAL rinit0= 2.63518E+02; // mm 
REAL z0= 0.0; // mm   
REAL pr0 = -1.61600E-04; //rad  
REAL pz0= 0.0; //rad  

REAL gap=8.97*25.4; //mm
REAL frequency = 23.05507; //MHz
REAL harm = 5;

REAL hmass=0.939277; // GeV
REAL gamma=(Edes+hmass)/hmass;
REAL beta=sqrt(1-(1/gamma^2));
REAL gambet=gamma*beta;
REAL Moment0 = gamma*beta*hmass;
REAL brho = (hmass*1.0e9*gambet) / CLIGHT;

REAL phi01=(128+180)/360.0*twopi; // TAU0 + 90deg  in radian

value,{gamma,Edes,frequency,pr0};

cycTest: Cyclotron, TYPE="BANDRF", CYHARMON= harm, PHIINIT= phi0, PRINIT= pr0,ZINIT=z0, PZINIT=pz0, RINIT= rinit0, SYMMETRY= 1.0, FMAPFN="../Test_Bz_5p5.dat",
           RFPHI= phi01, RFFREQ= frequency, RFMAPFN= "../Test_y15_rd2x.h5part", ESCALE= 200.0,
	         MAXZ=15.875, MINZ=-15.875;  //  0.625*25.4
l1: Line = (cycTest);

Dist1:DISTRIBUTION, DISTRIBUTION=fromfile,FNAME="../ParDataBaseCluster.dat"; 

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

beam1: BEAM, PARTICLE=HMINUS, pc=Moment0, SPACECHARGE=false, NPART=1000001, BCURRENT=0.000500, CHARGE=1.0, BFREQ= frequency;

Select, Line=l1;
TRACK,LINE=l1, BEAM=beam1, MAXSTEPS=steps*turns, STEPSPERTURN=steps;
 run, method = "CYCLOTRON-T", beam=beam1, fieldsolver=Fs1, distribution=Dist1;
endtrack;
Stop;