def gp_sims_totalerrors_overlay(version):
inDir, outDir = getWorkDirs()<EOL>inDir = os.path.join(inDir, version)<EOL>data = OrderedDict()<EOL>for energy in energies:<EOL><INDENT>fname = os.path.join(inDir, '<STR_LIT>'+str(energy)+'<STR_LIT>')<EOL>data[energy] = np.loadtxt(open(fname, '<STR_LIT:rb>'))<EOL>data[energy][:,<NUM_LIT:1>] = data[energy][:,<NUM_LIT:4>]/data[energy][:,<NUM_LIT:1>]<EOL>data[energy][:,<NUM_LIT:2>:] = <NUM_LIT:0><EOL><DEDENT>make_plot(<EOL>data = list(data.values()),<EOL>properties = [<EOL>'<STR_LIT>' % (default_colors[i])<EOL>for i in range(len(energies))<EOL>],<EOL>titles = [<EOL>'<STR_LIT:U+0020>'.join([getEnergy4Key(str(energy)), '<STR_LIT>'])<EOL>for energy in energies<EOL>],<EOL>xlabel = xlabel, ylabel = '<STR_LIT>',<EOL>name = os.path.join(outDir, '<STR_LIT>'),<EOL>xr = [<NUM_LIT:0.>,<NUM_LIT>], yr = [<NUM_LIT>,<NUM_LIT>], key = ['<STR_LIT>'],<EOL>tmargin = <NUM_LIT>, rmargin = <NUM_LIT>, bmargin = <NUM_LIT>,<EOL>size = '<STR_LIT>',<EOL>)<EOL>
single plot comparing syst. uncertainties on total cocktails at all energies

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