  We present a quantitative and homogeneous analysis of the broadband
multicolor photometric data sets gathered so far on rapidly pulsating hot B
subdwarf stars. This concerns seven distinct data sets related to six different
stars. Our analysis is carried out within the theoretical framework developed
by Randall et al., which includes full nonadiabatic effects. The goal of this
analysis is partial mode identification, i.e., the determination of the degree
index l of each of the observed pulsation modes. We assume possible values of l
from 0 to 5 in our calculations. For each target star, we compute a specific
model atmosphere and a specific pulsation model using estimates of the
atmospheric parameters coming from time-averaged optical spectroscopy. For
every assumed value of l, we use a formal chi-squared approach to model the
observed amplitude-wavelength distribution of each mode, and we compute a
quality-of-fit Q probability to quantify the derived fit and to discriminate
objectively between the various solutions. We find that no completely
convincing and unambiguous l identification is possible on the basis of the
available data, although partial mode discrimination has been reached for 25
out of the 41 modes studied. A brief statistical study of these results
suggests that a majority of the modes must have l values of 0, 1, and 2, but
also that modes with l = 4 could very well be present while modes with l = 3
appear to be rarer. This is in line with recent results showing that l = 4
modes in rapidly pulsating B subdwarfs have a higher visibility in the optical
domain than modes with l = 3. Although somewhat disappointing in terms of mode
discrimination, our results still suggest that the full potential of multicolor
photometry for l identification in pulsating subdwarfs is within reach.
