Injected FT gives quaternion isn't unit quaternion. Norm: nan
This input file works with DistOk and gives an error quaternion isn't unit quaternion. Norm: nan if one tries to inject: Dist
OPTION, ECHO=FALSE;
OPTION, INFO=FALSE;
OPTION, PSDUMPFREQ = 1000000;
OPTION, STATDUMPFREQ = 10;
OPTION, VERSION=10901;
Title, string="Gun and cavity for Euclid";
REAL BEAM_BUNCH_CHARGE = 10e-12; // Beam bunch charge (C) = I/f = (density * 2Pi * r^2 )/freq
REAL RF_FREQUENCY = 26000.0; // using 26 Ghz
REAL gamma = 1.39139;
REAL beta = sqrt(1-(1/gamma^2));
REAL gammabeta = gamma*beta;
REAL P0 = gamma*beta*EMASS;
DR1: DRIFT, L = 0.5, ELEMEDGE = 0.0;
L1: Line = (DR1);
DistOk: DISTRIBUTION, TYPE = FLATTOP,
SIGMAX = 1.0e-6, SIGMAPX = 0.0,
SIGMAY = 1.0e-6, SIGMAPY = 0.0,
SIGMAZ = 5.0e-5, SIGMAPZ = 0.0,
TRISE = 1.0745e-12, // (s)
TFALL = 1.0745e-12, // See notes above for equation of TRISE
TPULSEFWHM = 1.5e-12, // FWHM length in ti
NBIN = 9,
EMISSIONSTEPS = 100,
EMISSIONMODEL = ASTRA,
EKIN = 0.2,
ELASER = 5.0,
W = 3.2,
FE = 3.2,
CATHTEMP = 321.95,
EMITTED = True,
WRITETOFILE = True;
Dist: DISTRIBUTION, TYPE = FLATTOP,
SIGMAX = 1.0e-6, SIGMAPX = 0.0,
SIGMAY = 1.0e-6, SIGMAPY = 0.0,
SIGMAZ = 5.0e-5, SIGMAPZ = 0.0,
TRISE = 1.0745e-12, // (s)
TFALL = 1.0745e-12, // See notes above for equation of TRISE
TPULSEFWHM = 1.5e-12, // FWHM length in ti
EKIN = 0.2,
ELASER = 5.0,
W = 3.2,
FE = 3.2,
CATHTEMP = 321.95,
EMITTED = False,
WRITETOFILE = True;
FS_SC: Fieldsolver, FSTYPE = FFT,
MX = 8, MY = 8, MT = 16,
PARFFTX = false, PARFFTY = false, PARFFTT = true,
BCFFTX = open, BCFFTY = open, BCFFTT = open,
BBOXINCR = 1, GREENSF = INTEGRATED;
value,{gamma,beta,gammabeta};
beam1: BEAM, PARTICLE = ELECTRON, PC = P0, NPART = 15000,
BCURRENT = BEAM_BUNCH_CHARGE * 26e9, BFREQ = RF_FREQUENCY, CHARGE = -1;
TRACK, LINE = L1, BEAM = beam1, MAXSTEPS = 100000000, DT = 1.0e-13, ZSTOP=0.5;
RUN, METHOD = "PARALLEL-T", BEAM = beam1, FIELDSOLVER = FS_SC, DISTRIBUTION = DIST;
ENDTRACK;
STOP;