single plot comparing total cocktails at all energies

  :param version: plot version / input subdir name
  :type version: str

def gp_sims_total_overlay(version):
  """single plot comparing total cocktails at all energies

  :param version: plot version / input subdir name
  :type version: str
  """
  inDir, outDir = getWorkDirs()
  inDir = os.path.join(inDir, version)
  data = OrderedDict()
  for energy in energies:
      fname = os.path.join(inDir, 'cocktail'+str(energy)+'.dat')
      data[energy] = np.loadtxt(open(fname, 'rb'))
      data[energy][:,2:] = 0
  make_plot(
      data = data.values(),
      properties = [
          'with lines lc %s lw 4 lt 1' % (default_colors[i])
          for i in xrange(len(energies))
      ],
      titles = [
          ' '.join([getEnergy4Key(str(energy)), 'GeV'])
          for energy in energies
      ],
      xlabel = xlabel, ylabel = ylabel,
      name = os.path.join(outDir, 'sims_total_overlay'),
      ylog = True, xr = [0.,3.2], yr = [1e-6,9],
      tmargin = 0.98, rmargin = 0.99, bmargin = 0.14,
      size = '8.5in,8in',
  )