  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.
