Recent studies reveal that at high energies, collisions of small system like p + p give signatures similar to that widely observed in heavy ion collisions hinting towards a possibility of forming a medium with collective behaviour. With this motivation, in this work, we have used 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 ( N coll ) and mean number of participants ( Npart ) at different multiplicities are determined. Having estimated N coll , we have calculated nuclear modification-like factor (Rpp) in p+p collisions. We also estimated eccentricity and elliptic flow as a function of charged particle multiplicity using linear response to initial geometry.