The treatment of the pairing interaction in mean-field-based models is addressed. In particular, the possibility to use pair transfers as a tool to better constrain this interaction is discussed. First, pairing inter- actions with various density dependencies (surface/volume mixing) are used in the microscopic Hartree-Fock- Bogoliubov + quasiparticle random-phase approximation model to generate the form factors to be used in reac- tion calculations. Cross sections for (p,t) two-neutron transfer reactions are calculated in the one-step zero-range distorted-wave Born approximation for some Tin isotopes and for incident proton energies from 15 to 35 MeV. Three different surface/volume mixings of a zero-range density-dependent pairing interaction are employed in the microscopic calculations and the sensitivity of the cross sections to the different mixings is analyzed. Differences among the three different theoretical predictions are found espacially for the nucleus 136 Sn and they are more important at the incident proton energy of 15 MeV. We thus indicate (p,t) two-neutron transfer reactions with very neutron-rich Sn isotopes and at proton energies around 15 MeV as good experimental cases where the sur- face/volume mixing of the pairing interaction may be probed. In the second part of the manuscript, ground-state to ground-state transitions are investigated. Approximations made to estimate two-nucleon transfer probabilities in ground-state to ground-state transitions and the physical interpretation of these probabilities are discussed. Probabilities are often calculated by approximating both ground states of the initial nucleus A and of the final nucleus A ±2 by the same quasiparticle vacuum. We analyze two improvements of this approach. First, the effect of using two different ground states with average numbers of particles A and A ±2 is quantified. Second, by using projection techniques, the role of particle number restoration is analyzed. Our analysis shows that the improved treatment plays a role close to magicity, leading to an enhancement of the pair-transfer probability. In midshell regions, part of the error made by approximating the initial and final ground states by a single vacuum is com- pensated by projecting onto a good particle number. Surface effects are analyzed by using pairing interactions with a different volume/surface mixing. Finally, a simple expression of the pair-transfer probability is given in terms of occupation probabilities in the canonical basis. We show that, in the canonical basis formulation, surface effects that are visible in the transfer probability are related to the fragmentation of single-particle occupancies close to the Fermi energy. This provides a complementary interpretation with respect to the standard quasiparti- cle representation where surface effects are generated by the integrated radial profiles of the contributing wave functions.