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: "Adelmann Andreas (PSI)" <andreas.adelmann AT psi.ch>, "opal AT lists.psi.ch" <Opal AT lists.psi.ch>
  • Subject: Re: [Opal] Pencil beam
  • Date: Thu, 14 Dec 2017 07:46:16 +0000
  • Accept-language: en-US, de-CH

There seams to be a problem with the injected distribution (https://gitlab.psi.ch/OPAL/src/issues/194

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 08:35, Nicole Neveu <nneveu AT hawk.iit.edu> wrote:

I wanted an injected beam.
Will this distribution only work when emitted?

Nicole

On Thu, Dec 14, 2017 at 1:20 AM, Adelmann Andreas (PSI) <andreas.adelmann AT psi.ch> wrote:
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






Archive powered by MHonArc 2.6.19.

Top of Page