Difference between revisions of "HRRL 03-18-2011"

From New IAC Wiki
Jump to navigation Jump to search
 
(108 intermediate revisions by 2 users not shown)
Line 1: Line 1:
−
= Scan Qaud 1 =
+
=Tune parameters=
 +
 
 +
Previous Tune at the higher energy
 +
{| border="1"
 +
| 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
 +
 
 +
{| border="3"  cellpadding="20" cellspacing="0"
 +
|V1 (Gun Vertical)||H1 (Gun Horizontal)|| Sol1 (Gun Sol) || V2 (Output Vertical) || H2 (Output horizontal || Sol 2 (Output solenoid)|| OTR Pict || FC pict
 +
|-
 +
|(A)||(A) || (A) ||(A)||(A) || (A) || ||
 +
|-
 +
| (-0.3)||+0 || 0 ||0||-0 || -0 ||[[File:OTR_Q1Scan03182011_SterOff.png | 100 px]] ||[[File:FC_Q1Scan03182011_SterOff.png | 100 px]]
 +
|-
 +
| (-0.2)||+0.4 || 5.4 ||-0.5||-0.5 || 5.5 ||[[File:OTR_Q1Scan03182011_SterOn.png | 100 px]] ||[[File:FC_Q1Scan03182011_SterOon.png | 100 px]]
 +
|-
 +
| (-0.2)||+0.4 || 6.8 ||-0.5||-0.5 || 5.5 ||[[File:OTR_Q1Scan03182011_1.png | 100 px]] ||[[File:FC_Q1Scan03182011_1.png | 100 px]]
 +
|-
 +
| (-0.2)||+0.4 || 6.8 ||-0.5||-0.5 || 10.4 ||[[File:OTR_Q1Scan03182011_2.png | 100 px]] ||[[File:FC_Q1Scan03182011_2.png | 100 px]]
 +
|-
 +
|}
 +
 
 +
= Quad Scans =
 
To test reproducibility from yesterday. Scan at 40, 20, 10, 5 mA peak currents. Scan twice on each Current.
 
To test reproducibility from yesterday. Scan at 40, 20, 10, 5 mA peak currents. Scan twice on each Current.
−
= Do Energy Scan =
+
== Do Energy Scan ==
 
Bend the beam through 45 bending dipole, take images on the scope of FC and YAG images.
 
Bend the beam through 45 bending dipole, take images on the scope of FC and YAG images.
 +
 +
{| border="3"  cellpadding="20" cellspacing="0"
 +
|D1 current||Energy || FC peak Volt || Pic
 +
|-
 +
|(A)||MeV || mV ||
 +
|-
 +
|28 ||11.85 ||39 ||[[File:FC2_0318011_28A_11.85.png| 200 px]]
 +
|-
 +
|29 ||12.3 ||74 ||[[File:FC2_0318011_28A_12.3.png| 200 px]]
 +
|-
 +
|30 || 12.74 ||93 ||[[File:FC2_0318011_30A_12.74.png| 200 px]]
 +
|-
 +
|30.5 || 12.96 ||65 ||[[File:FC2_0318011_30.5A_12.96.png| 200 px]]
 +
|-
 +
|30.7 ||large uncertainty ||22 ||[[File:FC2_0318011_30.7A.png| 200 px]]
 +
|-
 +
|}
 +
 +
== Quad 1 Scan==
 +
 +
 +
{| border="3"  cellpadding="20" cellspacing="0"
 +
|Quad 1|| OTR Pict || FC pict
 +
|-
 +
|(A) || ||
 +
|-
 +
| +0||[[File:OTR_Q1Scan03182011_10.png | 100 px]] ||[[File:FC_Q1Scan03182011_10.png | 100 px]]
 +
|-
 +
| +0.2||[[File:OTR_Q1Scan03182011_11.png | 100 px]] ||[[File:FC_Q1Scan03182011_11.png | 100 px]]
 +
|-
 +
| +0.4||[[File:OTR_Q1Scan03182011_12.png | 100 px]] ||[[File:FC_Q1Scan03182011_12.png | 100 px]]
 +
|-
 +
| +0.6||[[File:OTR_Q1Scan03182011_13.png | 100 px]] ||[[File:FC_Q1Scan03182011_13.png | 100 px]]
 +
|-
 +
| +0.8||[[File:OTR_Q1Scan03182011_14.png | 100 px]] ||[[File:FC_Q1Scan03182011_14.png | 100 px]]
 +
|-
 +
| +1.0||[[File:OTR_Q1Scan03182011_15.png | 100 px]] ||[[File:FC_Q1Scan03182011_15.png | 100 px]]
 +
|-
 +
| +1.2||[[File:OTR_Q1Scan03182011_16.png | 100 px]] ||[[File:FC_Q1Scan03182011_16.png | 100 px]]
 +
|-
 +
| +1.4||[[File:OTR_Q1Scan03182011_17.png | 100 px]] ||[[File:FC_Q1Scan03182011_17.png | 100 px]]
 +
|-
 +
| +1.6||[[File:OTR_Q1Scan03182011_18.png | 100 px]] ||[[File:FC_Q1Scan03182011_18.png | 100 px]]
 +
|-
 +
| +1.8 ||[[File:OTR_Q1Scan03182011_19.png | 100 px]] ||[[File:FC_Q1Scan03182011_19.png | 100 px]]
 +
|-
 +
| +2.0 ||[[File:OTR_Q1Scan03182011_20.png | 100 px]] ||[[File:FC_Q1Scan03182011_20.png | 100 px]]
 +
|}
 +
 
= Scan Quad 4 =
 
= Scan Quad 4 =
 +
 +
= Calibration =
 +
 +
 +
 
 +
 +
Estimate the error in the above calibration
 +
 +
 +
 +
{| border="3"  cellpadding="20" cellspacing="0"
 +
| image || OTR image || Calibration
 +
|-
 +
|[[File:HRRL_03-18-2011_cal_1_source.png | 200 px]]  || [[File:HRRL_03-18-2011_cal_1.png | 200 px]]  ||cal_h = 31.75 mm / 737 px = 0.04308 mm/px; cal_v = 31.75 mm / sqrt(2) / 531 px = 0.04228 mm/px;
 +
|-
 +
|[[File:HRRL_03-18-2011_cal_2_source.png | 200 px]]  || [[File:HRRL_03-18-2011_cal_2.png | 200 px]]  || cal_h = 31.75 mm / 738 px = 0.04302 mm/px; cal_v = 31.75 mm / sqrt(2) / 535 px = 0.04196 mm/px;
 +
|-
 +
|[[File:HRRL_03-18-2011_cal_3_source.png | 200 px]]  || [[File:HRRL_03-18-2011_cal_3.png | 200 px]]  || cal_h = 31.75 mm / 735 px = 0.04320 mm/px; cal_v = 31.75 mm / sqrt(2) / 532 px = 0.04220 mm/px;
 +
|-
 +
|[[File:HRRL_03-18-2011_cal_4_source.png | 200 px]]  || [[File:HRRL_03-18-2011_cal_4.png | 200 px]]  || cal_h = 31.75 mm / 735 px = 0.04320 mm/px; cal_v = 31.75 mm / sqrt(2) / 536 px = 0.04189 mm/px;
 +
|-
 +
|[[File:HRRL_03-18-2011_cal_5_source.png | 200 px]]  || [[File:HRRL_03-18-2011_cal_5.png | 200 px]]  || cal_h = 31.75 mm / 749 px = 0.04239 mm/px; cal_v = 31.75 mm / sqrt(2)) / 539 px = 0.04165  mm/px
 +
|-
 +
|[[File:HRRL_03-18-2011_cal_6_source.png | 200 px]]  || [[File:HRRL_03-18-2011_cal_6.png | 200 px]]  || cal_h = 31.75 mm / 732 px = 0.04337 mm/px; cal_v = 31.75 mm / sqrt(2) / 536 px = 0.04189 mm/px
 +
|-
 +
|[[File:HRRL_03-18-2011_cal_7_source.png | 200 px]]  || [[File:HRRL_03-18-2011_cal_7.png | 200 px]]  || cal_h = 31.75 mm / 735 px = 0.04320 mm/px; cal_v = 31.75 mm / sqrt(2)) / 532 px = 0.04220 mm/px
 +
|-
 +
|}
 +
 +
 +
 +
{| border="3"  cellpadding="20" cellspacing="0"
 +
|- 
 +
|cal_h (mm/px) || cal_v (mm/px) ||
 +
|- 
 +
|0.04308|| 0.04228 ||
 +
|-
 +
|0.04302 || 0.04196 ||
 +
|- 
 +
|0.04320  || 0.04220 ||
 +
|- 
 +
|0.04320 || 0.04189 ||
 +
|- 
 +
|0.04239 ||  0.04165 ||
 +
|-
 +
|0.04337 || 0.04189 ||
 +
|-
 +
|0.04320 || 0.04220 ||
 +
|-
 +
|  mean (mm/px) || mean (mm/px) ||
 +
|-
 +
|  <math> 0.04307 \pm 0.00032 </math>||  <math> 0.04201 \pm  0.00023 </math> ||
 +
|}
 +
 +
= 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:
 +
 +
[[File:HRRL_03-18-2011_Data_Analysis_Preliminary_Results_non-Gauss_profile_1.png | 300 px]]
 +
 +
Fits for x=[362,404] and y=[241,301]:
 +
 +
[[File:HRRL_03-18-2011_Data_Analysis_Preliminary_Results_non-Gauss_profile_2.png | 300 px]]
 +
 +
 +
=== Gaussian fitted for whole image ===
 +
 +
 +
12.74 MeV
 +
 +
Fitted Gaussian for whole image
 +
x-projection:
 +
 +
[[File:HRRL_03-18-2011_Data_Analysis_Preliminary_Results_Gauss_Fitted_Whole_Image.png | 300 px]]
 +
 +
=== 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:
 +
 +
 +
[[File:HRRL_03-18-2011_Data_Analysis_Preliminary_Results_Gauss_Fitted_Part_1_Image.png | 300 px]]
 +
 +
=== 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:
 +
 +
 +
[[File:HRRL_03-18-2011_Data_Analysis_Preliminary_Results_Gauss_Fitted_Part_2_Image.png | 300 px]]
 +
 +
 +
=== 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:
 +
 +
 +
[[File:HRRL_03-18-2011_Data_Analysis_Preliminary_Results_Gauss_Fitted_Part_3_Image.png | 300 px]]
 +
 +
 +
 +
=== 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:
 +
 +
 +
[[File:HRRL_03-18-2011_Data_Analysis_Preliminary_Results_Gauss_Fitted_Part_4_Image.png | 300 px]]
 +
 +
 +
 +
 +
== Analysis with root ==
 +
 +
===  Q1_Scan, 42mA peak current, Scan Coil Current at positive 1.8_Amp ===
 +
 +
Image
 +
[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-signal.png | 200 px]]
 +
 +
 +
Back grounds
 +
 +
[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-bg-1.png | 100 px]]
 +
[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-bg-2.png | 100 px]]
 +
[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-bg-3.png | 100 px]]
 +
[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-bg-4.png | 100 px]]
 +
[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-bg-5.png | 100 px]]
 +
[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-bg-6.png | 100 px]]
 +
[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-bg-7.png | 100 px]]
 +
[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-bg-8.png | 100 px]]
 +
[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-bg-9.png | 100 px]]
 +
[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-bg-10.png | 100 px]]
 +
 +
 +
 +
{| border="1"
 +
| Signal || Background
 +
|-
 +
| [[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-signal.png | 300 px]]    || [[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-bg-1.png | 300 px]]
 +
|-
 +
|These MATLAB fits are set to fit at 0.7 times above height of the peak.
 +
|-
 +
|  Matlab fit to x profile || Matlab fit to y profile
 +
|-
 +
|[[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-signal-Matlab_Fit-x.png | 300 px]]  ||  [[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-signal-Matlab_Fit-y.png | 300 px]]
 +
|-
 +
|sigma_x = 24.97 px || sigma_y = 44.45 px
 +
|-
 +
| Blue: Signal, Red Background, Black: Subtraction || [[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-signal_bg_sg-bg.png | 300 px]]
 +
|-
 +
|signal || [[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-surf-sg.png | 300 px]]
 +
|
 +
|-
 +
|background || [[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-surf-bg.png | 300 px]]
 +
|
 +
|-
 +
| signal - background|| [[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-surf-sg-bg.png | 300 px]]
 +
|
 +
|-
 +
|signal || [[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-contour-sg.png | 300 px]]
 +
|-
 +
|background || [[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-contour-bg.png | 300 px]]
 +
|
 +
|-
 +
| signal - background|| [[File:HRRL-03-18-2011-Analysis-with-root-Q1_Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8_Amp-contour-sg-bg.png | 300 px]]
 +
|
 +
 +
|}
 +
 +
[[File:TF_Q1_1.8_03182011_NOCUTS.png| 200 px]][[File:TF_Q1_1.8_03182011.png| 200 px]]
 +
 +
The red line is the background above which seems to be larger than the signal.  Why?
 +
 +
Below I zoom in on channel 250 where you can see that for this channel the background is larger.  I then plot, on the right, all the elements along the x direction what are added up to get the values on the right.
 +
 +
[[File:TF_Q1_1.8_03182011Problem.png| 200 px]]
 +
 +
Why isn't the addition making sense?
 +
 +
===  Q6_Scan, 42mA peak current, Scan Coil Current at negative 2.5_Amp ===
 +
 +
{| border="1"
 +
| Signal || Background
 +
|-
 +
| [[File:HRRL-03-18-2011-Analysis-with-root-Q6_Scan-42mA-peak-current-Scan-Coil-Current-at-negative-2.5_Amp-signal.png | 300 px]]    || [[File:HRRL-03-18-2011-Analysis-with-root-Q6_Scan-42mA-peak-current-Scan-Coil-Current-at-negative-2.5_Amp-bg-1.png | 300 px]]
 +
|-
 +
|  Matlab fit to x profile || Matlab fit to y profile
 +
|-
 +
| set to fit at 0.2 times above height of the peak. || set to fit at 0.7 times above height of the peak.
 +
|-
 +
 +
|[[File:HRRL-03-18-2011-Analysis-with-root-Q6_Scan-42mA-peak-current-Scan-Coil-Current-at-negative-2.5_Amp-signal-Matlab_Fit-x.png | 300 px]]  ||  [[File:HRRL-03-18-2011-Analysis-with-root-Q6_Scan-42mA-peak-current-Scan-Coil-Current-at-negative-2.5_Amp-signal-Matlab_Fit-y.png | 300 px]]
 +
|-
 +
|sigma_x = 31.37 px || sigma_y = 18.63 px
 +
 +
|}
 +
 +
[[File:TF_Q6_2.5_03182011.png| 200px]]
 +
 +
== Lorentzian Fitting ==
 +
 +
Basic MATLAB Codes to fit Lorentzian:
 +
 +
[[File:Trial_my_Lorentzian_fit.txt]]
 +
 +
[[File:desum.txt]]
 +
 +
My MATLAB Fit code for beam
 +
 +
[[File:my_Lorentzian_fit.txt]]
 +
 +
[[File:desum.txt]]
 +
 +
Here is an example fitting:
 +
 +
{| border="1"
 +
| Signal || Background || fit
 +
|-
 +
| [[File:HRRL_03-18-2011_Lorentzian_Fitting_Example_sg.png| 300 px]] || [[File:HRRL_03-18-2011_Lorentzian_Fitting_Example_bg.png| 300 px]]  || [[File:HRRL_03-18-2011_Lorentzian_Fitting_Example_1.png| 300 px]]
 +
|-
 +
| 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:
 +
<math> \frac {2 \times 45.4280}{ 2.3548} = 38.5833192 </math>
 +
 +
 +
== Super Gaussian Fitting ==
 +
 +
Beam Distributions Beyond RMS:
 +
[[File:Beam_Distributions_Beyond_RMS.pdf]]
 +
 +
Bib: http://www.osti.gov/bridge/product.biblio.jsp?osti_id=82539#
 +
 +
 +
Basic MATLAB Codes to fit Super Gaussian:
 +
 +
[[File:Trial_my_Super_Gaussian_fit.txt]]
 +
 +
[[File:SupGau_desum.txt]]
 +
 +
My MATLAB Fit code for beam
 +
 +
[[File:my_Super_Gaussian_fit.txt]]
 +
 +
[[File:SupGau_desum.txt]]
 +
 +
Here is an example fitting:
 +
 +
{| border="1"
 +
| Signal || Background || fit
 +
|-
 +
| [[File:HRRL_03-18-2011_Lorentzian_Fitting_Example_sg.png| 300 px]] || [[File:HRRL_03-18-2011_Lorentzian_Fitting_Example_bg.png| 300 px]]  || [[File:HRRL_03-18-2011_SupGau_Fitting_Example_1.png| 300 px]]
 +
 +
 +
 +
[[HRRL_Emittance_Measurements_March14-18-2011]]
 +
 +
|-
 +
| 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:
 +
<math> \frac {2 \times 45.4280}{ 2.3548} = 38.5833192 </math>
 +
 +
This is very close to the sigma we got from Super Gaussian, which is 38.56
 +
 +
= Results =
 +
 +
== Q1 ==
 +
 +
=== Positive Current, X projection ===
 +
 +
==== 1st fit ====
 +
 +
 +
[[File:HRRL 03-18-2011_Results_Q1_42mA_Pos_X.png | 300 px]]
 +
 +
 +
==== 2nd fit after cut off some points ====
 +
 +
 +
[[File:HRRL 03-18-2011_Results_Q1_42mA_Pos_X_2.png | 300 px]]
 +
 +
== Q4 ==
 +
 +
Did not reach minimum.
 +
 +
 +
 +
  
  
 
[[HRRL_Emittance_Measurements_March14-18-2011]]
 
[[HRRL_Emittance_Measurements_March14-18-2011]]

Latest revision as of 19:13, 11 January 2012

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

V1 (Gun Vertical) H1 (Gun Horizontal) Sol1 (Gun Sol) V2 (Output Vertical) H2 (Output horizontal Sol 2 (Output solenoid) OTR Pict FC pict
(A) (A) (A) (A) (A) (A)
(-0.3) +0 0 0 -0 -0 OTR Q1Scan03182011 SterOff.png FC Q1Scan03182011 SterOff.png
(-0.2) +0.4 5.4 -0.5 -0.5 5.5 OTR Q1Scan03182011 SterOn.png FC Q1Scan03182011 SterOon.png
(-0.2) +0.4 6.8 -0.5 -0.5 5.5 OTR Q1Scan03182011 1.png FC Q1Scan03182011 1.png
(-0.2) +0.4 6.8 -0.5 -0.5 10.4 OTR Q1Scan03182011 2.png FC Q1Scan03182011 2.png

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 FC2 0318011 28A 11.85.png
29 12.3 74 FC2 0318011 28A 12.3.png
30 12.74 93 FC2 0318011 30A 12.74.png
30.5 12.96 65 FC2 0318011 30.5A 12.96.png
30.7 large uncertainty 22 FC2 0318011 30.7A.png

Quad 1 Scan

Quad 1 OTR Pict FC pict
(A)
+0 OTR Q1Scan03182011 10.png FC Q1Scan03182011 10.png
+0.2 OTR Q1Scan03182011 11.png FC Q1Scan03182011 11.png
+0.4 OTR Q1Scan03182011 12.png FC Q1Scan03182011 12.png
+0.6 OTR Q1Scan03182011 13.png FC Q1Scan03182011 13.png
+0.8 OTR Q1Scan03182011 14.png FC Q1Scan03182011 14.png
+1.0 OTR Q1Scan03182011 15.png FC Q1Scan03182011 15.png
+1.2 OTR Q1Scan03182011 16.png FC Q1Scan03182011 16.png
+1.4 OTR Q1Scan03182011 17.png FC Q1Scan03182011 17.png
+1.6 OTR Q1Scan03182011 18.png FC Q1Scan03182011 18.png
+1.8 OTR Q1Scan03182011 19.png FC Q1Scan03182011 19.png
+2.0 OTR Q1Scan03182011 20.png FC Q1Scan03182011 20.png

Scan Quad 4

Calibration

Estimate the error in the above calibration


image OTR image Calibration
HRRL 03-18-2011 cal 1 source.png HRRL 03-18-2011 cal 1.png cal_h = 31.75 mm / 737 px = 0.04308 mm/px; cal_v = 31.75 mm / sqrt(2) / 531 px = 0.04228 mm/px;
HRRL 03-18-2011 cal 2 source.png HRRL 03-18-2011 cal 2.png cal_h = 31.75 mm / 738 px = 0.04302 mm/px; cal_v = 31.75 mm / sqrt(2) / 535 px = 0.04196 mm/px;
HRRL 03-18-2011 cal 3 source.png HRRL 03-18-2011 cal 3.png cal_h = 31.75 mm / 735 px = 0.04320 mm/px; cal_v = 31.75 mm / sqrt(2) / 532 px = 0.04220 mm/px;
HRRL 03-18-2011 cal 4 source.png HRRL 03-18-2011 cal 4.png cal_h = 31.75 mm / 735 px = 0.04320 mm/px; cal_v = 31.75 mm / sqrt(2) / 536 px = 0.04189 mm/px;
HRRL 03-18-2011 cal 5 source.png HRRL 03-18-2011 cal 5.png cal_h = 31.75 mm / 749 px = 0.04239 mm/px; cal_v = 31.75 mm / sqrt(2)) / 539 px = 0.04165 mm/px
HRRL 03-18-2011 cal 6 source.png HRRL 03-18-2011 cal 6.png cal_h = 31.75 mm / 732 px = 0.04337 mm/px; cal_v = 31.75 mm / sqrt(2) / 536 px = 0.04189 mm/px
HRRL 03-18-2011 cal 7 source.png HRRL 03-18-2011 cal 7.png cal_h = 31.75 mm / 735 px = 0.04320 mm/px; cal_v = 31.75 mm / sqrt(2)) / 532 px = 0.04220 mm/px


cal_h (mm/px) cal_v (mm/px)
0.04308 0.04228
0.04302 0.04196
0.04320 0.04220
0.04320 0.04189
0.04239 0.04165
0.04337 0.04189
0.04320 0.04220
mean (mm/px) mean (mm/px)
[math] 0.04307 \pm 0.00032 [/math] [math] 0.04201 \pm 0.00023 [/math]

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:

HRRL 03-18-2011 Data Analysis Preliminary Results non-Gauss profile 1.png

Fits for x=[362,404] and y=[241,301]:

HRRL 03-18-2011 Data Analysis Preliminary Results non-Gauss profile 2.png


Gaussian fitted for whole image

12.74 MeV

Fitted Gaussian for whole image x-projection:

HRRL 03-18-2011 Data Analysis Preliminary Results Gauss Fitted Whole Image.png

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:


HRRL 03-18-2011 Data Analysis Preliminary Results Gauss Fitted Part 1 Image.png

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:


HRRL 03-18-2011 Data Analysis Preliminary Results Gauss Fitted Part 2 Image.png


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:


HRRL 03-18-2011 Data Analysis Preliminary Results Gauss Fitted Part 3 Image.png


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:


HRRL 03-18-2011 Data Analysis Preliminary Results Gauss Fitted Part 4 Image.png



Analysis with root

Q1_Scan, 42mA peak current, Scan Coil Current at positive 1.8_Amp

Image HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-signal.png


Back grounds

HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-bg-1.png HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-bg-2.png HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-bg-3.png HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-bg-4.png HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-bg-5.png HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-bg-6.png HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-bg-7.png HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-bg-8.png HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-bg-9.png HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-bg-10.png


Signal Background
HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-signal.png HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-bg-1.png
These MATLAB fits are set to fit at 0.7 times above height of the peak.
Matlab fit to x profile Matlab fit to y profile
HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-signal-Matlab Fit-x.png HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-signal-Matlab Fit-y.png
sigma_x = 24.97 px sigma_y = 44.45 px
Blue: Signal, Red Background, Black: Subtraction HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-signal bg sg-bg.png
signal HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-surf-sg.png
background HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-surf-bg.png
signal - background HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-surf-sg-bg.png
signal HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-contour-sg.png
background HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-contour-bg.png
signal - background HRRL-03-18-2011-Analysis-with-root-Q1 Scan-42mA-peak-current-Scan-Coil-Current-at-positive-1.8 Amp-contour-sg-bg.png

TF Q1 1.8 03182011 NOCUTS.pngTF Q1 1.8 03182011.png

The red line is the background above which seems to be larger than the signal. Why?

Below I zoom in on channel 250 where you can see that for this channel the background is larger. I then plot, on the right, all the elements along the x direction what are added up to get the values on the right.

TF Q1 1.8 03182011Problem.png

Why isn't the addition making sense?

Q6_Scan, 42mA peak current, Scan Coil Current at negative 2.5_Amp

Signal Background
HRRL-03-18-2011-Analysis-with-root-Q6 Scan-42mA-peak-current-Scan-Coil-Current-at-negative-2.5 Amp-signal.png HRRL-03-18-2011-Analysis-with-root-Q6 Scan-42mA-peak-current-Scan-Coil-Current-at-negative-2.5 Amp-bg-1.png
Matlab fit to x profile Matlab fit to y profile
set to fit at 0.2 times above height of the peak. set to fit at 0.7 times above height of the peak.
HRRL-03-18-2011-Analysis-with-root-Q6 Scan-42mA-peak-current-Scan-Coil-Current-at-negative-2.5 Amp-signal-Matlab Fit-x.png HRRL-03-18-2011-Analysis-with-root-Q6 Scan-42mA-peak-current-Scan-Coil-Current-at-negative-2.5 Amp-signal-Matlab Fit-y.png
sigma_x = 31.37 px sigma_y = 18.63 px

TF Q6 2.5 03182011.png

Lorentzian Fitting

Basic MATLAB Codes to fit Lorentzian:

File:Trial my Lorentzian fit.txt

File:Desum.txt

My MATLAB Fit code for beam

File:My Lorentzian fit.txt

File:Desum.txt

Here is an example fitting:

Signal Background fit
HRRL 03-18-2011 Lorentzian Fitting Example sg.png HRRL 03-18-2011 Lorentzian Fitting Example bg.png HRRL 03-18-2011 Lorentzian Fitting Example 1.png
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: [math] \frac {2 \times 45.4280}{ 2.3548} = 38.5833192 [/math]


Super Gaussian Fitting

Beam Distributions Beyond RMS: File:Beam Distributions Beyond RMS.pdf

Bib: http://www.osti.gov/bridge/product.biblio.jsp?osti_id=82539#


Basic MATLAB Codes to fit Super Gaussian:

File:Trial my Super Gaussian fit.txt

File:SupGau desum.txt

My MATLAB Fit code for beam

File:My Super Gaussian fit.txt

File:SupGau desum.txt

Here is an example fitting:

Signal Background fit
HRRL 03-18-2011 Lorentzian Fitting Example sg.png HRRL 03-18-2011 Lorentzian Fitting Example bg.png HRRL 03-18-2011 SupGau Fitting Example 1.png


HRRL_Emittance_Measurements_March14-18-2011

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: [math] \frac {2 \times 45.4280}{ 2.3548} = 38.5833192 [/math]

This is very close to the sigma we got from Super Gaussian, which is 38.56

Results

Q1

Positive Current, X projection

1st fit

HRRL 03-18-2011 Results Q1 42mA Pos X.png


2nd fit after cut off some points

HRRL 03-18-2011 Results Q1 42mA Pos X 2.png

Q4

Did not reach minimum.




HRRL_Emittance_Measurements_March14-18-2011

Retrieved from "https://wiki.iac.isu.edu/index.php?title=HRRL_03-18-2011&oldid=70454"