0001 clear all
0002 close all
0003 addpath('/mntdirect/_machfs/liuzzo/CODE/LatticeTuningFunctions/errors')
0004 
0005 
0006 load ../../ESRFLattice.mat
0007 
0008 
0009 indm=indBPM;
0010 indq=find(atgetcells(ring,'Class','Quadrupole')); 
0011 
0012 rerr=atsetrandomerrors(...
0013     ring,...
0014     indq,...
0015     indm,...
0016     1,...
0017     1e-1,...
0018     2.5,...
0019     'dpb2');
0020 
0021 figure('units','normalized','position',[0.1 0.4 0.65 0.35])
0022 atplot(rerr,[0,100],'comment',[],@plPolynomBComp);
0023 saveas(gca,'FieldIntegral.fig')
0024 export_fig('FieldIntegral.jpg','-r300')
0025