The charm quark mass to two-loop order

The truncation of the perturbative series at one loop order for the mass renormalization constants remain a significant systematic uncertainty in the determination of heavy quark masses in lattice QCD. We present here a high beta Monte Carlo calculation of the two loop mass renormalization constant for clover-improved fermions near the charm mass in the Fermilab heavy quark formalism. A preliminary value for the charm quark mass in the over MS scheme at two loop order is reported