Skip to content

GitLab

  • Projects
  • Groups
  • Snippets
  • Help
    • Loading...
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
  • Sign in
S src
  • Project overview
    • Project overview
    • Details
    • Activity
    • Releases
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
    • Locked Files
  • Issues 62
    • Issues 62
    • List
    • Boards
    • Labels
    • Service Desk
    • Milestones
    • Iterations
  • Merge requests 4
    • Merge requests 4
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Operations
    • Operations
    • Incidents
    • Environments
  • Analytics
    • Analytics
    • CI/CD
    • Code Review
    • Issue
    • Repository
    • Value Stream
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Members
    • Members
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • OPAL
  • src
  • Issues
  • #194

Closed
Open
Created Dec 14, 2017 by adelmann@adelmann🎗Owner

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;
Edited Dec 14, 2017 by snuverink_j
Assignee
Assign to
OPAL 1.9.x
Milestone
OPAL 1.9.x
Assign milestone
Time tracking