2024ApJ...970..134M

We use data from the Tibet ASγ experiment for 4 teraelectronvolt (TeV) cosmic rays as an example to perform a nonlinear interstellar distribution model regression according to the way the observed anisotropy is typically presented, from which we extract normalization factors that allow us to obtain a true relative intensity sky map from the measurements. By using various test statistics, we show that the nonlinear fit significantly outperforms the direct linear fit in its ability to model cosmic-ray anisotropy. The procedure also allows us to produce normalization constants that can trace minute latitudinal variations of experimental response to cosmic-ray intensity. Applying the correction of the latitudinal response function to the Tibet ASγ data, we generate a sky map of true relative intensity. As a result, we observe that the measured and corrected sky maps show significant differences in intensity and angular spectral power. Our full anisotropy sky map of true relative intensity contradicts the assumption that the latitudinal variation in longitudinally averaged flux is negligible. The result further confirms that TeV cosmic-ray anisotropy is dominated by a dipole (ℓ = 1) aligned with the interstellar magnetic field's direction. Our results also confirm the existence of much weaker middle-scale interstellar anisotropy between ℓ = 2 and ℓ = 13.