code: def gp_sims_total_overlay(version):
    
    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',
    )
docs: single plot comparing total cocktails at all energies
  
    :param version: plot version / input subdir name
    :type version: str