A generalization of the Hermite–Hadamard (HH) inequality for a positive convex stochastic process, by means of a newly proposed fractional integral operator, is hereby established. Results involving the Riemann– Liouville, Hadamard, Erdélyi–Kober, Katugampola, Weyl and Liouville fractional integrals are deduced as particular cases of our main result. In addition, we also apply some known HH results to obtain some estimates for the expectations of integrals of convex and