We characterize the spatial density of the Pan-STARRS1 (PS1) sample of RR Lyrae stars, to study the properties of the old Galactic stellar halo as traced by RRab stars. This sample of 44,403 sources spans Galactocentric radii of $0.55 ; mathrm{kpc} leq R_{mathrm{gc}} leq 141 ; mathrm{kpc}$ with a distance precision of 3% and thus is able to trace the halo out to larger distances than most previous studies. After excising stars that are attributed to dense regions such as stellar streams, the Galactic disc and bulge as well as halo globular clusters, the sample contains ${sim}11,000$ sources within $20 ; mathrm{kpc} leq R_{mathrm{gc}} leq 131 ; mathrm{kpc}$. We then apply forward modeling using ellipsoidal stellar density models $rho(l,b,R_{mathrm{gc}})$ both with a constant and a radius-dependent halo flattening $q(R_{mathrm{gc}})$. Assuming constant flattening $q$, the distribution of the sources is reasonably well fit from $20 ; mathrm{kpc}$ to $131 ; mathrm{kpc}$ by a single power law with $n=4.40^{+0.05}_{-0.04}$ and $q=0.918^{+0.016}_{-0.014}$. The distance distribution is fit comparably well by an Einasto profile with $n=9.53^{+0.27}_{-0.28}$, an effective radius $r_{mathrm{eff}}=1.07 pm 0.10 ; mathrm{kpc}$ and a halo flattening of $q=0.923 pm 0.007$. If we allow for a radius-dependent flattening $q(R_{mathrm{gc}})$, we find evidence for a distinct flattening of $q{sim}0.8$ of the inner halo at ${sim} 25 ; mathrm{kpc}$. Additionally, we find that the south Galactic hemisphere is more flattened than the north Galactic hemisphere. The results of our work are largely consistent with many earlier results, e.g. cite{Watkins2009}, cite{Iorio2017}. We find that the stellar halo, as traced in RR Lyrae stars, exhibits a substantial number of further significant over- and underdensities, even after all known overdensities have been masked.