Distribution of p-value for the scanned parameter combinations presented in Table 6, with different values of η parameter, using the H.E.S.S. and HAWC datasets in combination with E_c fixed to 1 PeV. This distribution is presented as violin plot, illustrating kernel probability density. The width of the shaded area represents the proportion of the number of fits located at these position. The upper and lower bar represent the maximum and minimum values. Models with p-value <0.003 are excluded.