Current status of MSSM Higgs sector with LHC 13 TeV data

ATLAS and CMS collaborations have reported the results on the Higgs search analyzing $\sim 36$ fb$^{-1}$ data from Run-II of LHC at 13 TeV. In this work, we study the Higgs sector of the phenomenological Minimal Supersymmetric Standard Model, in light of the recent Higgs data, by studying separately the impact of Run-I and Run-II data. One of the major impacts of the new data on the parameter space comes from the direct searches of neutral CP-even and CP-odd heavy Higgses ($H$ and $A$, respectively) in the $H/A \to \tau^{+} \tau^{-}$ channel which disfavours high $\tan\beta$ regions more efficiently than Run-I data. Secondly, we show that the latest result of the rare radiative decay of $B$ meson imposes a slightly stronger constraint on low $\tan \beta$ and low $M_A$ region of the parameter space, as compared to its previous measurement. Further, we find that in a global fit Run-II light Higgs signal strength data is almost comparable in strength with the corresponding Run-I data. Finally, we discuss scenarios with the Heavy Higgs boson decaying into electroweakinos and third generation squarks and sleptons.