Analytical formulae for the calculation of secondary particle yields in p-A
interactions are given. These formulae can be of great practical importanc
e for fast calculations of neutrino fluxes and for designing new neutrino b
eam-lines. The formulae are based on a parameterization of the inclusive in
variant cross sections for secondary particle production measured in p-Be i
nteractions. Data collected in different energy ranges and kinematic region
s are used. The accuracy of the fit to the data with the empirical formulae
adopted is within the experimental uncertainties. Prescriptions to extrapo
late this parameterization to finite targets and to targets of different ma
terials are given. The results obtained are then used as an input for the s
imulation of neutrino beams. We show that our approach describes well the m
ain characteristics of measured neutrino spectra at CERN. Thus it may be us
ed in fast simulations aiming at the optimisation of the long-baseline neut
rino beams at CERN and FNAL. In particular we will show our predictions for
the CNGS beam from CERN to Gran Sasso.