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: Huiwen Koay <koay AT rcnp.osaka-u.ac.jp>
  • To: "Gsell Achim (PSI)" <achim.gsell AT psi.ch>
  • 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: Fri, 19 Jul 2019 10:52:39 +0900

Dear Achim,

Oh that sounds great! 
Will check it soon.
Thank you very much!.
 
Best regards,
Huiwen


On Thu, Jul 18, 2019 at 10:39 PM Gsell Achim (PSI) <achim.gsell AT psi.ch> wrote:
Dear Huiwen,

there is a reviewed ascii2h5block in the master branch of the OPAL repository. Ippl is not longer required.

Best

Achim



> On 17 Jul 2019, at 06:50, Huiwen Koay <koay AT rcnp.osaka-u.ac.jp> wrote:
>
> Dear all,
>
> Hi, is there any recommendtion to access to Ippl.h library from the binary package of OPAL so that I can adapt the ascii2h5part.cpp?
>
> Thanks again for your precious time.
>
> Best regards,
> Huiwen
>
>
> On Wed, Jul 10, 2019 at 11:54 AM Huiwen Koay <koay AT rcnp.osaka-u.ac.jp> wrote:
> Dear Achim,
>
> Thank you for your reply.
> I checked through the link you gave. It could be a solution.
> However, as I am using binary package currently, I do not had have Ippl.h in the directory (or it is hidden in some directories which I don't know?)
> Do you have any suggestion to obtain this ?
>
> Thanks again.

> Best regards,
> Huiwen
>
>
> On Mon, Jul 8, 2019 at 10:55 PM Gsell Achim (PSI) <achim.gsell AT psi.ch> wrote:
>
>
> > 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