A new paradigm of multisensory observations joined multiband measurements of the radiation coming from celestial objects to develop and confirm models of the origin of high-energy cosmic rays (CR). The integral parameters of the cosmic ray flux, such as energy spectra and mass composition, mostly measured in the last century, deliver useful information on CR origination. Especially useful was an approach to disentangle the overall cosmic ray flux and obtain separate energy spectra of different mass groups. The MAKET-ANI surface array operated on Aragats Mt. in Armenia from 1997 to 2004 turned out to be very well suited for the energy and composition measurements at the “knee” of the cosmic ray spectrum. The problem of event-by-event classification of EAS has been solved by applying Bayesian and neural network techniques. The main results of MAKET-ANI data can be summarized as follows: The estimated energy spectrum of the light mass group of nuclei shows a sharp knee: ∆γ ~0.9, compared to ~0.3 for the all-particle energy spectra. The energy spectrum of the heavy mass group of cosmic rays shows no break in the energy interval of 10^15 – 10^16 eV. In the new era of EAS studies by HAWK, LHAASO, and other experiments aimed to detect point sources of ultra-high energy gamma radiation (PeVatrons) to solve the millennium problem of cosmic ray origin, it is interesting to present and analyze the full pattern of available energy spectra in the energy range from 10^13 to 10^16 eV.