function[] = lineprofile()
pix = 25*1000/1024; %nm
L = 1.5;
L_c = [0 0 0];
T1 = readmatrix("Plot Values_LEEM.csv");
T2 = readmatrix("Plot Values_pp.csv");
T3 = readmatrix("Plot Values_ip80.csv");
T4 = readmatrix("Plot Values_ip170.csv");
Len_x = T1(:,1)*pix;
L1 = T1(:,2);
L2 = T2(:,2);
L3 = T3(:,2);
L4 = T4(:,2);
% stripe position
x1 = 85*25/15;
x2 = 150*25/15;
x3 = 250*25/15;
x4 = 330*25/15;
x5 = 545*25/15;
x6 = 615*25/15;
x7 = 703*25/15;
x8 = 825*25/15;
x9 = 960*25/15;
x10 = 1030*25/15;
h1 = figure(1);
hold on
box on
plot([x1 x1],[min(L1)*0.5 max(L1)*1.1],'--','linewidth',L,'color',L_c)
plot([x2 x2],[min(L1)*0.5 max(L1)*1.1],'--','linewidth',L,'color',L_c)
plot([x3 x3],[min(L1)*0.5 max(L1)*1.1],'--','linewidth',L,'color',L_c)
plot([x4 x4],[min(L1)*0.5 max(L1)*1.1],'--','linewidth',L,'color',L_c)
plot([x5 x5],[min(L1)*0.5 max(L1)*