linac errors

Moderators: cyao, michael_borland

Post Reply
astecpete
Posts: 35
Joined: 24 Jul 2008, 04:01
Location: Daresbury Laboratory, UK

linac errors

Post by astecpete » 14 Dec 2010, 09:54

Hi Michael,

We're doing a simple error analysis and have a little problem. Attched are a couple of files. We run twice, the first is a fiducial run with one particle, the second we apply an error to the voltage of cavities and track 1000 particles. The trouble comes when we make this error small. In the file i have made the error zero. When we do this and plot out Cdelta, we see that it still diverges from the centroid and ends up at ~6e-4 fractional difference. When we have large errors, we see that in Cdelta but as we turn it down, Cdelta asymptotes to this 6e-4 value. Looking in the .error file, the values for VOLT are correct (i.e. the values from the second run are the same as the first). This effect seems to be undiminished when i increase particle number so doesn't seem to be statistical.
Can you help!

Pete (and Sara from MaxLab)
Attachments
wakefieldlong.sdds
(154.57 KiB) Downloaded 853 times
max4_mb.lte
(10.95 KiB) Downloaded 1057 times
max4_mb.ele
(1.83 KiB) Downloaded 987 times

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

Re: linac errors

Post by michael_borland » 14 Dec 2010, 10:22

Pete and Sara,

This happens because of the first_is_fiducial=1 flag in &bunched_beam. This causes elegant to fiducialize with a single-particle bunch. That bunch apparently sees very different wakes than a realistic bunch, so it isn't a good fiducial. I should probably put a warning about this in the manual.

I recommend setting first_is_fiducial=0 in &bunched_beam. Also, you might want to consider setting one_random_bunch=0, so that the same particles are used for all runs. That way the effects you'll see will only be due to random errors assigned via &error_element.

--Michael

astecpete
Posts: 35
Joined: 24 Jul 2008, 04:01
Location: Daresbury Laboratory, UK

Re: linac errors

Post by astecpete » 14 Dec 2010, 11:32

Ah, i didn't think of turning the wakes off.
Thanks very much Michael

Post Reply