Gaussian fits with sddsemitproc
Posted: 26 Jan 2011, 20:11
Hello,
I have a beam which has a mostly gaussian distribution but some tails. Although there is <10% of the beam in the tails, the sigmas (Sx, Sy) are correspondingly much larger than the values I get from a gaussian fit. To make a more realistic comparison with how I will actually measure the emittance using a quad scan, is there an intelligent way to use sddsemitproc with values from gaussian fits to a watch file?
I'll likely write something that steps the quad, dumps the watch file, fits the distributions, iterates then write the all the values to some input file 'quadScan.data' that I can feed into sddsemitproc. This just seems so much more complicated than the simple/quick combination of vary_element and sddsprocemit, so I was wondering if anyone else has already scripted this.
Update: I just realized (must have forgotten this) that the watch files generated by vary element have several pages of data (one for each step) so getting the fits should be much easier. I'll have a go and post if I'm successful.
Thanks,
Mark
I have a beam which has a mostly gaussian distribution but some tails. Although there is <10% of the beam in the tails, the sigmas (Sx, Sy) are correspondingly much larger than the values I get from a gaussian fit. To make a more realistic comparison with how I will actually measure the emittance using a quad scan, is there an intelligent way to use sddsemitproc with values from gaussian fits to a watch file?
I'll likely write something that steps the quad, dumps the watch file, fits the distributions, iterates then write the all the values to some input file 'quadScan.data' that I can feed into sddsemitproc. This just seems so much more complicated than the simple/quick combination of vary_element and sddsprocemit, so I was wondering if anyone else has already scripted this.
Update: I just realized (must have forgotten this) that the watch files generated by vary element have several pages of data (one for each step) so getting the fits should be much easier. I'll have a go and post if I'm successful.
Thanks,
Mark