We present the results of a 6.4 square degrees imaging survey of the Pleiades cluster in the I and Z-bands. The survey extends up to 3{deg} from the cluster center and is 90% complete down to I~=22. It covers a mass range from 0.03M_{sun}_ to 0.48M_{sun}_ and yields 40 brown dwarf candidates (BDCs) of which 29 are new. The spatial distribution of BDCs is fitted by a King profile in order to estimate the cluster substellar core radius. The Pleiades mass function is then derived across the stellar-substellar boundary and we find that, between 0.03M_{sun}_ and 0.48M_{sun}_, it is well represented by a single power-law, dN/dM{prop.to}M^-{alpha}^, with an index alpha=0.60+/-0.11. Over a larger mass domain, however, from 0.03M_{sun}_ to 10M_{sun}_, the mass function is better fitted by a log-normal function. We estimate that brown dwarfs represent about 25% of the cluster population which nevertheless makes up less than 1.5% of the cluster mass. The early dynamical evolution of the cluster appears to have had little effect on its present mass distribution at an age of 120Myr. Comparison between the Pleiades mass function and the Galactic field mass function suggests that apparent differences may be mostly due to unresolved binary systems.