[4] | 1 | function [CorFamRes,CorDevRes,CorElemRes,MaxEffRes,BpmFamRes,BpmDevRes]= mec(QuadFamily,QuadDev,Plane,CorBpmResp,CorNumb) |
---|
| 2 | % MEC |
---|
| 3 | % function [CorFam,CorDev,MaxEff]= mec(QuadFamily,QuadDev,Plane,CorBpmResp) |
---|
| 4 | % |
---|
| 5 | % INPUTS -> VALUES |
---|
| 6 | % |
---|
| 7 | % 1. QuadFamily: Family of the quadrupole studied -> 'Q2' |
---|
| 8 | % 2. QuadDev: [cell elt] of the quadrupole studied -> [1 1] |
---|
| 9 | % 3. Plane: Plane=1 horizontal only |
---|
| 10 | % =2 vertical only |
---|
| 11 | % 4. CorBpmResp: Response matrix of correctors obtained from "respcor(plane)" |
---|
| 12 | % |
---|
| 13 | % |
---|
| 14 | % OUPUTS |
---|
| 15 | % |
---|
| 16 | % 1. CorFamRes: Family of the Most Effective Corrector found |
---|
| 17 | % 2. CorDevRes: [cell elt] of the MEC found |
---|
| 18 | % 3. MaxEffRes: Efficacity of the MEC in mm/A |
---|
| 19 | % |
---|
| 20 | |
---|
| 21 | if ~iscellstr(QuadFamily) |
---|
| 22 | error('la famille entree n est pas au format cell') |
---|
| 23 | end |
---|
| 24 | |
---|
| 25 | tailleQDev=size(QuadDev(1,:)); |
---|
| 26 | tailleCorBpm=size(CorBpmResp); |
---|
| 27 | |
---|
| 28 | if isempty(QuadDev) |
---|
| 29 | error('Veuillez entrer le device du quadripole à étudier au format QuadDev=[cell elt]') |
---|
| 30 | elseif tailleQDev(2)~=2 |
---|
| 31 | error('Trop d''elements pour 1 device de QP: QuadDev=[cell elt]') |
---|
| 32 | end |
---|
| 33 | |
---|
| 34 | if Plane<0 || Plane>2 |
---|
| 35 | error('Veuillez entrer le plan étudié: Horizontal: Plane=1 ; Vertical:Plane=2 ; Both:Plane=0') |
---|
| 36 | end |
---|
| 37 | |
---|
| 38 | if tailleCorBpm(1)~=56 & tailleCorBpm(2)~=120 |
---|
| 39 | error('Veuillez vérifier la matrice réponse correcteurs/BPM') |
---|
| 40 | end |
---|
| 41 | |
---|
| 42 | |
---|
| 43 | |
---|
| 44 | |
---|
| 45 | |
---|
| 46 | load('respcor.mat') |
---|
| 47 | |
---|
| 48 | if Plane==1 |
---|
| 49 | BpmFam='BPMx'; |
---|
| 50 | CorFam='HCOR'; |
---|
| 51 | %CorBpmResp=efficacy.HPLane; %effic(Cor,Bpm) |
---|
| 52 | elseif Plane==2 |
---|
| 53 | BpmFam='BPMz'; |
---|
| 54 | CorFam='VCOR'; |
---|
| 55 | %CorBpmResp=efficacy.VPlane; |
---|
| 56 | else |
---|
| 57 | error('input must "1" for horizontal plane or "2" for vertical plane'); |
---|
| 58 | end |
---|
| 59 | |
---|
| 60 | |
---|
| 61 | |
---|
| 62 | %Most Effective Corrector |
---|
| 63 | |
---|
| 64 | %BPM le plus proche du Qp (QuadFamily,[QuadDev]); |
---|
| 65 | [BpmFam2,BpmDev,BpmPos]=proche2(BpmFam,QuadFamily,QuadDev,0); |
---|
| 66 | |
---|
| 67 | indBpm=dev2elem(BpmFam,BpmDev); |
---|
| 68 | % Faire un tri sur EffCor de maniÚre ᅵ pouvoir choisir par exemple le 2Úme correcteurs le plus |
---|
| 69 | % efficace si le premier ne convient pas... |
---|
| 70 | EffCor=CorBpmResp(:,indBpm); % vecteur colonne de longueur nb HCOR ou nb VCOR |
---|
| 71 | [posEff,indEff]=sort(abs(EffCor),'descend'); |
---|
| 72 | |
---|
| 73 | %MaxEff=max(abs(EffCor)); % extraction de l'efficacité max |
---|
| 74 | MaxEff=posEff(CorNumb); |
---|
| 75 | %indCor=find(abs(EffCor)==MaxEff); % recherche de l'indice de l'efficacité max dans "EffCor" |
---|
| 76 | indCor=indEff(CorNumb); |
---|
| 77 | CorDevList=getlist(CorFam); |
---|
| 78 | |
---|
| 79 | %Rï¿œsultats |
---|
| 80 | CorFamRes=CorFam; |
---|
| 81 | CorDevRes=CorDevList(indCor,:); |
---|
| 82 | CorElemRes=indCor; % Recherche du correcteur dᅵlivrant l'efficacitᅵ prᅵcᅵdente |
---|
| 83 | MaxEffRes=MaxEff; |
---|
| 84 | BpmFamRes=BpmFam; |
---|
| 85 | BpmDevRes=BpmDev; |
---|
| 86 | |
---|
| 87 | |
---|
| 88 | |
---|
| 89 | |
---|
| 90 | %************************************************************************** |
---|
| 91 | |
---|
| 92 | |
---|