You should retrieve the covariance matrix from your fit and then pass it to the TF1::IntegralError procedure.
If you cannot get the covariance matrix at all, try: Error on Integral of function not Resulting From Fit
BTW. Your function’s formula seems to be quite simple so, if you do not have any information about the covariance matrix, I think you can calculate both the integral and any errors analytically.