The OPAL Discussion Forum

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

Hi Nicol, please but me the input file at a readable spot on Bebop

Tnx, 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 CAB F37 or 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