glm(log(y) ~ a + b, family = gamma, data = ...)

glm(y ~ a + b, family = gamma(link="log"), data = ...)