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]
return
*
*
macro set_histo_limits izone plot_ratio
* 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.015
        ymax = 0.19
      endif
  (3) var = R2o0
      ymin = 0.0
      ymax = 49.9
  (4) var = R2o2
      if ([plot_ratio].eq.0) then
        ymin = -5.0
        ymax = 0.49
      else
        ymin = -0.19
        ymax = 0.015
      endif
  (5) var = R2s0
      ymin = 0.0
      ymax = 49.9
  (6) var = R2s2
      if ([plot_ratio].eq.0) then
        ymin = -0.49
        ymax = 5.0
      else
        ymin = -0.03
        ymax = 0.19
      endif
  (7) var = R2l0
      ymin = 0.0
      ymax = 149.0
  (8) var = R2l2
      if ([plot_ratio].eq.0) then
        ymin = -0.29
        ymax = 0.29
      else
        ymin = -0.018
        ymax = 0.018
      endif
  endcase
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 ------------
case [icalc] in
(1) cap='Central'
(2) cap='Mid-Central'
(3) cap='Peripheral'
endcase
key 2 [ykey] [lt] [cap] 1.4 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
*
*