We perform high-precision joint fitting of the ensemble of Auger extragalactic spectrum, composition, IceCube neutrino spectrum and upper limits, and Fermi-LAT diffuse gamma spectrum, to constrain the properties of UHECR sources, the source environments, source evolution and distance to nearest source. At present, the largest source of uncertainty comes from the LHC-tuned hadronic event generators, limiting the ability to derive decisive conclusions on all of the above topics. Nonetheless, interesting constraints on UHECR, neutrino and gamma ray sources can (and will) be presented.