GitLab

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;