  We propose trial wave functions for quasiparticle and exciton excitations of
the Moore-Read Pfaffian fractional quantum Hall states, both for bosons and for
fermions, and study these numerically. Our construction of trial wave functions
employs a picture of the bosonic Moore-Read state as a symmetrized double layer
composite fermion state. We obtain the number of independent angular momentum
multiplets of quasiparticle and exciton trial states for systems of up to 20
electrons. We find that the counting for quasielectrons at large angular
momentum on the sphere matches that expected from the CFT which describes the
Moore-Read state's boundary theory. In particular, the counting for
quasielectrons is the same as for quasiholes, in accordance with the idea that
the CFT describing both sides of the FQH plateau should be the same. We also
show that our trial wave functions have good overlaps with exact wave functions
obtained using various interactions, including second Landau level Coulomb
interactions and the 3-body delta interaction for which the Pfaffian states and
their quasiholes are exact ground states. We discuss how these results relate
to recent work by Sreejith et al. on a similar set of trial wave functions for
excitations over the Pfaffian state as well as to earlier work by Hansson et
al., which has produced trial wave functions for quasiparticles based on
conformal field theory methods and by Bernevig and Haldane, which produced
trial wave functions based on clustering properties and `squeezing'.
