Skip to Content.
Sympa Menu

opal - Re: [Opal] Inquiries regarding the sigma matrix of GaussMatched distribution

opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

Re: [Opal] Inquiries regarding the sigma matrix of GaussMatched distribution


Chronological Thread 
  • From: "Adelmann Andreas (PSI)" <andreas.adelmann AT psi.ch>
  • To: "Frey Matthias (PSI)" <matthias.frey AT psi.ch>
  • Cc: "Adelmann Andreas (PSI)" <andreas.adelmann AT psi.ch>, Huiwen Koay <koay AT rcnp.osaka-u.ac.jp>, opal <opal AT lists.psi.ch>
  • Subject: Re: [Opal] Inquiries regarding the sigma matrix of GaussMatched distribution
  • Date: Sun, 29 Sep 2019 18:39:19 +0000
  • Accept-language: en-US, de-CH
  • Authentication-results: localhost; iprev=pass (edge10.ethz.ch) smtp.remote-ip=82.130.75.186; spf=pass smtp.mailfrom=psi.ch; dmarc=skipped

Hi Matthias Huiwen will send us the input deck and I will have a first look at this
together with Jochem.

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
Phone Home: xx41 62 891 91 44
-------------------------------------------------------
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 25 Sep 2019, at 12:56, Frey Matthias (PSI) <matthias.frey AT psi.ch> wrote:

Hi Andreas

Unfortunately not really. This is however since I'm a bit short on time for other things.

Best,
Matthias


Von: Adelmann Andreas (PSI)
Gesendet: Mittwoch, 25. September 2019 12:48
An: Frey Matthias (PSI)
Cc: Huiwen Koay; opal
Betreff: Re: [Opal] Inquiries regarding the sigma matrix of GaussMatched distribution

Hi Matthias, I just talked to Huiwen about the matched distribution problem. 
Is there any progress on that?

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
Phone Home: xx41 62 891 91 44
-------------------------------------------------------
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 4 Sep 2019, at 16:31, Adelmann Andreas (PSI) <andreas.adelmann AT psi.ch> wrote:

When was the last commit of Chris?

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
Phone Home: xx41 62 891 91 44
-------------------------------------------------------
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 4 Sep 2019, at 14:11, Frey Matthias (PSI) <matthias.frey AT psi.ch> wrote:

Hi Huiwen

We have currently an open issue on the MatchedGauss feature. Please check also


@cortes_c seems to work on that.

Best,
Matthias


Von: opal-request AT lists.psi.ch [opal-request AT lists.psi.ch]" im Auftrag von "Huiwen Koay [koay AT rcnp.osaka-u.ac.jp]
Gesendet: Mittwoch, 4. September 2019 11:40
An: opal
Betreff: [Opal] Inquiries regarding the sigma matrix of GaussMatched distribution

Hi.
I have a question regarding the sigma-matrix shown in the output file of OPAL-cycl execution.
In the manual, it is stated that the sigma-matrix computation is based on (x,x′,y,y′,z,δ) instead of (x,px,y,py,z,pz). However, when I try to convert the x' to the sigmapx as shown in the output of MatchedGauss distribution, I couldnt reproduce the same value. 

For example, the output of OPAL shows:
OPAL> * I= 2.2 (mA)  E= 580 (MeV)
OPAL> * EX= 1e-06* EY= 1e-06* ET= 1e-06* FMAPFN ring590_bfld.dat
OPAL> * FMSYM= 8* HN= 6* PHIINIT= 0
OPAL> * ----------------------------------------------------
OPAL> * Converged (Ex, Ey, Ez) = (0.313794, 0.313794, 0.313794) pi mm mrad for E= 580 (MeV)
OPAL> * Sigma-Matrix
OPAL>       3.914 &     -0.6486 &           0 &           0 &      0.4542 &       1.362 \\
OPAL>     -0.6486 &      0.6396 &           0 &           0 &      0.7265 &     -0.2685 \\
OPAL>           0 &           0 &       6.239 &       1.198 &           0 &           0 \\
OPAL>           0 &           0 &       1.198 &      0.3859 &           0 &           0 \\
OPAL>      0.4542 &      0.7265 &           0 &           0 &       2.363 &      0.1752 \\
OPAL>       1.362 &     -0.2685 &           0 &           0 &      0.1752 &      0.8944 \\

OPAL> * ************* D I S T R I B U T I O N ********************************************
OPAL> *
OPAL> * Distribution type: MATCHEDGAUSS
OPAL> * SIGMAX   = 0.001978 [m]
OPAL> * SIGMAY   = 0.002498 [m]
OPAL> * SIGMAZ   = 0.001537 [m]
OPAL> * SIGMAPX  = 4.129e-05 [Beta Gamma]
OPAL> * SIGMAPY  = 3.639e-05 [Beta Gamma]
OPAL> * SIGMAPZ  = 4.49e-05 [Beta Gamma]
OPAL> * AVRGPZ   = 0 [Beta Gamma]
OPAL> * CORRX    = -0.4099
OPAL> * CORRY    = 0.7722
OPAL> * CORRZ    = 0.1205
OPAL> * R61      = 0.7277
OPAL> * R62      = -0.355
OPAL> * R51      = 0.1493
OPAL> * R52      = 0.5909
OPAL> * CUTOFFX  = 3 [units of SIGMAX]
OPAL> * CUTOFFY  = 3 [units of SIGMAY]
OPAL> * CUTOFFZ  = 3 [units of SIGMAZ]
OPAL> * CUTOFFPX = 3 [units of SIGMAPX]
OPAL> * CUTOFFPY = 3 [units of SIGMAPY]
OPAL> * CUTOFFPZ = 3 [units of SIGMAPY]
As the definition of rms-emittance is given by emittance=sqrt(sigma_11*sigma_22-sigma12*sigma21). However, if I used only the red coloured matrix, I can neither reproduce the matched emittance at 0.3137 pi mm mrad, nor the rms-px by taking sqrt(639)/1000*betagamma. I can only obtain the rms of X by taking the sqrt of 3.914.
I would appreciate if somebody could help and shed a light on this matter.
Thanks again.

Best regards,
Huiwen




Archive powered by MHonArc 2.6.19.

Top of Page