Skip to Content.
Sympa Menu

opal - [Opal] Pencil beam

opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

[Opal] Pencil beam


Chronological Thread 
  • From: Nicole Neveu <nneveu AT hawk.iit.edu>
  • To: opal <opal AT lists.psi.ch>
  • Subject: [Opal] Pencil beam
  • Date: Thu, 14 Dec 2017 00:29:55 -0600

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



Archive powered by MHonArc 2.6.19.

Top of Page