function varargout = Emit_Parabola_Fit_kl_XProjection(varargin)
%global n  k sigx  sigy  gaux_sigma_sqrd   gaux_sigma_sqrd 
%sig_st = fopen('sig_st.txt', 'r');
%for n = 1:21,
%fprintf(sig_st, '%12d %12.8f %12.8f %12.8f %12.8f %12.8f \n', n, k(n), sigx(n), sigy(n), gaux_sigma_sqrd(n), gaux_sigma_sqrd(n));    
%end 
%fclose(sig_st);
   
%load -ascii HRRL_Emit_Parbola_Fit_Xaxis_as_BfieldStrength;
%Sig = HRRL_Emit_Parbola_Fit_Xaxis_as_BfieldStrength
load -ascii sigx.txt;
Sig = sigx
n = Sig(:,1);
B = Sig(:,2);
Er_B = Sig(:,3);
gaux_sigma_sqrd = Sig(:,4);
Er_sigmax_sqrd = Sig(:,5);

r_quad = 0.0254; %distance from center of dipole to the pole face is 1 inche
g = B/r_quad; %quadrupole gradiant
En = 10; % e- energy in Mev
me = 0.511; % electron rest mass in MeV
p = sqrt(En*En - me*me)/1000;% e- momentum in GeV/c
k1=0.2998*(g/p);%quadruople strength

% K1*L, L = 0.15 m
% K1*L, L = 0.165 m
% K1*L, L = 0.18 m
% K1*L, L = 0.20 m
% K1*L, L = 0.22 m

kl = 0.15*Sig(:,2); %k1*L
%kl = 0.165*Sig(:,2);
%kl = 0.18*Sig(:,2);
%kl = 0.20*Sig(:,2);
%kl = 0.22*Sig(:,2);

%S12=2.07500;%m thin lense aproximation.
S12=2;%m (2.075 - 0.15/2)
%S12=1.9925;%m (2.075 - 0.165/2)
%S12=1.985;%m (2.075 - 0.18/2)
%S12=1.975;%m (2.075 - 0.20/2)
%S12=1.965;%m (2.075 - 0.22/2)

% kl, S^2 fitting
px = polyfit(kl,gaux_sigma_sqrd,2);
    
ax = px(1);
bx = px(2);
cx = px(3);
   
fitx = ax*(kl.*kl) + bx*kl + cx;

pardatx = plot(kl,gaux_sigma_sqrd,'.');
hold on

parfitx = plot(kl,fitx,'r');
set(parfitx,'LineWidth',1);
legend('Data Points','X Projetion Fit');
hx = xlabel('QM Field: k_{1}*L (T)');
hy = ylabel('\sigma^2 (m^2)');
fit=title('k_{1}*L vs \sigma^2');
set(fit,'FontSize',12,'FontWeight','bold');
drawnow
   
gamma = En/me;

Ax=ax;
Bx=bx/(-2*ax);
Cx=cx-(bx)^2/(4*ax);
emit_x = sqrt(Ax*Cx)/((S12)^2);
emit_nx = emit_x*gamma
beta_x = sqrt(Ax/Cx)
alpha_x = sqrt(Ax/Cx)*(Bx + 1/S12)

    
fit_result = fopen('Parabole_fit_result_XProjection.txt', 'w');
   %fprintf(fit_result, ' x-projection a=%12.9f, b=%12.9f, c=%12.9f,A=%12.9f, B=%12.9f, C=%12.9f, emit=%12.15f, emit_normx=%12.15f \n', ax, bx, cx, Ax, Bx, Cx, emit_x, emit_nx ); 
    fprintf(fit_result, ' x-projection: emit=%12.15f, emit_normx=%12.15f \n', emit_x, emit_nx ); 
   fclose(fit_result);
%return