Recently the operator product expansion coefficients of the 3D Ising model universality class have been calculated by studying via Monte Carlo simulation the two-point functions perturbed from the critical point with a relevant field. We show that this method can be applied also when the perturbation is performed with a relevant field coupled to a non uniform potential acting as a trap. This setting is described by the trap size scaling ansatz, that can be combined with the general framework of the conformal perturbation in order to write down the correlators < σ(r)σ(0) >, < σ(r) (0) > and < (r) (0) >, from which the operator product expansion coefficients can be estimated. We find C σ σ = 1.051(3) , in agreement with the results already known in the literature, and C = 1.32(15) , confirming and improving the previous estimate obtained in the uniform perturbation case.