The rapidity distributions of [Formula: see text] mesons produced in proton-lead ([Formula: see text]-Pb) collisions at center-of-mass energy per nucleon pair [Formula: see text][Formula: see text]TeV are studied by using a multisource thermal model and compared with the experimental data of the LHCb and ALICE Collaborations. Correspondingly, the pseudorapidity distributions are accurately obtained from the parameters extracted from the rapidity distributions. At the same time, the transverse momentum distributions in the same experiments are described by the simplest Erlang distribution which is the folding result of two exponential distributions which are contributed by the target and projectile partons, respectively.