The OPAL Discussion Forum

["Adelmann Andreas (PSI)" <andreas.adelmann AT psi.ch>]

Hi Nicole you need to specify

 EMITTED = True 

Here a Distribution that works:

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,

      CUTOFFPX = 5,

      CUTOFFPY = 5, CUTOFFLONG = 5, CUTOFFPZ = 5,

      TRISE = 1.0745e-12,       // (s)

      TFALL = 1.0745e-12,       // See notes above for equation of TRISE

      TPULSEFWHM = 1.5e-12,     // FWHM length in time (s)

      CUTOFFLONG = 4.0,

      NBIN = 9,

      EMISSIONSTEPS = 100,

      EMISSIONMODEL = ASTRA, //NONEQUIL,

      EKIN = 0.2,

      ELASER = 5.0,

      W = 3.2,

      FE = 3.2,

      CATHTEMP = 321.95,

      EMITTED = True,

      WRITETOFILE = True;

You maybe want to adjust some emission parameters according to your problem.

A

------
Dr. sc. math. Andreas (Andy) Adelmann
Staff Scientist
Paul Scherrer Institut WBBA/219 CH-5232 Villigen PSI
Phone Office: xx41 56 310 42 33 Fax: xx41 56 310 31 91
Phone Home: xx41 62 891 91 44
-------------------------------------------------------
Friday: ETH HPK G 28   +41 44 632 75 22
============================================
The more exotic, the more abstract the knowledge, 
the more profound will be its consequences.
Leon Lederman 
============================================

On 14 Dec 2017, at 07:29, Nicole Neveu <nneveu AT hawk.iit.edu> wrote:

Hi All,

I'm currently on OPAL-t git version b85df44 (v1.9).

I'm trying to make a mono-energetic injected pencil-like beam.

I defined P0, and used the following flattop distribution:

DIST1: DISTRIBUTION, TYPE = FLATTOP,
    SIGMAX = 1.0e-6, SIGMAPX = 0.0,
    SIGMAY = 1.0e-6, SIGMAPY = 0.0,
    SIGMAZ = 5.0e-5, SIGMAPZ = 0.0,
    WRITETOFILE = True;

There's not much else in the file besides a drift definition and a space charge solver.

I'm getting the following error:

[neveu@beboplogin3 euclid]$ more bdw.187328.bdw-0197.error                                    
Error{0}>             
Error{0}> *** User error detected by function "Quaternion::rotate()"                          
Error{0}> *** in line 108 of file "emgun.in" at end of statement:                             
Error{0}>     quaternion isn't unit quaternion. Norm: nan  

Am I missing a required input to the distribution?

I've tried a couple permutations of options, but same result.

Also, the following Gauss distribution works in the same file:

DIST2: DISTRIBUTION,
    SIGMAX = 1.0e-6, SIGMAPX = 0.0,
    SIGMAY = 1.0e-6, SIGMAPY = 0.0,
    SIGMAZ = 6.25e-5, SIGMAPZ = 0.0,
    CORRX = 0.0, CORRY = 0.0, CORRZ = 0,
    R51 = 0, R52 = 0, R61 = 0, R62 = 0,
    CUTOFFX = 5, CUTOFFPX = 5, CUTOFFY = 5,
    CUTOFFPY = 5, CUTOFFLONG = 5, CUTOFFPZ = 5,
    OFFSETX = 0.0, OFFSETY =0.0, OFFSETZ = 0.0,
    TYPE = GAUSS, INPUTMOUNITS = NONE,
    WRITETOFILE = True;

Thanks,

Nicole