Closed orbit for off-momentum particles
Posted: 09 Jun 2020, 03:21
Dear all,
I'm having troubles getting the closed_orbit command to work with off-momentum particles.
What I'm trying to do is calculate orbit and Twiss for an off-momentum particle. There is one way I can get what I want, by inserting a MALIGN element and using its DP item to shift the momentum. This is demonstrated by the script dporbit_malign.ele in the attached archive.
However, I actually wanted to use a different approach, shifting the momentum using the CENTROID item of the bunched_beam command and then setting start_from_dp_centroid=1 in the closed_orbit command to set up the calculation of the off-momentum closed orbit. However, when I do this, I always end up with a zero orbit, no matter what I try. My failed attempt is demonstrated by the script dporbit.ele in the attached archive.
In both cases, I created bunch and output file to check the validity of my input parameters for the momentum shift. I set p_central to 1 and intend to apply a 5e-3 momentum shift. From the output files, which both indicate the momentum of the single particle I'm using as 1.005, those input parameters seem to be okay. Also, in the second case, I really do see an off-momentum track in the centroid file, which does, however, not coincide with the off-momentum orbit since the computed orbit appears to be zero.
I'm pretty sure I'm just missing something here, so any help on how I could get the off-momentum orbit to work using the centroid shift (as suggested by the documentation of the closed_orbit command) would be highly appreciated.
Thanks in advance and best regards
David
I'm having troubles getting the closed_orbit command to work with off-momentum particles.
What I'm trying to do is calculate orbit and Twiss for an off-momentum particle. There is one way I can get what I want, by inserting a MALIGN element and using its DP item to shift the momentum. This is demonstrated by the script dporbit_malign.ele in the attached archive.
However, I actually wanted to use a different approach, shifting the momentum using the CENTROID item of the bunched_beam command and then setting start_from_dp_centroid=1 in the closed_orbit command to set up the calculation of the off-momentum closed orbit. However, when I do this, I always end up with a zero orbit, no matter what I try. My failed attempt is demonstrated by the script dporbit.ele in the attached archive.
In both cases, I created bunch and output file to check the validity of my input parameters for the momentum shift. I set p_central to 1 and intend to apply a 5e-3 momentum shift. From the output files, which both indicate the momentum of the single particle I'm using as 1.005, those input parameters seem to be okay. Also, in the second case, I really do see an off-momentum track in the centroid file, which does, however, not coincide with the off-momentum orbit since the computed orbit appears to be zero.
I'm pretty sure I'm just missing something here, so any help on how I could get the off-momentum orbit to work using the centroid shift (as suggested by the documentation of the closed_orbit command) would be highly appreciated.
Thanks in advance and best regards
David