macro set_param1
* this macro tells what are the number of calculations and max pT range
calcmax = 3
ptmax = 0.8
vec/in pass(1) [calcmax]
vec/in pass(2) [ptmax]
* following is for the little red arrows
vec/crea yarro([calcmax]) r 0 0.083 0.155
*
vec/crea calcs([calcmax]) r 1 23 31

return
*
*
macro set_histo_limits izone plot_ratio show_kolb
* set max and min limits for y-axis
  case [izone] in
  (1) | nominally R^2_{os,0} -- use it for key
      ykey=9
      ymin=0
      ymax=0 
  (2) var = R2os2
      if ([plot_ratio].eq.0) then
        ymin = -0.49
        ymax = 5.0
      else
        ymin = -0.09
        ymax = 0.29
      endif
  (3) var = R2o0
      ymin = 0.0
      ymax = 39.9
  (4) var = R2o2
      if ([plot_ratio].eq.0) then
        ymin = -5.0
        ymax = 0.49
      else
        ymin = -0.29
        ymax = 0.09
      endif
  (5) var = R2s0
      ymin = 0.0
      ymax = 39.9
  (6) var = R2s2
      if ([plot_ratio].eq.0) then
        ymin = -0.49
        ymax = 5.0
      else
        ymin = -0.09
        ymax = 0.29
      endif
  (7) var = R2l0
      ymin = 0.0
      ymax = 140.0
  (8) var = R2l2
      if ([plot_ratio].eq.0) then
        ymin = -0.19
        ymax = 0.19
      else
        ymin = -0.019
        ymax = 0.019
      endif
  endcase
*
if ([show_kolb].ne.0) then
  case [izone] in   
  (2) var = R2os2
      if ([plot_ratio].eq.0) then
        ymin = -0.49
        ymax = 5.0
      else
        ymin = -0.119
        ymax = 0.29
      endif
  (4) var = R2o2
      if ([plot_ratio].eq.0) then
        ymin = -5.0
        ymax = 0.49
      else
        ymin = -0.29
        ymax = 0.119
      endif
  (6) var = R2s2
      if ([plot_ratio].eq.0) then
        ymin = -0.49
        ymax = 5.0
      else
        ymin = -0.119
        ymax = 0.29
      endif
  (8) var = R2l2
      if ([plot_ratio].eq.0) then
        ymin = -0.19
        ymax = 0.19
      else
        ymin = -0.029
        ymax = 0.029
      endif
  endcase
endif
vec/in pass(1) [ymin]
vec/in pass(2) [ymax]
return
*
*
macro do_caption icalc lt
* Parameter numbers:
* 1- Ry;  2-Rx;  3-T;  4-rho0;  5-rhoa;  6-tau0;  7-a_s;  8-Delta_tau;
if ([icalc].eq.1) then
  vec/crea yk(1) r 9
endif
ykey = $SIGMA(yk(1))
* ----------- change part below this line ------------
pval = $SIGMA(params(1))
pval2 = $SIGMA(params(2))
key 2 [ykey] [lt] 'R?y! = '//[pval]//' fm' 1.4 L
ykey = [ykey] - 1.4
key 2 [ykey] [lt] 'R?x! = '//[pval2]//' fm' 0.0 L
* ----------- change part above this line ------------
ykey = [ykey]-2
vec/in yk(1) [ykey]   | save it in a vector for the next time 'round
return
*