  In more than four years of observation the Large Area Telescope on board the
Fermi satellite has identified pulsed $\gamma$-ray emission from more than 80
young pulsars, providing light curves with high statistics. Fitting the
observations with geometrical models can provide estimates of the magnetic
obliquity $\alpha$ and aspect angle $\zeta$, yielding estimates of the
radiation beaming factor and luminosity. Using $\gamma$-ray emission geometries
(Polar Cap, Slot Gap, Outer Gap, One Pole Caustic) and radio emission geometry,
we fit $\gamma$-ray light curves for 76 young pulsars and we jointly fit their
$\gamma$-ray plus radio light curves when possible. We find that a joint radio
plus $\gamma$-ray fit strategy is important to obtain ($\alpha$, $\zeta$)
estimates that can explain simultaneous radio and $\gamma$-ray emission. The
intermediate-to-high altitude magnetosphere models, Slot Gap, Outer Gap, and
One pole Caustic, are favoured in explaining the observations. We find no
evolution of $\alpha$ on a time scale of a million years. For all emission
geometries our derived $\gamma$-ray beaming factors are generally less than one
and do not significantly evolve with the spin-down power. A more pronounced
beaming factor vs. spin-down power correlation is observed for Slot Gap model
and radio-quiet pulsars and for the Outer Gap model and radio-loud pulsars. For
all models, the correlation between $\gamma$-ray luminosity and spin-down power
is consistent with a square root dependence. The $\gamma$-ray luminosities
obtained by using our beaming factors not exceed the spin-down power. This
suggests that assuming a beaming factor of one for all objects, as done in
other studies, likely overestimates the real values. The data show a relation
between the pulsar spectral characteristics and the width of the accelerator
gap that is consistent with the theoretical prediction for the Slot Gap model.
