×
近期发现有不法分子冒充我刊与作者联系,借此进行欺诈等不法行为,请广大作者加以鉴别,如遇诈骗行为,请第一时间与我刊编辑部联系确认(《中国物理C》(英文)编辑部电话:010-88235947,010-88236950),并作报警处理。
本刊再次郑重声明:
(1)本刊官方网址为cpc.ihep.ac.cn和https://iopscience.iop.org/journal/1674-1137
(2)本刊采编系统作者中心是投稿的唯一路径,该系统为ScholarOne远程稿件采编系统,仅在本刊投稿网网址(https://mc03.manuscriptcentral.com/cpc)设有登录入口。本刊不接受其他方式的投稿,如打印稿投稿、E-mail信箱投稿等,若以此种方式接收投稿均为假冒。
(3)所有投稿均需经过严格的同行评议、编辑加工后方可发表,本刊不存在所谓的“编辑部内部征稿”。如果有人以“编辑部内部人员”名义帮助作者发稿,并收取发表费用,均为假冒。
                  
《中国物理C》(英文)编辑部
2024年10月30日

Amplitude analysis of the decays ${\boldsymbol D^{\bf 0}\bf\to\pi^+\pi^-\pi^+\pi^-} $ and ${\boldsymbol D^{\bf 0}\bf\to\pi^+\pi^-\pi^0\pi^0} $

Figures(7) / Tables(20)

Get Citation
M. Ablikim, M. N. Achasov, P. Adlarson, O. Afedulidis, X. C. Ai, R. Aliberti, A. Amoroso, Q. An, Y. Bai, O. Bakina, I. Balossino, Y. Ban, H.-R. Bao, V. Batozskaya, K. Begzsuren, N. Berger, M. Berlowski, M. Bertani, D. Bettoni, F. Bianchi, E. Bianco, A. Bortone, I. Boyko, R. A. Briere, A. Brueggemann, H. Cai, X. Cai, A. Calcaterra, G. F. Cao, N. Cao, S. A. Cetin, J. F. Chang, W. L. Chang, G. R. Che, G. Chelkov, C. Chen, C. H. Chen, Chao Chen, G. Chen, H. S. Chen, M. L. Chen, S. J. Chen, S. L. Chen, S. M. Chen, T. Chen, X. R. Chen, X. T. Chen, Y. B. Chen, Y. Q. Chen, Z. J. Chen, Z. Y. Chen, S. K. Choi, X. Chu, G. Cibinetto, F. Cossio, J. J. Cui, H. L. Dai, J. P. Dai, A. Dbeyssi, R. E. de Boer, D. Dedovich, C. Q. Deng, Z. Y. Deng, A. Denig, I. Denysenko, M. Destefanis, F. De Mori, B. Ding, X. X. Ding, Y. Ding, Y. Ding, J. Dong, L. Y. Dong, M. Y. Dong, X. Dong, M. C. Du, S. X. Du, Z. H. Duan, P. Egorov, Y. H. Fan, J. Fang, J. Fang, S. S. Fang, W. X. Fang, Y. Fang, Y. Q. Fang, R. Farinelli, L. Fava, F. Feldbauer, G. Felici, C. Q. Feng, J. H. Feng, Y. T. Feng, K. Fischer, M. Fritsch, C. D. Fu, J. L. Fu, Y. W. Fu, H. Gao, Y. N. Gao, Yang Gao, S. Garbolino, I. Garzia, P. T. Ge, Z. W. Ge, C. Geng, E. M. Gersabeck, A. Gilman, K. Goetzen, L. Gong, W. X. Gong, W. Gradl, S. Gramigna, M. Greco, M. H. Gu, Y. T. Gu, C. Y. Guan, Z. L. Guan, A. Q. Guo, L. B. Guo, M. J. Guo, R. P. Guo, Y. P. Guo, A. Guskov, J. Gutierrez, K. L. Han, T. T. Han, X. Q. Hao, F. A. Harris, K. K. He, K. L. He, F. H. Heinsius, C. H. Heinz, Y. K. Heng, C. Herold, T. Holtmann, P. C. Hong, G. Y. Hou, X. T. Hou, Y. R. Hou, Z. L. Hou, B. Y. Hu, H. M. Hu, J. F. Hu, T. Hu, Y. Hu, G. S. Huang, K. X. Huang, L. Q. Huang, X. T. Huang, Y. P. Huang, T. Hussain, F. H\"olzken, N. H\"usken, N. in der Wiesche, M. Irshad, J. Jackson, S. Janchiv, J. H. Jeong, Q. Ji, Q. P. Ji, W. Ji, X. B. Ji, X. L. Ji, Y. Y. Ji, X. Q. Jia, Z. K. Jia, D. Jiang, H. B. Jiang, P. C. Jiang, S. S. Jiang, T. J. Jiang, X. S. Jiang, Y. Jiang, J. B. Jiao, J. K. Jiao, Z. Jiao, S. Jin, Y. Jin, M. Q. Jing, X. M. Jing, T. Johansson, S. Kabana, N. Kalantar-Nayestanaki, X. L. Kang, X. S. Kang, M. Kavatsyuk, B. C. Ke, V. Khachatryan, A. Khoukaz, R. Kiuchi, O. B. Kolcu, B. Kopf, M. Kuessner, X. Kui, A. Kupsc, W. K\"uhn, J. J. Lane, P. Larin, L. Lavezzi, T. T. Lei, Z. H. Lei, H. Leithoff, M. Lellmann, T. Lenz, C. Li, C. Li, C. H. Li, Cheng Li, D. M. Li, F. Li, G. Li, H. Li, H. B. Li, H. J. Li, H. N. Li, Hui Li, J. R. Li, J. S. Li, K. Li, L. J. Li, L. K. Li, Lei Li, M. H. Li, P. R. Li, Q. M. Li, Q. X. Li, R. Li, S. X. Li, T. Li, W. D. Li, W. G. Li, X. Li, X. H. Li, X. L. Li, X. Y. Li, Y. G. Li, Z. J. Li, Z. X. Li, C. Liang, H. Liang, H. Liang, Y. F. Liang, Y. T. Liang, G. R. Liao, L. Z. Liao, Y. P. Liao, J. Libby, A. Limphirat, D. X. Lin, T. Lin, B. J. Liu, B. X. Liu, C. Liu, C. X. Liu, F. Liu, F. H. Liu, Feng Liu, G. M. Liu, H. Liu, H. B. Liu, H. H. Liu, H. M. Liu, Huihui Liu, J. B. Liu, J. Y. Liu, K. Liu, K. Y. Liu, Ke Liu, L. Liu, L. C. Liu, Lu Liu, M. H. Liu, P. L. Liu, Q. Liu, S. B. Liu, T. Liu, W. K. Liu, W. M. Liu, X. Liu, X. Liu, Y. Liu, Y. Liu, Y. B. Liu, Z. A. Liu, Z. D. Liu, Z. Q. Liu, X. C. Lou, F. X. Lu, H. J. Lu, J. G. Lu, X. L. Lu, Y. Lu, Y. P. Lu, Z. H. Lu, C. L. Luo, M. X. Luo, T. Luo, X. L. Luo, X. R. Lyu, Y. F. Lyu, F. C. Ma, H. Ma, H. L. Ma, J. L. Ma, L. L. Ma, M. M. Ma, Q. M. Ma, R. Q. Ma, X. T. Ma, X. Y. Ma, Y. Ma, Y. M. Ma, F. E. Maas, M. Maggiora, S. Malde, A. Mangoni, Y. J. Mao, Z. P. Mao, S. Marcello, Z. X. Meng, J. G. Messchendorp, G. Mezzadri, H. Miao, T. J. Min, R. E. Mitchell, X. H. Mo, B. Moses, N. Yu. Muchnoi, J. Muskalla, Y. Nefedov, F. Nerling, I. B. Nikolaev, Z. Ning, S. Nisar, Q. L. Niu, W. D. Niu, Y. Niu, S. L. Olsen, Q. Ouyang, S. Pacetti, X. Pan, Y. Pan, A. Pathak, P. Patteri, Y. P. Pei, M. Pelizaeus, H. P. Peng, Y. Y. Peng, K. Peters, J. L. Ping, R. G. Ping, S. Plura, V. Prasad, F. Z. Qi, H. Qi, H. R. Qi, M. Qi, T. Y. Qi, S. Qian, W. B. Qian, C. F. Qiao, J. J. Qin, L. Q. Qin, X. S. Qin, Z. H. Qin, J. F. Qiu, S. Q. Qu, Z. H. Qu, C. F. Redmer, K. J. Ren, A. Rivetti, M. Rolo, G. Rong, Ch. Rosner, S. N. Ruan, N. Salone, A. Sarantsev, Y. Schelhaas, K. Schoenning, M. Scodeggio, K. Y. Shan, W. Shan, X. Y. Shan, J. F. Shangguan, L. G. Shao, M. Shao, C. P. Shen, H. F. Shen, W. H. Shen, X. Y. Shen, B. A. Shi, H. C. Shi, J. L. Shi, J. Y. Shi, Q. Q. Shi, R. S. Shi, S. Y. Shi, X. Shi, X. D. Shi, J. J. Song, T. Z. Song, W. M. Song, Y. J. Song, Y. X. Song, S. Sosio, S. Spataro, F. Stieler, Y. J. Su, G. B. Sun, G. X. Sun, H. Sun, H. K. Sun, J. F. Sun, K. Sun, L. Sun, S. S. Sun, T. Sun, W. Y. Sun, Y. Sun, Y. J. Sun, Y. Z. Sun, Z. Q. Sun, Z. T. Sun, C. J. Tang, G. Y. Tang, J. Tang, Y. A. Tang, L. Y. Tao, Q. T. Tao, M. Tat, J. X. Teng, V. Thoren, W. H. Tian, Y. Tian, Z. F. Tian, I. Uman, Y. Wan, S. J. Wang, B. Wang, B. L. Wang, Bo Wang, D. Y. Wang, F. Wang, H. J. Wang, J. P. Wang, K. Wang, L. L. Wang, M. Wang, Meng Wang, N. Y. Wang, S. Wang, S. Wang, T. Wang, T. J. Wang, W. Wang, W. Wang, W. P. Wang, X. Wang, X. F. Wang, X. J. Wang, X. L. Wang, X. N. Wang, Y. Wang, Y. D. Wang, Y. F. Wang, Y. L. Wang, Y. N. Wang, Y. Q. Wang, Yaqian Wang, Yi Wang, Z. Wang, Z. L. Wang, Z. Y. Wang, Ziyi Wang, D. Wei, D. H. Wei, F. Weidner, S. P. Wen, Y. R. Wen, U. Wiedner, G. Wilkinson, M. Wolke, L. Wollenberg, C. Wu, J. F. Wu, L. H. Wu, L. J. Wu, X. Wu, X. H. Wu, Y. Wu, Y. H. Wu, Y. J. Wu, Z. Wu, L. Xia, X. M. Xian, B. H. Xiang, T. Xiang, D. Xiao, G. Y. Xiao, S. Y. Xiao, Y. L. Xiao, Z. J. Xiao, C. Xie, X. H. Xie, Y. Xie, Y. G. Xie, Y. H. Xie, Z. P. Xie, T. Y. Xing, C. F. Xu, C. J. Xu, G. F. Xu, H. Y. Xu, Q. J. Xu, Q. N. Xu, W. Xu, W. L. Xu, X. P. Xu, Y. C. Xu, Z. P. Xu, Z. S. Xu, F. Yan, L. Yan, W. B. Yan, W. C. Yan, X. Q. Yan, H. J. Yang, H. L. Yang, H. X. Yang, T. Yang, Y. Yang, Y. F. Yang, Y. F. Yang, Y. X. Yang, Z. W. Yang, Z. P. Yao, M. Ye, M. H. Ye, J. H. Yin, Z. Y. You, B. X. Yu, C. X. Yu, G. Yu, J. S. Yu, T. Yu, X. D. Yu, C. Z. Yuan, J. Yuan, L. Yuan, S. C. Yuan, Y. Yuan, Z. Y. Yuan, C. X. Yue, A. A. Zafar, F. R. Zeng, S. H. Zeng, X. Zeng, Y. Zeng, Y. J. Zeng, Y. J. Zeng, X. Y. Zhai, Y. C. Zhai, Y. H. Zhan, A. Q. Zhang, B. L. Zhang, B. X. Zhang, D. H. Zhang, G. Y. Zhang, H. Zhang, H. C. Zhang, H. H. Zhang, H. H. Zhang, H. Q. Zhang, H. Y. Zhang, J. Zhang, J. Zhang, J. J. Zhang, J. L. Zhang, J. Q. Zhang, J. W. Zhang, J. X. Zhang, J. Y. Zhang, J. Z. Zhang, Jianyu Zhang, L. M. Zhang, Lei Zhang, P. Zhang, Q. Y. Zhang, S. H. Zhang, Shulei Zhang, X. D. Zhang, X. M. Zhang, X. Y. Zhang, Y. Zhang, Y. Zhang, Y. T. Zhang, Y. H. Zhang, Y. M. Zhang, Yan Zhang, Z. D. Zhang, Z. H. Zhang, Z. L. Zhang, Z. Y. Zhang, Z. Y. Zhang, G. Zhao, J. Y. Zhao, J. Z. Zhao, L. Zhao, Lei Zhao, M. G. Zhao, R. P. Zhao, S. J. Zhao, Y. B. Zhao, Y. X. Zhao, Z. G. Zhao, A. Zhemchugov, B. Zheng, J. P. Zheng, W. J. Zheng, Y. H. Zheng, B. Zhong, X. Zhong, H. Zhou, J. Y. Zhou, L. P. Zhou, X. Zhou, X. K. Zhou, X. R. Zhou, X. Y. Zhou, Y. Z. Zhou, J. Zhu, K. Zhu, K. J. Zhu, L. Zhu, L. X. Zhu, S. H. Zhu, S. Q. Zhu, T. J. Zhu, W. D. Zhu, Y. C. Zhu, Z. A. Zhu, J. H. Zou, J. Zu and (BESIII Collaboration). Amplitude analysis of the decays ${\boldsymbol D^{\bf 0}\bf\to\pi^+\pi^-\pi^+\pi^-} $ and ${\boldsymbol D^{\bf 0}\bf\to\pi^+\pi^-\pi^0\pi^0} $[J]. Chinese Physics C. doi: 10.1088/1674-1137/ad3d4d
M. Ablikim, M. N. Achasov, P. Adlarson, O. Afedulidis, X. C. Ai, R. Aliberti, A. Amoroso, Q. An, Y. Bai, O. Bakina, I. Balossino, Y. Ban, H.-R. Bao, V. Batozskaya, K. Begzsuren, N. Berger, M. Berlowski, M. Bertani, D. Bettoni, F. Bianchi, E. Bianco, A. Bortone, I. Boyko, R. A. Briere, A. Brueggemann, H. Cai, X. Cai, A. Calcaterra, G. F. Cao, N. Cao, S. A. Cetin, J. F. Chang, W. L. Chang, G. R. Che, G. Chelkov, C. Chen, C. H. Chen, Chao Chen, G. Chen, H. S. Chen, M. L. Chen, S. J. Chen, S. L. Chen, S. M. Chen, T. Chen, X. R. Chen, X. T. Chen, Y. B. Chen, Y. Q. Chen, Z. J. Chen, Z. Y. Chen, S. K. Choi, X. Chu, G. Cibinetto, F. Cossio, J. J. Cui, H. L. Dai, J. P. Dai, A. Dbeyssi, R. E. de Boer, D. Dedovich, C. Q. Deng, Z. Y. Deng, A. Denig, I. Denysenko, M. Destefanis, F. De Mori, B. Ding, X. X. Ding, Y. Ding, Y. Ding, J. Dong, L. Y. Dong, M. Y. Dong, X. Dong, M. C. Du, S. X. Du, Z. H. Duan, P. Egorov, Y. H. Fan, J. Fang, J. Fang, S. S. Fang, W. X. Fang, Y. Fang, Y. Q. Fang, R. Farinelli, L. Fava, F. Feldbauer, G. Felici, C. Q. Feng, J. H. Feng, Y. T. Feng, K. Fischer, M. Fritsch, C. D. Fu, J. L. Fu, Y. W. Fu, H. Gao, Y. N. Gao, Yang Gao, S. Garbolino, I. Garzia, P. T. Ge, Z. W. Ge, C. Geng, E. M. Gersabeck, A. Gilman, K. Goetzen, L. Gong, W. X. Gong, W. Gradl, S. Gramigna, M. Greco, M. H. Gu, Y. T. Gu, C. Y. Guan, Z. L. Guan, A. Q. Guo, L. B. Guo, M. J. Guo, R. P. Guo, Y. P. Guo, A. Guskov, J. Gutierrez, K. L. Han, T. T. Han, X. Q. Hao, F. A. Harris, K. K. He, K. L. He, F. H. Heinsius, C. H. Heinz, Y. K. Heng, C. Herold, T. Holtmann, P. C. Hong, G. Y. Hou, X. T. Hou, Y. R. Hou, Z. L. Hou, B. Y. Hu, H. M. Hu, J. F. Hu, T. Hu, Y. Hu, G. S. Huang, K. X. Huang, L. Q. Huang, X. T. Huang, Y. P. Huang, T. Hussain, F. H\"olzken, N. H\"usken, N. in der Wiesche, M. Irshad, J. Jackson, S. Janchiv, J. H. Jeong, Q. Ji, Q. P. Ji, W. Ji, X. B. Ji, X. L. Ji, Y. Y. Ji, X. Q. Jia, Z. K. Jia, D. Jiang, H. B. Jiang, P. C. Jiang, S. S. Jiang, T. J. Jiang, X. S. Jiang, Y. Jiang, J. B. Jiao, J. K. Jiao, Z. Jiao, S. Jin, Y. Jin, M. Q. Jing, X. M. Jing, T. Johansson, S. Kabana, N. Kalantar-Nayestanaki, X. L. Kang, X. S. Kang, M. Kavatsyuk, B. C. Ke, V. Khachatryan, A. Khoukaz, R. Kiuchi, O. B. Kolcu, B. Kopf, M. Kuessner, X. Kui, A. Kupsc, W. K\"uhn, J. J. Lane, P. Larin, L. Lavezzi, T. T. Lei, Z. H. Lei, H. Leithoff, M. Lellmann, T. Lenz, C. Li, C. Li, C. H. Li, Cheng Li, D. M. Li, F. Li, G. Li, H. Li, H. B. Li, H. J. Li, H. N. Li, Hui Li, J. R. Li, J. S. Li, K. Li, L. J. Li, L. K. Li, Lei Li, M. H. Li, P. R. Li, Q. M. Li, Q. X. Li, R. Li, S. X. Li, T. Li, W. D. Li, W. G. Li, X. Li, X. H. Li, X. L. Li, X. Y. Li, Y. G. Li, Z. J. Li, Z. X. Li, C. Liang, H. Liang, H. Liang, Y. F. Liang, Y. T. Liang, G. R. Liao, L. Z. Liao, Y. P. Liao, J. Libby, A. Limphirat, D. X. Lin, T. Lin, B. J. Liu, B. X. Liu, C. Liu, C. X. Liu, F. Liu, F. H. Liu, Feng Liu, G. M. Liu, H. Liu, H. B. Liu, H. H. Liu, H. M. Liu, Huihui Liu, J. B. Liu, J. Y. Liu, K. Liu, K. Y. Liu, Ke Liu, L. Liu, L. C. Liu, Lu Liu, M. H. Liu, P. L. Liu, Q. Liu, S. B. Liu, T. Liu, W. K. Liu, W. M. Liu, X. Liu, X. Liu, Y. Liu, Y. Liu, Y. B. Liu, Z. A. Liu, Z. D. Liu, Z. Q. Liu, X. C. Lou, F. X. Lu, H. J. Lu, J. G. Lu, X. L. Lu, Y. Lu, Y. P. Lu, Z. H. Lu, C. L. Luo, M. X. Luo, T. Luo, X. L. Luo, X. R. Lyu, Y. F. Lyu, F. C. Ma, H. Ma, H. L. Ma, J. L. Ma, L. L. Ma, M. M. Ma, Q. M. Ma, R. Q. Ma, X. T. Ma, X. Y. Ma, Y. Ma, Y. M. Ma, F. E. Maas, M. Maggiora, S. Malde, A. Mangoni, Y. J. Mao, Z. P. Mao, S. Marcello, Z. X. Meng, J. G. Messchendorp, G. Mezzadri, H. Miao, T. J. Min, R. E. Mitchell, X. H. Mo, B. Moses, N. Yu. Muchnoi, J. Muskalla, Y. Nefedov, F. Nerling, I. B. Nikolaev, Z. Ning, S. Nisar, Q. L. Niu, W. D. Niu, Y. Niu, S. L. Olsen, Q. Ouyang, S. Pacetti, X. Pan, Y. Pan, A. Pathak, P. Patteri, Y. P. Pei, M. Pelizaeus, H. P. Peng, Y. Y. Peng, K. Peters, J. L. Ping, R. G. Ping, S. Plura, V. Prasad, F. Z. Qi, H. Qi, H. R. Qi, M. Qi, T. Y. Qi, S. Qian, W. B. Qian, C. F. Qiao, J. J. Qin, L. Q. Qin, X. S. Qin, Z. H. Qin, J. F. Qiu, S. Q. Qu, Z. H. Qu, C. F. Redmer, K. J. Ren, A. Rivetti, M. Rolo, G. Rong, Ch. Rosner, S. N. Ruan, N. Salone, A. Sarantsev, Y. Schelhaas, K. Schoenning, M. Scodeggio, K. Y. Shan, W. Shan, X. Y. Shan, J. F. Shangguan, L. G. Shao, M. Shao, C. P. Shen, H. F. Shen, W. H. Shen, X. Y. Shen, B. A. Shi, H. C. Shi, J. L. Shi, J. Y. Shi, Q. Q. Shi, R. S. Shi, S. Y. Shi, X. Shi, X. D. Shi, J. J. Song, T. Z. Song, W. M. Song, Y. J. Song, Y. X. Song, S. Sosio, S. Spataro, F. Stieler, Y. J. Su, G. B. Sun, G. X. Sun, H. Sun, H. K. Sun, J. F. Sun, K. Sun, L. Sun, S. S. Sun, T. Sun, W. Y. Sun, Y. Sun, Y. J. Sun, Y. Z. Sun, Z. Q. Sun, Z. T. Sun, C. J. Tang, G. Y. Tang, J. Tang, Y. A. Tang, L. Y. Tao, Q. T. Tao, M. Tat, J. X. Teng, V. Thoren, W. H. Tian, Y. Tian, Z. F. Tian, I. Uman, Y. Wan, S. J. Wang, B. Wang, B. L. Wang, Bo Wang, D. Y. Wang, F. Wang, H. J. Wang, J. P. Wang, K. Wang, L. L. Wang, M. Wang, Meng Wang, N. Y. Wang, S. Wang, S. Wang, T. Wang, T. J. Wang, W. Wang, W. Wang, W. P. Wang, X. Wang, X. F. Wang, X. J. Wang, X. L. Wang, X. N. Wang, Y. Wang, Y. D. Wang, Y. F. Wang, Y. L. Wang, Y. N. Wang, Y. Q. Wang, Yaqian Wang, Yi Wang, Z. Wang, Z. L. Wang, Z. Y. Wang, Ziyi Wang, D. Wei, D. H. Wei, F. Weidner, S. P. Wen, Y. R. Wen, U. Wiedner, G. Wilkinson, M. Wolke, L. Wollenberg, C. Wu, J. F. Wu, L. H. Wu, L. J. Wu, X. Wu, X. H. Wu, Y. Wu, Y. H. Wu, Y. J. Wu, Z. Wu, L. Xia, X. M. Xian, B. H. Xiang, T. Xiang, D. Xiao, G. Y. Xiao, S. Y. Xiao, Y. L. Xiao, Z. J. Xiao, C. Xie, X. H. Xie, Y. Xie, Y. G. Xie, Y. H. Xie, Z. P. Xie, T. Y. Xing, C. F. Xu, C. J. Xu, G. F. Xu, H. Y. Xu, Q. J. Xu, Q. N. Xu, W. Xu, W. L. Xu, X. P. Xu, Y. C. Xu, Z. P. Xu, Z. S. Xu, F. Yan, L. Yan, W. B. Yan, W. C. Yan, X. Q. Yan, H. J. Yang, H. L. Yang, H. X. Yang, T. Yang, Y. Yang, Y. F. Yang, Y. F. Yang, Y. X. Yang, Z. W. Yang, Z. P. Yao, M. Ye, M. H. Ye, J. H. Yin, Z. Y. You, B. X. Yu, C. X. Yu, G. Yu, J. S. Yu, T. Yu, X. D. Yu, C. Z. Yuan, J. Yuan, L. Yuan, S. C. Yuan, Y. Yuan, Z. Y. Yuan, C. X. Yue, A. A. Zafar, F. R. Zeng, S. H. Zeng, X. Zeng, Y. Zeng, Y. J. Zeng, Y. J. Zeng, X. Y. Zhai, Y. C. Zhai, Y. H. Zhan, A. Q. Zhang, B. L. Zhang, B. X. Zhang, D. H. Zhang, G. Y. Zhang, H. Zhang, H. C. Zhang, H. H. Zhang, H. H. Zhang, H. Q. Zhang, H. Y. Zhang, J. Zhang, J. Zhang, J. J. Zhang, J. L. Zhang, J. Q. Zhang, J. W. Zhang, J. X. Zhang, J. Y. Zhang, J. Z. Zhang, Jianyu Zhang, L. M. Zhang, Lei Zhang, P. Zhang, Q. Y. Zhang, S. H. Zhang, Shulei Zhang, X. D. Zhang, X. M. Zhang, X. Y. Zhang, Y. Zhang, Y. Zhang, Y. T. Zhang, Y. H. Zhang, Y. M. Zhang, Yan Zhang, Z. D. Zhang, Z. H. Zhang, Z. L. Zhang, Z. Y. Zhang, Z. Y. Zhang, G. Zhao, J. Y. Zhao, J. Z. Zhao, L. Zhao, Lei Zhao, M. G. Zhao, R. P. Zhao, S. J. Zhao, Y. B. Zhao, Y. X. Zhao, Z. G. Zhao, A. Zhemchugov, B. Zheng, J. P. Zheng, W. J. Zheng, Y. H. Zheng, B. Zhong, X. Zhong, H. Zhou, J. Y. Zhou, L. P. Zhou, X. Zhou, X. K. Zhou, X. R. Zhou, X. Y. Zhou, Y. Z. Zhou, J. Zhu, K. Zhu, K. J. Zhu, L. Zhu, L. X. Zhu, S. H. Zhu, S. Q. Zhu, T. J. Zhu, W. D. Zhu, Y. C. Zhu, Z. A. Zhu, J. H. Zou, J. Zu and (BESIII Collaboration). Amplitude analysis of the decays ${\boldsymbol D^{\bf 0}\bf\to\pi^+\pi^-\pi^+\pi^-} $ and ${\boldsymbol D^{\bf 0}\bf\to\pi^+\pi^-\pi^0\pi^0} $[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ad3d4d shu
Milestone
Received: 2023-12-06
Article Metric

Article Views(567)
PDF Downloads(13)
Cited by(0)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Amplitude analysis of the decays ${\boldsymbol D^{\bf 0}\bf\to\pi^+\pi^-\pi^+\pi^-} $ and ${\boldsymbol D^{\bf 0}\bf\to\pi^+\pi^-\pi^0\pi^0} $

  • 1. Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China
  • 2. Beihang University, Beijing 100191, China
  • 3. Bochum Ruhr-University, D-44780 Bochum, Germany
  • 4. Budker Institute of Nuclear Physics SB RAS (BINP), Novosibirsk 630090, Russia
  • 5. Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
  • 6. Central China Normal University, Wuhan 430079, China
  • 7. Central South University, Changsha 410083, China
  • 8. China Center of Advanced Science and Technology, Beijing 100190, China
  • 9. China University of Geosciences, Wuhan 430074, China
  • 10. Chung-Ang University, Seoul, 06974, Republic of Korea
  • 11. COMSATS University Islamabad, Lahore Campus, Defence Road, Off Raiwind Road, 54000 Lahore, Pakistan
  • 12. Fudan University, Shanghai 200433, China
  • 13. GSI Helmholtzcentre for Heavy Ion Research GmbH, D-64291 Darmstadt, Germany
  • 14. Guangxi Normal University, Guilin 541004, China
  • 15. Guangxi University, Nanning 530004, China
  • 16. Hangzhou Normal University, Hangzhou 310036, China
  • 17. Hebei University, Baoding 071002, China
  • 18. Helmholtz Institute Mainz, Staudinger Weg 18, D-55099 Mainz, Germany
  • 19. Henan Normal University, Xinxiang 453007, China
  • 20. Henan University, Kaifeng 475004, China
  • 21. Henan University of Science and Technology, Luoyang 471003, China
  • 22. Henan University of Technology, Zhengzhou 450001, China
  • 23. Huangshan College, Huangshan 245000, China
  • 24. Hunan Normal University, Changsha 410081, China
  • 25. Hunan University, Changsha 410082, China
  • 26. Indian Institute of Technology Madras, Chennai 600036, India
  • 27. Indiana University, Bloomington, Indiana 47405, USA
  • 28A. INFN Laboratori Nazionali di Frascati, INFN Laboratori Nazionali di Frascati, I-00044, Frascati, Italy
  • 28B. INFN Laboratori Nazionali di Frascati, INFN Sezione di Perugia, I-06100, Perugia, Italy
  • 28C. INFN Laboratori Nazionali di Frascati, University of Perugia, I-06100, Perugia, Italy
  • 29A. INFN Sezione di Ferrara, INFN Sezione di Ferrara, I-44122, Ferrara, Italy
  • 29B. INFN Sezione di Ferrara, University of Ferrara, I-44122, Ferrara, Italy
  • 30. Inner Mongolia University, Hohhot 010021, China
  • 31. Institute of Modern Physics, Lanzhou 730000, China
  • 32. Institute of Physics and Technology, Peace Avenue 54B, Ulaanbaatar 13330, Mongolia
  • 33. Instituto de Alta Investigaci\'on, Universidad de Tarapac\'a, Casilla 7D, Arica 1000000, Chile
  • 34. Jilin University, Changchun 130012, China
  • 35. Johannes Gutenberg University of Mainz, Johann-Joachim-Becher-Weg 45, D-55099 Mainz, Germany
  • 36. Joint Institute for Nuclear Research, 141980 Dubna, Moscow region, Russia
  • 37. Justus-Liebig-Universitaet Giessen, II. Physikalisches Institut, Heinrich-Buff-Ring 16, D-35392 Giessen, Germany
  • 38. Lanzhou University, Lanzhou 730000, China
  • 39. Liaoning Normal University, Dalian 116029, China
  • 40. Liaoning University, Shenyang 110036, China
  • 41. Nanjing Normal University, Nanjing 210023, China
  • 42. Nanjing University, Nanjing 210093, China
  • 43. Nankai University, Tianjin 300071, China
  • 44. National Centre for Nuclear Research, Warsaw 02-093, Poland
  • 45. North China Electric Power University, Beijing 102206, China
  • 46. Peking University, Beijing 100871, China
  • 47. Qufu Normal University, Qufu 273165, China
  • 48. Renmin University of China, Beijing 100872, China
  • 49. Shandong Normal University, Jinan 250014, China
  • 50. Shandong University, Jinan 250100, China
  • 51. Shanghai Jiao Tong University, Shanghai 200240, China
  • 52. Shanxi Normal University, Linfen 041004, China
  • 53. Shanxi University, Taiyuan 030006, China
  • 54. Sichuan University, Chengdu 610064, China
  • 55. Soochow University, Suzhou 215006, China
  • 56. South China Normal University, Guangzhou 510006, China
  • 57. Southeast University, Nanjing 211100, China
  • 58. State Key Laboratory of Particle Detection and Electronics, Beijing 100049, Hefei 230026, China
  • 59. Sun Yat-Sen University, Guangzhou 510275, China
  • 60. Suranaree University of Technology, University Avenue 111, Nakhon Ratchasima 30000, Thailand
  • 61. Tsinghua University, Beijing 100084, China
  • 62A. Turkish Accelerator Center Particle Factory Group, Istinye University, 34010, Istanbul, Turkey
  • 62B. Turkish Accelerator Center Particle Factory Group, Near East University, Nicosia, North Cyprus, 99138, Mersin 10, Turkey
  • 63. University of Chinese Academy of Sciences, Beijing 100049, China
  • 64. University of Groningen, NL-9747 AA Groningen, The Netherlands
  • 65. University of Hawaii, Honolulu, Hawaii 96822, USA
  • 66. University of Jinan, Jinan 250022, China
  • 67. University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
  • 68. University of Muenster, Wilhelm-Klemm-Strasse 9, 48149 Muenster, Germany
  • 69. University of Oxford, Keble Road, Oxford OX13RH, United Kingdom
  • 70. University of Science and Technology Liaoning, Anshan 114051, China
  • 71. University of Science and Technology of China, Hefei 230026, China
  • 72. University of South China, Hengyang 421001, China
  • 73. University of the Punjab, Lahore-54590, Pakistan
  • 74A. University of Turin and INFN, University of Turin, I-10125, Turin, Italy
  • 74B. University of Turin and INFN, University of Eastern Piedmont, I-15121, Alessandria, Italy
  • 74C. University of Turin and INFN, INFN, I-10125, Turin, Italy
  • 75. Uppsala University, Box 516, SE-75120 Uppsala, Sweden
  • 76. Wuhan University, Wuhan 430072, China
  • 77. Yantai University, Yantai 264005, China
  • 78. Yunnan University, Kunming 650500, China
  • 79. Zhejiang University, Hangzhou 310027, China
  • 80. Zhengzhou University, Zhengzhou 450001, China
  • a. Also at the Moscow Institute of Physics and Technology, Moscow 141700, Russia
  • b. Also at the Novosibirsk State University, Novosibirsk, 630090, Russia
  • c. Also at the NRC "Kurchatov Institute", PNPI, 188300, Gatchina, Russia
  • d. Also at Goethe University Frankfurt, 60323 Frankfurt am Main, Germany
  • e. Also at Key Laboratory for Particle Physics, Astrophysics and Cosmology, Ministry of Education; Shanghai Key Laboratory for Particle Physics and Cosmology; Institute of Nuclear and Particle Physics, Shanghai 200240, China
  • f. Also at Key Laboratory of Nuclear Physics and Ion-beam Application (MOE) and Institute of Modern Physics, Fudan University, Shanghai 200443, China
  • g. Also at State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing 100871, China
  • h. Also at School of Physics and Electronics, Hunan University, Changsha 410082, China
  • i. Also at Guangdong Provincial Key Laboratory of Nuclear Science, Institute of Quantum Matter, South China Normal University, Guangzhou 510006, China
  • j. Also at MOE Frontiers Science Center for Rare Isotopes, Lanzhou University, Lanzhou 730000, China
  • k. Also at Lanzhou Center for Theoretical Physics, Lanzhou University, Lanzhou 730000, China
  • l. Also at the Department of Mathematical Sciences, IBA, Karachi 75270, Pakistan
  • m. Also at Ecole Polytechnique Federale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland

Abstract: Using $ e^+e^- $ annihilation data corresponding to an integrated luminosity of 2.93 $ \rm fb^{-1} $ taken at the center-of-mass energy $ \sqrt{s}=3.773 $ GeV with the BESIII detector, a joint amplitude analysis is performed on the decays $ D^0\to\pi^+\pi^-\pi^+\pi^- $ and $ D^0\to\pi^+\pi^-\pi^0\pi^0 $(non-η). The fit fractions of individual components are obtained, and large interferences among the dominant components of the decays $ D^{0}\to a_{1}(1260)\pi $, $ D^{0}\to\pi(1300)\pi $, $ D^{0}\to\rho(770)\rho(770) $, and $ D^{0}\to2(\pi\pi)_{S} $ are observed in both channels. With the obtained amplitude model, the $ CP $-even fractions of $ D^0\to \pi^+\pi^-\pi^+\pi^- $ and $ D^0\to\pi^+\pi^-\pi^0\pi^0 $(non-η) are determined to be $ (75.2\; \pm\; 1.1_{\rm stat.}\; \pm\; 1.5_{\rm syst.})\ $% and $(68.9\; \pm\; 1.5_{\rm stat.}\; \pm\; 2.4_{\rm syst.})$%, respectively. The branching fractions of $ D^0\to \pi^+\pi^-\pi^+\pi^- $ and $ D^0\to\pi^+\pi^-\pi^0\pi^0 $(non-η) are measured to be $(0.688\; \pm\; 0.010_{\rm stat.}\; \pm\; 0.010_{\rm syst.})$% and $(0.951\; \pm\; 0.025_{\rm stat.}\; \pm\; 0.021_{\rm syst.})$%, respectively. The amplitude analysis provides an important model for the binning strategy in measuring the strong phase parameters of $ D^0 \to 4\pi $ when used to determine the CKM angle $ \gamma (\phi_{3}) $ via the $ B^{-}\to D K^{-} $ decay.

    HTML

    I.   INTRODUCTION
    • Precision measurements of the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix and the test of the unitarity of the CKM triangle [1, 2] are essential goals in the field of flavor physics. One of the three angles of the unitarity triangle, $ \gamma(\phi_{3})\equiv \mathrm{arg}(-V_{ud}V_{ub}^{*}/ V_{cd}V_{cb}^{*}) $, can be measured with the tree-level decay $ B^{\pm}\to DK^{\pm} $ through the interference between the $ B^{-}\to D^0 K^{-} (b\to c\bar{u}s) $ and $ B^{-}\to \bar{D}^0 K^{-} \;\; (b\to u\bar{c}s) $ amplitudes, which is one of the most important measurements of LHCb and Belle II. Several approaches have been proposed to measure the CKM angle γ via the decay $ B^{\pm}\to DK^{\pm} $ [35]. Here, the relative magnitude and phase between the $ D^0 $ and $ \bar{D}^0 $ mesons decay into the same final states, and the $ D^0 $ decay parameters ($ e.g. \;\; CP $-even fraction $ F_+ $) are the critical inputs. With more data collected by the LHCb and Belle II experiments in the future, the decay parameters of $ D^0 $ ($ \bar{D}^0 $) will become the dominant source of uncertainty in the γ measurement. Therefore, precision measurements of the $ D^0 $ decay parameters are urgently required to improve the precision of the γ measurement.

      The decay $ D^{0}\to4\pi $ is regarded as a sensitive mode [6, 7] to extract the CKM angle γ via $ B^{-}\to D K^{-} $. Their $ CP $-even fractions ($ F_+ $) and relatively strong phase parameters in the different phase space (PHSP) bins ($ c_i $/$ s_i $) serve as the direct inputs in the GLW [3] and BPGGSZ [5, 8] methods, respectively. A reliable decay amplitude model of $ D^{0}\to4\pi $ is critical to precisely extract $ F_+ $ and a model-independent $ c_i $/$ s_i $ [6], and to search for $ CP $ violation in $ D^0\to 4\pi $ [9]. Moreover, the four-body $ D^0 $ hadronic decays provide an excellent platform to study the two-body decays $ D^0\to VV $ and $ D^0\to AP $, where V, A, and P denote vector, axial-vector, and pseudo-scalar mesons, respectively. These decays thus enhance the understanding of the decay dynamics of the $ D^0 $ meson [10, 11].

      Currently, the experimental studies of $ D^0\to 4\pi $ are limited. The FOCUS experiment performed an amplitude analysis of $ D^0\to\pi^+\pi^-\pi^+\pi^- $ based on $ \sim6000 $ candidate events with a background fraction of $ \sim10 $% [12]. An amplitude analysis of $ D^0\to\pi^+\pi^-\pi^+\pi^- $ was also conducted with the CLEO-c data containing $ \sim7000 $ candidate events with a background fraction of $ \sim20 $% [13]. However, no amplitude analysis of $ D^0\to\pi^+\pi^-\pi^0\pi^0 $ has been performed yet. BESIII has collected 2.93 $ \rm{fb}^{-1} $ of $ e^+e^- $ collision data at the center-of-mass energy $\sqrt{s}= 3.773$ GeV, where the $ D\bar{D} $ pair is produced without any additional hadrons. This data sample provides an ideal environment for studying D meson decays with the double tag (DT) technique [14, 15]. In this method, a single tag (ST) candidate requires that only one D meson is reconstructed via the ST mode. A DT candidate requires that D and $ \bar{D} $ are reconstructed via the signal and ST modes, respectively. Based on this data sample and the DT method with three ST modes $ \bar{D}^{0}\to K^{+}\pi^- $, $ \bar{D}^{0}\to K^{+}\pi^-\pi^0 $, and $ \bar{D}^{0}\to K^{+}\pi^-\pi^+\pi^- $, we report a joint amplitude analysis of $ D^0\to\pi^+\pi^-\pi^+\pi^- $ and $D^0\to \pi^+\pi^-\pi^0\pi^0$(non-η). Furthermore, we determine the model-dependent $ CP $-even, absolute branching, and individual component fractions. Throughout this study, the charge-conjugated processes are always implied.

    II.   BESIII DETECTOR AND MONTE CARLO SIMULATION
    • The BESIII detector [16] records symmetric $ e^+e^- $ collisions provided by the BEPCII storage ring [17] in the center-of-mass energy range from 2.0 to 4.95 GeV, with a peak luminosity of $ 1 \times 10^{33}\;\text{cm}^{-2}\text{s}^{-1} $ achieved at $\sqrt{s} = 3.77\;\text{GeV}$. BESIII has collected large data samples in this energy region [18]. The cylindrical core of the $ \rm{BESIII} $ detector covers 93% of the full solid angle. It comprises a helium-based multilayer drift chamber (MDC), a plastic scintillator time-of-flight system (TOF), and a CsI(Tl) electromagnetic calorimeter (EMC),all enclosed in a superconducting solenoidal magnet providing a 1.0 T magnetic field. The solenoid is supported by an octagonal flux-return yoke with resistive plate counter muon identification modules interleaved with steel. The charged-particle momentum resolution at $ 1\; {\rm GeV}/c $ is $ 0.5 $%, and the $ {\rm d}E/{\rm d}x $ resolution is $ 6 $% for resulting Bhabha scattering electrons. The EMC measures photon energies with a resolution of $ 2.5 $% ($ 5 $%) at $ 1 $ GeV in the barrel (end cap) region. The time resolution in the TOF barrel region is 68 ps, while that in the end cap region is 110 ps.

      Simulated data samples produced with a geant4-based [19] Monte Carlo (MC) package, which includes the geometric description [20] of the BESIII detector and the detector response, are used to determine detection efficiencies and estimate backgrounds. The simulation models the beam energy spread and initial state radiation (ISR) in the $ e^+e^- $ annihilations with the generator kkmc [21, 22]. The inclusive MC sample including the production of $ D\bar{D} $ pairs (including quantum coherence for the neutral D channels), the non-$ D\bar{D} $ decays of the $ \psi(3770) $, the ISR production of the $ J/\psi $ and $ \psi(3686) $ states, and the continuum processes incorporated in kkmc are generated to estimate the background and ST efficiencies. All particle decays are modeled with evtgen [23, 24] using branching fractions either taken from the Particle Data Group (PDG) [25], when available, or otherwise estimated with lundcharm [26, 27]. Final state radiation from charged final state particles is incorporated using the photos package [28]. Signal MC samples of $ D^0\to\pi^+\pi^-\pi^+\pi^- $ and $ D^0\to\pi^+\pi^-\pi^0\pi^0 $ generated uniformly in PHSP are used to normalize probability density functions (PDFs) in the amplitude analysis, while those generated according to the amplitude analysis results are used to estimate the DT efficiencies.

    III.   EVENT SELECTION
    • Charged tracks detected in the MDC are required to be within a polar angle (θ) range of $ |\rm{cos\theta}|<0.93 $, where θ is defined with respect to the z axis, which is the symmetry axis of the MDC. The distance of the closest approach of these charged tracks to the interaction point must be less than 10 cm along the z axis and less than 1 cm in the transverse plane. Particle identification (PID) for charged tracks combines measurements of the specific ionization energy loss in the MDC (dE/dx) and the flight time in the TOF to form likelihoods $ \mathcal{L}(h)\; (h=K,\pi) $ for each hadron h hypothesis. The charged kaons and pions are identified by comparing the likelihoods for the kaon and pion hypotheses, $ \mathcal{L}(K)>\mathcal{L}(\pi) $ and $ \mathcal{L}(\pi)>\mathcal{L}(K) $, respectively.

      Photon candidates are identified using showers in the EMC. The deposited energy of each shower must be $ >25 $ MeV in the barrel region ($ |\cos \theta|< 0.80 $) and $ >50 $ MeV in the end cap region ($ 0.86 <|\cos \theta|< 0.92 $). To exclude showers originating from charged tracks, the angle subtended by the EMC shower and the position of the closest charged track at the EMC must be $ >10^{\circ} $, as measured from the interaction point. To suppress electronic noises and reject showers unrelated to the event, the difference between the EMC time and the event start time must be within [0, 700] ns. The $ \pi^0 $ candidates are reconstructed from pairs of photon candidates with the invariant mass being in the interval (0.115, 0.150) GeV/$ c^2 $. To improve the momentum resolution, a kinematic fit constraining the two-photon invariant mass to the known $ \pi^0 $ mass [25] is performed, and the four-momenta updated by this kinematic fit are used in the subsequent analysis.

      The signal candidates of $ D^0\to\pi^+\pi^-\pi^+\pi^- $ and $D^0\to \pi^+\pi^-\pi^0\pi^0$ are selected with the DT method [14, 15]. First, the ST $ \bar{D}^{0} $ mesons are reconstructed with the three hadronic decay modes $ \bar{D}^{0}\to K^{+}\pi^- $, $ \bar{D}^{0}\to K^{+}\pi^-\pi^0 $, and $\bar{D}^{0}\to K^{+}\pi^-\pi^+\pi^-$. Two kinematic variables, the energy difference with respect to the beam energy $ \Delta E $ and the beam energy constrained mass $ M_{\rm{bc}} $, are defined as

      $ \Delta E=E_{\bar{D}^{0}} - E_{\rm{beam}}, $

      (1)

      $ M_{\rm{bc}} =\sqrt{E_{\rm{beam}}^{2}-\vec{p}_{\bar{D}^{0}}^2}\; , $

      (2)

      where $ E_{\rm{beam}} $ is the beam energy and $ E_{\bar{D}^{0}} $ and $ \vec{p}_{\bar{D}^{0}} $ are the energy and momentum of the ST $ \bar{D}^{0} $ candidate in the $ e^+e^- $ center-of-mass frame. For multiple $ \bar{D}^{0} $ candidates in each ST mode, only the one with the smallest $ |\Delta E| $ is kept for further analysis. To reject the backgrounds from cosmic rays and Bhabha events in the ST mode $ \bar{D}^{0}\to K^{+}\pi^- $, the requirements described in Ref. [29] are applied. To reject the peaking background from the decay $ \bar{D}^{0}\to K^{+}K_{S}^{0}\pi^- $ in the ST mode $ \bar{D}^{0}\to K^{+}\pi^-\pi^+\pi^- $, the events with $|M(\pi^+\pi^-) -0.4976| < 0.03$ GeV/$ c^{2} $ are vetoed. To further reject combinatorial backgrounds, the $ \bar{D}^{0} $ candidates are required to have $ \Delta E $ within a given interval defined in Table 1, approximately three times the resolution for each ST mode. In the sub-sample containing ST candidates, the signal candidates of $ D^0\to\pi^+\pi^-\pi^+\pi^- $ and $ D^0\to\pi^+\pi^-\pi^0\pi^0 $ are reconstructed with the $ \pi^{\pm} $ and $ \pi^0 $ candidates, which have not been used in the ST side (namely DT thereafter). Similar kinematic variables $ \Delta E $ and $ M_{\rm{bc}} $ are formed for signal $ D^0 $ candidates, and the corresponding $ \Delta E $ requirements are listed in Table 1. For multiple signal $ D^0 $ candidates, only the one giving the smallest $ |\Delta E| $ is kept.

      Decay mode$ \Delta E $ /GeV
      $ \bar{D}^{0}\to K^{+}\pi^{-} $$ (-0.027,0.026) $
      $ \bar{D}^{0}\to K^{+}\pi^{-}\pi^{0} $$ (-0.057,0.043) $
      $ \bar{D}^{0}\to K^{+}\pi^{-}\pi^{+}\pi^{-} $$ (-0.020,0.018) $
      $ D^{0}\to \pi^{+}\pi^{-}\pi^{+}\pi^{-} $$ (-0.032,0.028) $
      $ D^{0}\to \pi^{+}\pi^{-}\pi^{0}\pi^{0} $$ (-0.066,0.041) $

      Table 1.  The $ \Delta E $ requirements for different decay modes.

      To improve the purity of signal candidates in the amplitude analysis, some further selection criteria are applied. The studies based on the inclusive MC sample indicate that most backgrounds from the continuum process $ e^+e^- \to q\bar{q} $ include a $ K_{S}^{0} $ in the final state. Therefore, common and secondary vertex fits are performed on all $ \pi^{+}\pi^{-} $ pairs in the event to reconstruct the $ K_{S}^{0} $ candidate. The candidate events are rejected if any $ K_{S}^{0} $ candidate with $ \pi^+\pi^- $ invariant mass exists within the interval $ |M(\pi^{+}\pi^{-})-0.4976|<0.03\; {\rm GeV}/c^{2} $ with a decay length greater than twice its resolution. The MC studies also show that most backgrounds from $e^+e^- \to D^+D^-$ contain decay $ D^{-}\to K^{+}\pi^{-}\pi^{-} $ owing to its large branching fraction and a similar topology as the signal. The candidate events are rejected if any $ K^{+}\pi^{-}\pi^{-} $ combinations exist with $ 1.863<M_{\rm{bc}}(K^{+}\pi^{-}\pi^{-})<1.878\; {\rm GeV}/c^{2} $ and $|\Delta E(K^{+}\pi^{-}\pi^{-})| < 0.03 \; \rm{GeV}$. The background from the process $ e^+e^- \to D^0\bar{D}^0 $ is also considered. Events with any $ \pi^{+}\pi^{-}\pi^{0} $ combinations satisfying $1.859 < M_{\rm{bc}}(\pi^{+}\pi^{-}\pi^{0})< 1.873 {\rm GeV}/c^{2}$ and $ -0.057<\Delta E(\pi^{+}\pi^{-}\pi^{0})<0.043 \; \rm{GeV} $ are rejected to eliminate the background of the decay $D^{0}\to \pi^{+}\pi^{-}\pi^{0}$ in the signal process $ D^{0}\to\pi^{+}\pi^{-}\pi^{0}\pi^{0} $. All the above backgrounds do not form peaks in the $ M_{\rm bc} $ distribution of the signal side. Additionally, some backgrounds have the same final states as those in the signal mode and produce peaks in the $ M_{\rm bc} $ distributions of both the ST and signal sides. To reject the peaking backgrounds from the decays $ D^{0}\to K_{S}^{0}(\rightarrow \pi^{+}\pi^{-})\pi^{+}\pi^{-} $ and $ D^{0}\to K^{-}(\rightarrow \pi^{+}\pi^{-}\pi^{-})\pi^{+} $ in the signal process $ D^0\to \pi^+\pi^-\pi^+\pi^- $, events with $ |M(\pi^{+}\pi^{-})-0.4976|< 0.03\; {\rm GeV}/c^{2} $ or $M(\pi^{+}\pi^{-}\pi^{-}) < 0.51 {\rm GeV}/c^{2}$ are vetoed. To reject the peaking backgrounds from the decays $ D^{0}\to K_{S}^{0} (\rightarrow \pi^{0}\pi^{0})\pi^{+}\pi^{-} $, $D^{0}\to K_{S}^{0}(\rightarrow \pi^{+}\pi^{-})\pi^{0}\pi^{0}$ and $ D^{0}\to K^{-}(\rightarrow \pi^{-}\pi^{0})\pi^{+}\pi^{0} $ in the signal process $ D^0\to \pi^{+}\pi^{-}\pi^0\pi^{0} $, events with $0.4376 < M(\pi^{0}\pi^{0}) <0.5276\; {\rm GeV}/c^{2}$ or $ |M(\pi^{+}\pi^{-})-0.4976|<0.03\; {\rm GeV}/c^{2} $ or $ 0.4677<M(\pi^{-}\pi^{0})< 0.5067\; {\rm GeV}/c^{2} $ are vetoed. Since the interference between the decays $ D^{0}\to \pi^{0}\eta(\rightarrow \pi^{+}\pi^{-}\pi^{0}) $ and $ D^{0}\to \pi^{+}\pi^{-}\pi^{0}\pi^{0} $ is negligible, events with $M(\pi^{+}\pi^{-}\pi^{0}) < 0.57\; {\rm GeV}/c^{2}$ are also vetoed in the signal process $ D^0\to \pi^{+}\pi^{-}\pi^0\pi^{0} $. Moreover, we reject events if the photon energy from $ \pi^{0} $ decays is below 50 MeV, as this indicates an incorrect choice of a soft photon.

    IV.   BACKGROUND STUDY AND SIGNAL EXTRACTION
    • Based on the above selection criteria, the two-dimensional (2D) distributions of $ M_{\rm bc}^{\rm sig} $ versus $ M_{\rm bc}^{\rm tag} $ of the survived events are shown in the left columns of Figs. 2 and 3, where $ M_{\rm bc}^{\rm sig} $ and $ M_{\rm bc}^{\rm tag} $ are for the signal and tag sides, respectively. Typically, the signal events are accumulated around the intersection of $ M_{\rm bc}^{\rm sig} $ and $ M_{\rm bc}^{\rm tag} $ at the $ D^0 $ nominal mass. The background, arising from miscombinations in the signal and tag sides (namely BKGI thereafter), is located at the diagonal band. Additionally, the backgrounds from $ e^+e^-\to D^0\bar{D}^0 $ but with the wrong reconstruction of the signal side exist (namely BKGII) or tag side (namely BKGIII), distributed as the vertical and horizontal bands, respectively. Detailed MC studies indicate that most of the backgrounds from $ e^+e^-\to D^0\bar{D}^0 $ do not form individual peaks in the distribution of $ M_{\rm bc}^{\rm sig} $ or $ M_{\rm bc}^{\rm tag} $. However, backgrounds from $ e^+e^-\to D^0\bar{D}^0 $ with $\bar{D}^0 \to K^+\pi^-\pi^0$ (namely BKGIV) or $ D^0\to \pi^+\pi^-\pi^0\pi^0 $ (namely BKGV) with the wrong reconstruction of $ \pi^0 $ exist, which are exactly consistent with the signal but with the wrong reconstruction of $ \pi^0 $ and produce a relatively broad peak in the $ M_{\rm bc}^{\rm tag} $ or $ M_{\rm bc}^{\rm sig} $ distribution, respectively. TBackgrounds from the decays $ D^0\to K_S^0\pi^+\pi^- $, $D^0\to K_S^0\pi^0\pi^0$, $ D^0\to K^-\pi^+\pi^0 $, and $ D^0\to K_S^0\eta^{\prime} $ (namely BKGVI), which have analogous final states as the signal with the specific decay modes of $ K_S^0\to\pi^+\pi^- $, $ K_S^0\to\pi^0\pi^0 $, $ K^-\to\pi^-\pi^0 $, and $ \eta^{\prime}\to\gamma\pi^+\pi^- $. These backgrounds are not directly distinguished from the signal in the $ M_{\rm bc}^{\rm sig} $ and $ M_{\rm bc}^{\rm tag} $ distributions and are estimated by MC simulation. The corresponding yields are summarized in Table 2.

      Signal mode$ \pi^{+}\pi^{-}\pi^{+}\pi^{-} $$ \pi^{+}\pi^{-}\pi^{0}\pi^{0} $
      Tag mode$ K^{+}\pi^{-} $$ K^{+}\pi^{-}\pi^{0} $$ K^{+}\pi^{-}\pi^{+}\pi^{-} $$ K^{+}\pi^{-} $$ K^{+}\pi^{-}\pi^{0} $$ K^{+}\pi^{-}\pi^{+}\pi^{-} $
      $ N(K_{S}^{0}\pi^{+}\pi^{-}) $$ 24.7 \pm 2.3 $$ 42.4 \pm 4.3 $$ 22.9 \pm 2.6 $$ 12.6 \pm 3.1 $$ 27.7 \pm 6.2 $$ 29.2 \pm 6.2 $
      $ N(K_{S}^{0}\pi^{0}\pi^{0}) $$ 1.7 \pm 0.4 $$ 4.3 \pm 0.9 $$ 2.4 \pm 0.6 $
      $ N(K^{-}\pi^{+}\pi^{0}) $$ 12.6 \pm 1.2 $$ 23.1 \pm 2.2 $$ 24.0 \pm 1.8 $
      $ N(K_{S}^{0}\eta') $$ 3.0 \pm 1.2 $$ 3.9 \pm 1.7 $$ 2.8 \pm 1.2 $

      Table 2.  The estimated numbers of peaking background events. The uncertainties include the statistical uncertainties of the estimated background yields and the uncertainties of the quoted branching fractions of different background processes.

      Figure 2.  (color online) The 2D distributions (left) of $ M_{\rm bc}^{\rm tag} $ versus $ M_{\rm bc}^{\rm sig} $ and the projections on $ M_{\rm bc}^{\rm tag} $ (middle) and $ M_{\rm bc}^{\rm sig} $ (right) of the 2D fits on the DT events tagged by $ \bar{D}^{0}\to K^{+}\pi^- $ (first row), $ \bar{D}^{0}\to K^{+}\pi^-\pi^0 $ (second row), and $ \bar{D}^{0}\to K^{+}\pi^-\pi^+\pi^- $ (third row) in the $ D^0\to\pi^+\pi^-\pi^+\pi^- $ decay. The dots with error bars are the data. The blue solid curves are the total fit results, and the red dashed curves show the signal. The green, blue, orange, and magenta dashed curves are BKGI, BKGIII, BKGVI, and the combination of BKGII and BKGV, respectively.

      Figure 3.  (color online) The 2D distributions (left) of $ M_{\rm bc}^{\rm tag} $ versus $ M_{\rm bc}^{\rm sig} $ and the projections on $ M_{\rm bc}^{\rm tag} $ (middle) and $ M_{\rm bc}^{\rm sig} $ (right) of the 2D fits on the DT events tagged by $ \bar{D}^{0}\to K^{+}\pi^- $ (first row), $ \bar{D}^{0}\to K^{+}\pi^-\pi^0 $ (second row) and $ \bar{D}^{0}\to K^{+}\pi^-\pi^+\pi^- $ (third row) in the $ D^0\to\pi^+\pi^-\pi^0\pi^0 $ decay. The dots with error bars are the data. The blue solid curves are the total fit results, and the red dashed curves show the signal. The green dashed curves are the combination of BKGI and BKGIV, the magenta dashed curves are the combination of BKGII and BKGV, the blue and orange dashed curves are BKGIII and BKGVI, respectively.

      To extract the ST and DT yields, unbinned maximum likelihood fits are performed on the $ M_{\rm bc}^{\rm tag} $ distribution of the remaining ST candidates and the 2D distributions of $ M_{\rm bc}^{\rm sig} $ versus $ M_{\rm bc}^{\rm tag} $ of the remaining DT events in individual tag modes. In the fit to the $ M_{\rm bc}^{\rm tag} $ distribution, the signal shape is described by the MC simulated shape of the truth-matched ST events, and the background shape is described by an ARGUS function with the cut-off parameter fixed at $ 1.8865 $ GeV/$ c^2 $ [30]. For the $ \bar{D}^0\to K^+\pi^-\pi^0 $ tag mode, an additional peaking background is considered, representing the wrong reconstruction of $ \pi^0 $, as discussed above. This peaking background shape is described by the MC-simulated shape of the truth-matched ST events convolved with a bifurcated Gaussian function with fixed parameters obtained from the fit to the corresponding simulated background events. The results of the fits to the $ M_{\rm bc}^{\rm tag} $ distributions are shown in Fig. 1.

      Figure 1.  (color online) The $ M_{\rm bc}^{\rm tag} $ distributions of the ST candidates for $ \bar{D}^{0}\to K^+\pi^- $ (left), $ \bar{D}^{0}\to K^+\pi^-\pi^0 $ (middle) and $ \bar{D}^{0}\to K^+\pi^-\pi^+\pi^- $ (right). The dots with error bars are the data, the blue solid curves are the total fit results, and the red and blue dashed curves are the signal and background, respectively.

      In the fit to the 2D distribution of $ M_{\rm bc}^{\rm sig} $ (labelled as $ y_1 $) versus $ M_{\rm bc}^{\rm tag} $ (labelled as $ y_2 $), the signal is described by

      $ \begin{array}{*{20}{l}} S(y_1,y_2) \otimes G(y_1;\mu_{y_1},\sigma_{y_1}) \otimes G(y_2;\mu_{y_2},\sigma_{y_2}), \end{array} $

      where $ S(y_1,y_2) $ is the signal MC shape derived from the truth-matched events by using a smoothed 2D histogram, and $ G(\mu,\sigma) $ is the Gaussian function describing the resolution difference between data and MC simulation by fixing the parameters $ \mu_{y_{1(2)}} $ and $ \sigma_{y_{1(2)}} $ to those from the fits on the one-dimensional $ M_{\rm bc} $ distributions.

      The BKGI is described by

      $ \begin{array}{*{20}{l}} &T(y_1-y_2; \mu, \sigma(y_1+y_2), n) \times A(y_1; m_{y_1}, z'_{y_1}, \rho'_{y_1}) \\ &\times A(y_2; m_{y_2}, z'_{y_2}, \rho'_{y_2}), \end{array} $

      where T is the student's function defined as $T(y; \mu, \sigma, n) = \dfrac{\Gamma(n/2+0.5)}{\sigma\sqrt{n\pi}\Gamma(n/2)}\left[1 + \dfrac{1}{n}\left(\dfrac{y-\mu}{\sigma}\right)^{2}\right]^{-\frac{n+1}{2}}$, $\sigma(y_1 + y_2) = \sigma_0 + \sigma_{1}(y_1 +y_2 -m_{y_1} -m_{y_2})$ and A is the ARGUS function defined as $ A(y;m,z,\rho) = y\left(1-\dfrac{y^2}{m^2}\right)^{\rho}{\rm e}^{z(1-\frac{y^2}{m^2})} $. In the fit, only the cut-off parameter $ m_{y_{1(2)}} $ is fixed to $ 1.8865\; {\rm GeV}/c^{2} $, the other parameters are free to vary.

      The BKGII and BKGIII are described by

      $ \begin{array}{*{20}{l}} &[S(y_{2(1)}) \otimes G(y_{2(1)};\mu_{y_{2(1)}},\sigma_{y_{2(1)}})] \\ &\times A(y_{1(2)}; m_{y_{1(2)}}, z_{y_{1(2)}}, \rho_{y_{1(2)}}), \end{array} $

      where A is the ARGUS function as described above, $ S(y_{2(1)}) $ is the projection of $ S(y_1,y_2) $ on $ y_{2(1)} $, and G is the same Gaussian function, as described above. In the fit, the parameter $ m_{y_{1(2)}} $ is fixed to $ 1.8865\; {\rm GeV}/c^{2} $, $ \rho_{y_{1(2)}} $ is fixed to the values obtained from the fits to the inclusive MC sample, and $ z_{y_{1(2)}} $ is a free parameter.

      The BKGIV is described by

      $ \begin{array}{*{20}{l}} B_1(y_1,y_2) \otimes G(y_1;\mu_{y_1},\sigma_{y_1}) \otimes G(y_2;\mu_{y_2},\sigma_{y_2}), \end{array} $

      where $ B_1(y_1,y_2) $ is the truth-matched signal MC shape $ S(y_1,y_2) $ convolved with a bifurcated Gaussian function in $ y_2 $ and G is the same Gaussian function, as described above. In the fit, the parameters of the bifurcated Gaussian function are fixed to the values obtained from the fit to the $ M_{\rm bc}^{\rm tag} $ distribution of the corresponding simulated background events.

      The BKGV in $ D^0\to\pi^+\pi^-\pi^0\pi^0 $ is described by

      $ \begin{array}{*{20}{l}} B_2(y_1,y_2) \otimes G(y_1;\mu_{y_1},\sigma_{y_1}) \otimes G(y_2;\mu_{y_2},\sigma_{y_2}), \end{array} $

      where $ B_2(y_1,y_2) $ is the truth-matched signal MC shape $ S(y_1,y_2) $ convolved with a bifurcated Gaussian function in $ y_1 $ and G is the same Gaussian function, as described above. In the fit, the parameters of the bifurcated Gaussian function are fixed to the values obtained from the fit to the $ M_{\rm bc}^{\rm sig} $ distribution of the corresponding simulated background events.

      The BKGVI is described by

      $ \begin{array}{*{20}{l}} B_3(y_1,y_2) \otimes G(y_1;\mu_{y_1},\sigma_{y_1}) \otimes G(y_2;\mu_{y_2},\sigma_{y_2}), \end{array} $

      where $ B_3(y_1,y_2) $ is the truth-matched MC shape $ S(y_1,y_2) $ convolved with a bifurcated Gaussian function in $ y_1 $ and G are the Gaussian functions. In the fit, the parameters of the bifurcated Gaussian functions are fixed to the value obtained from the fit to the $ M_{\rm bc}^{\rm sig} $ distribution of the corresponding simulated background events from the inclusive MC sample, and the yields are fixed to those summarized in Table 2.

      The projections of the 2D fits on $ M_{\rm bc}^{\rm tag} $ and $ M_{\rm bc}^{\rm sig} $ are shown in Figs. 2 and 3, respectively. The obtained ST and DT yields are summarized in Table 3.

      Tag mode$ \bar{D}^{0}\to K^+\pi^- $$ \bar{D}^{0}\to K^+\pi^-\pi^0 $$ \bar{D}^{0}\to K^+\pi^-\pi^+\pi^- $
      $ N^{\rm ST} $$ 549586 \pm 778 $$ 914531 \pm 1321 $$ 600316 \pm 856 $
      $ \epsilon^{\rm ST} $$ 0.6762 \pm 0.0004 $$ 0.2953 \pm 0.0001 $$ 0.3365 \pm 0.0002 $
      $ N^{\rm DT}_{\pi^+\pi^-\pi^+\pi^-} $$ 1719 \pm 42 $$ 2560 \pm 56 $$ 1520 \pm 43 $
      $ \epsilon^{\rm DT}_{\pi^+\pi^-\pi^+\pi^-} $$ 0.2792 \pm 0.0006 $$ 0.1187 \pm 0.0002 $$ 0.1221 \pm 0.0003 $
      $ N^{\rm DT}_{\pi^+\pi^-\pi^0\pi^0}\;(\rm non{\text -}\eta) $$ 721 \pm 32 $$ 917 \pm 39 $$ 562 \pm 29 $
      $ \epsilon^{\rm DT}_{\pi^+\pi^-\pi^0\pi^0}\;(\rm non{\text -}\eta) $$ 0.0818 \pm 0.0003 $$ 0.0319 \pm 0.0001 $$ 0.0334 \pm 0.0001 $
      $ \dfrac{2rR\rm{cos}\delta}{1+r^2} $$ -0.114 \pm 0.004 $ [31]$ -0.067 \pm 0.005 $ [32]$ -0.046^{+0.013}_{-0.011} $ [32]

      Table 3.  The ST and DT yields, efficiencies, and quantum correlated correction factors for three tag modes.

    V.   AMPLITUDE FORMULA
    • To improve the resolutions of kinematic variables, a one-constraint kinematic fit with the hypothesis of $ D^0\to 4 \pi $ by constraining the $ 4\pi $ invariant mass to the known $ D^0 $ mass [25] is performed on the candidate events, and the updated kinematic variables are used in the amplitude analysis. Using the GPUPWA framework [33], a joint amplitude analysis is performed on the candidate events of $ D^0\to\pi^+\pi^-\pi^+\pi^- $ and $ D^0\to\pi^+\pi^-\pi^0\pi^0 $(non-η). The general amplitude of the $ D^{0}\to f $ decay is given by

      $ A_{f}(p) = \sum\limits_{i} \Lambda_iU_i(p)\; , $

      (3)

      where $ U_i $ is the amplitude of the i-th intermediate process with a complex coupling factor $ \Lambda_i $ and p is a set of four-momenta of final states. For the $ \bar{D}^{0} $ decay amplitude, assuming $ CP $ conservation, the $ CP $-conjugate PHSP $ \bar{p} $ is defined by interchanging the charge of the final state and the three-momenta reversal. Then, the amplitude of the $ \bar{D}^{0}\to \bar{f} $ decay is given as

      $ \bar{A}_{\bar{f}}(p)= \sum\limits_{i}\Lambda_i\bar{U}_i(p)=\sum\limits_{i}\Lambda_iU_i(\bar{p})\; . $

      (4)

      Since the $ D^{0}\bar{D}^0 $ pair is produced in the $ \psi(3770) $ decay, quantum correlation between $ D^{0} $ and $ \bar{D}^0 $ needs to be considered. By ignoring the effects of $ CP $ violation and $ D^{0} $-$ \bar{D}^0 $ mixing, the observed differential cross section $ |M_f(p)|^{2} $ of $ D^0\to f $ and the tag mode $ \bar{D}^0\to \bar{g} $ is

      $ \begin{array}{*{20}{l}} |A_f(p)-r_{g}R_{g}{\rm e}^{-{\rm i} \delta_{g}}\bar{A}_f(p)|^2+r_{g}^2(1-R_{g}^2)|\bar{A}_f(p)|^2 \; , \end{array} $

      (5)

      where $ r_{g}^2 = \dfrac{\int |\bar{A}_{g}|^2{\mathrm{d}\Phi_{g}}}{\int |{A}_{g}|^2{\mathrm{d}\Phi_{g}}} $, $R_g {\rm e}^{-{\rm i}\delta_g} = \dfrac{ \int A_{g}^*\bar{A}_{g} {\mathrm{d}\Phi_{g}}}{\sqrt {\int |A_{g}|^2{\mathrm{d}\Phi_{g}}\int |\bar{A}_{g}|^2{\mathrm{d}\Phi_{g}}}}$. The values of r, R, and δ for the three tag modes are summarized in Table 4. The second term in Eq. (5) is ignored in practice owing to the relatively small value of $ r_{g}^2(1-R_{g}^2) $.

      Mode$ K^-\pi^+ $$ K^-\pi^+\pi^0 $$ K^-\pi^+\pi^-\pi^+ $
      $ r(\%) $$ 5.86 \pm 0.02 $ [31]$ 4.41 \pm 0.11 $ [32]$ 5.50 \pm 0.07 $ [32]
      R1$ 0.79 \pm 0.04 $ [32]$ 0.44_{-0.10}^{+0.09} $ [32]
      $ \delta(^{\circ}) $$ 192.1^{+8.6}_{-10.2} $ [31]$ 196 \pm 11 $ [32]$ 161^{+28}_{-18} $ [32]

      Table 4.  The input values of r, R, and δ for the three tag modes.

      The amplitude $ U_i $ is constructed with the spin factor, Blatt-Weisskopf barrier factors [34], and resonance propagators. To construct $ U_i $ of the four-body $ D^0 $ decay, the isobar model is applied, factorizing the decay into subsequent two-body decay amplitudes [3537]. The general amplitude of the i-th intermediate process in a four-body decay is given by

      $ \begin{array}{*{20}{l}} U_i(p)=S_i(p)B_{L_D}(p)P_{R_1}(p)B_{L_{R_1}}(p)P_{R_2}(p)B_{L_{R_2}}(p)\; , \end{array} $

      (6)

      where $ S_i $ is the spin factor of the i-th decay amplitude, $ B_{L_X} $ ($ X=D $, $ R_1 $, $ R_2 $) are the Blatt-Weisskopf barrier factors for the D meson and the resonances $ R_1 $ and $ R_2 $, and $ P_{R_1} $ and $ P_{R_2} $ are the propagators of $ R_{1} $ and $ R_{2} $, respectively. The amplitude is constructed with the exchange symmetry for indistinguishable pions.

      The spin factor is constructed with the covariant Zemach (Rarita-Schwinger) tensor formalism [3841] by combining pure-orbital-angular-momentum covariant tensors $ \tilde{t}^{(L)}_{\mu_1\cdots \mu_l} $ and the momenta of parent particles together with Minkowski metric $ g_{\mu\nu} $ and Levi-Civita symbol $ \epsilon_{\mu\nu\lambda\sigma} $. For a process $ a\to bc $, the covariant tensors $ \tilde{t}^{(L)}_{\mu_1\cdots \mu_l} $ for the final states of pure orbital angular momentum L are

      $ \begin{array}{*{20}{l}} \tilde{t}^{(L)}_{\mu_1\cdots \mu_l} = (-1)^{L}P^{(L)}_{\mu_1\cdots \mu_L\mu'_1\cdots \mu'_L}(p_a)r^{\mu'_1}\cdots r^{\mu'_L}\; , \end{array} $

      (7)

      where $ r=p_b-p_c $ and $ P^{(L)}_{\mu_1\cdots \mu_L\mu'_1\cdots \mu'_L}(p_a) $ is the spin projection operator of the particle a,

      $ P^{(0)}(p_a) =1\; , $

      (8)

      $ P^{(1)}_{\mu\mu'}(p_a) =-g_{\mu\mu'}+\frac{p_{a\mu}p_{a\mu'}}{p_a^2}\; , $

      (9)

      $ \begin{aligned}[b] P^{(2)}_{\mu\nu\mu'\nu'}(p_a)=\;&\frac{1}{2}[P^{(1)}_{\mu\mu'}(p_a)P^{(1)}_{\nu\nu'}(p_a)+P^{(1)}_{\mu\nu'}(p_a)P^{(1)}_{\nu\mu'}(p_a)] \\ &-\frac{1}{3}P^{(1)}_{\mu\nu}(p_a)P^{(1)}_{\mu'\nu'}(p_a)\; . \end{aligned} $

      (10)

      Following the isobar model, the spin factors of the four-body decay $ D^{0}\to P_1P_2P_3P_4 $ are summarized in Table 5.

      Decay chainSpin factor
      $ D[S]\to P P_1, P[S]\to S P_2, S[S]\to P_3 P_4 $1
      $ D[S]\to P P_1, P[P]\to V P_2, V[P]\to P_3 P_4 $$ \tilde{t}_{\mu}(P)\tilde{t}^{\mu}(V) $
      $ D[S]\to P P_1, P[D]\to T P_2, T[D]\to P_3 P_4 $$ \tilde{t}_{\mu\nu}(P)\tilde{t}^{\mu\nu}(T) $
      $ D[P]\to A P_1, A[S]\to V P_2, V[P]\to P_3 P_4 $$ \tilde{t}_{\mu}(D)P^{\mu\nu}(A)\tilde{t}_{\nu}(V) $
      $ D[P]\to A P_1, A[D]\to V P_2, V[P]\to P_3 P_4 $$ \tilde{t}_{\mu}(D)\tilde{t}^{\mu\nu}(A)\tilde{t}_{\nu}(V) $
      $ D[P]\to A P_1, A[P]\to S P_2, S[S]\to P_3 P_4 $$ \tilde{t}_{\mu}(D)\tilde{t}^{\mu}(A) $
      $ D[P]\to A P_1, A[P]\to T P_2, T[D]\to P_3 P_4 $$ \tilde{t}_{\mu}(D)\tilde{t}_{\nu}(A)P^{\mu\nu\rho\sigma}(A)\tilde{t}_{\rho\sigma}(T) $
      $ D[P]\to V_1 P_1, V_1[P]\to V_2 P_2, V_2[P]\to P_3 P_4 $$ \tilde{t}_{\mu}(D)P^{\mu\nu}(V_1)\epsilon_{\nu\lambda\sigma\alpha}p_{V_1}^{\alpha}\tilde{t}^{\lambda}(V_1)\tilde{t}^{\sigma}(V_2) $
      $ D[D]\to PT P_1, PT[S]\to T P_2, T[D]\to P_3 P_4 $$ \tilde{t}_{\mu\nu}(D)P^{\mu\nu\alpha\beta}(PT)\tilde{t}_{\alpha\beta}(T) $
      $ D[D]\to PT P_1, PT[P]\to V P_2, V[P]\to P_3 P_4 $$ \tilde{t}_{\mu\nu}(D)P^{\mu\nu\alpha\beta}(PT)\tilde{t}_{\alpha}(PT)\tilde{t}_{\beta}(V) $
      $ D[D]\to PT P_1, PT[D]\to S P_2, S[S]\to P_3 P_4 $$ \tilde{t}_{\mu\nu}(D)\tilde{t}^{\mu\nu}(PT) $
      $ D[D]\to T P_1, T[D]\to V P_2, V[P]\to P_3 P_4 $$ \tilde{t}_{\mu\nu}(D)P^{\mu\nu\alpha\beta}(T)\epsilon_{\alpha\lambda\sigma\rho}p_{T}^{\rho}\tilde{t}_{\beta}^{\lambda}(T)P^{\sigma\gamma}(T)\tilde{t}_{\gamma}(V) $
      $ D[D]\to T_1 P_1, T_1[P]\to T_2 P_2, T_2[D]\to P_3 P_4 $$ \tilde{t}_{\mu\nu}(D)P^{\mu\nu\alpha\beta}(T_1)\epsilon_{\alpha\lambda\sigma\rho}p_{T_1}^{\rho}\tilde{t}^{\lambda}(T_1)\tilde{t}^{\sigma}_{\beta}(T_2) $
      $ D[S]\to S_1 S_2, S_1[S]\to P_1 P_2, S_2[S]\to P_3 P_4 $1
      $ D[P]\to V S, V[P]\to P_1 P_2, S[S]\to P_3 P_4 $$ \tilde{t}_{\mu}(D)\tilde{t}^{\mu}(V) $
      $ D[S]\to V_1 V_2, V_1[P]\to P_1 P_2, V_2[P]\to P_3 P_4 $$ \tilde{t}_{\mu}(V_1)\tilde{t}^{\mu}(V_2) $
      $ D[P]\to V_1 V_2, V_1[P]\to P_1 P_2, V_2[P]\to P_3 P_4 $$ \epsilon_{\mu\nu\alpha\beta}p_{D}^{\beta}\tilde{t}^{\mu}(D)\tilde{t}^{\nu}(V_1)\tilde{t}^{\alpha}(V_2) $
      $ D[D]\to V_1 V_2, V_1[P]\to P_1 P_2, V_2[P]\to P_3 P_4 $$ \tilde{t}_{\mu\nu}(D)\tilde{t}^{\mu}(V_1)\tilde{t}^{\nu}(V_2) $
      $ D[D]\to T S, T[D]\to P_1 P_2, S[S]\to P_3 P_4 $$ \tilde{t}_{\mu\nu}(D)\tilde{t}^{\mu\nu}(T) $
      $ D[P]\to T V, T[D]\to P_1 P_2, V[P]\to P_3 P_4 $$ \tilde{t}_{\mu}(D)\tilde{t}^{\mu\nu}(T)\tilde{t}_{\nu}(V) $
      $ D[D]\to T V, T[D]\to P_1 P_2, V[P]\to P_3 P_4 $$ \epsilon_{\mu\nu\alpha\beta}p_{D}^{\beta}\tilde{t}^{\mu\rho}(D)\tilde{t}^{\nu}_{\rho}(T)\tilde{t}^{\alpha}(V) $
      $ D[S]\to T_1 T_2, T_1[D]\to P_1 P_2, T_2[D]\to P_3 P_4 $$ \tilde{t}_{\mu\nu}(T_1)\tilde{t}^{\mu\nu}(T_2) $
      $ D[P]\to T_1 T_2, T_1[D]\to P_1 P_2, T_2[D]\to P_3 P_4 $$ \epsilon_{\mu\nu\alpha\beta}p_{D}^{\beta}\tilde{t}^{\mu}(D)\tilde{t}^{\nu\rho}(T_1)\tilde{t}_{\rho}^{\alpha}(T_2) $
      $ D[D]\to T_1 T_2, T_1[D]\to P_1 P_2, T_2[D]\to P_3 P_4 $$ \tilde{t}_{\mu\nu}(D)\tilde{t}^{\mu\rho}(T_1)\tilde{t}^{\nu}_{\rho}(T_2) $

      Table 5.  Summary of the spin factors in this analysis, where S, P, V, A, T, and $ PT $ denote scalar, pseudo-scalar, vector, axial-vector, tensor, and pseudo-tensor mesons, respectively. $ [S] $, $ [P] $, and $ [D] $ represent orbital angular momenta L = 0, 1, and 2 in the decays, respectively.

      The Blatt-Weisskopf barrier factors $ B_L(q) $ are derived by assuming a square well interaction potential as

      $ B_{L=0}(q)=1\; , $

      (11)

      $ B_{L=1}(q)=\sqrt{\frac{2}{q^2+q_R^2}}\; , $

      (12)

      $ B_{L=2}(q)=\sqrt{\frac{13}{q^4+3q^2q_R^2+9q_R^4}}\; , $

      (13)

      where q is the momentum of the daughter particle in the rest frame of the mother particle, L is the orbital angular momentum, and $ q_R=1/R $ is a hadron "scale" parameter (R denotes the radius of the centrifugal barrier). In this analysis, the radius R is taken to be $ 5.0\; {\rm{GeV}^{-1}}c $ for $ D^{0} $ mesons, and $ 3.0\; {\rm{GeV}^{-1}}c $ for other intermediate resonances.

      The resonance propagators are generally described by a relativistic Breit-Wigner function

      $ P(s)=\frac{1}{m_0^2-s- {\rm i}\sqrt{s}\Gamma(s)}\; , $

      (14)

      where $ m_0 $ is the nominal mass of the resonance. For a resonance decaying into two scalar particles $ a\to bc $, $ \Gamma_{a\to bc}(s) $ is given by

      $ \Gamma_{a\to bc}(s)=\Gamma_0^{a\to bc}\left(\frac{q}{q_0}\right)^{2L+1}\left(\frac{m_0}{\sqrt{s}}\right)^{2}\left(\frac{B_L(q)}{B_L(q_0)}\right)^2\; , $

      (15)

      where $ \Gamma_0^{a\to bc} $ is the nominal width when $ s=m_0^2 $ and $ q_0 $ is the corresponding momentum of the daughter particle in the rest frame of the mother particle.

      In this analysis, the propagators of $ \rho(770) $ and $ \rho(1450) $ are described by the Gounaris-Sakurai parametrization [42]. The propagator of $ f_0(980) $ in the decay $a_1(1420)\to f_0(980)\pi$ is described by a Flatt$ \rm{\acute{e}} $ parametrization of $ \pi\pi $ and $ KK $ coupled channels with parameters from Ref. [43].

      The K-matrix parametrization [44, 45] instead of the Breit-Wigner formula is adopted for the $ \pi\pi $ S-wave. The P-vector parametrization of the K-matrix for the $ \pi\pi $ S-wave with five coupled channels $ \pi\pi $, $ KK $, $ \pi\pi\pi\pi $, $ \eta\eta $, $ \eta\eta' $ and five poles is written as

      $ \begin{array}{*{20}{l}} F_{\mu}(s) = [I-{\rm i}K(s)\rho(s)]^{-1}_{\mu\nu}P_{\nu}(s)\; , \end{array} $

      (16)

      where I is the identity matrix, K is the K-matrix describing the scattering process, ρ is the PHSP matrix, and P is the initial production vector (P-vector). The indices μ and ν denote the coupled channels ($ \pi\pi $, $ KK $, $ \pi\pi\pi\pi $, $ \eta\eta $, and $ \eta\eta' $), and only $ F_{\mu}(s) $ of the $ \pi\pi $ component is used as the propagator of the $ \pi\pi $ S-wave. The K-matrix is given by

      $ K_{\mu\nu}(s) = \left(\sum\limits_{\alpha}\frac{g^{\alpha}_{\mu}g^{\alpha}_{\nu}}{m_{\alpha}^2-s}+f_{\mu\nu}^{\rm scatt}\frac{1-s_{0}^{\rm scatt}}{s-s_{0}^{\rm scatt}}\right)f_{A0}(s)\; , $

      (17)

      containing five poles ($ \alpha = 1 $$ 5 $). The K-matrix parmeters are fixed to those in Ref. [45]. The P-vector is written as

      $ P_{\nu}(s) = \left(\sum\limits_{\alpha}\frac{\beta_{\alpha}g^{\alpha}_{\nu}}{m_{\alpha}^2-s}+f_{\nu}^{\rm prod}\frac{1-s_{0}^{\rm prod}}{s-s_{0}^{\rm prod}}\right)\; , $

      (18)

      where the poles are the same as those for the K-matrix. The parameters $ \beta_{\alpha} $ and $ f_{1\nu}^{\rm prod} $ are free in the fit. For the parameter $ s_0^{\rm prod} $ in the P-vector, we assume that they take the same value in all $ \pi\pi $ S-waves and fix them to $ -5\; {\rm GeV}^{2}/c^{4} $ since they are insensitive to any choice if $ s_0^{\rm prod}\leq-5\; {\rm GeV}^{2}/c^{4} $ in our fit.

      In practice, two $ \pi\pi $ S-waves exist in the decay $ D^0\to S_1S_2 $, and the corresponding propagator is

      $ \begin{aligned}[b] F'_{\mu\nu}(s_1,s_2)=\;& [I-{\rm i}K(s_1)\rho(s_1)]^{-1}_{\mu\rho}[I-{\rm i}K(s_2)\rho(s_2)]^{-1}_{\nu\sigma} \\ & \times P_{\rho\sigma}(s_1,s_2)\; , \end{aligned} $

      (19)

      where $ s_1 $ and $ s_2 $ are the invariant masses squared of $ \pi\pi $, $ P_{\rho\sigma}(s_1,s_2) $ is the expansion of the product of two P-vectors $ P_{\rho}(s_1) $ and $ P_{\sigma}(s_2) $, and the corresponding coefficients of each term are taken to be independent parameters in the fit. Benefiting from the exchange symmetry of two $ \pi\pi $ S-waves in the amplitude, $ P_{\rho\sigma}(s_1,s_2) $ can be written as

      $ \begin{aligned}[b] P_{\rho\sigma}(s_1,s_2)=&\sum\limits_{\alpha,\beta}^{\alpha\le\beta}a_{\alpha,\beta}\bigg[\frac{g^{\alpha}_{\rho}g^{\beta}_{\sigma}}{(m_{\alpha}^2-s_1)(m_{\beta}^2-s_2)} + \frac{g^{\beta}_{\rho}g^{\alpha}_{\sigma}}{(m_{\beta}^2-s_1)(m_{\alpha}^2-s_2)}\bigg] \\ &+\sum\limits_{\alpha}b_{\alpha,\rho}\frac{g^{\alpha}_{\sigma}(1-s_{0}^{\rm prod})}{(s_1-s_{0}^{\rm prod})(m_{\alpha}^2-s_2)} \\& + \sum\limits_{\alpha}b_{\alpha,\sigma}\frac{g^{\alpha}_{\rho}(1-s_{0}^{\rm prod})}{(s_2-s_{0}^{\rm prod})(m_{\alpha}^2-s_1)} \\ &+c_{[\rho,\sigma]}\frac{(1-s_{0}^{\rm prod})^2}{(s_1-s_{0}^{\rm prod})(s_2-s_{0}^{\rm prod})}, \end{aligned} $

      (20)

      where $ c_{[\rho,\sigma]}=c_{\rho,\sigma} + c_{\sigma,\rho} $, and $ c_{[\rho,\sigma]} $ with $ \rho>\sigma $ is fixed to zero in the fit. Similar to one $ \pi\pi $ S-wave, only the ($ \pi\pi,\pi\pi $) component of $ F'_{\mu\nu} $ is used.

      The propagators of the resonances $ a_1(1260) $, $ a_{1}(1640) $, $ h_1(1170) $ and $ \pi(1300) $ decaying into $ 3\pi $ are described by relativistic Breit-Wigner functions with the coupled channel $ 3\pi $, where $ \Gamma(s) $ are obtained by integrating the amplitude squared over PHSP

      $ \Gamma_{a\to bcd}(s) \propto \frac{m_0}{\sqrt{s}}\int \sum\limits_{\rm spin}|A_{a \to bcd}|^2 \mathrm{d}\Phi_3\; \; , $

      (21)

      based on the amplitudes $ A_{R\to3\pi} $ obtained in this analysis. Figure 4 shows the $ \Gamma(s) $ obtained in the amplitude analysis. For other resonances decaying into $ 3\pi $, a relativistic Breit-Wigner function with a constant $ \Gamma(s) $ is used. The parameters of $ a_1(1260) $ and $ \pi(1300) $ are determined in the fit and others are fixed to their respective PDG values [25]. Owing to its small contribution, the $ a_1(1420) $ is described as a relativistic Breit-Wigner function with a constant $ \Gamma(s) $, even though it has been regarded as the effect of the Triangle Singularity [46]. Only the decay of $ a_1(1420)\to f_0(980)\pi $ is considered, and the corresponding resonant parameters are taken from Ref. [47].

      Figure 4.  (color online) The $ \Gamma(s) $ of $ a_1(1260) $, $ a_{1}(1640) $, $ h_1(1170) $, and $ \pi(1300) $ resonances obtained in the amplitude analysis.

      To account for the resolution effect for the narrow resonances $ \phi(1020) $ and $ \omega(782) $, the shapes of the corresponding relativistic Breit-Wigner functions in the $ M(\pi^+\pi^-\pi^0) $ distribution are convolved with a Gaussian function. The parameters of this Gaussian function are obtained by the fit to the η peak in the $ M(\pi^+\pi^-\pi^0) $ distribution with the Breit-Wigner function of η convolved with the same Gaussian function.

      In this analysis, two decay topologies, $ D\to R_1 R_2 $ and $ D\to R_1\pi \to R_2\pi\pi $, are considered. Meanwhile, the isospin symmetry is considered by applying the Clebsch-Gordan (CG) coefficients between the isospin processes. The isospin states for the u and d quarks are

      $ |u\rangle:|\frac{1}{2}, \frac{1}{2}\rangle,\; |\bar{u}\rangle:-|\frac{1}{2},-\frac{1}{2}\rangle,\; |d\rangle:|\frac{1}{2}, -\frac{1}{2}\rangle,\; |\bar{d}\rangle:|\frac{1}{2},\frac{1}{2}\rangle. $

      (22)

      Based on the quark components of hadrons, the isospin states of the hadrons are shown in Table 6. The isospin amplitudes for typical processes are written as

      $ \pi^+,b^+,\rho^+,a^+ $$ \pi^0,b^0,\rho^0,a^0 $$ \pi^-,b^-,\rho^-,a^- $
      $ I=1 $$ |u\bar{d}\rangle:|1,1\rangle $$ |\dfrac{u\bar{u}-d\bar{d}}{\sqrt{2}}\rangle:-|1,0\rangle $$ |d\bar{u}\rangle:-|1,-1\rangle $
      $ h,\omega,\phi,f $
      $ I=0 $$ |(u\bar{u}+d\bar{d})+g(s\bar{s})\rangle:-|0,0\rangle $

      Table 6.  The quark components and isospin states for different hadrons used in the analysis.

      $ \begin{aligned}[b] A(D^{0}\to a_1^+\pi^-)=\;& A(\rho^0_{\pi^+\pi^-}\pi^+,\pi^-) + A(\rho^+_{\pi^+\pi^{0}}\pi^{0},\pi^-) \\&- c[2A(f_{\pi^+\pi^-}\pi^+,\pi^-) + A(f_{\pi^{0}\pi^{0}}\pi^+,\pi^-)], \end{aligned} $

      (23)

      $ \begin{aligned}[b] A(D^{0}\to a_1^0\pi^{0})\; =\;& -A(\rho^+_{\pi^+\pi^{0}}\pi^-,\pi^{0}) + A(\rho^-_{\pi^{0}\pi^-}\pi^+,\pi^{0}) \\&- c[2A(f_{\pi^+\pi^-}\pi^{0},\pi^{0}) + A(f_{\pi^{0}\pi^{0}}\pi^{0},\pi^{0})], \end{aligned} $

      (24)

      $ \begin{aligned}[b] A(D^{0}\to a_1^-\pi^+)=\;& -A(\rho^0_{\pi^+\pi^-}\pi^-,\pi^+) - A(\rho^-_{\pi^{0}\pi^-}\pi^{0},\pi^+) \\&- c[2A(f_{\pi^+\pi^-}\pi^-,\pi^+) + A(f_{\pi^{0}\pi^{0}}\pi^-,\pi^+)], \end{aligned} $

      (25)

      $ \begin{aligned}[b] A(D^{0}\to h_1^0\pi^{0})\; =\;& -A(\rho^+_{\pi^+\pi^{0}}\pi^-,\pi^{0}) - A(\rho^-_{\pi^{0}\pi^-}\pi^+,\pi^{0}) \\&+ A(\rho^0_{\pi^+\pi^-}\pi^{0},\pi^{0}), \end{aligned} $

      (26)

      $ \begin{aligned}[b] A(D^{0}\to f f')=\;& 4A(f_{\pi^+\pi^-},f'_{\pi^+\pi^-}) + 2[A(f_{\pi^+\pi^-},f'_{\pi^{0}\pi^{0}})\\&+ A(f_{\pi^{0}\pi^{0}},f'_{\pi^+\pi^-})] + A(f_{\pi^{0}\pi^{0}},f'_{\pi^{0}\pi^{0}}), \end{aligned} $

      (27)

      $ A(D^{0}\to \rho^0 f)= 2A(\rho^0_{\pi^+\pi^-},f_{\pi^+\pi^-}) + A(\rho^0_{\pi^+\pi^-},f_{\pi^{0}\pi^{0}}). $

      (28)

      Here, $ a_1 $ and ρ are isospin vectors, whereas $ h_1 $ and $ f^{(')} $ are isospin scalars. The normalization factors are dropped here and c is the relative difference (including magnitude and phase) between two isospin processes determined in the amplitude analysis.

    VI.   AMPLITUDE FIT PROCESS AND RESULT
    • In the amplitude fit, the sPlot technique [48] is applied to deal with the background effect. In the sPlot method, the four-momenta of the final states are chosen as control variables x, $ M_{\rm bc}^{\rm sig} $ and $ M_{\rm bc}^{\rm tag} $ are chosen as discriminating variables y, and the weight for each candidate event is obtained by an unbinned maximum likelihood fit on the 2D distribution of $ M_{\rm bc}^{\rm sig} $ versus $ M_{\rm bc}^{\rm tag} $, as described in Sec. IV.

      According to the sPlot technique and based on the above DT fit, the weight $ W_s(y) $ of each candidate event is calculated as

      $W_s(y) = \frac{\sum _{i}V_{si}F_{i}(y)}{\sum _{i}N_jF_j(y)+\sum _{k}N'_kF'_k(y)}, $

      (29)

      where the index s denotes the signal, $ F_{i(j)}(y) $ are the PDFs of the signal ($ i(j)=s $) and backgrounds ($ i(j)\neq s $) from BKGI to BGKV with floating yields, $ F'_{k}(y) $ are the PDFs of BKGVI with fixed yields, and $ N_{i(j)} $ and $ N_k^{'} $ are the corresponding yields for each component. The inverse of covariance matrix $ V_{ij} $ is calculated with

      $ V_{ij}^{-1}=\sum\limits_{n\in {\rm Data}}\frac{F_i(y_n)F_j(y_n)}{(\sum _{k}N_kF_k(y_n)+\sum\limits_{l}N'_lF'_l(y_n))^{2}}, $

      (30)

      where the summations for indices k and l run over all the components of $ F_k(y) $ and $ F'_l(y) $, as described above, respectively. The summation for index n runs over all the events in the DT fit.

      Based on the sPlot technique, $ W_s(y) $, defined in Eq. (29), includes the effects of the backgrounds with a fixed number of events. To account for these effects, the coefficients $ c_{i} $ are calculated as

      $ c_{i} = \sum\limits_{j}V_{sj}\nu_{ij}, $

      (31)

      where the summation runs over the signal and backgrounds from BKGI to BKGV with floating yields in the fit, i represents the different components in the BKGVI, and $ \nu_{ij} $ is

      $ \nu_{ij}=\sum\limits_{n\in {\rm Data}}\frac{F'_i(y_n)F_j(y_n)}{(\sum\limits_{k}N_kF_k(y_n)+\sum\limits_{l}N'_lF'_l(y_n))^{2}}. $

      (32)

      With the obtained $ W_s(y) $ and $ c_i $, the x distribution of all candidate events with the weight $ W_s(y) $ is the real signal together with the contributions of the peaking background in BKGVI,

      $ N_{s}P_s(x) + \sum\limits_{j} c_{j}N'_{j}P'_{j}(x), $

      (33)

      where $ P_s(x) $ and $ P'_{j}(x) $ are the x distributions of the signal and the peaking background j in the BKGVI and $ N_{s} $ and $ N'_{j} $ are the corresponding yields obtained from the above 2D fit and MC simulation, respectively.

      In practice, the PDF of observing a signal event with the given final kinematic p is written as

      $ P_{s}(p) = \frac{\epsilon_{s}(p)|M_f(p)|^{2}\phi_4(p)}{\int\epsilon_{s}(p)|M_f(p)|^{2}\mathrm{d}\Phi_{4}}\; , $

      (34)

      where $ \epsilon_s(p) $ is the signal efficiency, $ |M_f(p)|^2 $ is the differential cross section, as discussed in Sec. V, and $ \phi_4(p) $ is the PHSP density. The normalization factor is calculated by MC integration with the PHSP signal MC sample after event selection,

      $ \int \epsilon_{s}(p)|M_f(p)|^{2}\mathrm{d}\Phi_{4} \propto \frac{1}{N_{\rm MC}}\sum\limits_{i=1}^{N_{\rm MC}} |M_f(p_i)|^{2}. $

      (35)

      where the $ N_{\rm MC} $ is the number of events of the PHSP signal MC. According to Eq. (33), the weighted likelihood $ \mathrm{ln}L $ is given by

      $ f\left[\sum\limits_{i \in {\rm Data}} W_s(y_i)\mathrm{ln}P_{s}(p_i) - \sum\limits_{j}\sum\limits_{i\in {\rm BMC}_j} \omega_{j} \mathrm{ln}P_{s}(p_i)\right], $

      (36)

      where the second term in the bracket is the contribution from the peaking background, estimated using the corresponding simulated background MC (BMC). The normalization factor $ \omega_j $ is given by

      $ \omega_j = c_{j}N_{{\rm BKG}_j}^{\rm Data}/N_{{\rm BKG}_j}^{\rm MC}. $

      (37)

      Here, $ c_j $ is obtained from Eq. (31), and $ N_{{\rm BKG}_j}^{\rm Data} $ and $ N_{{\rm BKG}_j}^{\rm MC} $ are the background yields in the data and simulated background events, respectively. The factor f in Eq. (36),

      $ f = \frac{\sum _{i\in {\rm Data}}W_s(y_i) + \sum _{j}\omega_jN^{\rm MC}_{{\rm BKG}_j}}{\sum _{i\in {\rm Data}}W^2_s(y_i) + \sum _{j}\omega_j^2N^{\rm MC}_{{\rm BKG}_j}}\; , $

      (38)

      is the global factor to correct the statistical bias in the weighted maximum likelihood fit.

      The total likelihood function in this analysis is summed over the two signal channels and three tag modes,

      $ \mathrm{ln}L_{\rm total} = \sum\limits_{i\in {\rm tag}} \left(\mathrm{ln}L_{i}^{\pi^+\pi^-\pi^+\pi^-} + \mathrm{ln}L_{i}^{\pi^+\pi^-\pi^0\pi^0}\right), $

      (39)

      and the free parameters are optimized via a maximum likelihood fit using the MINUIT [49] package.

      Generally, all the possible intermediate processes including the resonances listed in Table 7 and based on $ J^{PC} $ conservation are considered in the fit. The only exception are $ D^0 $ decays, where parity conservation is not required. Only processes with a significance greater than $ 5\sigma $ are kept during the fit unless otherwise noted. Here, the significance of a specific amplitude is calculated according to the Wilks's Theorem by comparing the change of log-likelihood (2ΔlnL) to the expected values from the chi-square distribution ($ \chi^{2}_{\Delta N_{\rm para}} $) with the number of degrees of freedom (NDF) equal to the change of the numbers of fit parameters ($ \Delta N_{\rm para} $). If the significance of the amplitudes containing the isospin vector in the $ 3\pi $ invariant mass spectrum is greater than 5σ, the corresponding isospin partners with significance greater than 3σ are also kept. For the P-vector of the $ \pi\pi $ S-wave, only the parameters before the 1st and 2nd poles and $ \pi\pi $, $ KK $ non-resonant terms are considered, and the others are fixed to zero in the nominal fit. Meanwhile, only those with significance greater than $ 3\sigma $ are kept.

      $ J^{P}=0^+ $$ J^{P}=0^- $$ J^{P}=1^+ $$ J^{P}=1^- $$ J^{P}=2^+ $$ J^{P}=2^{-} $
      $ \pi\pi $$ (\pi\pi)_{S} $$ \rho(770) $$ f_{2}(1270) $
      $ \rho(1450) $
      $ \pi\pi\pi $$ \pi(1300) $$ a_{1}(1260) $$ \omega(782) $$ a_{2}(1320) $$ \pi_{2}(1670) $
      $ a_{1}(1420) $$ \phi(1020) $
      $ a_{1}(1640) $$ \pi_{1}(1400) $
      $ h_{1}(1170) $$ \pi_{1}(1600) $

      Table 7.  Resonances considered in this analysis.

      To find the optimal solution, the baseline model containing the processes from Refs. [12, 13] is built up first. Next, starting from the baseline model, the significance of each possible process is tested, and the most significant among those satisfying the significance requirement is added to the current model. This step is repeated until no additional processes can be added. After this step, the significances of individual processes in the existing model are tested again, and those that do not satisfy the significance requirement are removed. The above steps are repeated until all the processes in the model satisfy the significance requirement, and no further ones can be added. The nominal amplitude model is obtained with this strategy and no multiple solutions are found.

      Using the nominal amplitude model, the fit fraction (FF) of a specific amplitude i is calculated as

      $ {\rm FF}_i = \frac{\int |\Lambda_iU_i(p)|^2\mathrm{d}\Phi}{\int |\sum _j \Lambda_j U_j(p)|^2\mathrm{d}\Phi}. $

      (40)

      The FF of the interference between two amplitudes i and j ($ j\neq i $) is calculated as

      $ {\rm FF}_{ij} = \frac{\int 2{\rm Re}[\Lambda_iU_i(p)\Lambda^{*}_jU^{*}_j(p)]\mathrm{d}\Phi}{\int |\sum _k \Lambda_k U_k(p)|^2\mathrm{d}\Phi}. $

      (41)

      The $ CP $-even fraction $ F_+ $ is calculated as

      $ F_+^{f} = \frac{\int |A^+_{f}(p)|^2\mathrm{d}\Phi}{\int |A^+_{f}(p)|^2+|A^-_{f}(p)|^2 \mathrm{d}\Phi}, $

      (42)

      where $ A^{\pm}_{f}(p)=\dfrac{1}{\sqrt{2}}[A_{f}(p) \; \pm\; \bar{A}_{f}(p)] $ is the amplitude of the $ D^0\to f $ decay in a $ CP $-even or $ CP $-odd state. The FF results of different amplitudes are summarized in Tables 8, 9 and 10. The resonant parameters, masses, and widths of $ a_1(1260) $ and $ \pi(1300) $ are determined by the parameter scans, as shown in Fig. 5, which are compatible with the values from PDG. The fit fractions of the interference terms are summarized in Table 11. The fit results show large interferences among the dominant intermediate processes $ D^{0}\to a_{1}(1260)\pi $, $ D^{0}\to\pi(1300)\pi $, $ D^{0}\to\rho(770)\rho(770) $, and $ D^{0}\to2(\pi\pi)_{S} $. The results for $ F_+ $ obtained in this analysis are presented in Table 12, which show good agreement with other measurements. In all the above tables, the mean values are obtained based on the output from the MINUIT fit, while the corresponding statistical uncertainties are estimated by the bootstrap method [50] based on the data since the weighted maximum likelihood fit is used, and the statistical uncertainties given by the inverse second derivative of the negative logarithmic likelihood are no longer asymptotically correct [51]. In the bootstrap method, the bootstrap samples are generated by repeatedly resampling the data set with replacement a thousand times, and both sPlot and amplitude analysis are performed for these samples. The width of the distribution of the estimated parameter values is used to estimate the parameter uncertainty. The comparisons between data and MC projection based on the nominal model for various invariant mass and angle distributions are shown in Figs. 6 and 7.

      AmplitudeMagnitudePhase (rad)FF (%)Significance (σ)
      $ \pi^{+}\pi^{-}\pi^{+}\pi^{-} $$ \pi^{+}\pi^{-}\pi^{0}\pi^{0} $
      $ a_{1}(1260)^{+}\pi^{-} $100(fixed)0(fixed)82.2 $ \pm $ 3.3 $ \pm $ 2.3 $ \pm $ 15.857.4 $ \pm $ 2.7 $ \pm $ 3.0 $ \pm $ 6.5$> $10
      $ a_{1}(1260)^{-}\pi^{+} $$ 35.3 \pm 2.7 \pm 0.8 \pm 4.6 $$ \phantom{+}0.23 \pm 0.07 \pm 0.02 \pm 0.16 $10.3 $ \pm $ 1.5 $ \pm $ 0.3 $ \pm $ 2.57.2 $ \pm $ 1.1 $ \pm $ 0.2 $ \pm $ 2.2$> $10
      $ a_{1}(1260)^{0}\pi^{0} $$ 50.9 \pm 3.1 \pm 0.4 \pm 4.7 $$ -2.99 \pm 0.06 \pm 0.08 \pm 0.14 $32.9 $ \pm $ 3.2 $ \pm $ 1.6 $ \pm $ 8.3$> $10
      $ a_{1}(1420)^{+}\pi^{-} $$ 19.0 \pm 3.6 \pm 1.3 \pm 3.9 $$ \phantom{+}2.70 \pm 0.18 \pm 0.05 \pm 1.09 $0.6 $ \pm $ 0.2 $ \pm $ 0.0 $ \pm $ 0.20.3 $ \pm $ 0.1 $ \pm $ 0.0 $ \pm $ 0.16.0
      $ a_{1}(1640)^{+}\pi^{-} $$ 20.1 \pm 3.0 \pm 2.6 \pm 5.8 $$ -2.07 \pm 0.16 \pm 0.02 \pm 0.28 $1.7 $ \pm $ 0.5 $ \pm $ 0.4 $ \pm $ 0.81.1 $ \pm $ 0.3 $ \pm $ 0.2 $ \pm $ 0.67.3
      $ a_{1}(1640)^{-}\pi^{+} $$ 10.5 \pm 2.8 \pm 0.6 \pm 3.8 $$ -1.26 \pm 0.29 \pm 0.23 \pm 0.50 $0.5 $ \pm $ 0.3 $ \pm $ 0.0 $ \pm $ 0.40.3 $ \pm $ 0.2 $ \pm $ 0.0 $ \pm $ 0.25.2
      $ a_{2}(1320)^{+}\pi^{-} $$ 0.23 \pm 0.07 \pm 0.03 \pm 0.05 $$ -2.92 \pm 0.30 \pm 0.14 \pm 0.23 $0.2 $ \pm $ 0.1 $ \pm $ 0.0 $ \pm $ 0.10.2 $ \pm $ 0.1 $ \pm $ 0.0 $ \pm $ 0.14.6
      $ a_{2}(1320)^{-}\pi^{+} $$ 0.30 \pm 0.05 \pm 0.01 \pm 0.04 $$ -0.47 \pm 0.21 \pm 0.06 \pm 0.15 $0.4 $ \pm $ 0.1 $ \pm $ 0.0 $ \pm $ 0.10.3 $ \pm $ 0.1 $ \pm $ 0.0 $ \pm $ 0.16.4
      $ h_{1}(1170)^{0}\pi^{0} $$ 9.7 \pm 2.2 \pm 1.5 \pm 3.6 $$ -0.59 \pm 0.27 \pm 0.09 \pm 0.35 $1.3 $ \pm $ 0.6 $ \pm $ 0.4 $ \pm $ 1.06.5
      $ \pi(1300)^{+}\pi^{-} $$ 76.3 \pm 3.6 \pm 4.6 \pm 4.3 $$ -2.325 \pm 0.044 \pm 0.038 \pm 0.297 $32.3 $ \pm $ 2.6 $ \pm $ 1.6 $ \pm $ 4.215.6 $ \pm $ 1.4 $ \pm $ 1.3 $ \pm $ 2.2$> $10
      $ \pi(1300)^{-}\pi^{+} $$ 65.1 \pm 3.4 \pm 3.3 \pm 3.9 $$ -2.631 \pm 0.045 \pm 0.083 \pm 0.208 $23.5 $ \pm $ 2.3 $ \pm $ 0.5 $ \pm $ 3.911.4 $ \pm $ 1.1 $ \pm $ 0.6 $ \pm $ 2.3$> $10
      $ \pi(1300)^{0}\pi^{0} $$ 61.1 \pm 3.2 \pm 3.3 \pm 4.0 $$ \phantom{+}0.61 \pm 0.05 \pm 0.08 \pm 0.29 $23.2 $ \pm $ 2.8 $ \pm $ 1.4 $ \pm $ 3.1$> $10
      $ \pi_{2}(1670)^{0}\pi^{0} $$ 12.2 \pm 1.5 \pm 1.5 \pm 2.1 $$ -1.11 \pm 0.14 \pm 0.13 \pm 0.36 $1.1 $ \pm $ 0.2 $ \pm $ 0.2 $ \pm $ 0.36.9
      $ \rho(770)^{0}\rho(770)^{0} $28.0 $ \pm $ 1.9 $ \pm $ 0.6 $ \pm $ 3.0$> $10
      $ [S] $$ 6.1 \pm 1.1 \pm 0.3 \pm 1.4 $$ -3.10 \pm 0.17 \pm 0.11 \pm 0.50 $1.7 $ \pm $ 0.6 $ \pm $ 0.1 $ \pm $ 0.56.5
      $ [P] $$ 6.17 \pm 0.36 \pm 0.13 \pm 0.58 $$ \phantom{+}1.62 \pm 0.07 \pm 0.02 \pm 0.09 $9.8 $ \pm $ 1.0 $ \pm $ 0.4 $ \pm $ 0.8$> $10
      $ [D] $$ 4.54 \pm 0.22 \pm 0.06 \pm 0.34 $$ -3.06 \pm 0.05 \pm 0.01 \pm 0.20 $23.1 $ \pm $ 2.1 $ \pm $ 0.8 $ \pm $ 2.3$> $10
      $ \rho(770)^{0}\rho(1450)^{0} $2.5 $ \pm $ 0.9 $ \pm $ 0.0 $ \pm $ 1.28.0
      $ [P] $$ 13.9 \pm 2.5 \pm 0.7 \pm 1.5 $$ \phantom{+}0.68 \pm 0.20 \pm 0.09 \pm 0.15 $1.0 $ \pm $ 0.4 $ \pm $ 0.1 $ \pm $ 0.66.4
      $ [D] $$ 5.6 \pm 1.3 \pm 0.3 \pm 1.2 $$ \phantom{+}3.08 \pm 0.20 \pm 0.07 \pm 0.43 $1.5 $ \pm $ 0.9 $ \pm $ 0.2 $ \pm $ 1.15.0
      $ \rho(770)^{+}\rho(770)^{-} $90.9 $ \pm $ 3.9 $ \pm $ 3.5 $ \pm $ 6.6$> $10
      $ [S] $$ 13.7 \pm 1.2 \pm 0.6 \pm 1.3 $$ \phantom{+}3.03 \pm 0.09 \pm 0.09 \pm 0.17 $13.0 $ \pm $ 2.0 $ \pm $ 1.2 $ \pm $ 3.2$> $10
      $ [P] $$ 7.10 \pm 0.36 \pm 0.23 \pm 0.35 $$ -1.69 \pm 0.07 \pm 0.02 \pm 0.14 $19.6 $ \pm $ 1.3 $ \pm $ 1.3 $ \pm $ 1.3$> $10
      $ [D] $$ 4.59 \pm 0.22 \pm 0.05 \pm 0.24 $$ \phantom{+}0.06 \pm 0.05 \pm 0.02 \pm 0.09 $36.0 $ \pm $ 3.0 $ \pm $ 0.6 $ \pm $ 2.4$> $10
      $ \rho(770)^{+}\rho(1450)^{-}[D] $$ 8.1 \pm 1.7 \pm 1.4 \pm 2.7 $$ -1.01 \pm 0.18 \pm 0.11 \pm 0.26 $1.7 $ \pm $ 0.8 $ \pm $ 0.6 $ \pm $ 1.76.3
      $ \rho(770)^{0}(\pi\pi)_{S} $2.7 $ \pm $ 0.6 $ \pm $ 0.3 $ \pm $ 1.71.0 $ \pm $ 0.2 $ \pm $ 0.1 $ \pm $ 0.4$> $10
      $ \beta_{1} $$ 8.4 \pm 3.6 \pm 1.1 \pm 1.6 $$ -1.68 \pm 0.50 \pm 0.26 \pm 0.29 $
      $ f_{\pi\pi}^{\rm prod} $$ 40.7 \pm 5.0 \pm 3.7 \pm 3.9 $$ -0.50 \pm 0.14 \pm 0.09 \pm 0.14 $
      $ f_{KK}^{\rm prod} $$ 121 \pm 25 \pm 11 \pm 16 $$ \phantom{+}1.73 \pm 0.23 \pm 0.03 \pm 0.16 $
      $ (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $62.8 $ \pm $ 4.6 $ \pm $ 0.5 $ \pm $ 9.737.4 $ \pm $ 3.0 $ \pm $ 1.8 $ \pm $ 4.8$> $10
      $ a_{1,1} $$ 2224 \pm 35 \pm 26 \pm 33 $$ -1.044 \pm 0.019 \pm 0.008 \pm 0.074 $
      $ a_{1,2} $$ 7287 \pm 62 \pm 46 \pm 93 $$ \phantom{+}1.727 \pm 0.009 \pm 0.004 \pm 0.062 $
      $ b_{2,\pi\pi} $$ 8816 \pm 120 \pm 181\pm 90 $$ -1.107 \pm 0.014 \pm 0.002 \pm 0.082 $
      $ c_{[\pi\pi,\pi\pi]} $$ 2433 \pm 96 \pm 67 \pm 90 $$ \phantom{+}1.796 \pm 0.043 \pm 0.042 \pm 0.044 $
      $ c_{[\pi\pi,KK]} $$ 5417 \pm 477 \pm 47 \pm 462 $$ \phantom{+}2.68 \pm 0.10 \pm 0.06 \pm 0.09 $
      $ f_{2}(1270)^{0}(\pi\pi)_{S} $1.8 $ \pm $ 0.4 $ \pm $ 0.0 $ \pm $ 1.31.1 $ \pm $ 0.2 $ \pm $ 0.0 $ \pm $ 0.79.1
      $ f_{\pi\pi}^{\rm prod} $$ 18.3 \pm 1.8 \pm 1.1 \pm 3.5 $$ -1.39 \pm 0.10 \pm 0.04 \pm 0.20 $
      $ f_{KK}^{\rm prod} $$ 56 \pm 10 \pm 14 \pm 8 $$ \phantom{+}2.29 \pm 0.20 \pm 0.05 \pm 0.39 $
      $ \omega(782)\pi^{0} $$ 1.58 \pm 0.30 \pm 0.05 \pm 0.13 $$ -0.50 \pm 0.44 \pm 0.15 \pm 0.23 $0.9 $ \pm $ 0.4 $ \pm $ 0.0 $ \pm $ 0.26.1
      $ \phi(1020)\pi^{0} $$ 0.44 \pm 0.06 \pm 0.03 \pm 0.05 $$ \phantom{+}2.51 \pm 0.41 \pm 0.10 \pm 0.22 $1.5 $ \pm $ 0.4 $ \pm $ 0.2 $ \pm $ 0.27.4

      Table 8.  The fit parameters, FFs, and significances of individual amplitudes, where the first uncertainties are statistical, the second are the experimental systematic uncertainties and the third are the model-dependent systematic uncertainties.

      AmplitudeMagnitudePhase (rad)Relative FF (%)Significance
      (σ)
      $\pi^{+}\pi^{-}\pi^{+}\pi^{-} $$ \pi^{+}\pi^{-}\pi^{0}\pi^{0} $
      charge $ = \pm 1 $charge $ = \pm 1 $charge $ =0 $
      ${ a_{1}(1260)\to \rho(770)\pi[S]} $1(fixed)0(fixed)79.7±2.2±1.8±3.181.2±2.0±1.6±2.878.9±2.2±1.5±3.5>10
      ${ a_{1}(1260)\to \rho(770)\pi[D] }$${ 0.060 \pm 0.009 \pm 0.002 \pm 0.022 }$${ \phantom{+}0.01 \pm 0.17 \pm 0.09 \pm 0.09 }$1.3±0.3±0.1±0.71.2±0.3±0.1±0.61.3±0.3±0.1±0.77.5
      $ a_{1}(1260)\to f_{2}(1270)\pi[P] $$ 0.311 \pm 0.033 \pm 0.019 \pm 0.021 $$ -1.58 \pm 0.11 \pm 0.08 \pm 0.25 $1.8±0.4±0.3±0.50.8±0.2±0.1±0.31.2±0.2±0.1±0.3>10
      $ a_{1}(1260)\to (\pi^{+}\pi^{-})_{S}\pi[P] $5.4±0.6±0.5±0.83.1±0.4±0.4±0.73.9±0.5±0.5±0.9>10
      $ \beta_{1} $$ 0.83 \pm 0.13 \pm 0.11 \pm 0.34 $$ -2.18 \pm 0.14 \pm 0.19 \pm 0.16 $
      $ f_{\pi\pi}^{\rm prod} $$ 2.47 \pm 0.16 \pm 0.12 \pm 0.12 $$ \phantom{+}0.34 \pm 0.08 \pm 0.06 \pm 0.12 $
      $ f_{KK}^{\rm prod} $$ 6.59 \pm 0.84 \pm 0.85 \pm 1.24 $$ \phantom{+}1.76 \pm 0.12 \pm 0.08 \pm 0.13 $
      $ \pi(1300)\to \rho(770)\pi $1(fixed)0(fixed)53.8±2.9±0.9±8.779.0±2.0±0.8±6.873.6±2.6±0.8±9.1>10
      $ \pi(1300)\to (\pi^{+}\pi^{-})_{S}\pi $51.1±2.9±0.9±8.426.2±2.1±0.8±6.536.7±2.6±1.0±8.3>10
      $ \beta_{1} $$ 5.0 \pm 0.5 \pm 0.2 \pm 0.2 $$ -0.64 \pm 0.10 \pm 0.04 \pm 0.25 $
      $ f_{KK}^{\rm prod} $$ 43.8 \pm 2.8 \pm 0.5 \pm 2.7 $$ \phantom{+}0.34 \pm 0.06 \pm 0.02 \pm 0.06 $
      $ a_{1}(1420)\to f_{0}(980)\pi[P] $1(fixed)0(fixed)1001006.0
      $ a_{1}(1640)\to \rho(770)\pi[S] $1(fixed)0(fixed)1001009.1
      $ a_{2}(1320)\to \rho(770)\pi[D] $1(fixed)0(fixed)1001007.3
      $ h_{1}(1170)\to \rho(770)\pi[S] $1(fixed)0(fixed)1006.5
      $ \pi_{2}(1670)\to f_2(1270)\pi[S] $1(fixed)0(fixed)1006.9

      Table 9.  The fit parameters, FFs and significance for the three-body decays of $ a_{1}(1260) $, $ a_1(1420) $, $ a_1(1640) $, $ a_2(1320) $, $ h_1(1170) $, $ \pi(1300) $, and $ \pi_2(1670) $. The first uncertainties are the statistical uncertainties, the second are the experimental systematic uncertainties, and the third are the model-dependent systematic uncertainties.

      This workPDG
      ${\rm Mass}/ ({\rm GeV}/c^{2})$${\rm Width} / {\rm GeV}$${\rm Mass} /({\rm GeV}/c^{2})$${\rm Width}/{\rm GeV}$
      $ a_1(1260) $$ 1.193 \pm 0.005 \pm 0.003 \pm 0.023 $$ 0.487 \pm 0.009 \pm 0.015 \pm 0.039 $$ 1.230 \pm 0.040 $$ 0.250-0.600 $
      $ \pi(1300) $$ 1.534 \pm 0.011 \pm 0.009 \pm 0.020 $$ 0.610 \pm 0.030 \pm 0.021 \pm 0.090 $$ 1.300 \pm 0.100 $$ 0.200-0.600 $

      Table 10.  Masses and widths of $ a_1(1260) $ and $ \pi(1300) $ obtained via parameter scans and presented in the PDG [25]. The first uncertainties in this fit are the statistical uncertainties, the second are the experimental systematic uncertainties, and the third are the model-dependent systematic uncertainties.

      Figure 5.  Likelihood scans for masses and widths of $ a_{1}(1260) $ and $ \pi(1300) $.

      FF(%)1234567891011121314
      1. $D^{0}\to a_{1}(1260)^{+}\pi^{-}$82.2 $\pm$ 3.39.9 $\pm$ 0.80.5 $\pm$ 0.3–15.3 $\pm$ 2.5–0.4 $\pm$ 0.70.0 $\pm$ 0.00.0 $\pm$ 0.01.7 $\pm$ 0.2–27.3 $\pm$ 2.1–4.6 $\pm$ 1.9–5.2 $\pm$ 1.5–1.3 $\pm$ 2.35.2 $\pm$ 1.1–1.4 $\pm$ 0.2
      2. $D^{0}\to a_{1}(1260)^{-}\pi^{+}$10.3 $\pm$ 1.5–0.1 $\pm$ 0.1–1.4 $\pm$ 0.3–2.5 $\pm$ 0.80.0 $\pm$ 0.00.0 $\pm$ 0.0–10.8 $\pm$ 1.2–0.1 $\pm$ 0.1–1.8 $\pm$ 0.7–1.5 $\pm$ 0.5–0.3 $\pm$ 0.83.4 $\pm$ 0.5–0.4 $\pm$ 0.1
      3. $D^{0}\to a_{1}(1420)^{+}\pi^{-}$0.6 $\pm$ 0.20.1 $\pm$ 0.1–0.0 $\pm$ 0.00.0 $\pm$ 0.00.0 $\pm$ 0.0–0.2 $\pm$ 0.00.1 $\pm$ 0.1–0.3 $\pm$ 0.1–0.1 $\pm$ 0.1–0.2 $\pm$ 0.11.0 $\pm$ 0.2–0.2 $\pm$ 0.0
      4. $D^{0}\to a_{1}(1640)^{+}\pi^{-}$1.7 $\pm$ 0.50.2 $\pm$ 0.1–0.0 $\pm$ 0.0–0.0 $\pm$ 0.00.2 $\pm$ 0.14.9 $\pm$ 0.7–1.1 $\pm$ 0.30.3 $\pm$ 0.10.1 $\pm$ 0.3–1.6 $\pm$ 0.60.1 $\pm$ 0.0
      5. $D^{0}\to a_{1}(1640)^{-}\pi^{+}$0.5 $\pm$ 0.3–0.0 $\pm$ 0.0–0.0 $\pm$ 0.02.7 $\pm$ 0.80.3 $\pm$ 0.1–0.2 $\pm$ 0.20.1 $\pm$ 0.10.4 $\pm$ 0.3–1.9 $\pm$ 0.5–0.0 $\pm$ 0.0
      6. $D^{0}\to a_{2}(1320)^{+}\pi^{-}$0.2 $\pm$ 0.1–0.1 $\pm$ 0.00.0 $\pm$ 0.00.0 $\pm$ 0.0–0.8 $\pm$ 0.30.2 $\pm$ 0.1–0.0 $\pm$ 0.0–0.0 $\pm$ 0.0–0.0 $\pm$ 0.0
      7. $D^{0}\to a_{2}(1320)^{-}\pi^{+}$0.3 $\pm$ 0.1–0.0 $\pm$ 0.0–0.0 $\pm$ 0.01.0 $\pm$ 0.2–0.2 $\pm$ 0.10.0 $\pm$ 0.00.0 $\pm$ 0.00.0 $\pm$ 0.0
      8. $D^{0}\to \pi(1300)^{+}\pi^{-}$32.3 $\pm$ 2.69.2 $\pm$ 1.5–9.3 $\pm$ 0.7–2.8 $\pm$ 0.84.0 $\pm$ 1.6–49.2 $\pm$ 4.10.3 $\pm$ 0.3
      9. $D^{0}\to \pi(1300)^{-}\pi^{+}$23.5 $\pm$ 2.3–6.0 $\pm$ 0.7–2.6 $\pm$ 0.7–2.8 $\pm$ 1.5–39.8 $\pm$ 3.9–0.1 $\pm$ 0.3
      10. $D^{0}\to \rho(770)^{0}\rho(770)^{0}$28.0 $\pm$ 1.90.0 $\pm$ 2.10.0 $\pm$ 0.02.0 $\pm$ 0.5–1.2 $\pm$ 0.2
      11. $D^{0}\to \rho(770)^{0}\rho(1450)^{0}$2.5 $\pm$ 0.90.0 $\pm$ 0.0–0.6 $\pm$ 0.3–0.3 $\pm$ 0.1
      12. $D^{0}\to \rho(770)^{0}(\pi\pi)_{S}$2.7 $\pm$ 0.60.0 $\pm$ 0.0–0.0 $\pm$ 0.0
      13. $D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S}$62.8 $\pm$ 4.6–1.3 $\pm$ 0.4
      14. $D^{0}\to f_{2}(1270)^{0}(\pi\pi)_{S}$1.8 $\pm$ 0.4

      Table 11.  The fit fractions of interference terms in $D^0\to\pi^+\pi^-\pi^+\pi^-$ (first) and $D^0\to\pi^+\pi^-\pi^0\pi^0$ (second). The uncertainties are statistical only.

      FF(%)1234567891011121314151617181920
      1. $D^{0}\to a_{1}(1260)^{+}\pi^{-}$57.5±2.74.6±0.48.5±0.60.7±0.2–10.2±1.7–0.3±0.30.0±0.0–0.0±0.0–1.6±0.4–0.0±0.00.8±0.1–18.6±1.80.6±0.2–33.5±1.8–1.4±1.2–1.8±0.32.4±0.8–0.4±0.10.0±0.0–0.0±0.0
      2. $D^{0}\to a_{1}(1260)^{-}\pi^{+}$7.2±1.13.1±0.3–0.0±0.0–0.6±0.1–1.8±0.5–0.0±0.00.0±0.00.6±0.10.3±0.0–0.0±0.0–7.6±1.00.1±0.1–12.1±1.10.4±0.10.6±0.12.0±0.4–0.1±0.00.0±0.0–0.0±0.0
      3. $D^{0}\to a_{1}(1260)^{0}\pi^{0}$32.9±3.2–0.1±0.0–1.5±0.3–0.2±0.30.0±0.00.0±0.00.0±0.0–11.4±1.1–12.0±1.00.7±0.10.3±0.1–34.3±2.60.3±0.80.0±0.04.8±1.1–0.3±0.1–0.0±0.0–0.0±0.0
      4. $D^{0}\to a_{1}(1420)^{+}\pi^{-}$0.3±0.10.0±0.1–0.0±0.0–0.0±0.00.0±0.0–0.0±0.00.0±0.00.1±0.0–0.3±0.1–0.0±0.0–0.4±0.10.0±0.0–0.1±0.00.7±0.10.0±0.0–0.0±0.00.0±0.0
      5. $D^{0}\to a_{1}(1640)^{+}\pi^{-}$1.1±0.30.1±0.0–0.0±0.00.0±0.00.3±0.1–0.0±0.0–0.1±0.03.7±0.60.0±0.03.0±0.50.2±0.10.1±0.0–1.3±0.40.0±0.00.0±0.00.0±0.0
      6. $D^{0}\to a_{1}(1640)^{-}\pi^{+}$0.3±0.20.0±0.0–0.0±0.0–0.0±0.1–0.1±0.00.0±0.01.7±0.50.0±0.01.1±0.4–0.0±0.0–0.0±0.0–1.4±0.4–0.0±0.0–0.0±0.0–0.0±0.0
      7. $D^{0}\to a_{2}(1320)^{+}\pi^{-}$0.2±0.1–0.0±0.00.0±0.00.0±0.0–0.0±0.00.0±0.00.0±0.0–0.6±0.20.0±0.00.0±0.0–0.0±0.0–0.0±0.00.0±0.0–0.0±0.0
      8. $D^{0}\to a_{2}(1320)^{-}\pi^{+}$0.3±0.1–0.0±0.00.0±0.0–0.0±0.0–0.0±0.00.0±0.01.0±0.20.0±0.0–0.0±0.0–0.0±0.0–0.0±0.00.0±0.0–0.0±0.0
      9. $D^{0}\to h_{1}(1170)^{0}\pi^{0}$1.3±0.61.5±0.5–1.7±0.4–0.0±0.00.0±0.00.0±0.00.1±0.10.3±0.20.0±0.0–0.0±0.00.0±0.00.0±0.0
      10. $D^{0}\to \pi(1300)^{+}\pi^{-}$15.6±1.44.2±0.4–3.7±0.4–0.2±0.0–5.3±0.7–1.4±0.71.1±0.2–14.6±1.20.2±0.10.0±0.0–0.0±0.0
      11. $D^{0}\to \pi(1300)^{-}\pi^{+}$11.4±1.1–3.1±0.4–0.1±0.0–1.9±0.60.3±0.2–1.0±0.2–11.1±1.10.1±0.10.0±0.0–0.0±0.0
      12. $D^{0}\to \pi(1300)^{0}\pi^{0}$23.2±2.80.4±0.1–5.3±1.2–1.1±0.60.0±0.0–21.5±2.20.2±0.2–0.0±0.00.0±0.0
      13. $D^{0}\to \pi_{2}(1670)^{0}\pi^{0}$1.1±0.2–0.6±0.2–0.0±0.00.0±0.0–0.0±0.0–0.5±0.20.0±0.0–0.0±0.0
      14. $D^{0}\to \rho(770)^{+}\rho(770)^{-}$90.9±3.9–5.1±3.40.0±0.0–2.1±1.0–2.4±0.30.0±0.0–0.0±0.0
      15. $D^{0}\to \rho(770)^{+}\rho(1450)^{-}[D]$1.7±0.80.0±0.0–1.0±0.30.2±0.1–0.0±0.00.0±0.0
      16. $D^{0}\to \rho(770)^{0}(\pi\pi)_{S}$1.0±0.20.0±0.00.0±0.00.0±0.0–0.0±0.0
      17. $D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S}$37.4±3.0–0.0±0.0–0.0±0.00.0±0.0
      18. $D^{0}\to f_{2}(1270)^{0}(\pi\pi)_{S}$1.1±0.20.0±0.0–0.0±0.0
      19. $D^{0}\to \omega(782)\pi^{0}$0.9±0.4–0.0±0.0
      20. $D^{0}\to \phi(1020)\pi^{0}$1.5±0.4

      Figure 6.  (color online) The mass and angular distributions for $ D^0\to \pi^+\pi^-\pi^+\pi^- $, where $ \theta_{a}^{ab} $ is the helicity angle of a in the $ ab $ system, $ \phi_{ab,cd} $ is the angle between the decay planes of the $ ab $ and $ cd $ systems in the $ D^0 $ rest frame.

      Figure 7.  (color online) The mass and angular distributions in $ D^0\to \pi^{+}\pi^{-}\pi^{0}\pi^{0} $, where $ \theta_{a}^{ab} $ is the helicity angle of a in the $ ab $ system and $ \phi_{ab,cd} $ is the angle between the decay planes of the $ ab $ and $ cd $ systems in the $ D^0 $ rest frame.

    VII.   SYSTEMATIC UNCERTAINTY OF THE AMPLITUDE ANALYSIS
    • The systematic uncertainties of the amplitude analysis result from two aspects. One is the experimental systematic uncertainty, including those from background estimation, detection efficiency over PHSP, and fit bias. Another is the model-dependent systematic uncertainty, including those from the resonance line shape, the radii of Blatt-Weisskopf barrier factors, quantum correlation parameters, and extra amplitudes. The fit is performed with alternative conditions to estimate these uncertainties for each source, and the deviations from the nominal results are taken as the corresponding uncertainties. Table 12 and 13 summarize the systematic uncertainties on the magnitude and phase of the fit parameters in the unit of the statistical uncertainty, respectively. Table 15 summarized the systematic uncertainties on the FFs, resonance parameters, and $ CP $-even fractions in the unit of the statistical uncertainty. The total systematic uncertainties are the square roots of the quadrature sums of the individual contributions. The individual uncertainties are obtained as follows:

      $ F_{+}^{\pi^{+}\pi^{-}\pi^{+}\pi^{-}} $$ F_{+}^{\pi^{+}\pi^{-}\pi^{0}\pi^{0}} $(non-η)
      This work (model-dependent)$ (75.2 \pm 1.1_{\rm stat.} \pm 1.5_{\rm syst.}) $%$ (68.9 \pm 1.5_{\rm stat.} \pm 2.4_{\rm syst.}) $%
      CLEO-c (model-dependent)$ (72.9 \pm 0.9_{\rm stat.} \pm 1.5_{\rm syst.} \pm 1.0_{\rm model}) $% [13]
      CLEO-c (model-independent, global)$ (73.7 \pm 2.8) $% [52]
      CLEO-c (model-independent, binned)$ (76.9 \pm 2.1_{\rm stat.} \pm 1.0_{\rm syst.} \pm 0.2_{K_{S} {\rm veto}}) $% [6]
      BESIII (model-independent, global)$ (73.4 \pm 1.5_{\rm stat.} \pm 0.8_{\rm syst.}) $% [53]$ (68.2 \pm 7.7) $% [54]

      Table 12.  The $ CP $-even fractions obtained in this work and comparisons with the CLEO-c and prior BESIII measurements.

      1234567Total
      magnitude($ D^{0}\to a_{1}(1260)^{-}\pi^{+} $)0.20.10.21.60.10.70.3 $ \pm $ 1.7
      magnitude($ D^{0}\to a_{1}(1260)^{0}\pi^{0} $)0.10.00.11.30.00.80.1 $ \pm $ 1.5
      magnitude($ D^{0}\to a_{1}(1420)^{+}\pi^{-} $)0.20.00.31.00.10.50.4 $ \pm $ 1.1
      magnitude($ D^{0}\to a_{1}(1640)^{+}\pi^{-} $)0.80.20.21.40.21.30.9 $ \pm $ 1.9
      magnitude($ D^{0}\to a_{1}(1640)^{-}\pi^{+} $)0.10.10.20.60.21.20.2 $ \pm $ 1.4
      magnitude($ D^{0}\to a_{2}(1320)^{+}\pi^{-} $)0.20.10.30.70.10.20.4 $ \pm $ 0.7
      magnitude($ D^{0}\to a_{2}(1320)^{-}\pi^{+} $)0.00.00.20.60.10.50.2 $ \pm $ 0.8
      magnitude($ D^{0}\to h_{1}(1170)^{0}\pi^{0} $)0.60.10.11.30.01.00.7 $ \pm $ 1.7
      magnitude($ D^{0}\to \pi(1300)^{+}\pi^{-} $)1.30.00.01.10.10.41.3 $ \pm $ 1.2
      magnitude($ D^{0}\to \pi(1300)^{-}\pi^{+} $)0.80.10.61.10.10.11.0 $ \pm $ 1.1
      magnitude($ D^{0}\to \pi(1300)^{0}\pi^{0} $)1.00.20.11.20.20.31.0 $ \pm $ 1.2
      magnitude($ D^{0}\to \pi_{2}(1670)^{0}\pi^{0} $)0.70.10.60.30.11.41.0 $ \pm $ 1.4
      magnitude($ D^{0}\to \rho(770)^{0}\rho(770)^{0}[S] $)0.00.00.31.20.00.40.3 $ \pm $ 1.3
      magnitude($ D^{0}\to \rho(770)^{0}\rho(770)^{0}[P] $)0.30.00.21.50.40.40.4 $ \pm $ 1.6
      magnitude($ D^{0}\to \rho(770)^{0}\rho(770)^{0}[D] $)0.10.00.21.50.00.20.3 $ \pm $ 1.5
      magnitude($ D^{0}\to \rho(770)^{0}\rho(1450)^{0}[P] $)0.00.20.20.50.30.20.3 $ \pm $ 0.6
      magnitude($ D^{0}\to \rho(770)^{0}\rho(1450)^{0}[D] $)0.20.10.10.30.10.90.2 $ \pm $ 0.9
      magnitude($ D^{0}\to \rho(770)^{+}\rho(770)^{-}[S] $)0.50.10.10.60.10.90.5 $ \pm $ 1.1
      magnitude($ D^{0}\to \rho(770)^{+}\rho(770)^{-}[P] $)0.60.10.00.70.00.60.6 $ \pm $ 1.0
      magnitude($ D^{0}\to \rho(770)^{+}\rho(770)^{-}[D] $)0.20.10.10.90.00.60.2 $ \pm $ 1.1
      magnitude($ D^{0}\to \rho(770)^{+}\rho(1450)^{-}[D] $)0.80.10.11.40.10.70.8 $ \pm $ 1.6
      magnitude($ D^{0}\to \rho(770)^{0}(\pi\pi)_{S} $, $ \beta_1 $)0.20.00.20.20.00.40.3 $ \pm $ 0.4
      magnitude($ D^{0}\to \rho(770)^{0}(\pi\pi)_{S} $, $ f_{\pi\pi}^{\rm prod} $)0.60.00.40.40.10.70.7 $ \pm $ 0.8
      magnitude($ D^{0}\to \rho(770)^{0}(\pi\pi)_{S} $, $ f_{KK}^{\rm prod} $)0.40.00.00.20.10.60.4 $ \pm $ 0.7
      magnitude($ D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $, $ a_{1,1} $)0.70.10.00.50.10.80.7 $ \pm $ 0.9
      magnitude($ D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $, $ a_{1,2} $)0.60.20.31.20.20.90.7 $ \pm $ 1.5
      magnitude($ D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $, $ b_{2,\pi\pi} $)1.50.00.00.30.10.71.5 $ \pm $ 0.7
      magnitude($ D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $, $ c_{[\pi\pi,\pi\pi]} $)0.60.10.20.40.10.80.7 $ \pm $ 0.9
      magnitude($ D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $, $ c_{[\pi\pi,KK]} $)0.10.00.10.60.20.80.1 $ \pm $ 1.0
      magnitude($ D^{0}\to f_{2}(1270)^{0}(\pi\pi)_{S} $, $ f_{\pi\pi}^{\rm prod} $)0.50.00.30.40.11.90.6 $ \pm $ 1.9
      magnitude($ D^{0}\to f_{2}(1270)^{0}(\pi\pi)_{S} $, $ f_{KK}^{\rm prod} $)1.20.10.70.30.10.71.4 $ \pm $ 0.8
      magnitude($ D^{0}\to \omega(782)\pi^{0} $)0.10.00.10.20.00.40.2 $ \pm $ 0.4
      magnitude($ D^{0}\to \phi(1020)\pi^{0} $)0.20.10.40.60.10.50.5 $ \pm $ 0.8
      magnitude($ a_{1}(1260)\to \rho(770)\pi[D] $)0.10.00.12.40.00.20.2 $ \pm $ 2.4
      magnitude($ a_{1}(1260)\to f_{2}(1270)\pi[P] $)0.60.00.10.50.10.30.6 $ \pm $ 0.6
      magnitude($ a_{1}(1260)\to (\pi^{+}\pi^{-})_{S}\pi[P] $, $ \beta_1 $)0.90.00.00.60.02.50.9 $ \pm $ 2.6
      magnitude($ a_{1}(1260)\to (\pi^{+}\pi^{-})_{S}\pi[P] $, $ f_{\pi\pi}^{\rm prod} $)0.70.00.00.70.10.20.7 $ \pm $ 0.7
      magnitude($ a_{1}(1260)\to (\pi^{+}\pi^{-})_{S}\pi[P] $, $ f_{KK}^{\rm prod} $)1.00.00.20.80.11.31.0 $ \pm $ 1.5
      magnitude($ \pi(1300)\to (\pi^{+}\pi^{-})_{S}\pi $, $ \beta_1 $)0.20.10.30.40.00.10.3 $ \pm $ 0.5
      magnitude($ \pi(1300)\to (\pi^{+}\pi^{-})_{S}\pi $, $ f_{KK}^{\rm prod} $)0.20.00.10.80.10.50.2 $ \pm $ 1.0

      Table 13.  Systematic uncertainties in the magnitude of fit parameters in units of statistical standard deviations. 1: Background estimation. 2: Detection efficiency. 3: Fit bias. 4: Resonance parameters. 5: Radii of Blatt-Weisskopf barrier factors. 6: Quantum correlation correction. 7: Extra amplitudes. In the "Total" column, the term before "$ \pm $" is the total experimental systematic uncertainty (1, 2, and 3), and the term after "$ \pm $" is the total model-dependent systematic uncertainty (4, 5, 6, and 7). The systematic uncertainties from "5"are not considered here.

      1234567Total
      phase($ D^{0}\to a_{1}(1260)^{-}\pi^{+} $)0.20.00.02.20.20.10.20.2 $ \pm $ 2.2
      phase($ D^{0}\to a_{1}(1260)^{0}\pi^{0} $)1.40.10.10.71.80.01.31.4 $ \pm $ 2.3
      phase($ D^{0}\to a_{1}(1420)^{+}\pi^{-} $)0.00.00.36.00.60.00.60.3 $ \pm $ 6.1
      phase($ D^{0}\to a_{1}(1640)^{+}\pi^{-} $)0.00.00.11.60.60.00.20.1 $ \pm $ 1.7
      phase($ D^{0}\to a_{1}(1640)^{-}\pi^{+} $)0.80.00.21.11.20.20.60.8 $ \pm $ 1.7
      phase($ D^{0}\to a_{2}(1320)^{+}\pi^{-} $)0.50.10.10.60.20.00.50.5 $ \pm $ 0.8
      phase($ D^{0}\to a_{2}(1320)^{-}\pi^{+} $)0.30.00.00.50.30.10.40.3 $ \pm $ 0.7
      phase($ D^{0}\to h_{1}(1170)^{0}\pi^{0} $)0.30.10.20.40.50.11.10.3 $ \pm $ 1.3
      phase($ D^{0}\to \pi(1300)^{+}\pi^{-} $)0.90.00.06.22.40.11.30.9 $ \pm $ 6.8
      phase($ D^{0}\to \pi(1300)^{-}\pi^{+} $)1.80.10.13.62.50.01.41.9 $ \pm $ 4.6
      phase($ D^{0}\to \pi(1300)^{0}\pi^{0} $)1.60.10.05.12.20.11.61.6 $ \pm $ 5.8
      phase($ D^{0}\to \pi_{2}(1670)^{0}\pi^{0} $)0.90.00.01.70.60.11.80.9 $ \pm $ 2.6
      phase($ D^{0}\to \rho(770)^{0}\rho(770)^{0}[S] $)0.60.00.12.90.60.20.40.6 $ \pm $ 2.9
      phase($ D^{0}\to \rho(770)^{0}\rho(770)^{0}[P] $)0.10.20.11.20.00.20.50.2 $ \pm $ 1.3
      phase($ D^{0}\to \rho(770)^{0}\rho(770)^{0}[D] $)0.20.00.13.70.70.11.00.2 $ \pm $ 3.9
      phase($ D^{0}\to \rho(770)^{0}\rho(1450)^{0}[P] $)0.40.20.10.70.10.10.10.4 $ \pm $ 0.8
      phase($ D^{0}\to \rho(770)^{0}\rho(1450)^{0}[D] $)0.40.00.01.71.00.20.60.4 $ \pm $ 2.1
      phase($ D^{0}\to \rho(770)^{+}\rho(770)^{-}[S] $)1.00.00.10.30.80.11.71.0 $ \pm $ 1.9
      phase($ D^{0}\to \rho(770)^{+}\rho(770)^{-}[P] $)0.10.10.31.80.80.20.40.3 $ \pm $ 2.1
      phase($ D^{0}\to \rho(770)^{+}\rho(770)^{-}[D] $)0.40.00.11.70.10.10.40.4 $ \pm $ 1.8
      phase($ D^{0}\to \rho(770)^{+}\rho(1450)^{-}[D] $)0.60.10.01.00.90.00.50.6 $ \pm $ 1.4
      phase($ D^{0}\to \rho(770)^{0}(\pi\pi)_{S} $, $ \beta_1 $)0.50.10.20.10.50.00.30.5 $ \pm $ 0.6
      phase($ D^{0}\to \rho(770)^{0}(\pi\pi)_{S} $, $ f_{\pi\pi}^{\rm prod} $)0.60.00.20.30.30.00.90.6 $ \pm $ 1.0
      phase($ D^{0}\to \rho(770)^{0}(\pi\pi)_{S} $, $ f_{KK}^{\rm prod} $)0.10.00.00.20.50.00.30.1 $ \pm $ 0.7
      phase($ D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $, $ a_{1,1} $)0.40.00.10.43.80.00.50.4 $ \pm $ 3.9
      phase($ D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $, $ a_{1,2} $)0.40.00.01.06.60.11.50.4 $ \pm $ 6.9
      phase($ D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $, $ b_{2,\pi\pi} $)0.10.00.00.45.60.11.50.2 $ \pm $ 5.8
      phase($ D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $, $ c_{[\pi\pi,\pi\pi]} $)1.00.00.10.30.70.00.71.0 $ \pm $ 1.0
      phase($ D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $, $ c_{[\pi\pi,KK]} $)0.40.00.30.50.50.10.60.6 $ \pm $ 0.9
      phase($ D^{0}\to f_{2}(1270)^{0}(\pi\pi)_{S} $, $ f_{\pi\pi}^{\rm prod} $)0.40.10.10.41.40.01.30.4 $ \pm $ 2.0
      phase($ D^{0}\to f_{2}(1270)^{0}(\pi\pi)_{S} $, $ f_{KK}^{\rm prod} $)0.20.00.20.31.60.11.10.2 $ \pm $ 2.0
      phase($ D^{0}\to \omega(782)\pi^{0} $)0.30.00.10.50.30.00.00.3 $ \pm $ 0.5
      phase($ D^{0}\to \phi(1020)\pi^{0} $)0.10.00.20.30.40.00.10.2 $ \pm $ 0.5
      phase($ a_{1}(1260)\to \rho(770)\pi[D] $)0.50.10.10.40.30.10.00.5 $ \pm $ 0.5
      phase($ a_{1}(1260)\to f_{2}(1270)\pi[P] $)0.70.00.12.20.20.00.70.7 $ \pm $ 2.3
      phase($ a_{1}(1260)\to (\pi^{+}\pi^{-})_{S}\pi[P] $, $ \beta_1 $)1.30.20.30.60.90.20.31.4 $ \pm $ 1.1
      phase($ a_{1}(1260)\to (\pi^{+}\pi^{-})_{S}\pi[P] $, $ f_{\pi\pi}^{\rm prod} $)0.70.10.10.50.50.11.20.7 $ \pm $ 1.5
      phase($ a_{1}(1260)\to (\pi^{+}\pi^{-})_{S}\pi[P] $, $ f_{KK}^{\rm prod} $)0.60.00.00.30.40.10.90.6 $ \pm $ 1.0
      phase($ \pi(1300)\to (\pi^{+}\pi^{-})_{S}\pi $, $ \beta_1 $)0.30.10.20.32.40.10.80.4 $ \pm $ 2.5
      phase($ \pi(1300)\to (\pi^{+}\pi^{-})_{S}\pi $, $ f_{KK}^{\rm prod} $)0.20.00.20.50.30.10.80.3 $ \pm $ 1.0

      Table 14.  Systematic uncertainties in the phase of fit parameters in units of statistical standard deviations. 1: Background estimation. 2: Detection efficiency. 3: Fit bias. 4: Resonance parameters. 5: Radii of Blatt-Weisskopf barrier factors. 6: Quantum correlation correction. 7: Extra amplitudes. In the "Total" column, the term before "$ \pm $" is the total experimental systematic uncertainty (1, 2, and 3), and the term after "$ \pm $" is the total model-dependent systematic uncertainty (4, 5, 6, and 7).

      ● Background estimation

      The sPlot technique is used to estimate the signal weight and background; the corresponding sources of systematic uncertainty are all the PDFs and the magnitudes of peaking backgrounds in BKGVI in the 2D fit on the $ M_{\rm bc}^{\rm tag} $ versus $ M_{\rm bc}^{\rm sig} $ distribution. To estimate the corresponding uncertainties, an alternative fit on the $ M_{\rm bc}^{\rm tag} $ versus $ M_{\rm bc}^{\rm sig} $ distribution is performed by varying the means and widths of two Gaussian functions by $ \; \pm\; 1\sigma $ for the signal, varying the fixed parameters of the ARGUS function by $ \; \pm\; 1\sigma $ for BKGI and BKGII, changing the Student's function to the bifurcated Student's function with different n and σ on the left and right sides of the maximum value and two ARGUS functions to one ARGUS function on the ($ M_{\rm bc}^{\rm sig}+M_{\rm bc}^{\rm tag} $) dimension for BKGIII, varying the means and widths of the bifurcated Gaussian functions by $ \; \pm\; 1\sigma $ for BKGIV and BKGV, and varying the fixed peaking background yields in BKGVI by $ \; \pm\; 1\sigma $ according to the estimated values in Table 2. The output results are used in the sPlot technique and amplitude analysis.

      ● Detection efficiency

      The systematic uncertainties due to detection efficiencies are from the $ \pi^{\pm} $ tracking/PID and $ \pi^0 $ reconstruction. The corresponding uncertainties are obtained by weighting the PHSP signal MC sample with a factor $ \epsilon_{\rm Data}/\epsilon_{\rm MC} $ when normalizing the PDF in the amplitude analysis, where $ \epsilon_{\rm Data} $ and $ \epsilon_{\rm MC} $ are the efficiencies of the data and MC simulation, respectively. The $ \pi^{\pm} $ tracking/PID efficiencies are quoted from Ref. [55], while the efficiency of the $ \pi^0 $ reconstruction is studied with the control samples of $ D^0\to K^-\pi^+\pi^0 $ and $ D^0\to\pi^+\pi^-\pi^0 $ versus the three tag modes used in this analysis.

      ● Fit bias

      The sPlot technique, MC integration, detection resolution, and many other potential problems from fit tools may lead to fit bias. To estimate this fit bias, 100 sets of samples with the same size and signal purity as the data are generated, where the signal events are generated according to the the nominal model and the background events are from the inclusive MC sample. The amplitude fit is performed on these samples, and the average deviations of the fitted parameters from their input values are taken as uncertainties due to fit bias.

      ● The resonance parameters and models

      The systematic uncertainties associated with the masses and widths of resonances in the amplitude analysis are estimated by shifting the corresponding fixed PDG [25] values or optimized values by $ \; \pm\; 1\sigma $. The systematic uncertainties associated with the model of $ \pi\pi $ S-wave are estimated by replacing the K-matrix formula with the sum of three independent resonances of $ f_0(500) $, $ f_0(980) $, and $ f_0(1370) $, described by the formula in Ref. [56], a Flatt$ \rm{\acute{e}} $ parametrization, and a relative Breit-Wigner function. For the magnitude and phase of the fit parameters of $ \pi\pi $ S-wave, the systematic uncertainties associated with the $ \pi\pi $ S-wave model are not considered. The overall uncertainties are the square roots of the quadrature sums of the individual contributions.

      ● Radii of Blatt-Weisskopf barrier factors

      The systematic uncertainties due to the radii of Blatt-Weisskopf barrier factors are obtained by varying the radii of the $ D^{0} $ meson and other resonances by $ \; \pm\; 1\; {\rm{GeV}^{-1}}c $ in the amplitude analysis, and the square root of the quadrature sum of the individual effects is taken as the uncertainty. The radii choice significantly affects the magnitudes of the fit parameters, especially for the components with high orbital angular momentum. Therefore, this part is not included in the systematic uncertainties of the magnitudes of the fit parameters.

      ● Quantum correlation correction

      The systematic uncertainty due to the quantum correlation correction is estimated by varying the input quantum correlation parameters by $ \; \pm\; 1\sigma $ in the amplitude analysis.

      ● Extra amplitudes

      The systematic uncertainties due to the extra amplitudes are estimated with the alternative fits with the models, including an additional amplitude of $ f_2(1270)f_2(1270)[D] $ or $ 1^{++}[\rho(770)^+\pi^{-} [S]]_{\rm NR}\pi^{0} $ ("NR" represents non-resonant contribution), which are of 3σ to 5σ significance based on the nominal model. The square root of the quadrature sum of the individual effects is taken as the uncertainty.

    VIII.   MEASUREMENT OF BRANCHING FRACTIONS
    • The absolute branching fractions of $ D^0\to \pi^+\pi^-\pi^+\pi^- $ and $ D^0\to\pi^+\pi^-\pi^0\pi^0 $(non-η) are measured using the DT method. The numbers of ST events ($ N_{g}^{\rm ST} $) for the tag mode g and DT events ($ N_{fg}^{\rm DT} $) for the self-conjugated signal model f with the tag mode g are given by

      $ N_{g}^{\rm ST} =2N_{D^{0}\bar{D}^{0}}\mathcal{B}_{g}\epsilon_{g}^{\rm ST}(1+y_D^2)(1+r_{g}^{2} -2r_{g}R_{g}y_D\mathrm{cos}\delta_{g}), $

      (43)

      $ N_{fg}^{\rm DT} = 2N_{D^{0}\bar{D}^{0}}\mathcal{B}_{f}\mathcal{B}_{g}\epsilon_{fg}^{\rm DT}(1+y_D^2)[1 + r_{g}^{2} - 2r_{g}R_{g}\mathrm{cos}\delta_{g}(2F_+^{f} - 1)]. $

      (44)

      Here, $ N_{D^{0}\bar{D}^{0}} $ is the total number of $ D^{0}\bar{D}^{0} $ pairs in the data, $ \mathcal{B}_{f} $ and $ \mathcal{B}_{g} $ are the branching fractions of the signal mode f and tag mode g, respectively, $ \epsilon_{f}^{\rm ST} $ and $ \epsilon_{fg}^{\rm ST} $ are the corresponding ST and DT efficiencies, and $ y_D $ is the $ D^0 $-$ \bar{D}^0 $ mixing parameter. By combining Eqs. (43) and (44) and ignoring the term $ 2r_{g}R_{g}y_D\mathrm{cos}\delta_{g} $ in Eq. (43), the branching fraction of $ D^0\to f $ is given as

      $ \mathcal{B}_{f} = \frac{\sum _{g}N_{fg}^{\rm DT}}{\sum _{g}N_{g}^{\rm ST}(\epsilon_{fg}^{\rm DT}/\epsilon_{g}^{\rm ST})\left[1-\dfrac{2r_{g}R_{g}\mathrm{cos}\delta_{g}}{1+r_{g}^{2}}(2F_+^{f}-1)\right]}. $

      (45)

      As described in Sec. IV, the ST and DT yields are extracted by performing an unbinned maximum likelihood fit on the $ M_{\rm bc}^{\rm tag} $ distribution and the 2D distribution of $ M_{\rm bc}^{\rm tag} $ versus $ M_{\rm bc}^{\rm sig} $, respectively. The corresponding ST and DT efficiencies are obtained with the similar fit processes on the signal MC sample, as described in Sec. II. The corresponding ST and DT yields and the ST and DT efficiencies are summarized in Table 3. According to Eq. (45), the values in Table 3 and the $ CP $-even fractions obtained in the amplitude analysis, the branching fractions of $ D^0\to \pi^+\pi^-\pi^+\pi^- $ and $ D^0\to \pi^+\pi^-\pi^0\pi^0 $(non-η) are calculated and summarized in Table 16. According to the FFs obtained in the amplitude analysis, the branching fractions of the intermediate processes are also determined. The obtained results are summarized in Table 17.

      1234567Total
      FF($ D^{0}\to a_{1}(1260)^{+}\pi^{-} $)0.6/1.10.0/0.10.3/0.13.8/0.53.0/2.10.0/0.10.2/1.10.7 $ \pm $ 4.8/1.1 $ \pm $ 2.4
      FF($ D^{0}\to a_{1}(1260)^{-}\pi^{+} $)0.1/0.10.1/0.10.2/0.20.7/1.51.4/1.20.1/0.10.7/0.70.2 $ \pm $ 1.7/0.2 $ \pm $ 2.0
      FF($ D^{0}\to a_{1}(1260)^{0}\pi^{0} $)–/0.4–/0.1–/0.2–/1.8–/1.7–/0.0–/0.9–/0.5 $ \pm $ 2.6
      FF($ D^{0}\to a_{1}(1420)^{+}\pi^{-} $)0.1/0.20.0/0.00.1/0.00.4/0.70.4/0.50.0/0.10.6/0.50.1 $ \pm $ 0.8/0.2 $ \pm $ 1.0
      FF($ D^{0}\to a_{1}(1640)^{+}\pi^{-} $)0.7/0.80.2/0.20.2/0.10.6/0.91.1/1.10.2/0.21.2/1.30.8 $ \pm $ 1.7/0.8 $ \pm $ 1.9
      Continued on next page

      Table 15.  Systematic uncertainties in the FFs, resonance parameters, and $ CP $–even fractions in units of statistical standard deviations. 1: Background estimation. 2: Detection efficiency. 3: Fit bias. 4: Resonance parameters. 5: Radii of Blatt-Weisskopf barrier factors. 6: Quantum correlation correction. 7: Extra amplitudes. For the items with "/", the quantities before and after "/" are for $ D^0\to\pi^{+}\pi^{-}\pi^{+}\pi^{-} $ and $ D^0\to\pi^{+}\pi^{-}\pi^{0}\pi^{0} $, respectively. In the "Total" column, the term before "$ \pm $" is the total experimental systematic uncertainty (1, 2, and 3), and the term after "$ \pm $" is the total model-dependent systematic uncertainty (4, 5, 6, and 7).

      Table 15-continued from previous
      1234567Total
      FF($ D^{0}\to a_{1}(1640)^{-}\pi^{+} $)0.1/0.10.1/0.10.0/0.00.4/0.20.6/0.60.2/0.20.9/1.00.1 $ \pm $ 1.2/0.1 $ \pm $ 1.2
      FF($ D^{0}\to a_{2}(1320)^{+}\pi^{-} $)0.1/0.20.2/0.20.3/0.30.5/0.40.2/0.10.1/0.10.3/0.20.4 $ \pm $ 0.6/0.4 $ \pm $ 0.5
      FF($ D^{0}\to a_{2}(1320)^{-}\pi^{+} $)0.1/0.00.0/0.00.1/0.10.7/0.30.5/0.60.1/0.10.4/0.50.1 $ \pm $ 1.0/0.1 $ \pm $ 0.8
      FF($ D^{0}\to h_{1}(1170)^{0}\pi^{0} $)–/0.7–/0.1–/0.2–/1.1–/0.7–/0.0–/1.0–/0.7 $ \pm $ 1.6
      FF($ D^{0}\to \pi(1300)^{+}\pi^{-} $)0.6/0.90.1/0.10.0/0.00.9/1.01.1/1.20.2/0.10.6/0.20.6 $ \pm $ 1.6/0.9 $ \pm $ 1.6
      FF($ D^{0}\to \pi(1300)^{-}\pi^{+} $)0.1/0.40.2/0.20.0/0.11.1/1.41.2/1.50.1/0.10.5/0.10.2 $ \pm $ 1.7/0.5 $ \pm $ 2.1
      FF($ D^{0}\to \pi(1300)^{0}\pi^{0} $)–/0.5–/0.2–/0.0–/0.9–/0.6–/0.0–/0.2–/0.5 $ \pm $ 1.1
      FF($ D^{0}\to \pi_{2}(1670)^{0}\pi^{0} $)–/0.8–/0.1–/0.1–/0.3–/0.4–/0.1–/1.4–/0.8 $ \pm $ 1.5
      FF($ D^{0}\to \rho(770)^{0}\rho(770)^{0} $)0.3/–0.0/–0.0/–1.5/–0.2/–0.1/–0.6/–0.3 $ \pm $ 1.6/–
      FF($ D^{0}\to \rho(770)^{0}\rho(770)^{0}[S] $)0.1/–0.0/–0.1/–0.7/–0.3/–0.1/–0.3/–0.1 $ \pm $ 0.8/–
      FF($ D^{0}\to \rho(770)^{0}\rho(770)^{0}[P] $)0.0/–0.0/–0.4/–0.5/–0.5/–0.3/–0.1/–0.4 $ \pm $ 0.8/–
      FF($ D^{0}\to \rho(770)^{0}\rho(770)^{0}[D] $)0.3/–0.0/–0.2/–0.9/–0.3/–0.1/–0.5/–0.4 $ \pm $ 1.1/–
      FF($ D^{0}\to \rho(770)^{0}\rho(1450)^{0} $)0.0/–0.0/–0.0/–0.7/–0.7/–0.1/–0.8/–0.0 $ \pm $ 1.3/–
      FF($ D^{0}\to \rho(770)^{0}\rho(1450)^{0}[P] $)0.1/–0.2/–0.2/–1.3/–0.6/–0.3/–0.1/–0.3 $ \pm $ 1.5/–
      FF($ D^{0}\to \rho(770)^{0}\rho(1450)^{0}[D] $)0.1/–0.1/–0.1/–0.2/–0.9/–0.1/–0.8/–0.2 $ \pm $ 1.2/–
      FF($ D^{0}\to \rho(770)^{+}\rho(770)^{-} $)–/0.9–/0.0–/0.3–/1.3–/0.5–/0.1–/1.0–/0.9 $ \pm $ 1.7
      FF($ D^{0}\to \rho(770)^{+}\rho(770)^{-}[S] $)–/0.6–/0.1–/0.1–/1.1–/0.6–/0.0–/1.0–/0.6 $ \pm $ 1.6
      FF($ D^{0}\to \rho(770)^{+}\rho(770)^{-}[P] $)–/1.0–/0.2–/0.1–/0.5–/0.4–/0.1–/0.8–/1.0 $ \pm $ 1.0
      FF($ D^{0}\to \rho(770)^{+}\rho(770)^{-}[D] $)–/0.1–/0.1–/0.2–/0.4–/0.3–/0.1–/0.6–/0.2 $ \pm $ 0.8
      FF($ D^{0}\to \rho(770)^{+}\rho(1450)^{-}[D] $)–/0.7–/0.1–/0.4–/1.9–/0.8–/0.1–/0.7–/0.8 $ \pm $ 2.2
      FF($ D^{0}\to \rho(770)^{0}(\pi\pi)_{S} $)0.5/0.30.0/0.00.2/0.22.7/1.80.1/0.10.0/0.01.0/0.70.5 $ \pm $ 2.9/0.4 $ \pm $ 1.9
      FF($ D^{0}\to (\pi^{+}\pi^{-})_{S}(\pi\pi)_{S} $)0.1/0.60.0/0.00.1/0.22.0/1.30.6/0.80.1/0.00.3/0.50.1 $ \pm $ 2.1/0.6 $ \pm $ 1.6
      FF($ D^{0}\to f_{2}(1270)^{0}(\pi\pi)_{S} $)0.0/0.10.0/0.00.1/0.12.6/2.91.1/1.00.0/0.01.4/1.60.1 $ \pm $ 3.2/0.1 $ \pm $ 3.5
      FF($ D^{0}\to \omega(782)\pi^{0} $)–/0.1–/0.0–/0.1–/0.3–/0.1–/0.0–/0.4–/0.1 $ \pm $ 0.5
      FF($ D^{0}\to \phi(1020)\pi^{0} $)–/0.2–/0.1–/0.4–/0.3–/0.2–/0.1–/0.5–/0.5 $ \pm $ 0.6
      FF($ a_{1}(1260)^{\pm}\to \rho(770)\pi[S] $)0.8/0.80.1/0.10.1/0.01.2/1.20.4/0.40.1/0.10.6/0.60.8 $ \pm $ 1.4/0.8 $ \pm $ 1.4
      FF($ a_{1}(1260)^{\pm}\to \rho(770)\pi[D] $)0.1/0.10.0/0.00.2/0.22.2/2.10.5/0.40.0/0.00.1/0.10.2 $ \pm $ 2.3/0.2 $ \pm $ 2.1
      FF($ a_{1}(1260)^{\pm}\to f_{2}(1270)\pi[P] $)0.7/0.70.0/0.00.1/0.10.4/0.31.2/1.20.1/0.10.3/0.30.7 $ \pm $ 1.3/0.7 $ \pm $ 1.3
      FF($ a_{1}(1260)^{\pm}\to (\pi^{+}\pi^{-})_{S}\pi[P] $)0.8/0.90.1/0.00.2/0.21.0/1.60.3/0.10.2/0.10.9/0.80.8 $ \pm $ 1.4/0.9 $ \pm $ 1.8
      FF($ a_{1}(1260)^{0}\to \rho(770)\pi[S] $)–/0.7–/0.1–/0.0–/1.4–/0.3–/0.1–/0.6–/0.7 $ \pm $ 1.6
      FF($ a_{1}(1260)^{0}\to \rho(770)\pi[D] $)–/0.1–/0.0–/0.2–/2.2–/0.4–/0.0–/0.1–/0.2 $ \pm $ 2.2
      FF($ a_{1}(1260)^{0}\to f_{2}(1270)\pi[P] $)–/0.7–/0.0–/0.1–/0.4–/1.4–/0.1–/0.3–/0.7 $ \pm $ 1.5
      FF($ a_{1}(1260)^{0}\to (\pi^{+}\pi^{-})_{S}\pi[P] $)–/0.8–/0.0–/0.3–/1.6–/0.1–/0.1–/0.8–/0.9 $ \pm $ 1.8
      FF($ \pi(1300)^{\pm}\to \rho(770)\pi $)0.0/0.20.1/0.10.3/0.32.0/0.82.2/3.30.1/0.10.4/0.30.3 $ \pm $ 3.0/0.4 $ \pm $ 3.4
      FF($ \pi(1300)^{\pm}\to (\pi^{+}\pi^{-})_{S}\pi $)0.1/0.10.1/0.10.3/0.41.8/1.22.3/2.90.1/0.10.3/0.20.3 $ \pm $ 2.9/0.4 $ \pm $ 3.1
      FF($ \pi(1300)^{0}\to \rho(770)\pi $)–/0.2–/0.1–/0.2–/1.0–/3.3–/0.1–/0.4–/0.3 $ \pm $ 3.5
      FF($ \pi(1300)^{0}\to (\pi^{+}\pi^{-})_{S}\pi $)–/0.0–/0.1–/0.4–/1.9–/2.6–/0.1–/0.1–/0.4 $ \pm $ 3.2
      $ M_{a_{1}(1260)} $0.50.10.04.01.20.11.60.5 $ \pm $ 4.5
      $ \Gamma_{a_{1}(1260)} $1.70.20.03.22.30.41.71.7 $ \pm $ 4.3
      $ M_{\pi(1300)} $0.80.10.01.11.10.01.00.8 $ \pm $ 1.8
      $ \Gamma_{\pi(1300)} $0.70.10.02.90.40.20.50.7 $ \pm $ 3.0
      $ F^{+} $0.7/1.00.0/0.20.1/0.10.7/0.80.3/0.80.1/0.11.0/0.60.7 $ \pm $ 1.3/1.0 $ \pm $ 1.3
      This workPDG
      $ D^0 \to \pi^+\pi^-\pi^+\pi^- $$ (0.688 \pm 0.010 \pm 0.010) $%$ (0.756 \pm 0.020) $%
      $ D^0 \to \pi^+\pi^-\pi^0\pi^0 $(non-η)$ (0.951 \pm 0.025 \pm 0.021) $%$ (1.005 \pm 0.090) $%

      Table 16.  Branching fractions of $ D^0\to \pi^+\pi^-\pi^+\pi^- $ and $ D^0\to\pi^+\pi^-\pi^0\pi^0 $(non-η) from this work compared to the PDG [25] values, where the first uncertainties are statistical and the second are systematic.

      ComponentBranching fraction (%)
      $ \pi^+\pi^-\pi^+\pi^- $$ \pi^+\pi^-\pi^0\pi^0 $
      $ D^0\to a_{1}(1260)^{+}\pi^- $$ 0.566 \pm 0.024 \pm 0.008 \pm 0.110 $$ 0.546 \pm 0.027 \pm 0.011 \pm 0.069 $
      $ D^0\to a_{1}(1260)^{-}\pi^+ $$ 0.071 \pm 0.010 \pm 0.001 \pm 0.017 $$ 0.068 \pm 0.011 \pm 0.001 \pm 0.021 $
      $ D^0\to a_{1}(1260)^{0}\pi^0 $$ 0.313 \pm 0.031 \pm 0.007 \pm 0.082 $
      $ D^0\to a_{1}(1640)^{+}\pi^- $$ 0.012 \pm 0.003 \pm 0.000 \pm 0.006 $$ 0.010 \pm 0.003 \pm 0.000 \pm 0.007 $
      $ D^0\to h_{1}(1170)^{0}\pi^0 $$ 0.012 \pm 0.006 \pm 0.000 \pm 0.010 $
      $ D^0\to \pi(1300)^{+}\pi^- $$ 0.222 \pm 0.018 \pm 0.003 \pm 0.031 $$ 0.148 \pm 0.014 \pm 0.003 \pm 0.025 $
      $ D^0\to \pi(1300)^{-}\pi^+ $$ 0.162 \pm 0.016 \pm 0.002 \pm 0.028 $$ 0.108 \pm 0.011 \pm 0.002 \pm 0.021 $
      $ D^0\to \pi(1300)^{0}\pi^0 $$ 0.221 \pm 0.027 \pm 0.005 \pm 0.033 $
      $ D^0\to \pi_{2}(1670)^{0}\pi^0 $$ 0.010 \pm 0.002 \pm 0.000 \pm 0.004 $
      $ D^0\to \rho(770)^{0}\rho(770)^{0} $$ 0.193 \pm 0.013 \pm 0.003 \pm 0.022 $
      $ D^0\to \rho(770)^{0}\rho(770)^{0}[S] $$ 0.012 \pm 0.004 \pm 0.000 \pm 0.003 $
      $ D^0\to \rho(770)^{0}\rho(770)^{0}[P] $$ 0.067 \pm 0.007 \pm 0.001 \pm 0.006 $
      $ D^0\to \rho(770)^{0}\rho(770)^{0}[D] $$ 0.159 \pm 0.015 \pm 0.002 \pm 0.017 $
      $ D^0\to \rho(770)^{0}\rho(1450)^{0} $$ 0.017 \pm 0.006 \pm 0.000 \pm 0.008 $
      $ D^0\to \rho(770)^{0}\rho(1450)^{0}[P] $$ 0.007 \pm 0.003 \pm 0.000 \pm 0.003 $
      $ D^0\to \rho(770)^{0}\rho(1450)^{0}[D] $$ 0.010 \pm 0.006 \pm 0.000 \pm 0.008 $
      $ D^0\to \rho(770)^{+}\rho(770)^{-} $$ 0.864 \pm 0.040 \pm 0.018 \pm 0.075 $
      $ D^0\to \rho(770)^{+}\rho(770)^{-}[S] $$ 0.124 \pm 0.019 \pm 0.003 \pm 0.033 $
      $ D^0\to \rho(770)^{+}\rho(770)^{-}[P] $$ 0.186 \pm 0.013 \pm 0.004 \pm 0.019 $
      $ D^0\to \rho(770)^{+}\rho(770)^{-}[D] $$ 0.342 \pm 0.029 \pm 0.007 \pm 0.024 $
      $ D^0\to \rho(770)^{+}\rho(1450)^{-}[D] $$ 0.016 \pm 0.008 \pm 0.000 \pm 0.016 $
      $ D^0\to \rho(770)^{0}(\pi\pi)_{S} $$ 0.019 \pm 0.004 \pm 0.000 \pm 0.012 $$ 0.010 \pm 0.002 \pm 0.000 \pm 0.004 $
      $ D^0\to (\pi^+\pi^-)_{S}(\pi\pi)_{S} $$ 0.432 \pm 0.032 \pm 0.006 \pm 0.066 $$ 0.356 \pm 0.029 \pm 0.007 \pm 0.049 $
      $ D^0\to f_{2}(1270)^{0}(\pi\pi)_{S} $$ 0.012 \pm 0.003 \pm 0.000 \pm 0.008 $$ 0.010 \pm 0.002 \pm 0.000 \pm 0.008 $
      $ D^0\to \omega(782)\pi^0 $$ 0.009 \pm 0.004 \pm 0.000 \pm 0.002 $
      $ D^0\to \phi(1020)\pi^0 $$ 0.014 \pm 0.004 \pm 0.000 \pm 0.003 $

      Table 17.  Branching fractions of the intermediate processes with FFs$ > $1% in $ D^0\to \pi^+\pi^-\pi^+\pi^- $ and $ D^0\to \pi^+\pi^-\pi^0\pi^0 $. The first uncertainties are the statistical uncertainties, the second are the systematic uncertainties from the total branching fraction measurements, and the third are the systematic uncertainties of the FFs from the amplitude analysis.

      Benefiting from the DT method, most systematic uncertainties associated with the ST selection are canceled. The relative systematic uncertainties in the branching fraction measurement are summarized in Table 18, where the total uncertainties are the square roots of the quadrature sums of the individual ones. Details of the systematic uncertainties in the branching fraction measurements are described below.

      Source$ \mathcal{B}(D^0\to\pi^{+}\pi^{-}\pi^{+}\pi^{-}) $$ \mathcal{B}(D^0\to\pi^{+}\pi^{-}\pi^{0}\pi^{0}) $
      $ \pi^{\pm} $ tracking0.40.2
      $ \pi^{\pm} $ PID1.20.6
      $ \pi^0 $ reconstruction-1.3
      $ \Delta E $ cut in signal side0.10.1
      ST fit0.50.5
      DT fit0.61.4
      Amplitude model0.30.6
      Quantum correlation correction0.30.4
      MC statistics0.10.2
      Total1.52.2

      Table 18.  Relative systematic uncertainties (%) of the measured branching fractions.

      ● Tracking efficiency

      The tracking efficiency of $ \pi^{\pm} $ has been studied in Ref. [55]. The systematic uncertainties are assigned by re-weighting the data-MC difference in $ \pi^{\pm} $ tracking efficiencies according to the momentum distribution of $ \pi^{\pm} $ in the signal MC sample. They are assigned as 0.4% and 0.2% for $ D^0\to\pi^{+}\pi^{-}\pi^{+}\pi^{-} $ and $ D^0\to\pi^{+}\pi^{-}\pi^{0}\pi^{0} $, respectively.

      ● PID efficiency

      The PID efficiency of $ \pi^{\pm} $ has been studied in Ref. [55]. The systematic uncertainties are assigned by re-weighting the data-MC difference in $ \pi^{\pm} $ PID efficiencies according to the momentum distribution of $ \pi^{\pm} $ in the signal MC sample. They are assigned as 1.2% and 0.6% for $ D^0\to\pi^{+}\pi^{-}\pi^{+}\pi^{-} $ and $ D^0\to\pi^{+}\pi^{-}\pi^{0}\pi^{0} $, respectively.

      $ \pi^0 $ reconstruction efficiency

      The $ \pi^0 $ reconstruction efficiency is studied by the control samples $ \bar{D}^{0}\to K^+\pi^- $, $ K^+\pi^-\pi^0 $, $ K^+\pi^-\pi^+\pi^- $ versus $ D^{0}\to K^-\pi^+\pi^0 $, $ \pi^+\pi^-\pi^0 $. The corresponding systematic uncertainty for $ D^0\to \pi^+\pi^-\pi^0\pi^0 $ is assigned as 1.3%, by re-weighting the data-MC difference in $ \pi^0 $ reconstruction efficiencies according to the momentum distribution of $ \pi^0 $ in the signal MC sample.

      $ \Delta E $ requirement of the signal side

      To study the systematic uncertainty from the $ \Delta E $ requirement of the signal side, the $ \Delta E $ of the signal MC sample is smeared by a Gaussian function, while the mean and width are obtained by performing a fit to the $ \Delta E $ distribution of the data with the signal MC shape convolved with a Gaussian function. The resultant changes of the efficiencies with respect to the nominal values, 0.1% for both $ D^0\to\pi^{+}\pi^{-}\pi^{+}\pi^{-} $ and $ D^0\to\pi^{+}\pi^{-}\pi^{0}\pi^{0} $, are taken as the systematic uncertainties.

      ● ST fit

      The systematic uncertainty of the ST fit is studied by changing the signal and background shapes. The signal shape is changed by replacing the Gaussian resolution with a Crystal Ball function [57]. The background shape is changed with a floating cut-off parameter for the ARGUS function. The square root of the quadrature sum of the relative change of the ST yields, which gives 0.5%, is taken as the systematic uncertainty.

      ● DT fit

      As described in Sec. VII, the systematic uncertainties of the DT fit are obtained by varying the signal and background shapes and the yields of the peaking backgrounds in BKGVI in the fit. The square roots of the quadrature sums of the relative changes of the DT yields in the individual changes are 0.6% and 1.4% for $ D^0\to \pi^{+}\pi^{-}\pi^{+}\pi^{-} $ and $ D^0\to\pi^{+}\pi^{-}\pi^{0}\pi^{0} $, respectively, which are taken as the systematic uncertainties.

      ● Amplitude model

      The DT efficiencies are obtained by varying the parameters of the amplitude model within their uncertainties. The resultant standard deviations of the DT efficiencies, which are 0.3% and 0.6% for $ D^0\to \pi^{+}\pi^{-}\pi^{+}\pi^{-} $ and $ D^0\to\pi^{+}\pi^{-}\pi^{0}\pi^{0} $, respectively, are taken as the systematic uncertainties.

      ● Quantum correlation corrections

      The uncertainties associated with the input parameters of quantum correlation (r, R, δ, and $ F_{+} $) are 0.3% and 0.4% for $ D^0\to \pi^{+}\pi^{-}\pi^{+}\pi^{-} $ and $ D^0\to\pi^{+}\pi^{-}\pi^{0}\pi^{0} $, respectively, which are assigned according to uncertainty propagation.

      ● MC statistics

      The systematic uncertainties related with the limited statistics of the signal MC samples are 0.1% and 0.2% for $ D^0\to \pi^{+}\pi^{-}\pi^{+}\pi^{-} $ and $ D^0\to\pi^{+}\pi^{-}\pi^{0}\pi^{0} $, respectively.

    IX.   SUMMARY
    • Using 2.93 $ \rm fb^{-1} $ of the $ e^+e^- $ collision data taken at $ \sqrt{s}=3.773\; \rm{GeV} $ with the BESIII detector, a joint amplitude analysis of $ D^{0}\to\pi^{+}\pi^{-}\pi^{+}\pi^{-} $ and $ D^{0}\to\pi^{+}\pi^{-}\pi^{0}\pi^{0} $(non-η) is performed. Large interferences between the dominant amplitudes of $ D^{0}\to a_{1}(1260)\pi $, $ D^0\to\pi(1300)\pi $, $ D^0\to\rho(770)\rho(770) $, and $ D^0\to2(\pi\pi)_{S} $ are observed. Based on the amplitude model, the model dependent $ CP $-even fractions of $ D^{0}\to\pi^{+}\pi^{-}\pi^{+}\pi^{-} $ and $ D^{0}\to\pi^{+}\pi^{-}\pi^{0}\pi^{0} $(non-η) are determined to be $ (75.2\; \pm\; 1.1_{\rm stat.}\; \pm\; 1.5_{\rm syst.})\ $% and $(68.9\; \pm\; 1.5_{\rm stat.}\; \pm\; 2.4_{\rm syst.})$%, respectively, consistent with the previous measurements performed by the CLEO [6, 13, 52] and BESIII [53, 54] collaboration. The branching fractions of $ D^{0}\to\pi^{+}\pi^{-}\pi^{+}\pi^{-} $ and $ D^{0}\to \pi^{+}\pi^{-}\pi^{0}\pi^{0} $(non-η) are measured to be $(0.688\; \pm\; 0.010_{\rm stat.}\; \pm\; 0.010_{\rm syst.})$% and $(0.951\; \pm\; 0.025_{\rm stat.}\; \pm\; 0.021_{\rm syst.})$%, respectively, where the former is $ 3\sigma $ lower than the PDG [25] value. This $ 3\sigma $ deviation may be due to amplitude model variations, which affect the global reconstruction efficiency. These results provide essential information to the search for $ CP $ violation [9] and measurements of the binned strong phase parameter [6], which are important inputs in the γ measurement via the $ B^{-}\to D K^{-} $ decay.

    ACKNOWLEDGEMENTS
    • The BESIII Collaboration thanks the staff of BEPCII, the IHEP Computing Center, and the Supercomputing Center of the University of Science and Technology of China (USTC) for their immense support.

Reference (57)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return