dist-ghtFit function

GHT distribution fit