Difference between revisions of "HRRL 03-18-2011"
(104 intermediate revisions by 2 users not shown) | |||
Line 1: | Line 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" | {| 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 || | + | |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
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
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) | |
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:
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:
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:
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:
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:
Analysis with root
Q1_Scan, 42mA peak current, Scan Coil Current at positive 1.8_Amp
Back grounds
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.
Why isn't the addition making sense?
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
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
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
1st fit
2nd fit after cut off some points
Q4
Did not reach minimum.