Skip to Content.
Sympa Menu

opal - Re: [Opal] Pencil beam

opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

Re: [Opal] Pencil beam


Chronological Thread 
  • From: "Adelmann Andreas (PSI)" <andreas.adelmann AT psi.ch>
  • To: Nicole Neveu <nneveu AT hawk.iit.edu>
  • Cc: opal <opal AT lists.psi.ch>
  • Subject: Re: [Opal] Pencil beam
  • Date: Thu, 14 Dec 2017 06:42:42 +0000
  • Accept-language: en-US, de-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




Archive powered by MHonArc 2.6.19.

Top of Page