import numpy as np
    import matplotlib.pyplot as plt

    data = ((2970, 110.6/10.1, 'Saturn V'),
            (773,  50.5/5.4,   'Ariane 5 ES'),
            (549,  70./3.7,    'Falcon 9 FT'),
            (283,  62.5/5.0,   'Delta IV medium'),
            (12.5, 17./1.2,    'Electron'),
            (2.6,  14.5/0.44,  'Black Brandt 10'))
    data1 = [x for x in data if not 'Falcon' in x[2]]


    if True:
        fig = plt.figure()
        ax  = fig.add_subplot(1, 1, 1)
        for (x, y, name) in data:
            if any([n in name for n in ('Brandt', 'Falcon')]):
                ax.plot(x, y, 'or', markersize=10)
            else:
                ax.plot(x, y, 'ok', markersize=10)
        ax.text(445, 20.8, 'F9', fontsize=16)
        ax.text( 1.9,    35, 'BB 10', fontsize=16)
        ax.set_xscale('log')
        x0, x1 = ax.get_xlim()
        ax.plot([x0, x1], [9, 9], '-k')
        ax.plot([x0, x1], [15, 15], '-k')
        ax.set_ylabel('height/diameter', fontsize=18)
        ax.set_xlabel('mass (metric tons)', fontsize=18)
        ax.set_ylim(0, 40)
        plt.setp(ax.get_xticklabels(), fontsize=14)
        plt.setp(ax.get_yticklabels(), fontsize=14)

        plt.show()