Difference between revisions of "HRRL 03-18-2011"
Line 311: | Line 311: | ||
|} | |} | ||
+ | |||
+ | [[File:TF_Q6_2.5_03182011.png| 200px]] | ||
== Lorentzian Fitting == | == Lorentzian Fitting == |
Revision as of 02:18, 9 April 2011
Tune parameters
Previous Tune at the higher energy
Solenoid 1 | 5.4 A |
Solenoid 2 | 5.5 A |
Gun Ver | -0.2 A |
Gun Hor | +0.4 A |
Output Hor | -0.5 A |
Output Ver | -0.5 A |
Gun HV | +9.75 (Knob Setting) |
Gun Grid Voltage | 5.25 (Knob Setting) |
RF frequency | 2855.813 MHz |
Modulator HV Power Supply | 4.11 (Knob Setting) |
RF macro Pulse Length (FWHM) | 200 ns |
Solenoid comparison with yesterdays higher energy tune
Quad Scans
To test reproducibility from yesterday. Scan at 40, 20, 10, 5 mA peak currents. Scan twice on each Current.
Do Energy Scan
Bend the beam through 45 bending dipole, take images on the scope of FC and YAG images.
D1 current | Energy | FC peak Volt | Pic |
(A) | MeV | mV | |
28 | 11.85 | 39 | |
29 | 12.3 | 74 | |
30 | 12.74 | 93 | |
30.5 | 12.96 | 65 | |
30.7 | large uncertainty | 22 |
Quad 1 Scan
Quad 1 | OTR Pict | FC pict |
(A) | ||
+0 | ||
+0.2 | ||
+0.4 | ||
+0.6 | ||
+0.8 | ||
+1.0 | ||
+1.2 | ||
+1.4 | ||
+1.6 | ||
+1.8 | ||
+2.0 |
Scan Quad 4
Calibration
Estimate the error in the above calibration
Data Analysis
Preliminary results
Here is the study on how to fit non_Gaussian curve. I fitted for whole image as well as parts of the image.
Here is beam profile, which is non-Gaussian.
Fits for whole image:
Fits for x=[362,404] and y=[241,301]:
Gaussian fitted for whole image
12.74 MeV
Fitted Gaussian for whole image x-projection:
emit=0.521 +- 0.013 mm*mrad, emit_norm=13.002 +- 0.330 mm*mrad
beta=0.884 +- 0.025, alpha=0.469 +-0.027
Gaussian fitted for part 1 of the image
Fitted Gaussian for part 1 of the image
xl = 362;% center at 382 xr = 404; yl = 241;% center at 271 yr = 301;
xl: x left, xr: x right. yl: y left, yr: y right.
x-projection:
emit=0.322 +- 0.001 mm*mrad, emit_norm=8.022 +- 0.033 mm*mrad
beta=3.110 +- 0.003, alpha=0.906 +-0.002
Gaussian fitted for part 2 of the image
Fitted Gaussian for part 2 of the image
xl = 342;% center at 382 xr = 424; yl = 241;% center at 271 yr = 301;
x-projection:
emit=0.360 +- 0.001 mm*mrad, emit_norm=8.970 +- 0.031 mm*mrad
beta=1.993 +- 0.001, alpha=0.672 +-0.001
Gaussian fitted for part 3 of the image
Fitted Gaussian for part 3 of the image
xl = 342;% center at 382 xr = 424; yl = 221;% center at 271 yr = 321;
x-projection:
emit=0.360 +- 0.001 mm*mrad, emit_norm=8.970 +- 0.031 mm*mrad
beta=1.993 +- 0.001, alpha=0.672 +-0.001
Gaussian fitted for part 4 of the image
Fitted Gaussian for part 4 of the image
xl = 322;% center at 382 xr = 444; yl = 221;% center at 271 yr = 321;
x-projection:
emit=0.350 +- 0.001 mm*mrad, emit_norm=8.738 +- 0.032 mm*mrad
beta=1.344 +- 0.002, alpha=0.464 +-0.002
Analysis with root
Q1_Scan, 42mA peak current, Scan Coil Current at positive 1.8_Amp
Back grounds
Q6_Scan, 42mA peak current, Scan Coil Current at negative 2.5_Amp
Lorentzian Fitting
Basic MATLAB Codes to fit Lorentzian:
File:Trial my Lorentzian fit.txt
My MATLAB Fit code for beam
Here is an example fitting:
Signal | Background | fit |
height of the peak | location of the peak | half-width at half-maximum (HWHM) |
1.6935e+006 | 358.9910 | 45.4280 |
If we were to apply the relation of the Gaussian rms to its FWHM to Lorentzian, we will get the sigma (or rms) of Lorentzian to be:
Super Gaussian Fitting
Beam Distributions Beyond RMS: File:Beam Distributions Beyond RMS.pdf
Basic MATLAB Codes to fit Super Gaussian:
File:Trial my Super Gaussian fit.txt
My MATLAB Fit code for beam
File:My Super Gaussian fit.txt
Here is an example fitting:
Signal | Background | fit |
| ||
base | Amplitude | center |
202 | 8.0480e+003 | 359.6730 |
sigma_0 | N | sigma |
20.92 | 0.8494 | 38.56 |
If we were to apply the relation of the Gaussian rms to its FWHM to Lorentzian, we will get the sigma (or rms) of Lorentzian to be:
This is very close to the sigma we got from Super Gaussian, which is 38.56
Results
Q1
Positive Current, X projection
emit=0.370 +- 0.001 mm*mrad, emit_norm=9.225 +- 0.030 mm*mrad
beta=1.329 +- 0.001, alpha=0.754 +-0.001