//-----------------
// Simple quad test
//
//-----------------
TITLE, STRING="Simple quad test";


OPTION, ECHO=TRUE; //ECHO Controls printing of an echo of input lines on the standard error file.
//OPTION, TFS =FALSE; //TFS If true, the format is in TFS format, otherwise in human-readable format.
OPTION, INFO=TRUE; //INFO If this option is turned off, OPAL suppresses all information messages.
OPTION, PSDUMPFREQ=1000; //Defines after how many time steps the phase space is dumped into the H5Part file. Default value is 10.
OPTION, STATDUMPFREQ=1000;
OPTION, AUTOPHASE = 1;
OPTION, EBDUMP=TRUE;;

QB               = 100e-12;
BFREQ            = 1.3e9;
BCURRENT         = QB*BFREQ;

call, file = "lat.opal";

file_distribution:DISTRIBUTION, DISTRIBUTION=fromfile, FNAME="lat.opal_particles";


sigT = 1.50828e-12;
tlen = sqrt(12) * sigT;
xyrms  = 0.000456094;
//radius = xyrms * 2;
//KT = 0.120 //eV;
//E0 = 3.0/2.0*KT;


thermal_distribution: DISTRIBUTION, DISTRIBUTION = GUNGAUSSFLATTOPTH,
    sigmax = xyrms,
	sigmay = xyrms,
  pt = 1.0, 
	tRise = 0,
	tFall = 0, 
	tPulseFWHM = tlen, 
	ekin = .180, // = 3/2*kT when kT = 0.120 eV 
	NBIN = 50, DEBIN = 80;

astra_distribution:DISTRIBUTION, DISTRIBUTION = "ASTRAFLATTOPTH",
  SIGMAX = 2*0.000275, SIGMAPX = 0.0, CORRX = 0.0,
  SIGMAY = 2*0.000275, SIGMAPY = 0.0, CORRY = 0.0,
  SIGMAT = 0.0, PT = 0.0, SIGMAPT = 0.0, CORRT = 0.0,
  TRISE = 7.82e-13, TFALL = 7.82e-13, TPULSEFWHM = 9.9e-12,
  EKIN = 0.63, NBIN = 20, DEBIN = 80;

no_fieldsolver:FIELDSOLVER,
	FSTYPE = NONE, MX = 16, MY = 16, MT = 32,
	PARFFTX = true, PARFFTY = true, PARFFTT = true,
	BCFFTX = open, BCFFTY = open, BCFFTT = open,
	BBOXINCR = 1, GREENSF = INTEGRATED;

fft_fieldsolver:FIELDSOLVER, FSTYPE=FFT, MX=32, MY=32, MT=32, 
       PARFFTX=true, PARFFTY=true, PARFFTT=true,
       BCFFTX=open, BCFFTY=open, BCFFTT=open, 
       BBOXINCR=1, GREENSF=INTEGRATED;


my_beam: BEAM,
	PARTICLE=ELECTRON, 
	//pc = 10,          // momentum per particle in GeV/c
	NPART = 100,     //number of test particles simulated
	BCURRENT = BCURRENT, // beam current in A
	BFREQ = BFREQ;
//  CHARGE = -1;

SELECT, LINE=lattice;

// Simulation
TRACK, 	LINE = lattice, 
	BEAM = my_beam, 
	MAXSTEPS = 1000000,
	DT=10.0e-13, //Timestep size in s
	ZSTOP=6.0; //Stopping position z in meters
	
RUN,	METHOD = "PARALLEL-T",  
	BEAM = my_beam, 
	FIELDSOLVER = no_fieldsolver, 
    DISTRIBUTION = file_distribution;
	//DISTRIBUTION = thermal_distribution;
	//DISTRIBUTION = astra_distribution;
ENDTRACK;

STOP;