  Since the detection of very high energy (VHE) $\gamma$-rays from Mrk 501, its
broad band emission of radiation was mostly and quite effectively modeled using
one zone emission scenario. However, broadband spectral and flux variability
studies enabled by the multiwavelength campaigns carried out during the recent
years have revealed rather complex behavior of Mrk 501. The observed emission
from Mrk 501 could be due to a complex superposition of multiple emission
zones. Moreover new evidences of detection of very hard intrinsic $\gamma$-ray
spectra obtained from {\it Fermi}--LAT observations have challenged the
theories about origin of VHE $\gamma$-rays. Our studies based on {\it
Fermi}--LAT data indicate the existence of two separate components in the
spectrum, one for low energy $\gamma$-rays and the other for high energy
$\gamma$-rays. Using multiwaveband data from several ground and space based
instruments, in addition to HAGAR data, the spectral energy distribution of
Mrk~501 is obtained for various flux states observed during 2011. In the
present work, this observed broadband spectral energy distribution is
reproduced with a leptonic, multi-zone Synchrotron Self-Compton model.
