Touschek lifetime and +/- momentum aperture

Moderators: cyao, michael_borland

Post Reply
NicolaCarmignani
Posts: 8
Joined: 20 Dec 2012, 08:33

Touschek lifetime and +/- momentum aperture

Post by NicolaCarmignani » 04 Jan 2013, 10:07

At ESRF we are doing some lattice studies for our current and upgraded storage rings, using Elegant to compute the momentum aperture and resulting Touschek lifetime.
Elegant computes the positive and negative momentum aperture to be used in the lifetime calculation. The corresponding scattering rates should then be added to get the total scattering rate (1/tau = 1/tau_p + 1/tau_n).
However, looking at the touschekLifetime C code, it seems that the minimum of the momentum aperture is used.
We have cases with big enough asymmetry between delta_p and delta_n that this difference in evaluation of lifetime has significant impact.
Did we read the source code correctly?
Are we right in thinking the lifetime should be computed including the positive and negative terms separately?

Nicola Carmignani
Boaz Nash

xiaoam
Posts: 17
Joined: 12 Jun 2008, 09:27

Re: Touschek lifetime and +/- momentum aperture

Post by xiaoam » 04 Jan 2013, 19:16

Nicola,

You are right. We did assume symmetry (using the minimum value) momentum aperture in our code for reducing calculation work. Thanks for point out the need to treat negative and positive part separately. The code has been updated and should be realeased next time. Anybody need the code before the new realease date can send me an email.

Thank you very much,

aimin

hywel
Posts: 6
Joined: 07 Jan 2013, 08:21

Re: Touschek lifetime and +/- momentum aperture

Post by hywel » 07 Jan 2013, 09:47

That's really helpful. We're about to do some +/- Touschek calculations with asymmetric momentum aperture too, so this is a handy update when it comes.

michael_borland
Posts: 1797
Joined: 19 May 2008, 09:33
Location: Argonne National Laboratory
Contact:

Re: Touschek lifetime and +/- momentum aperture

Post by michael_borland » 07 Jan 2013, 11:22

Since it seems this is of some interest and urgency, we'll try to push out the update in the next week or so.

--Michael

xiaoam
Posts: 17
Joined: 12 Jun 2008, 09:27

Re: Touschek lifetime and +/- momentum aperture

Post by xiaoam » 07 Jan 2013, 11:39

A quick solution to this problem is as following:
Assume the input file twiss: a.twi, aperture: a.mmap
sddsprocess a.mmap pos.mmap -redef=col,deltaNegative,"deltaPositive -1 *"
sddsprocess a.mmap neg.mmap -redef=col,deltaPositive,"deltaNegative -1 *"
touschekLifetime neg.out -twiss=a.twi -aperture=neg.mmap ...
touschekLifetime pos.out -twiss=a.twi -aperture=pos.mmap ...
1/lifetime = 1/(2*lifetime_neg) + 1/(2*lifetime_pos)

This only works for command line calculation. For Monte Carlo simulation inside elegant, you have to wait for the new release.

Thanks,

aimin

hywel
Posts: 6
Joined: 07 Jan 2013, 08:21

Re: Touschek lifetime and +/- momentum aperture

Post by hywel » 07 Jan 2013, 11:58

That's a great fix, thanks!

Post Reply