/ATF2/FlightSim/coreApps/changeP.m

http://atf2flightsim.googlecode.com/ · MATLAB · 94 lines · 73 code · 6 blank · 15 comment · 14 complexity · cd6d1067026bee46f53151940c482b2d MD5 · raw file 

    

  1. function stat=changeP(dp,df)

    2. 

  2. %

    3. 

  3. % stat=changeP(0.01); % P=P0*(1+0.01)

    4. 

  4. %

    5. 

  5. % Change beam momentum by dp*(design momentum) for dispersion measurement

    6. 

  6. % (change relative to beam first time this function is run; reset reference P

    7. 

  7. % for this by issuing "clear changeP" command)

    8. 

  8. %

    9. 

  9. % stat=changeP(~,'on');  % turn DR RF-ramp ON

    10. 

  10. % stat=changeP(~,-2);    % set DR RF-ramp dF= -2 kHz

    11. 

  11. % stat=changeP(~,'off'); % turn DR RF-ramp OFF

    12. 

  12. 

    13. 

  13. 

    14. 

  14. global FL BEAMLINE

    15. 

  15. persistent initX

    16. 

  16. 

    17. 

  17. debug=false;

    18. 

  18. if (debug)

    19. 

  19.   fprintf(1,'NOTE: "debug" flag set in changeP\n');

    20. 

  20. end

    21. 

  21. 

    22. 

  22. stat{1}=1;

    23. 

  23. if (~isnumeric(dp))

    24. 

  24.   stat{1}=-1;

    25. 

  25.   stat{2}='changeP: dp arg must be numeric';

    26. 

  26.   return

    27. 

  27. end

    28. 

  28. twoargs=exist('df','var');

    29. 

  29. 

    30. 

  30. if (isequal(FL.mode,'trusted')&&(~FL.isthread)) % server

    31. 

  31.   if (FL.SimMode)

    32. 

  32.     % save initial beam

    33. 

  33.     if (isempty(initX)||isempty(initX{FL.simBeamID}))

    34. 

  34.       initX{FL.simBeamID}=FL.SimBeam{FL.simBeamID}.Bunch.x;

    35. 

  35.     end

    36. 

  36.     if (twoargs&&ischar(df))

    37. 

  37.       return % 'on' and 'off' are no-ops in simulation mode

    38. 

  38.     end

    39. 

  39.     % for simulation case, change input beam momentum (dp required)

    40. 

  40.     if (debug)

    41. 

  41.       if (twoargs)

    42. 

  42.         fprintf(1,'changeP: dp=%s , df=%s\n',num2str(dp),num2str(df));

    43. 

  43.       else

    44. 

  44.         fprintf(1,'changeP: dp=%s\n',num2str(dp));

    45. 

  45.       end

    46. 

  46.     else

    47. 

  47.       FL.SimBeam{FL.simBeamID}.Bunch.x(6,:)= ...

    48. 

  48.         initX{FL.simBeamID}(6,:)+BEAMLINE{1}.P*dp ;

    49. 

  49.       % also change initial beam conditions

    50. 

  50.       FL.SimBeam{FL.simBeamID}.Bunch.x(1,:)= ...

    51. 

  51.         initX{FL.simBeamID}(1,:)+FL.SimModel.Initial.x.Twiss.eta*dp;

    52. 

  52.       FL.SimBeam{FL.simBeamID}.Bunch.x(2,:)= ...

    53. 

  53.         initX{FL.simBeamID}(2,:)+FL.SimModel.Initial.x.Twiss.etap*dp;

    54. 

  54.       FL.SimBeam{FL.simBeamID}.Bunch.x(3,:)= ...

    55. 

  55.         initX{FL.simBeamID}(3,:)+FL.SimModel.Initial.y.Twiss.eta*dp;

    56. 

  56.       FL.SimBeam{FL.simBeamID}.Bunch.x(4,:)= ...

    57. 

  57.         initX{FL.simBeamID}(4,:)+FL.SimModel.Initial.y.Twiss.etap*dp;

    58. 

  58.     end % if debug

    59. 

  59.   else

    60. 

  60.     % use control system interface (df required)

    61. 

  61.     if (~twoargs)

    62. 

  62.       stat{1}=-1;

    63. 

  63.       stat{2}='changeP: df arg required';

    64. 

  64.       return

    65. 

  65.     end

    66. 

  66.     if (debug)

    67. 

  67.       if (ischar(df))

    68. 

  68.         fprintf(1,'changeP: dp=%s , df=''%s''\n',num2str(dp),df);

    69. 

  69.       else

    70. 

  70.         fprintf(1,'changeP: dp=%s , df=%s\n',num2str(dp),num2str(df));

    71. 

  71.       end

    72. 

  72.       stat{1}=1;

    73. 

  73.     else

    74. 

  74.       stat=DRRF_RampControl(df); % df= 'on', 'off', or dF (kHz)

    75. 

  75.     end

    76. 

  76.   end % if FL.SimModel

    77. 

  77. else % client

    78. 

  78.   % forward request to server

    79. 

  79.   if (twoargs)

    80. 

  80.     cmd=['changeP(',num2str(dp,FL.comPrec),','];

    81. 

  81.     if (ischar(df))

    82. 

  82.       cmd=[cmd,'''',df,''')'];

    83. 

  83.     else

    84. 

  84.       cmd=[cmd,num2str(df,FL.comPrec),')'];

    85. 

  85.     end

    86. 

  86.   else

    87. 

  87.     cmd=['changeP(',num2str(dp,FL.comPrec),')'];

    88. 

  88.   end

    89. 

  89.   stat=sockrw('writechar','cas',cmd);

    90. 

  90.   if (stat{1}~=1),return,end

    91. 

  91.   sockrw('readchar','cas',10);

    92. 

  92. end % if trusted

    93. 

  93. 

    94. 

  94. end
    95.