## 3 Results
The angular distribution _tx-ty_ of the number of particles which hit the detector in the coordinate system of emulsion at different values of wall thickness L and density \(\rho\): (a) L = 0.5 m, \(\rho\) = 2.65 g/cm\({}^{3}\); (b) L = 1.5 m, \(\rho\) = 2.65 g/cm\({}^{3}\); (c) L = 0.5 m, \(\rho\) = 2.00 g/cm\({}^{3}\); (d) L = 1.5 m, \(\rho\) = 2.00 g/cm\({}^{3}\). The deficit of muons near _ty_\(\sim-0.5\) has to do with the geometry of the building: at this angle, the muon path in the wall is greater (compare cases (**a**) and (**b**) for different values of the wall thickness). The simulation was carried out for two values of wall material density (2.00 g/cm\({}^{3}\) and 2.65 g/cm\({}^{3}\)) and two values of wall thickness (0.5 m and 1.5 m). The muon decays and the angular dependence of their flux were taken into account, and calculations were made for particles with incident angles \(0^{\mathrm{o}}<\theta<90^{\mathrm{o}}\). The maximal value of the initial muon flux obviously corresponds to the vertical direction (with increasing zenith angle \(\theta\) the flux drops). In Figure 5, at a wall density of 2.00 g/cm\({}^{3}\), the peak at ty = -0.8 corresponding to the vertical direction is clearly seen. At a wall density of 2.65 g/cm\({}^{3}\), the peak is weaker due to a significant increase of muon absorption. The distributions in Figure 5 show that the change of wall thickness in the analyzed range has the largest impact on the number of passed particles, since the threefold increase in the wall thickness (from 0.5 m up to 1.5 m) changes the degree of absorption significantly more than the change in the density of material by 25%, from 2.00 g/cm\({}^{3}\) up to 2.65 g/cm\({}^{3}\). The test irradiation of emulsion detectors in the Naryn-Kala citadel lasted for two months. After exposition the detectors were moved to the surface and disassembled; the nuclear emulsions were developed and processed. Scanning of the emulsion layers and online processing of the obtained images was carried out on the automated scanning complex PAVICOM [20]. The most resource-consuming stage of the online processing is the reconstruction of microtracks - segments of tracks between pairs of track-forming grains, for which the track length and slope are within a certain range. Further offline data processing is performed with the FEDRA (Framework for Emulsion Data Reconstruction and Analysis) package [21] providing for reconstruction of tracks, the relative position of emulsion plates in space and the track reconstruction in the entire volume of the scanned data. Figure 6 shows examples of image reconstruction in two adjacent films of one element of five exposed emulsion detectors. Each black dot on the Figure corresponds to one so-called base track passing through the upper and lower surface of one emulsion plate. The presented distributions demonstrate irregularity in the density of the reconstructed tracks in the XY plane. Irregularity in the density of the reconstructed tracks can be associated, in particular, with characteristics of the emulsion, including a large fog (spontaneously developed silver grains) in one of the layers of the double-sided emulsion film. As a result of different number of the reconstructed grains forming a track on each of the emulsion layers, the number of tracks passing through the entire emulsion volume may include a number of background tracks. Another reason for the observed irregularity in the density of the reconstructed tracks may be the forced arrangement of detectors under conditions that are unfavorable for the exposure, when the water level at the detector location was  
Figure 6: Two-dimensional distributions of reconstructed base tracks in the XY plane in two adjacent films ((**a**)—plate 1, (**b**)—plate 2) of the same detector. The scanned emulsion area: from 10 mm to 110 mm on the X axis; from 10 mm to 90 mm on the Y axis. about 30 cm from the floor.