Bayesian derivation of plasma equilibrium distribution function for tokamak scenarios and the associated Landau collision operator