The heritage Andean churches of Arica and Parinacota are representative of the majority of Christian churches in the region. Judging by their architecture, materiality and decoration, they have many features in common and are faithful to the intentions of their original builders. However, their geographical locations and the landscape that surrounds these churches show specific features for each of them. In particular, the orientation of the main axes (towards the altar of the churches) shows some diversity. In this extended area, with little attention from parish priests, the local Aymara culture, with its own worldview, surely engaged in a dialogue with the Western tradition in order to design and build the Indian Reductions. In this paper we present the results obtained from the analysis of the precise spatial orientations of all the cataloged churches, employing the tools of archaeoastronomy, as a complement to the architectural, historical and cultural research of the temples. We also comment briefly on some progress that we hope to make in neighboring areas to the region here studied, which include a large number of very old temples and are part, along with the churches here considered, of a common and widespread cultural koine.

Deep neural networks are a powerful technique that have found ample applications in several branches of Physics. In this work, we apply machine learning algorithms to a specific problem of Cosmic Ray Physics: the estimation of the muon content of extensive air showers when measured at the ground. As a working case, we explore the performance of a deep neural network applied to the signals recorded by the water-Cherenkov detectors of the Surface Detector Array of the Pierre Auger Observatory. We apply deep learning architectures to large sets of simulated data. The inner structure of the neural network is optimized through the use of genetic algorithms. To obtain a prediction of the recorded muon signal in each individual detector, we train neural networks with a mixed sample of light, intermediate and heavy nuclei. When true and predicted signals are compared at detector level, the primary values of the Pearson correlation coefficients are above 95\%. The relative errors of the predicted muon signals are below 10\% and do not depend on the event energy, zenith angle, total signal size, distance range or the hadronic model used to generate the events.

After a historical and cultural overview of the "Indian republic", the Andean religiosity and the churches of the highlands of Arica, we detail our archaeoastronomical studies on the orientations of the Andean Christian churches in the region of Arica and Parinacota, and analyze some preliminary results. We conclude with a brief discussion on the patterns of orientation found in our measurements and on their possible causes.