code: def plotnoise(noisepkl, mergepkl, plot_width=950, plot_height=400):
      
  
      d = pickle.load(open(mergepkl))
      ndist, imstd, flagfrac = plotnoisedist(noisepkl, plot_width=plot_width/2, plot_height=plot_height)
      fluxscale = calcfluxscale(d, imstd, flagfrac)
      logger.info('Median image noise is {0:.3} Jy.'.format(fluxscale*imstd))
      ncum, imnoise = plotnoisecum(noisepkl, fluxscale=fluxscale, plot_width=plot_width/2, plot_height=plot_height)
  
      hndle = show(Row(ndist, ncum, width=plot_width, height=plot_height))
      return imnoise
docs: Make two panel plot to summary noise analysis with estimated flux scale