macro examine_ratios
*
* this is a modified version of make_figure, written to look at the
* ratios of the HBT radii
*
*
calcmax = 3
*
op ndat
op utit
set lwid 8
set fwid 8
set hwid 8
set csiz 0.5
set vsiz 0.6
set xval 0.45
set xwin 0.001
set ywin 0.001
set bcol 1
set ylab 0.9
set xlab 1.6
set gsiz 0.6
*
size 20 30
zone 2 4
*
title '  n=0                              n=2'
ptmax = 0.8
vec/crea par(2) r
do izone=1,8
  case [izone] in
  (1) | nominally R^2_{os,0} -- use it for key
      ykey=9
      lab = 'R^2?os,n! (fm^2!)'
  (2) var = R2os2
      ymin = 0.01
      ymax = 0.75
  (3) var = R2o0
      ymin = 0.01
      ymax = 0.75
      lab = 'R^2?o,n! (fm^2!)'
  (4) var = R2o2
      ymin = 0.01
      ymax = 0.75
  (5) var = R2s0
      ymin = 0.01
      ymax = 0.75
      lab = 'R^2?s,n! (fm^2!)'
  (6) var = R2s2
      ymin = 0.01
      ymax = 0.75
  (7) var = R2l0
      ymin = 0.01
      ymax = 1.1
      lab = 'R^2?l,n! (fm^2!)'
  (8) var = R2l2
      ymin = 0.01
      ymax = 1.1
  endcase
  if ([izone].eq.1) then
    zone 2 4
    hplo/null 0 10 0 10 a
  else
    zone 2 4 [izone] s
    if ([izone].eq.7) then
      ptmaxplot = [ptmax]-0.0001
    else
      ptmaxplot = [ptmax]
    endif
    2dhisto 100 ' ' 40 0 [ptmaxplot] 40 [ymin] [ymax]
  endif
  if ([izone].lt.7) then
    set yval 200
  else
      set yval 0.2
  endif
  if (([izone].eq.1).or.([izone].eq.3).or.([izone].eq.5).or.([izone].eq.7)) then
    set xval 0.1
  else
    set xval -9.1

  endif
  if ([izone].ne.1) then
    hi/plo 100
    graph/hplo/tick
  endif
  if ([izone].ge.7) then
      set asiz 0.6
      atit 'p?T! (GeV/c)'
  endif
  if (([izone].eq.1).or.([izone].eq.3).or.([izone].eq.5).or.([izone].eq.7)) then
    set asiz 0.65
    atit ' ' [lab]
  endif
*
  do icalc=1,[calcmax]
    save_vec = 'var'//[icalc]
    exec examine_ratios#read_data [icalc]
    vec/cop [var] [save_vec]
  enddo
*
  do irat=1,3
    set basl 0.01
    case [irat] in
    (1) lt = 1
        vdiv var1 var2 plot_this
        keytext = '9.75/13'
    (2) lt = 12
        vdiv var1 var3 plot_this
        keytext = '9.75/19.5'
    (3) set basl 0.025
        lt = 12
        vdiv var2 var3 plot_this
        keytext = '13/19.5'
    (4) lt = 13
    endcase
    set ltyp [lt]
    igset ltyp [lt]
    if ([izone].eq.1) then
      pval = $SIGMA(par(1))
      pval2 = $SIGMA(par(2))
      key 2 [ykey] [lt] 'R?y! = '//[keytext] 1.4 L
      ykey = [ykey] - 1.4
*      key 2 [ykey] [lt] 'R?x! = '//[pval2]//' fm' 0.0 L
      ykey = [ykey]-2
    else
      pline $VDIM(pt) pt plot_this
*
      vec/pri plot_this
      wait -10
*
    endif
  enddo  | irat
enddo    | izone
*


return


macro read_data icalc
*
vec/rea pt,npi,v2pi,R2o0,R2o2,r3,R2s0,R2s2,r6,R2os0,R2os2,r9,R2l0,R2l2,r12 'Calc'//[icalc]//'/pions-vs-pT.BW' ! ! -/%/
vec/rea params 'Calc'//[icalc]//'/input.inp'
vec/in par(1) $SIGMA(params(1))    | NOTE params(1) is Ry
vec/in par(2) $SIGMA(params(2))    | NOTE params(2) is Rx
*
* normalize to first point (the one at pT=0)
*
vscale npi $SIGMA(1.0/npi(1)) npi
*
return
*