Towards the determination of the photon parton distribution function constrained by LHC data