opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

[Opal] OPAL-T: Distribution SIGMAPX and SIGMAPY

From: "Dr. Simon Friederich" <s.friederich AT uni-mainz.de>
To: "opal AT lists.psi.ch" <opal AT lists.psi.ch>
Subject: [Opal] OPAL-T: Distribution SIGMAPX and SIGMAPY
Date: Wed, 24 Feb 2021 13:47:38 +0100
Authentication-results: localhost; iprev=pass (ironport-1.zdv.net) smtp.remote-ip=134.93.178.241; spf=pass smtp.mailfrom=uni-mainz.de; dkim=pass header.d=uni-mainz.de header.s=20190702 header.a=rsa-sha256; dmarc=pass header.from=uni-mainz.de
Ironport-phdr: 9a23:QX6q1xJmHfjfemApZ9mcpTZWNBhigK39O0sv0rFitYgXKv3yrarrMEGX3/hxlliBBdydt6sVzbqJ+Pi6EUU7or+5+EgYd5JNUxJXwe43pCcHRPC/NEvgMfTxZDY7FskRHHVs/nW8LFQHUJ2mPw6arXK99yMdFQviPgRpOOv1BpTSj8Oq3Oyu5pHfeQpFiCe5bL9oLxi7rwrdutUIjYB/Nqs/1xzFr2dSde9L321oP1WTnxj95se04pFu9jlbtuwi+cBdT6j0Zrw0QrNEAjsoNWA1/9DrugLYTQST/HscU34ZnQRODgPY8Rz1RJbxsi/9tupgxCmXOND9QL4oVTi+6apgVRnlgzoFOTEk6mHaktJ+gqJGrhKmpBJxzYDbb46JO/p/Za7deM8WSHBbU8pNWSFMAIWxZJYPAeobOuZYqpHwq1UAoxukHwWnGf7hyj5Nhn/x260xzuQvERvH3AM+At0Brm7UrNTpO6cSX+y+0a7FzS7ZY/NYwzf96ZbHfQ0/rvGWXLJ/bcrRxFIyFwPDj1WcsInlMCmU1uQLsmmX9fdvWvy2hmMhtgp+rSShyN02hYnVmoIa1ErE9SNhzYsrJdC1R012bMC4HJZeuCyWK4V7TMM/T2x1uys31rILt5C5cSUK1Zgqxh3SZvyJfoSU/h/uVOWcLDl3iX54fL+ygQu5/0u4yuDkWMm7zVlHojRAn9TDuH0ByRLe582ER/Rh40iuxDOC2gHP5uxFJE05kbDUJ4Ivz7M1kpccqkrDHiHzlUX4kqCaaF4r9+2z5Oj6f7nqu4eQOot0hw3kNqkih9azDfomPgQTRWSX5/iw2KPl8EbkWrtFlOc2nbPcsJ3CIMQbobO2DBFN34Y47ha/Ey+m0NMFnXkbNF5FeAyIj4zuO1zWPv/5AuqzjlqjnTt23vzGO6bhD5XNLnTajrvsfqxx6k9Bxwot19xQ6ZRUCrAYLPz0QEDxtdrYDx4jPwK6zOjrENR91oUAVmKTGqKZPr7dvFCW6u8qLOSBapEZtTj4JvUm+/HhkGI2lFoccKW13JYabXS1Ee56L0iXZHfhhswOEWkEvgc7QuzllkOOXiVWanmpQq4w+yw3BJqjDYjZR4CthLKB3D28Hp1Tfm1GEVGMEXnyd4WFRfgDdTydL9R8kjwCVLitUooh2guvtA/8yrtrNPfb9TMAuZL5yNd1//HTlQ019TFsCMSdy3uCT2Vvk2MMRj82x7x/oVdnxlebyqh4guZUGsBT56ABbgBvf5nExuhSD9HpRhmHZsyIDl2mCJ3yBik4Svo1wsQSeABlANTkhRCVjASwBLpA3ZCRQtQP86fb2TDfPcdwx3vc1aIhxRFyWY1BL2yriqNl9gPTL5PPkgCcnqesZ+IQ0Xiepy+40WOSsRQAA0ZLWqLfUCVHaw==
Ironport-sdr: 87NR7U32XxKsckqkL+H7MA66czGXfxNKGS0iD1rgJgGTeCvtADdIdaTey2iDAiNm5kAGUeA26m CtrjNbzCBTJI/fPoUaR12jkc44vs7Pvg5oi+DNdUsVRyHjEAY3gYwoQnR1SgSQZ3Ek7Km1N22y fcDHXnesU8QHWotMrU9HKjkozugENtC3FVSQsAcT38wtSZDedg+RAptxLrSJQMQcOvXNybroAU cBkO0XFjE6PFgNZbgANqGmQ1lJuP/RqzGrw0AKUmo7bbHGJkxdNfGxwIZ3yt4MenAeScK3f5ol nAQ=
Organization: Institut für Kernphysik Mainz

Dear Opal-Team,

I simulated only an injected GAUSS-type distribution with a non-zero SIGMAPX and SIGMAPY:
distr_gauss_injected: DISTRIBUTION,
    TYPE = GAUSS,
    SIGMAR = 1.0e-3, CUTOFFR = 3,
    SIGMAPX = 5.0e-3*beta*gamma, SIGMAPY=5.0e-3,
    SIGMAZ = 1e-6, CUTOFFLONG = 4,
    EMITTED = false;

I thought, I could use these paramters to have the distribution start with a transverse divergence of 5 mrad,
i.e. after Δs=1 m the rms size increases by Δσ=5mm (maybe this is where I am already mistaken).

I did not change INPUTMOUNITS (default = NONE) but neither
    SIGMAPX = 5.0e-3*beta*gamma // this I would expect from the explanation in the manual, right?
    nor
    SIGMAPY=5.0e-3
show the expected divergence (see screenshot in the attachement). At least SIGMAPX leads to Δσx=4mm though.

Beta and gamma are defined by:
REAL Ekin_GeV = 1e-4; // 100 keV
REAL gamma = (EMASS + Ekin_GeV)/EMASS;
REAL beta = sqrt(1-1/gamma^2);

Did I missinterpret SIGMAPX and SIGMAPY?

Best regards
Simon

-- 
Dr. Simon Friederich
Helmholtz-Institut Mainz | Insitut für Kernphysik
Kollaboration ACID       | B
Johannes Gutenberg-Universität Mainz
Staudingerweg 18         | Johann-Joachim-Becherweg 45
55128 Mainz, Deutschland

Tel.: +49 (0)6131 39-23160
E-Mail: s.friederich AT uni-mainz.de

Attachment: divergence.png
Description: PNG image

[Opal] OPAL-T: Distribution SIGMAPX and SIGMAPY, Dr. Simon Friederich, 02/24/2021
- Re: [Opal] OPAL-T: Distribution SIGMAPX and SIGMAPY [nosign], Snuverink Jochem (PSI), 02/24/2021