Glauber model for small system using anisotropic and inhomogeneous density profile of proton

Recent studies reveal that at high energies, collisions of small system like $p+p$ gives signatures similar to that widely observed in heavy ion collisions hinting towards a possibility of forming a medium with collective behaviour. With this motivation, we have used the Glauber model, which is traditionally applied to heavy ion collisions, in small system using anisotropic and inhomogeneous density profile of proton and found that the proposed model reproduces the charged particle multiplicity distribution of $p+p$ collisions at LHC energies very well. Collision geometric properties like mean impact parameter, mean number of binary collisions ($\langle N_{coll} \rangle$) and mean number of participants ($\langle N_{part} \rangle$) at different multiplicities are determined. Having estimated $\langle N_{coll} \rangle$, we have calculated nuclear modification-like factor ($R_{HL}$) in $p+p$ collisions. We also estimated eccentricity and elliptic flow as a function of charged particle multiplicity using the linear response to initial geometry.