In this paper we perform an amplitude analysis of essentially all published pion and kaon pair production data from two-photon collisions below 1.5 GeV. This includes all the high statistics results from Belle, as well as older data from Mark II at SLAC, CELLO at DESY, and Crystal Ball at SLAC. The purpose of this analysis is to provide as close to a model-independent determination of the γγ to meson pair amplitudes as possible. Having data with limited angular coverage, typically j cos θj < 0.6-0.8, and no polarization information for reactions in which spin is an essential complication, the determination of the underlying amplitudes might appear an intractable problem. However, imposing the basic constraints required by analyticity, unitarity, and crossing symmetry makes up for the experimentally missing information. Above 1.5 GeV multimeson production channels become important, and we have too little information to resolve the amplitudes. Nevertheless, below 1.5 GeV the two-photon production of hadron pairs serves as a paradigm for the application of S-matrix techniques. Final state interactions among the meson pairs are critical to this analysis. To fix these, we include the latest ππ → ππ, KK scattering amplitudes given by dispersive analyses, supplemented in the KK threshold region by the recent precision Dalitz plot analysis from BABAR. With these hadronic amplitudes built into unitarity, we can constrain the overall description of γγ → ππ and KK data sets, both integrated and differential cross sections, including the high statistics charged and neutral pion, as well as K s K s , data from Belle. Since this analysis invokes coupled hadronic channels, having data on both ππ and KK reduces the solution space to essentially a single form in the region where these channels saturate unitarity. For the ππ channel, the separation of isospin-0 and -2 and helicity-0 and -2 components is complete. We present the partial wave amplitudes, show how well they fit all the available data, and give the two-photon couplings of scalar and tensor resonances that appear. These partial waves are important inputs into forthcoming dispersive calculations of hadronic light-by-light scattering. PHYSICAL REVIEW D 90, 036004 (2014) 036004-2 LING-YUN DAI AND M. R. PENNINGTON PHYSICAL REVIEW D 90, 036004 (2014) 036004-4 FIG. 8 (color online). Fit to the γγ → π þ π − differential cross section of CELLO experiment. Here, Cello1 is from Harjes [3] and Cello2 from Behrend et al. [4]. The numbers give the central energy in GeV of each angular distribution listed in order of the cross section at z ¼ 0, where z ¼ cos θ. COMPREHENSIVE AMPLITUDE ANALYSIS OF … PHYSICAL REVIEW D 90, 036004 (2014) 036004-15 FIG. 9 (color online). Fit to the γγ → π þ π − differential cross section of the Belle experiment [13]. The numbers give the central energy in GeV of each angular distribution listed in order of the cross section at z ¼ 0, where z ¼ cos θ. The data are normalized so that the integrated cross section is just a sum of the differential cross sect...