model in its default configuration (i.e. coupled mode) returns a fit similar to the `Torus` one: $\chi^2$/$\nu$ = 217/116, $\Gamma \sim$ 1.4 and $N_{\rm H}$ $\sim$ 6 $\times$ 10$^{23}$ cm$^{-2}$. A common feature of these models is the underprediction of both the flux of the line component at 6-7 keV, and the emission below 5 keV. A `MYTorus` model in its decoupled mode is then applied ($\chi^2$/$\nu$ = 161/115). The front-scattered component vanishes, the photon index is $\sim$ 1.6 and the column density is $\sim$ 10$^{24}$ cm$^{-2}$. These values are consistent with what is found by applying phenomenological models in the first part of the analysis. We note that scattering $\sim$ 3% of the primary continuum into a scattered power law and adding a line feature at (6.8 $\pm$ 0.1) keV do not change the fundamental fit parameters, but improve it at more than 99% confidence limit ($\chi^2$/$\nu$ =113/112, $\Delta\chi^2$/$\Delta\nu$ = 48/3). Even if the best fit model is made up of a combination of `plcabs`, `pexrav` and `zgauss` models ($\chi^2$/$\nu$ = 85/111), we choose to rely on the best fit among the self-consistent ones, as for the other sources, which is the decoupled `MYTorus` model (see Figure \[fig:subfig1\]). Summarizing, 