Abstract: Photochemical hazes are expected to form in hot Jupiter atmospheres and may explain the strong scattering slopes and muted spectral features observed in the transmission spectra of many hot Jupiters. Absorption and scattering by photochemical hazes have the potential to drastically alter temperature structure and atmospheric circulation of these planets but have previously been neglected in general circulation models (GCMs). We present GCM simulations of hot Jupiter HD 189733b that include photochemical hazes as a radiatively active tracer fully coupled to atmospheric dynamics. The influence of haze radiative feedback strongly depends on the assumed haze optical properties. For soot hazes, two distinct thermal inversions form, separated by a local temperature minimum around 10$^{-5}$ bar caused by upwelling on the dayside mixing air with low haze abundance upwards. The equatorial jet broadens and slows down. The horizontal distribution of hazes remains relatively similar to simulations with radiatively passive tracers. For Titan-type hazes, the equatorial jet accelerates and extends to much lower pressures, resulting in a dramatically different 3D distribution of hazes compared to radiatively passive or soot hazes. Further experimental and observational studies to constrain the optical properties of photochemical hazes will therefore be crucial for understanding the role of hazes in exoplanet atmospheres. In the dayside emission spectrum, for both types of hazes the amplitude of near-infrared features is reduced, while the emitted flux at longer wavelengths ($>$4 $\mu$m) increases. Haze radiative feedback leads to increased phase curve amplitudes in many infrared wavelength regions, mostly due to stronger dayside emission.