Fig. 6. The matrix is determined using the TAUOLA decay model MC simulation. Clear correlations between the measured and generated values are seen along the diagonal parts of ˆ A . For the unfolding algorithm, we follow the method of the ALEPH Collaboration, which is based on the Singular Value Decomposition technique for matrix inversion [21] and advanced 14 FIG. 6: Response matrix ˆ A obtained by merging 16 sub-response matrices that represent the efficiency or fake rates among the four three-prong modes. The bin numbers correspond to the following ranges of the invariant masses: (1) from 0 to 138 , 0.41 GeV/ c 2 to 1.8 GeV/ c 2 of M ( πππ ); (2) from 139 to 241 , 0.77 GeV/ c 2 to 1.8 GeV/ c 2 of M ( Kππ ); (3) from 242 to 309 , 1.12 GeV/ c 2 to 1.8 GeV/ c 2 of M ( KKπ ); and (4) from 310 to 341 , 1.48 GeV/ c 2 to 1.8 GeV/ c 2 of M ( KKK ). All bin sizes are 10 MeV/ c $^{2}$. -7 10$^{-6}$ 10$^{-5}$ 10$^{-4}$ 10$^{-3}$ 10$^{-2}$ 10$^{-1}$ 10 0.41 1.80 0.77 1.80 1.12 1.80 1.48 1.80 0.41 1.80 0.77 1.80 1.12 1.80 1.48 1.80 ) ( gen M ) (K gen M ) (KK gen M (KKK) gen M ) ( rec M ) (K rec M ) (KK rec M (KKK) rec M ] 2 c [GeV/ ] 2 c [GeV/ Bin numbers of M(generated) Bin numbers of M(reconstructed) 0 100 200 300 0 100 200 300