[17] | 1 | function [mCHE, mCVE, mCHS, mCVS] = idCalcFeedForwardCorTables(idName, idParams, fileNamesMeasCOD, dirMeasCOD, dirMeasCorResp, arIndsBPMsToSkip) |
---|
| 2 | %Calculate Feed-Forward Correction Tables vs Undulator Parameters for 2 Horizontal and 2 Vertical correctors. |
---|
| 3 | %NOTE: Horizontal correctors are those displacing electron orbit in Horizontal direction (by means of Vertical magnetic field) |
---|
| 4 | |
---|
| 5 | mCHE = []; mCVE = []; mCHS = []; mCVS = []; |
---|
| 6 | |
---|
| 7 | numUndParams = length(idParams); |
---|
| 8 | if((numUndParams < 0) || (numUndParams > 2)) |
---|
| 9 | sprintf('Number of ID parameters for the correction tables should be 1 or 2'); |
---|
| 10 | return; |
---|
| 11 | end |
---|
| 12 | |
---|
| 13 | numValsUndParam1 = length(idParams{1}{2}); |
---|
| 14 | undParamName1 = idParams{1}{1}; |
---|
| 15 | |
---|
| 16 | numValsUndParam2 = 1; |
---|
| 17 | undParamName2 = ''; |
---|
| 18 | if(numUndParams > 1) |
---|
| 19 | numValsUndParam2 = length(idParams{2}{2}); |
---|
| 20 | undParamName2 = idParams{2}{1}; |
---|
| 21 | end |
---|
| 22 | |
---|
| 23 | if((numValsUndParam2 <= 1) && (strcmp(undParamName2, '') ~= 0)) |
---|
| 24 | numValsUndParam2 = 1; |
---|
| 25 | undParamName2 = 'dummy'; |
---|
| 26 | end |
---|
| 27 | |
---|
| 28 | mCHE = zeros(numValsUndParam2, numValsUndParam1); |
---|
| 29 | mCVE = zeros(numValsUndParam2, numValsUndParam1); |
---|
| 30 | mCHS = zeros(numValsUndParam2, numValsUndParam1); |
---|
| 31 | mCVS = zeros(numValsUndParam2, numValsUndParam1); |
---|
| 32 | |
---|
| 33 | mCorRespX = zeros(120, 2); |
---|
| 34 | mCorRespZ = zeros(120, 2); |
---|
| 35 | |
---|
| 36 | corEfficiencyDependsOnUndParam = 1; |
---|
| 37 | if((strcmp(undParamName1, 'gap') == 0) && (strcmp(undParamName2, 'gap') == 0)) |
---|
| 38 | [mCorRespX, mCorRespZ] = idCalcCorRespMatr(idName, 0, dirMeasCorResp); |
---|
| 39 | corEfficiencyDependsOnUndParam = 0; |
---|
| 40 | elseif((strcmp(undParamName1, 'gap') ~= 0) && (strcmp(undParamName2, '') ~= 0)) |
---|
| 41 | [mCorRespX, mCorRespZ] = idCalcCorRespMatr(idName, 0, dirMeasCorResp); |
---|
| 42 | mCorRespXgl = zeros(120, 2); |
---|
| 43 | mCorRespZgl = zeros(120, 2); |
---|
| 44 | [mCorRespXgl, mCorRespZgl] = idCalcCorRespMatr(idName, 250, dirMeasCorResp); |
---|
| 45 | corEfficiencyDependsOnUndParam = 0; |
---|
| 46 | for i = 1:length(mCorRespX) |
---|
| 47 | if((mCorRespXgl(i,1) ~= mCorRespX(i,1)) || (mCorRespXgl(i,2) ~= mCorRespX(i,2))) |
---|
| 48 | corEfficiencyDependsOnUndParam = 1; |
---|
| 49 | break; |
---|
| 50 | end |
---|
| 51 | if((mCorRespZgl(i,1) ~= mCorRespZ(i,1)) || (mCorRespZgl(i,2) ~= mCorRespZ(i,2))) |
---|
| 52 | corEfficiencyDependsOnUndParam = 1; |
---|
| 53 | break; |
---|
| 54 | end |
---|
| 55 | end |
---|
| 56 | end |
---|
| 57 | |
---|
| 58 | indMeasCOD = 1; |
---|
| 59 | for i = 1:numValsUndParam1 |
---|
| 60 | undParam1 = idParams{1}{2}(i); |
---|
| 61 | if(corEfficiencyDependsOnUndParam) |
---|
| 62 | if(strcmp(undParamName1, 'gap')) |
---|
| 63 | [mCorRespX, mCorRespZ] = idCalcCorRespMatr(idName, undParam1, dirMeasCorResp); |
---|
| 64 | end |
---|
| 65 | end |
---|
| 66 | |
---|
| 67 | for j = 1:numValsUndParam2 |
---|
| 68 | undParam2 = 0; |
---|
| 69 | if(numUndParams > 1) |
---|
| 70 | undParam2 = idParams{2}{2}(j); |
---|
| 71 | end |
---|
| 72 | if(corEfficiencyDependsOnUndParam) |
---|
| 73 | if(strcmp(undParamName2, 'gap')) |
---|
| 74 | [mCorRespX, mCorRespZ] = idCalcCorRespMatr(idName, undParam2, dirMeasCorResp); |
---|
| 75 | end |
---|
| 76 | end |
---|
| 77 | |
---|
| 78 | fileNameMeasMain = char(fileNamesMeasCOD{indMeasCOD}{1}); |
---|
| 79 | fileNameMeasBkg = char(fileNamesMeasCOD{indMeasCOD}{2}); |
---|
| 80 | |
---|
| 81 | dirStart = pwd; |
---|
| 82 | if strcmp(dirMeasCOD, '') |
---|
| 83 | dirMeasCOD = getfamilydata('Directory',idName); |
---|
| 84 | end |
---|
| 85 | cd(dirMeasCOD); |
---|
| 86 | stMeasMain = load(fileNameMeasMain); |
---|
| 87 | stMeasBkg = load(fileNameMeasBkg); |
---|
| 88 | cd(dirStart); |
---|
| 89 | |
---|
| 90 | vCODx = stMeasMain.X - stMeasBkg.X; |
---|
| 91 | curCorX = idLeastSqLinFit(mCorRespX, vCODx, arIndsBPMsToSkip); |
---|
| 92 | |
---|
| 93 | vCODz = stMeasMain.Z - stMeasBkg.Z; |
---|
| 94 | curCorZ = idLeastSqLinFit(mCorRespZ, vCODz, arIndsBPMsToSkip); |
---|
| 95 | |
---|
| 96 | mCHE(j, i) = -curCorX(1); |
---|
| 97 | mCVE(j, i) = -curCorZ(1); |
---|
| 98 | mCHS(j, i) = -curCorX(2); |
---|
| 99 | mCVS(j, i) = -curCorZ(2); |
---|
| 100 | indMeasCOD = indMeasCOD + 1; |
---|
| 101 | end |
---|
| 102 | end |
---|
| 103 | |
---|
| 104 | if(numUndParams > 1) |
---|
| 105 | fprintf('VERTICAL dimension in tables (i.e. row no. is variable; column no. is const) corresponds to variation of %s\n', undParamName2); |
---|
| 106 | fprintf('HORIZONTAL dimension in tables (i.e. row no. is const; column no. is variable) corresponds to variation of %s\n', undParamName1); |
---|
| 107 | end |
---|