def cell_type_specificity(ax):
    '''make an imshow of the intranetwork connectivity'''
    masked_array = np.ma.array(params.T_yX, mask=params.T_yX==0)
    # cmap = plt.get_cmap('hot', 20)
    # cmap.set_bad('k', 0.5)
    # im = ax.imshow(masked_array, cmap=cmap, vmin=0, interpolation='nearest')
    im = ax.pcolormesh(masked_array, cmap=cmap, vmin=0, ) #interpolation='nearest')
    ax.axis(ax.axis('tight'))
    ax.invert_yaxis()
    ax.xaxis.set_ticks_position('top')
    ax.set_xticks(np.arange(9)+0.5)
    ax.set_yticks(np.arange(16)+0.5)
    ax.set_xticklabels(params.X, rotation=270)
    ax.set_yticklabels(params.y, )
    ax.xaxis.set_label_position('top')
    ax.set_xlabel(r'$X$', labelpad=-1,fontsize=8)
    ax.set_ylabel(r'$y$', labelpad=0, rotation=0,fontsize=8)

    rect = np.array(ax.get_position().bounds)
    rect[0] += rect[2] + 0.01
    rect[2] = 0.01
    fig = plt.gcf()
    cax = fig.add_axes(rect)

    cbar = plt.colorbar(im, cax=cax)
    #cbar.set_label(r'$\mathcal{T}_{yX}$', ha='center')
    cbar.set_label(r'$\mathcal{T}_{yX}$', labelpad=0)
make an imshow of the intranetwork connectivity