The OPAL Discussion Forum

["Adelmann Andreas (PSI)" <andreas.adelmann AT psi.ch>]

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

https://gitlab.psi.ch/OPAL/src/issues/285

@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