Manually calculating robust standard errors for pglm using gradient and hessian matrix

Question

I want to manually calculate robust standard error for a fixed effect poisson model produced using the pglm function that, unlike the plm function, does support sandwich error m