function [] = plotvlines(Calx, MeanData, Points, LW)
for i=1:length(Points)
    indx = find(Calx==Points(i),1);
    plot([Points(i) Points(i)], [0 1.3*max(MeanData(:,indx))], 'k--', 'LineWidth', LW)
    text(Points(i), 1.35*max(MeanData(:,indx)), num2str(Points(i)), 'Rotation', 90, 'FontSize', 20)
end
end