Commit da650abf authored by Simone Liuzzo's avatar Simone Liuzzo

pinhole D09 instead of D01

parent 177a704b
......@@ -4,7 +4,7 @@ function qemb = qemoptics(qemres,qemb,qemb0)
rph = ebs.model.addpinholes(qemb.at);
bpmidx = find(atgetcells(rph,'Class','Monitor'))';
nbpm=length(bpmidx);
[idx,j,kdx]=unique([bpmidx qemres.id07idx qemres.d17idx qemres.id25idx qemres.d27idx qemres.d01idx]); %#ok<ASGLU>
[idx,j,kdx]=unique([bpmidx qemres.id07idx qemres.d09idx qemres.d17idx qemres.id25idx qemres.d27idx]); %#ok<ASGLU>
orbit0=findsyncorbit(rph,qemres.ct,bpmidx);
dpp=orbit0(5,1);
......@@ -14,10 +14,10 @@ lindata=lindata(kdx);
qemb.lindata=lindata(1:nbpm); % BPM results
qemb.id07data=lindata(nbpm+1); % ID07 results
qemb.d17data=lindata(nbpm+2); % D17 results
qemb.id25data=lindata(nbpm+3); % ID25 results
qemb.d27data=lindata(nbpm+4); % D27 results
qemb.d01data=lindata(nbpm+5); % D01 results
qemb.d09data=lindata(nbpm+2); % D09 results
qemb.d17data=lindata(nbpm+3); % D17 results
qemb.id25data=lindata(nbpm+4); % ID25 results
qemb.d27data=lindata(nbpm+5); % D27 results
dispersion=cat(2,qemb.lindata.Dispersion)';
qemb.tunes=qemb.pm.fulltunes;
qemb.beta=cat(1,qemb.lindata.beta);
......@@ -49,7 +49,7 @@ bpmidx = qemres.bpmidx;
qemres.brot,qemres.bhgain,qemres.bvgain,orbit0);
if nargin >= 3
lindata0=[qemb0.lindata qemb0.id07data qemb0.d17data qemb0.id25data qemb0.d27data qemb0.d01data];
lindata0=[qemb0.lindata qemb0.id07data qemb0.d09data qemb0.d17data qemb0.id25data qemb0.d27data];
beta0=cat(1,lindata0.beta);
beam66=cat(3,lindata.beam66);
sigma2=[squeeze(beam66(1,1,:)) squeeze(beam66(3,3,:))];
......
......@@ -237,7 +237,7 @@ qemres.skewl=atgetfieldvalues(qemres.at(qemres.skewidx),'Length');
% pinhole data from the model with pinholes
qemres.id07idx=atmodpinholes.get(0,'pinholeID07')';
qemres.id25idx=atmodpinholes.get(0,'pinholeID25')';
qemres.d01idx=atmodpinholes.get(0,'pinholeD01')';
qemres.d09idx=atmodpinholes.get(0,'pinholeD09')';
qemres.d17idx=atmodpinholes.get(0,'pinholeD17')';
qemres.d27idx=atmodpinholes.get(0,'pinholeD27')';
......
......@@ -7,17 +7,17 @@ alllist=1:nbpm; % all BPMS
% ebs.pinholenames(1:5);, the LATTICE order starting frommjnection (cell
% 04)
id07=nbpm+1; % ID07
d17=nbpm+2; % D17
id25=nbpm+3; % ID25
d27=nbpm+4; % D27
d01=nbpm+5; % D01 pinhole
pinholelist =[id07 d17 id25 d27 d01] ; % this order should be identical to the one in qemoptics.m
pinholename={'ID07';'D17';'ID25';'D27';'D01'};
d09=nbpm+2; % D09 pinhole
d17=nbpm+3; % D17
id25=nbpm+4; % ID25
d27=nbpm+5; % D27
pinholelist =[id07 d09 d17 id25 d27 ] ; % this order should be identical to the one in qemoptics.m
pinholename={'ID07';'D09';'D17';'ID25';'D27'};
emittances=qemb.emittances;
sbpm=cat(1,qemb.lindata.SPos);
sskew=findspos(qemb.at(:),qemres.skewidx);
spinhole=[qemb.id07data.SPos; qemb.d17data.SPos; qemb.id25data.SPos; qemb.d27data.SPos;qemb.d01data.SPos];
spinhole=[qemb.id07data.SPos; qemb.d09data.SPos; qemb.d17data.SPos; qemb.id25data.SPos; qemb.d27data.SPos];
mess={...
sprintf('em. H [nm]: %7.3f %7.3f %7.3f',1.e9*mean(emittances(pinholelist,1:2:end)));...
......@@ -25,10 +25,10 @@ mess={...
sprintf('V. dispersion [m]: %g',qemb.pm.alpha*std(qemb.frespz,1));...
'pinhole V [pm]:';...
sprintf('ID07:%8.4f',1.e12*emittances(id07,6));...
sprintf(' D09:%8.4f',1.e12*emittances( d09,6));...
sprintf(' D17:%8.4f',1.e12*emittances( d17,6));...
sprintf('ID25:%8.4f',1.e12*emittances(id25,6));...
sprintf(' D27:%8.4f',1.e12*emittances( d27,6));...
sprintf(' D01:%8.4f',1.e12*emittances( d01,6))};
sprintf(' D27:%8.4f',1.e12*emittances( d27,6))};
figure(1); % emittances
plot(sbpm,1.e12*emittances(alllist,2:2:end));
......
No preview for this file type
......@@ -6,12 +6,12 @@ sbpm=cat(1,qemb.lindata.SPos);
sskew=findspos(qemb.at(:),qemres.skewidx);
figure(fig);
switch plotcode
case 1 % D01 spot
semplotemit(1.0e12*qemb.d01data.beam66([1 3],[1 3]))
case 1 % ID07 spot
semplotemit(1.0e12*qemb.id07data.beam66([1 3],[1 3]))
xlabel('x [\mum]');
ylabel('z [\mum]');
case 2 % ID07 spot
semplotemit(1.0e12*qemb.id07data.beam66([1 3],[1 3]))
case 2 % D09 spot
semplotemit(1.0e12*qemb.d09data.beam66([1 3],[1 3]))
xlabel('x [\mum]');
ylabel('z [\mum]');
case 3 % D17 spot
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment