© 2021 World Scientific Publishing Company.111In, a known gamma emitter, is used for many medical purposes such as imaging of myocardial metastases. It can be produced by using different nuclear reactions. In this study, the reactions of 109Ag(?,2n)111In, 111Cd(p,n)111In, 112Cd(p,2n)111In, 113Cd(p,3n)111In and 114Cd(p,4n)111In, which are the production routes of 111In, were investigated. Production cross-section calculations were performed by using equilibrium and pre-equilibrium models of TALYS 1.95 and EMPIRE 3.2 nuclear reaction codes. Hauser-Feshbach Model was appointed in both codes for calculations of equilibrium approximations. Exciton and Hybrid Monte Carlo Simulation (HMS) models were used in the EMPIRE 3.2, whereas Two-Component Exciton and Geometry Dependent Hybrid Model, which is implemented to TALYS code, has been used in the TALYS 1.95 for pre-equilibrium reactions. Also, a weighting matrix of the nuclear models was obtained by using statistical variance analysis. The optimum beam energy to obtain 111In has been determined by using the results obtained from this weighting matrix.