  4C+28.07 is a $\gamma$-ray Flat Spectrum Radio Quasar (FSRQ) type source. It
is often monitored at different frequencies, though long-term multi-wavelength
data of this source have not been modelled in detail before. We have analyzed
$\sim 12$ years (Aug, 2008 - May, 2020) of Fermi-LAT data with a binning of 10
days time scale and observed three distinctive flaring states. Each flaring
state consists of different phases of activity, namely, pre-flare, flare \&
post-flare regions. $\gamma$-ray spectral analysis of these different activity
phases has been performed and the best fit model for its spectra is found to be
a Log-parabola model. We have also studied the correlation of simultaneous
$\gamma$-ray light curves with the optical \& radio counterparts in these
flaring states and report the DCF with a 95\% significance level. A large time
delay is found between radio and gamma-ray data for two flares, indicating two
zones of emission. We have fitted the multi-wavelength data with a two-zone
leptonic model. In our two-zone leptonic model, the maximum required power in
the jet is 9.64 $\times$ 10$^{46}$ erg sec$^{-1}$, which is lower than its
Eddington luminosity $2.29\times 10^{47}$ erg sec$^{-1}$.
