Float_t combined_error(Float_t error1, Float_t error2) {
  return 1.0/sqrt(1.0/(error1*error1) + 1.0/(error2*error2));
}