Skip to Content.
Sympa Menu

opal - Re: [Opal] Reference particle in OPAL

opal AT

Subject: The OPAL Discussion Forum

List archive

Re: [Opal] Reference particle in OPAL

Chronological Thread  
  • From: Nicole R Neveu <nneveu AT>
  • To: Wei Hou Tan <tanweihou AT>, opal <opal AT>
  • Subject: Re: [Opal] Reference particle in OPAL
  • Date: Fri, 22 Jan 2021 22:25:37 +0000
  • Accept-language: en-US
  • Arc-authentication-results: i=1; 1; spf=pass; dmarc=pass action=none; dkim=pass; arc=none
  • Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed;; s=arcselector9901; h=From:Date:Subject:Message-ID:Content-Type:MIME-Version:X-MS-Exchange-SenderADCheck; bh=ZtsWGJrfpR48vBsH7Kf5JoXfkK/bs6ytIx/zh+NDZ6w=; b=AbTiffxMIxuS6u4NkYywr8EVe0Gy6QQ8BpZ+scCL4SohTKSnuUeLZRDWDwxBekGmyBbpw2BxS1W8PAe9c3AVJ5nnFCI82b1hLPUMIWN1tLVZlZvva30uxMNz0rmqo9WL0inaoGtwriLiZCVZzqtVarFoa1oL0Q4oKig4rHT42FB5GuebGvYCYjDNKNlwEXhR2c7zRL456WwhzqNAGo7AUG5f9qNht1vQQoOlU0e/HxbRJ1MXEG87lgTAw+IrkejCJXJm1a0lZ7R5Rfu0+lfR3EvLj7Q+CWnw2UXpkGQIyFo3YU/HZOslpe8iGqAFpu6hyiF/PggIbcfaD7chmtmuHA==
  • Arc-seal: i=1; a=rsa-sha256; s=arcselector9901;; cv=none; b=c2uAWpGPAosYpgrr/BvNJ8MNeg/vbvG09YtP3QQ+KSQgjMKVeXk6bprzoPE+l0j+Irqh+nHJlaGGh3ZyRi5YbNfmS6RntBucklbTvtespu6ZAgkSr2yMZIMOdzaX67wqyBi9BWAQ3aM5ZdDJFi+3ngefpp6xm61RSv+Ppkdtr0P96rVEaYL9FD9h2UE/fSLS+kakY1WcNUJ9+lllY6u+Brq5tkz0hV/Z3NfoL/sYt4UNRPi6J8gb1w4+9OpC4nU3oC9f1bpYJy2SDmrfPrzuNqLEyYfL8nRYSzaLh/21WQGeSpwR0NFv6XfnkOJjgWa5+FM8Wl7dkkc8bnc61SBShQ==
  • Authentication-results:; dkim=none (message not signed) header.d=none;; dmarc=none action=none;

Have you tried more emission steps?

I think the default is only ‘1’, and this value was too low for my case.

Here’s my standard distribution:



                   FNAME = "opal_emitted.txt",

                   EMITTED = TRUE,

                   EMISSIONMODEL = NONE,

                   EMISSIONSTEPS = 100,



Other things that I needed to adjust were:

  • Grid size / # particles ratio (maybe you’ve already done a convergence study?)
  • Number of Fourier samples in the field maps (originally I had this number too high)






From: Wei Hou Tan <tanweihou AT>
Date: Friday, January 22, 2021 at 2:16 PM
To: Nicole R Neveu <nneveu AT>, opal <opal AT>
Subject: Re: [Opal] Reference particle in OPAL


Dear Nicole,


The distribution was generated by external scripts, which then converted into ASTRA and OPAL inputs.

The distribution comes with initial momentum pz, I did not add extra Ekin as shown in the code below.


Here's the distribution command look like



        FNAME = "opal.ini",



        EMITTED = True,


        NBIN = 100;





Wei Hou

From: Nicole R Neveu <nneveu AT>
Sent: Saturday, January 23, 2021 6:07 AM
To: Wei Hou Tan <tanweihou AT>; opal <opal AT>
Subject: Re: [Opal] Reference particle in OPAL


Hi Wei Hou,


When I did a similar exercise for the LCLS-II injector, one cause of discrepancy for me was the definition of Ekin.

Are you suppling the same particles distribution to ASTRA and OPAL?

Or letting the codes generate the distribution?

What does your distribution command look like?







From: <opal-request AT> on behalf of Wei Hou Tan <tanweihou AT>
Reply-To: Wei Hou Tan <tanweihou AT>
Date: Friday, January 22, 2021 at 1:59 PM
To: opal <opal AT>
Subject: [Opal] Reference particle in OPAL


Dear OPAL developers,


This is Wei Hou TAN, graduate student working with Prof. Philippe Piot in the Northern Illinois University.


I would like to learn OPAL, starting by cross-checking its results with ASTRA.


In my toy example, I simulated a rf gun with both ASTRA and OPAL.

However, I am having problems getting same results in OPAL.

I have checked that both ASTRA and OPAL have the same distributions, same RF phase.


I suspected that OPAL has a different definition of reference particle but I could not find any documentations on it.


What is OPAL convention of reference particle ?


Thank you.



Wei Hou

Archive powered by MHonArc 2.6.19.

Top of Page