int phaseShiftNumFrames = 16;

      local StFFT fft;
      fft.setWindowSize(fftBits);
      fft.setEnableWindow(2/*tri*/);
      fft.setEnableMagnitudeScaling(true);
      fft.setInputSamples(inputBuf, 1/*numCh*/, null/*img*/, 0/*imgSkip*/);
      int fftOffset = fftWinSize/2;
      int numBins = fftWinSize/2 +1;
      local FloatArray magnitudes; magnitudes.alloc(numBins);
      local IntArray magnitudesSortIndices; magnitudesSortIndices.allocAndFill(numBins, 0);
      local FloatArray phases; phases.alloc(numBins);

      fft.fft(fftOffset, 0/*inChOffset*/);
      int magNumBins = fft.calcMagnitudes(magnitudes);
      magnitudes.fixNaNToZero();
      magnitudes.sortByValue(magnitudesSortIndices);
      fft.calcPhases(phases);

      int sortedIdx = magnitudesSortIndices.last;

      // Repeat FFT to find true frequency
      local FloatArray inputBuf2;
      inputBuf2.visit(inputBuf, phaseShiftNumFrames, inputBuf.numElements-phaseShiftNumFrames);
      local FloatArray phases2; phases2.alloc(numBins);
      fft.setInputSamples(inputBuf2, 1/*numCh*/, null/*img*/, 0/*imgSkip*/);
      fft.fft(fftOffset, 0/*inChOffset*/);
      fft.calcPhases(phases2);

      float phase1 = phases[sortedIdx];
      float phase2 = phases2[sortedIdx];
      if(phase2 < phase1)
         phase2 += 2PI;

      float phaseDelta = phase2 - phase1;
   
      float binFreq = (sortedIdx * _sampleRate * 0.5) / numBins;
      float binPeriodSz = (_sampleRate / binFreq);
      float expectedPhaseDelta = (2PI * phaseShiftNumFrames) / binPeriodSz;
      float actualBinFreq = (binFreq * phaseDelta) / expectedPhaseDelta;