code: def newsnr_sgveto(snr, bchisq, sgchisq):
      
      nsnr = numpy.array(newsnr(snr, bchisq), ndmin=1)
      sgchisq = numpy.array(sgchisq, ndmin=1)
      t = numpy.array(sgchisq > 4, ndmin=1)
      if len(t):
          nsnr[t] = nsnr[t] / (sgchisq[t] / 4.0) ** 0.5
  
      # If snr input is float, return a float. Otherwise return numpy array.
      if hasattr(snr, '__len__'):
          return nsnr
      else:
          return nsnr[0]
docs: Combined SNR derived from NewSNR and Sine-Gaussian Chisq