The TeV blazar Ton 599 has exhibited a peculiar flare in 2017 November. The temporal variation of the source is studied using simultaneous $\gamma$-ray data from $\textit{Fermi}$ Large Area Telescope and radio data from Owens Valley Radio Observatory's 40 m telescope, over the period of nine years. Four major flaring periods are observed in the $\gamma$-ray energy band of 0.1-300 GeV. These periods are studied on a shorter timescale and modeled with a time-dependent function containing exponential rising and decaying components. The physical parameters of the jet are estimated numerically and compared with those reported in the literature. During the fourth flare a bunch of high energy photons ($>$10 GeV) were detected. The two highest energy photons having an energy of 76.9 GeV and 61.9 GeV are detected on MJD 58059.0 and 58073.3, respectively. This observation possibly constrains the $\gamma$-ray emission region to lie near outer edge or outside the broad line region of size $\sim$0.08 pc. The variation of equivalent width of a Mg-II line is studied using the spectroscopic data from Steward observatory. It was observed that the equivalent width of the line varies inversely with the underlying power-law continuum.

Aim : The nearby TeV blazar 1ES 1959+650 (z=0.047) was reported to be in flaring state during June - July 2016 by Fermi-LAT, FACT, MAGIC and VERITAS collaborations. We studied the spectral energy distributions (SEDs) in different states of the flare during MJD 57530 - 57589 using simultaneous multiwaveband data to understand the possible broadband emission scenario during the flare. Methods : The UV/optical and X-ray data from UVOT and XRT respectively on board Swift and high energy $\gamma$-ray data from Fermi-LAT are used to generate multiwaveband lightcurves as well as to obtain high flux states and quiescent state SEDs. The correlation and lag between different energy bands is quantified using discrete correlation function. The synchrotron self Compton (SSC) model was used to reproduce the observed SEDs during flaring and quiescent states of the source. Results : A decent correlation is seen between X-ray and high energy $\gamma$-ray fluxes. The spectral hardening with increase in the flux is seen in X-ray band. The powerlaw index vs flux plot in $\gamma$-ray band indicates the different emission regions for 0.1 - 3 GeV and 3-300 GeV energy photons. Two zone SSC model satisfactorily fits the observed broadband SEDs. The inner zone is mainly responsible for producing synchrotron peak and high energy $\gamma$-ray part of the SED in all states. The second zone is mainly required to produce less variable optical/UV and low energy $\gamma$-ray emission. Conclusions : Conventional single zone SSC model does not satisfactorily explain broadband emission during observation period considered. There is an indication of two emission zones in the jet which are responsible for producing broadband emission from optical to high energy $\gamma$-rays.