The OPAL Discussion Forum

[Wei Hou Tan <tanweihou AT outlook.com>]

Dear Andreas,

Thank you, unfortunately it still does not work.

Here is how I read the final output distribution from HDF5 file. Is it correct ?

Thanks.

Regards

Wei Hou

---

From: Adelmann Andreas (PSI) <andreas.adelmann AT psi.ch>
Sent: Tuesday, January 26, 2021 5:07 AM
To: Wei Hou Tan <tanweihou AT outlook.com>
Cc: Nicole R Neveu <nneveu AT stanford.edu>; opal AT lists.psi.ch <opal AT lists.psi.ch>
Subject: Re: [Opal] Reference particle in OPAL

 

Hello Wei Hou can you change EMISSIONSTEPS = 200, 

to lower values and see if the distribution changes? 

These small changes in dE are maybe caused by slightly 

different emission models. 

Cheers A

------
Dr. sc. math. Andreas (Andy) Adelmann
Head a.i. Labor for Scientific Computing and Modelling 
Paul Scherrer Institut OHSA/ CH-5232 Villigen PSI
Phone Office: xx41 56 310 42 33 Fax: xx41 56 310 31 91
Zoom ID: 470-582-4086 Password: AdA

Zoom Link: https://ethz.zoom.us/j/4705824086?pwd=dFcvT1pMMGY0bHg0dTNncUNZZTJkZz09

-------------------------------------------------------
Friday: ETH HPK G 28   +41 44 633 3076
============================================
The more exotic, the more abstract the knowledge, 
the more profound will be its consequences.
Leon Lederman 
============================================

On 23 Jan 2021, at 15:57, Wei Hou Tan <tanweihou AT outlook.com> wrote:

Hi Andreas,

Thanks for the suggestion, my apology for misunderstanding "energy" in <output>.stat.

It seems like the OPAL and ASTRA agree for the reference particle, but

My input bunch is quite long, about 90ps full length, and here is my distribution command.

Dist: DISTRIBUTION, TYPE = FROMFILE,

        FNAME = "opal.ini",

        EMISSIONMODEL = NONE,

        EKIN=0,

        EMITTED = True,

        WRITETOFILE=True,

        EMISSIONSTEPS = 200,

        NBIN = 200;

Is there anything wrong with it ?

<Screenshot 2021-01-23 082555.jpg>

<Screenshot 2021-01-23 085210.jpg>

Thanks.

Regards

Wei Hou

---

From: Adelmann Andreas (PSI) <andreas.adelmann AT psi.ch>
Sent: Saturday, January 23, 2021 4:59 PM
To: Wei Hou Tan <tanweihou AT outlook.com>
Cc: Nicole R Neveu <nneveu AT stanford.edu>; opal AT lists.psi.ch <opal AT lists.psi.ch>
Subject: Re: [Opal] Reference particle in OPAL

 

Hi Wei Hou check 6.7 in http://amas.web.psi.ch/opal/Documentation/master/OPAL_Manual.html#chp.opalt

here you see what type of output you have. 

Consider to plot column 23 (Z-component of normalized momentum of reference particle in floor coordinate system)

Cheers A

------
Dr. sc. math. Andreas (Andy) Adelmann
Head a.i. Labor for Scientific Computing and Modelling 
Paul Scherrer Institut OHSA/ CH-5232 Villigen PSI
Phone Office: xx41 56 310 42 33 Fax: xx41 56 310 31 91
Zoom ID: 470-582-4086 Password: AdA

Zoom Link: https://ethz.zoom.us/j/4705824086?pwd=dFcvT1pMMGY0bHg0dTNncUNZZTJkZz09

-------------------------------------------------------
Friday: ETH HPK G 28   +41 44 633 3076
============================================
The more exotic, the more abstract the knowledge, 
the more profound will be its consequences.
Leon Lederman 
============================================

On 23 Jan 2021, at 00:30, Wei Hou Tan <tanweihou AT outlook.com> wrote:

Dear Nicole,

Thanks for your suggestion.

Unfortunately, it did not help much.

For my example, I did not turn on space charge calculation.

From the plot below, it seems like there is a mismatch in RF phase, as seen by the reference particle.

Is there a way to check which reference particle that OPAL uses for its calculation ? (For ASTRA, it is the bunch centroid)

<dasdasdas.jpg>

Thanks.

Regards

Wei Hou

---

From: Nicole R Neveu <nneveu AT stanford.edu>
Sent: Saturday, January 23, 2021 6:25 AM
To: Wei Hou Tan <tanweihou AT outlook.com>; opal <opal AT lists.psi.ch>
Subject: Re: [Opal] Reference particle in OPAL

 

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:

 

Dist:DISTRIBUTION, TYPE = FROMFILE,

                   FNAME = "opal_emitted.txt",

                   EMITTED = TRUE,

                   EMISSIONMODEL = NONE,

                   EMISSIONSTEPS = 100,

                   EKIN=0;

 

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)

 

Thanks,

 

Nicole

 

From: Wei Hou Tan <tanweihou AT outlook.com>
Date: Friday, January 22, 2021 at 2:16 PM
To: Nicole R Neveu <nneveu AT stanford.edu>, opal <opal AT lists.psi.ch>
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

 

Dist: DISTRIBUTION, TYPE = FROMFILE,

        FNAME = "opal.ini",

        EMISSIONMODEL = NONE,

        EKIN=0,

        EMITTED = True,

        WRITETOFILE=True,

        NBIN = 100;

 

Thanks.

 

Regards

Wei Hou

---

From: Nicole R Neveu <nneveu AT stanford.edu>
Sent: Saturday, January 23, 2021 6:07 AM
To: Wei Hou Tan <tanweihou AT outlook.com>; opal <opal AT lists.psi.ch>
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?

 

Thanks,

 

Nicole

 

 

From: <opal-request AT lists.psi.ch> on behalf of Wei Hou Tan <tanweihou AT outlook.com>
Reply-To: Wei Hou Tan <tanweihou AT outlook.com>
Date: Friday, January 22, 2021 at 1:59 PM
To: opal <opal AT lists.psi.ch>
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.

 

Regards

Wei Hou