docs: def gp_sims_panel(version):
    """panel plot of cocktail simulations at all energies, includ. total
  
    :param version: plot version / input subdir name
    :type version: str
    """
    inDir, outDir = getWorkDirs()
    inDir = os.path.join(inDir, version)
    mesons = ['pion', 'eta', 'etap', 'rho', 'omega', 'phi', 'jpsi']
    fstems = ['cocktail_contribs/'+m for m in mesons] + ['cocktail', 'cocktail_contribs/ccbar']
    data = OrderedDict((energy, [
        np.loadtxt(open(os.path.join(inDir, fstem+str(energy)+'.dat'), 'rb'))
        for fstem in fstems
    ]) for energy in energies)
    for v in data.values():
        # keep syserrs for total cocktail
        v[-2][:,(2,3)] = 0
        # all errors zero for cocktail contribs
        v[-1][:,2:] = 0
        for d in v[:-2]: d[:,2:] = 0
    make_panel(
        dpt_dict = OrderedDict((
            ' '.join([getEnergy4Key(str(energy)), 'GeV']),
            [ data[energy], [
                'with %s lc %s lw 5 lt %d' % (
                    'lines' if i!=len(fstems)-2 else 'filledcurves pt 0',
                    default_colors[(-2*i-2) if i!=len(fstems)-2 else 0] \
                    if i!=len(fstems)-1 else default_colors[1], int(i==3)+1
                ) for i in xrange(len(fstems))
            ], [particleLabel4Key(m) for m in mesons] + [
                'Cocktail (w/o {/Symbol \162})', particleLabel4Key('ccbar')
            ]]
        ) for energy in energies),
        name = os.path.join(outDir, 'sims_panel'),
        ylog = True, xr = [0.,3.2], yr = [1e-6,9],
        xlabel = xlabel, ylabel = ylabel,
        layout = '2x2', size = '7in,9in',
        key = ['width -4', 'spacing 1.5', 'nobox', 'at graph 0.9,0.95']
    )