from statsforecast.arima import ARIMASummary

...

sf = StatsForecast(
    models = [
        AutoARIMA(),
    ],
    ...
)

sf.fit(df)

print(ARIMASummary(sf.fitted_[0][0].model_))