Skip to Content.
Sympa Menu

opal - Re: [Opal] Questions regarding .stat output and rfmap input

opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

Re: [Opal] Questions regarding .stat output and rfmap input


Chronological Thread 
  • From: "Gsell Achim (PSI)" <achim.gsell AT psi.ch>
  • To: Huiwen Koay <koay AT rcnp.osaka-u.ac.jp>
  • Cc: "Snuverink Jochem (PSI)" <jochem.snuverink AT psi.ch>, "Frey Matthias (PSI)" <matthias.frey AT psi.ch>, opal <opal AT lists.psi.ch>
  • Subject: Re: [Opal] Questions regarding .stat output and rfmap input
  • Date: Mon, 8 Jul 2019 13:53:33 +0000
  • Accept-language: en-US, de-CH
  • Authentication-results: localhost; iprev=pass (edge20.ethz.ch) smtp.remote-ip=82.130.99.26; spf=pass smtp.mailfrom=psi.ch; dmarc=skipped



> On 8 Jul 2019, at 10:57, Huiwen Koay <koay AT rcnp.osaka-u.ac.jp> wrote:
>
> Dear Jochem and Matthias,
>
> I finally managed get the estimated bunch size. Thanks for your help.
> Anyway, I am still looking for methods to convert Efield input from ascii
> to h5part.
> Do you have any suggestions?

Unfortunately we don’t have a general purpose program for this. But you can
adapt

https://gitlab.psi.ch/OPAL/src/blob/master/tools/BandRF/ascii2h5block_asgic.cpp

from the OPAL sources to your needs.

Best

Achim


>
> Thanks again.
>
> Best regards,
> Huiwen
>
>
> On Wed, Jun 26, 2019 at 10:29 PM Snuverink Jochem (PSI)
> <jochem.snuverink AT psi.ch> wrote:
> As Matthias pointed out to me, it should be:
>
> y - longitudinal
>
> Best regards,
> Jochem
>
>> On 26 Jun 2019, at 15:19, Snuverink Jochem (PSI) <jochem.snuverink AT psi.ch>
>> wrote:
>>
>> Dear Huiwen,
>>
>> Matthias and I can reproduce your finding.
>> The local frame (OPTION, PSDUMPFRAME=bunch_mean) in the stat file is
>> always:
>>
>> x - horizontal
>> y - radial
>> z - vertical
>>
>> I suppose your radial momentum (PRINIT attribute of the CYCLOTRON) is not
>> 0?
>>
>> In that case the local frame is rotated with respect to the initial
>> gaussian input distribution, so RMSx (in the rotated local frame) will no
>> longer be 0.
>>
>> Best regards,
>> Jochem
>>
>>> On 26 Jun 2019, at 04:16, Huiwen Koay <koay AT rcnp.osaka-u.ac.jp> wrote:
>>>
>>> Dear Matthias,
>>>
>>> Thank you for your prompt response.
>>> In your reply, RMSx = sqrt( <x^2> - <x>^2 ). So if I set the initial
>>> sigmax=0, I should have the initial RMSx=0.
>>> However, it is not so. As long as sigmay is not zero, RMSx will not be
>>> zero.
>>> It seems to me that
>>>
>>> RMSx = sqrt( <x^2> - <x>^2 ) + f(y)
>>>
>>> I am quite confused by this. Hope that you can help in clarifying this
>>> matter.
>>>
>>>
>>>
>>> Best regards,
>>> Huiwen
>>>
>>>
>>> On Wed, Jun 26, 2019 at 3:01 AM Frey Matthias (PSI)
>>> <matthias.frey AT psi.ch> wrote:
>>>
>>>
>>> Dear Huiwen,
>>>
>>> Regarding RMSx and RMSy. In OPAL-cycl, x is horizontal and y is
>>> longitudinal direction. Sure, when you change sigmax of the initial
>>> distribution you can have a different RMSy due to
>>> transversal-longitudinal couplings.
>>>
>>> RMS is the root-mean square, e.g. in x-direction it is computed by
>>>
>>> RMSx = sqrt( <x^2> - <x>^2 )
>>>
>>> where <x> is the centroid and <x^2> the second order moment.
>>>
>>>
>>> Best,
>>> Matthias
>>>
>>>
>>> Von: opal-request AT lists.psi.ch [opal-request AT lists.psi.ch]" im Auftrag
>>> von "Snuverink Jochem (PSI) [jochem.snuverink AT psi.ch]
>>> Gesendet: Dienstag, 25. Juni 2019 17:00
>>> An: Huiwen Koay
>>> Cc: opal
>>> Betreff: Re: [Opal] Questions regarding .stat output and rfmap input
>>>
>>> Dear Huiwen,
>>>
>>> Concerning questions about OPAL, it is best to use the email list as more
>>> users and members of the developer's team will see your message, and
>>> you'll get a quicker response.
>>>
>>> 1) @Matthias, can you answer that question perhaps?
>>>
>>> 2) I agree there is no reference in the Manual to the h5part format. That
>>> should be added and thanks for pointing that out. There is some
>>> information here: https://gitlab.psi.ch/H5hut/src/wikis/home
>>>
>>> @Achim, can you point Huiwen to how he can convert his ASCII field map to
>>> a h5part file?
>>>
>>> Best regards,
>>> Jochem
>>>
>>>> On 25 Jun 2019, at 15:41, Huiwen Koay <koay AT rcnp.osaka-u.ac.jp> wrote:
>>>>
>>>> Dear Jochem Snuverink,
>>>>
>>>> Hi. I have previously sent email to inquire about OPAL input for B
>>>> fieldmap. Currently I managed to simulate the single particle tracking
>>>> using opal-cycl. However, when I move onto deeper calculations, such as
>>>> bunch calculation and central region calculations, I face some troubles.
>>>> I did send inquiries to some of the members in the mailing lists
>>>> directly but I receive no reply. Thus, I would like to ask for help
>>>> again here. I sincerely appreciate if anyone could help in clearing my
>>>> confusion.
>>>>
>>>> 1) In manual 2.1, the coordinate system of .stat output in opal-t is
>>>> given. From my interpretation, the local frame in .stat file is similar
>>>> to that defined in opal-cycl. However, when I changed sigmax in
>>>> distribution input of type=GAUSS, both RMSx and RMSy changed. The RMSx
>>>> and RMSy in .stat file seems not to be consistent with the X and Y
>>>> defined in local frame of opal-cycl. So, in that case, I would like to
>>>> know the definition of RMSx and RMSy given in the .stat file.
>>>>
>>>> 2) I would like to use an eletric field map in OPAL for central region
>>>> calculations. However, in type=BANDRF, the input format for rfmap is in
>>>> .h5part. Regarding this, I would like to ask if you had any reference
>>>> format of .h5part file so that I can translate my ascii field map to
>>>> H5part?
>>>>
>>>> Thank you very much for your help.
>>>> I am looking forward to your soon reply.
>>>>
>>>>
>>>>
>>>> Best regards,
>>>> Huiwen
>>>
>>
>




Archive powered by MHonArc 2.6.19.

Top of Page