It is worthy to note further that
nuclear reaction calculations involve several (sub-)models to compute the key input quantities like level density, \(\gamma \)  strength function, optical potential, nuclear masses, pre-equilibrium and the compound process.
Their interplay makes the model evaluation almost impossible.
Instead of developing a consistent reaction model, we attempt to find the best model parameters of TALYS, which can well describe the
known experimental data. We are aware that the parameters might be mass-range dependent or only hold to a certain precision in applying to nuclei nearby, therefore we concentrate on the mass range of 160.
This strategy is practically useful to reduce the nuclear physics
inputs in \(\gamma \)  process simulations [1]}, [2]}.
We thus made a Monte Carlo calculation by randomly varying the key physical inputs parameters [3]}, and employed the least square method to find the calculation that matches best both the \(^{160}\) Dy(\(p,\gamma \) )\(^{161}\) Ho and \(^{162}\) Er(\(p,\gamma \) )\(^{163}\) Tm data. The results are shown as well in Figure REF .
