Abstract. We calculate O(c~s) corrections to large Pr double photon production in hadronic collisions. We find that these corrections to the basic qq ~ 77 subprocess are important, preventing to describe double photon production on the basis of the leading logarithm approximation only. We give a phenomenological discussion of the results obtained at SPS, ISR and SpOS energies. In particular, we investigate how the effect of the intrinsic parton < ~c r) may be disentangled from the perturbative contribution. We also calculate the aplanarity distribution characteristic of 3 jet events.