BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
D0To2pip2pim Class Reference

#include <D0To2pip2pim.h>

Public Member Functions

 D0To2pip2pim ()
virtual ~D0To2pip2pim ()
void init ()
complex< double > Amp (vector< double > Pip1, vector< double > Pim1, vector< double > Pip2, vector< double > Pim2)

Detailed Description

Definition at line 9 of file D0To2pip2pim.h.

Constructor & Destructor Documentation

◆ D0To2pip2pim()

D0To2pip2pim::D0To2pip2pim ( )
inline

Definition at line 12 of file D0To2pip2pim.h.

12{}

◆ ~D0To2pip2pim()

D0To2pip2pim::~D0To2pip2pim ( )
virtual

Definition at line 11 of file D0To2pip2pim.cxx.

11{}

Member Function Documentation

◆ Amp()

complex< double > D0To2pip2pim::Amp ( vector< double > Pip1,
vector< double > Pim1,
vector< double > Pip2,
vector< double > Pim2 )

Definition at line 1658 of file D0To2pip2pim.cxx.

1659 {
1660
1661 vector<double> Pip1Pim1;
1662 Pip1Pim1.clear();
1663 vector<double> Pip1Pim2;
1664 Pip1Pim2.clear();
1665 vector<double> Pip2Pim1;
1666 Pip2Pim1.clear();
1667 vector<double> Pip2Pim2;
1668 Pip2Pim2.clear();
1669
1670 Pip1Pim1 = sum_tensor( Pip1, Pim1 );
1671 Pip1Pim2 = sum_tensor( Pip1, Pim2 );
1672 Pip2Pim1 = sum_tensor( Pip2, Pim1 );
1673 Pip2Pim2 = sum_tensor( Pip2, Pim2 );
1674
1675 vector<double> Pip1Pip2Pim1;
1676 Pip1Pip2Pim1.clear();
1677 vector<double> Pip1Pip2Pim2;
1678 Pip1Pip2Pim2.clear();
1679 vector<double> Pim1Pim2Pip1;
1680 Pim1Pim2Pip1.clear();
1681 vector<double> Pim1Pim2Pip2;
1682 Pim1Pim2Pip2.clear();
1683
1684 Pip1Pip2Pim1 = sum_tensor( Pip1Pim1, Pip2 );
1685 Pip1Pip2Pim2 = sum_tensor( Pip1Pim2, Pip2 );
1686 Pim1Pim2Pip1 = sum_tensor( Pip1Pim1, Pim2 );
1687 Pim1Pim2Pip2 = sum_tensor( Pip2Pim1, Pim2 );
1688
1689 vector<double> D0;
1690 D0.clear();
1691 D0 = sum_tensor( Pip1Pip2Pim1, Pim2 );
1692
1693 double M2_Pip1Pim1 = contract_11_0( Pip1Pim1, Pip1Pim1 );
1694 double M2_Pip1Pim2 = contract_11_0( Pip1Pim2, Pip1Pim2 );
1695 double M2_Pip2Pim1 = contract_11_0( Pip2Pim1, Pip2Pim1 );
1696 double M2_Pip2Pim2 = contract_11_0( Pip2Pim2, Pip2Pim2 );
1697
1698 double M2_Pip1Pip2Pim1 = contract_11_0( Pip1Pip2Pim1, Pip1Pip2Pim1 );
1699 double M2_Pip1Pip2Pim2 = contract_11_0( Pip1Pip2Pim2, Pip1Pip2Pim2 );
1700 double M2_Pim1Pim2Pip1 = contract_11_0( Pim1Pim2Pip1, Pim1Pim2Pip1 );
1701 double M2_Pim1Pim2Pip2 = contract_11_0( Pim1Pim2Pip2, Pim1Pim2Pip2 );
1702 double M2_D0 = contract_11_0( D0, D0 );
1703
1704 complex<double> GS_rho770_11 =
1705 GS( M2_Pip1Pim1, m0_rho770, w0_rho770, m2_Pi, m2_Pi, rRes, 1 );
1706 complex<double> GS_rho770_12 =
1707 GS( M2_Pip1Pim2, m0_rho770, w0_rho770, m2_Pi, m2_Pi, rRes, 1 );
1708 complex<double> GS_rho770_21 =
1709 GS( M2_Pip2Pim1, m0_rho770, w0_rho770, m2_Pi, m2_Pi, rRes, 1 );
1710 complex<double> GS_rho770_22 =
1711 GS( M2_Pip2Pim2, m0_rho770, w0_rho770, m2_Pi, m2_Pi, rRes, 1 );
1712
1713 complex<double> GS_rho1450_11 =
1714 GS( M2_Pip1Pim1, m0_rho1450, w0_rho1450, m2_Pi, m2_Pi, rRes, 1 );
1715 complex<double> GS_rho1450_12 =
1716 GS( M2_Pip1Pim2, m0_rho1450, w0_rho1450, m2_Pi, m2_Pi, rRes, 1 );
1717 complex<double> GS_rho1450_21 =
1718 GS( M2_Pip2Pim1, m0_rho1450, w0_rho1450, m2_Pi, m2_Pi, rRes, 1 );
1719 complex<double> GS_rho1450_22 =
1720 GS( M2_Pip2Pim2, m0_rho1450, w0_rho1450, m2_Pi, m2_Pi, rRes, 1 );
1721
1722 complex<double> RBW_f21270_11 =
1723 RBW( M2_Pip1Pim1, m0_f21270, w0_f21270, m2_Pi, m2_Pi, rRes, 2 );
1724 complex<double> RBW_f21270_12 =
1725 RBW( M2_Pip1Pim2, m0_f21270, w0_f21270, m2_Pi, m2_Pi, rRes, 2 );
1726 complex<double> RBW_f21270_21 =
1727 RBW( M2_Pip2Pim1, m0_f21270, w0_f21270, m2_Pi, m2_Pi, rRes, 2 );
1728 complex<double> RBW_f21270_22 =
1729 RBW( M2_Pip2Pim2, m0_f21270, w0_f21270, m2_Pi, m2_Pi, rRes, 2 );
1730
1731 complex<double> PiPiS_11_0 = Fvector( M2_Pip1Pim1, s0_prod, 0 );
1732 complex<double> PiPiS_12_0 = Fvector( M2_Pip1Pim2, s0_prod, 0 );
1733 complex<double> PiPiS_21_0 = Fvector( M2_Pip2Pim1, s0_prod, 0 );
1734 complex<double> PiPiS_22_0 = Fvector( M2_Pip2Pim2, s0_prod, 0 );
1735
1736 complex<double> PiPiS_11_1 = Fvector( M2_Pip1Pim1, s0_prod, 1 );
1737 complex<double> PiPiS_12_1 = Fvector( M2_Pip1Pim2, s0_prod, 1 );
1738 complex<double> PiPiS_21_1 = Fvector( M2_Pip2Pim1, s0_prod, 1 );
1739 complex<double> PiPiS_22_1 = Fvector( M2_Pip2Pim2, s0_prod, 1 );
1740
1741 complex<double> PiPiS_11_5 = Fvector( M2_Pip1Pim1, s0_prod, 5 );
1742 complex<double> PiPiS_12_5 = Fvector( M2_Pip1Pim2, s0_prod, 5 );
1743 complex<double> PiPiS_21_5 = Fvector( M2_Pip2Pim1, s0_prod, 5 );
1744 complex<double> PiPiS_22_5 = Fvector( M2_Pip2Pim2, s0_prod, 5 );
1745
1746 complex<double> PiPiS_11_6 = Fvector( M2_Pip1Pim1, s0_prod, 6 );
1747 complex<double> PiPiS_12_6 = Fvector( M2_Pip1Pim2, s0_prod, 6 );
1748 complex<double> PiPiS_21_6 = Fvector( M2_Pip2Pim1, s0_prod, 6 );
1749 complex<double> PiPiS_22_6 = Fvector( M2_Pip2Pim2, s0_prod, 6 );
1750
1751 complex<double> RBW_a11260p_1 = RBWa1260( M2_Pip1Pip2Pim1, m0_a11260, g1_a11260, g2_a11260 );
1752 complex<double> RBW_a11260p_2 = RBWa1260( M2_Pip1Pip2Pim2, m0_a11260, g1_a11260, g2_a11260 );
1753 complex<double> RBW_a11260m_1 = RBWa1260( M2_Pim1Pim2Pip1, m0_a11260, g1_a11260, g2_a11260 );
1754 complex<double> RBW_a11260m_2 = RBWa1260( M2_Pim1Pim2Pip2, m0_a11260, g1_a11260, g2_a11260 );
1755
1756 complex<double> RBW_a21320p_1 =
1757 RBW( M2_Pip1Pip2Pim1, m0_a21320, w0_a21320, -1.0, -1.0, -1, -1 );
1758 complex<double> RBW_a21320p_2 =
1759 RBW( M2_Pip1Pip2Pim2, m0_a21320, w0_a21320, -1.0, -1.0, -1, -1 );
1760 complex<double> RBW_a21320m_1 =
1761 RBW( M2_Pim1Pim2Pip1, m0_a21320, w0_a21320, -1.0, -1.0, -1, -1 );
1762 complex<double> RBW_a21320m_2 =
1763 RBW( M2_Pim1Pim2Pip2, m0_a21320, w0_a21320, -1.0, -1.0, -1, -1 );
1764
1765 complex<double> RBW_pi1300p_1 = RBWpi1300( M2_Pip1Pip2Pim1, m0_pi1300, w0_pi1300 );
1766 complex<double> RBW_pi1300p_2 = RBWpi1300( M2_Pip1Pip2Pim2, m0_pi1300, w0_pi1300 );
1767 complex<double> RBW_pi1300m_1 = RBWpi1300( M2_Pim1Pim2Pip1, m0_pi1300, w0_pi1300 );
1768 complex<double> RBW_pi1300m_2 = RBWpi1300( M2_Pim1Pim2Pip2, m0_pi1300, w0_pi1300 );
1769
1770 complex<double> RBW_a11420p_1 =
1771 RBW( M2_Pip1Pip2Pim1, m0_a11420, w0_a11420, -1.0, -1.0, -1, -1 );
1772 complex<double> RBW_a11420p_2 =
1773 RBW( M2_Pip1Pip2Pim2, m0_a11420, w0_a11420, -1.0, -1.0, -1, -1 );
1774 complex<double> RBW_a11420m_1 =
1775 RBW( M2_Pim1Pim2Pip1, m0_a11420, w0_a11420, -1.0, -1.0, -1, -1 );
1776 complex<double> RBW_a11420m_2 =
1777 RBW( M2_Pim1Pim2Pip2, m0_a11420, w0_a11420, -1.0, -1.0, -1, -1 );
1778
1779 // D->XP 1-rank Projection
1780 vector<double> Proj1_3p1;
1781 Proj1_3p1.clear();
1782 vector<double> Proj1_3p2;
1783 Proj1_3p2.clear();
1784 vector<double> Proj1_3m1;
1785 Proj1_3m1.clear();
1786 vector<double> Proj1_3m2;
1787 Proj1_3m2.clear();
1788
1789 Proj1_3p1 = ProjectionTensors( Pip1Pip2Pim1, 1 );
1790 Proj1_3p2 = ProjectionTensors( Pip1Pip2Pim2, 1 );
1791 Proj1_3m1 = ProjectionTensors( Pim1Pim2Pip1, 1 );
1792 Proj1_3m2 = ProjectionTensors( Pim1Pim2Pip2, 1 );
1793
1794 // D->XP 2-rank Projection
1795 vector<double> Proj2_3p1;
1796 Proj2_3p1.clear();
1797 vector<double> Proj2_3p2;
1798 Proj2_3p2.clear();
1799 vector<double> Proj2_3m1;
1800 Proj2_3m1.clear();
1801 vector<double> Proj2_3m2;
1802 Proj2_3m2.clear();
1803
1804 Proj2_3p1 = ProjectionTensors( Pip1Pip2Pim1, 2 );
1805 Proj2_3p2 = ProjectionTensors( Pip1Pip2Pim2, 2 );
1806 Proj2_3m1 = ProjectionTensors( Pim1Pim2Pip1, 2 );
1807 Proj2_3m2 = ProjectionTensors( Pim1Pim2Pip2, 2 );
1808
1809 // X->PP Orbital
1810 vector<double> T1_Pip1Pim1;
1811 T1_Pip1Pim1.clear();
1812 vector<double> T1_Pip1Pim2;
1813 T1_Pip1Pim2.clear();
1814 vector<double> T1_Pip2Pim1;
1815 T1_Pip2Pim1.clear();
1816 vector<double> T1_Pip2Pim2;
1817 T1_Pip2Pim2.clear();
1818
1819 T1_Pip1Pim1 = OrbitalTensors( Pip1Pim1, Pip1, Pim1, rRes, 1 );
1820 T1_Pip1Pim2 = OrbitalTensors( Pip1Pim2, Pip1, Pim2, rRes, 1 );
1821 T1_Pip2Pim1 = OrbitalTensors( Pip2Pim1, Pip2, Pim1, rRes, 1 );
1822 T1_Pip2Pim2 = OrbitalTensors( Pip2Pim2, Pip2, Pim2, rRes, 1 );
1823
1824 vector<double> T1_Pim1Pip1;
1825 T1_Pim1Pip1.clear();
1826 vector<double> T1_Pim1Pip2;
1827 T1_Pim1Pip2.clear();
1828 vector<double> T1_Pim2Pip1;
1829 T1_Pim2Pip1.clear();
1830 vector<double> T1_Pim2Pip2;
1831 T1_Pim2Pip2.clear();
1832
1833 T1_Pim1Pip1 = OrbitalTensors( Pip1Pim1, Pim1, Pip1, rRes, 1 );
1834 T1_Pim1Pip2 = OrbitalTensors( Pip2Pim1, Pim1, Pip2, rRes, 1 );
1835 T1_Pim2Pip1 = OrbitalTensors( Pip1Pim2, Pim2, Pip1, rRes, 1 );
1836 T1_Pim2Pip2 = OrbitalTensors( Pip2Pim2, Pim2, Pip2, rRes, 1 );
1837
1838 vector<double> T2_Pip1Pim1;
1839 T2_Pip1Pim1.clear();
1840 vector<double> T2_Pip1Pim2;
1841 T2_Pip1Pim2.clear();
1842 vector<double> T2_Pip2Pim1;
1843 T2_Pip2Pim1.clear();
1844 vector<double> T2_Pip2Pim2;
1845 T2_Pip2Pim1.clear();
1846
1847 T2_Pip1Pim1 = OrbitalTensors( Pip1Pim1, Pip1, Pim1, rRes, 2 );
1848 T2_Pip1Pim2 = OrbitalTensors( Pip1Pim2, Pip1, Pim2, rRes, 2 );
1849 T2_Pip2Pim1 = OrbitalTensors( Pip2Pim1, Pip2, Pim1, rRes, 2 );
1850 T2_Pip2Pim2 = OrbitalTensors( Pip2Pim2, Pip2, Pim2, rRes, 2 );
1851
1852 // X->YP Orbital
1853 vector<double> T1_Pip1Pim1Pip2;
1854 T1_Pip1Pim1Pip2.clear();
1855 vector<double> T1_Pip2Pim1Pip1;
1856 T1_Pip2Pim1Pip1.clear();
1857 vector<double> T1_Pip1Pim2Pip2;
1858 T1_Pip1Pim2Pip2.clear();
1859 vector<double> T1_Pip2Pim2Pip1;
1860 T1_Pip2Pim2Pip1.clear();
1861 vector<double> T1_Pip1Pim1Pim2;
1862 T1_Pip1Pim1Pim2.clear();
1863 vector<double> T1_Pip1Pim2Pim1;
1864 T1_Pip1Pim2Pim1.clear();
1865 vector<double> T1_Pip2Pim1Pim2;
1866 T1_Pip2Pim1Pim2.clear();
1867 vector<double> T1_Pip2Pim2Pim1;
1868 T1_Pip2Pim2Pim1.clear();
1869
1870 T1_Pip1Pim1Pip2 = OrbitalTensors( Pip1Pip2Pim1, Pip1Pim1, Pip2, rRes, 1 );
1871 T1_Pip2Pim1Pip1 = OrbitalTensors( Pip1Pip2Pim1, Pip2Pim1, Pip1, rRes, 1 );
1872 T1_Pip1Pim2Pip2 = OrbitalTensors( Pip1Pip2Pim2, Pip1Pim2, Pip2, rRes, 1 );
1873 T1_Pip2Pim2Pip1 = OrbitalTensors( Pip1Pip2Pim2, Pip2Pim2, Pip1, rRes, 1 );
1874 T1_Pip1Pim1Pim2 = OrbitalTensors( Pim1Pim2Pip1, Pip1Pim1, Pim2, rRes, 1 );
1875 T1_Pip2Pim1Pim2 = OrbitalTensors( Pim1Pim2Pip2, Pip2Pim1, Pim2, rRes, 1 );
1876 T1_Pip1Pim2Pim1 = OrbitalTensors( Pim1Pim2Pip1, Pip1Pim2, Pim1, rRes, 1 );
1877 T1_Pip2Pim2Pim1 = OrbitalTensors( Pim1Pim2Pip2, Pip2Pim2, Pim1, rRes, 1 );
1878
1879 vector<double> T2_Pip1Pim1Pip2;
1880 T2_Pip1Pim1Pip2.clear();
1881 vector<double> T2_Pip2Pim1Pip1;
1882 T2_Pip2Pim1Pip1.clear();
1883 vector<double> T2_Pip1Pim2Pip2;
1884 T2_Pip1Pim2Pip2.clear();
1885 vector<double> T2_Pip2Pim2Pip1;
1886 T2_Pip2Pim2Pip1.clear();
1887 vector<double> T2_Pip1Pim1Pim2;
1888 T2_Pip1Pim1Pim2.clear();
1889 vector<double> T2_Pip2Pim1Pim2;
1890 T2_Pip2Pim1Pim2.clear();
1891 vector<double> T2_Pip1Pim2Pim1;
1892 T2_Pip1Pim2Pim1.clear();
1893 vector<double> T2_Pip2Pim2Pim1;
1894 T2_Pip2Pim2Pim1.clear();
1895
1896 T2_Pip1Pim1Pip2 = OrbitalTensors( Pip1Pip2Pim1, Pip1Pim1, Pip2, rRes, 2 );
1897 T2_Pip2Pim1Pip1 = OrbitalTensors( Pip1Pip2Pim1, Pip2Pim1, Pip1, rRes, 2 );
1898 T2_Pip1Pim2Pip2 = OrbitalTensors( Pip1Pip2Pim2, Pip1Pim2, Pip2, rRes, 2 );
1899 T2_Pip2Pim2Pip1 = OrbitalTensors( Pip1Pip2Pim2, Pip2Pim2, Pip1, rRes, 2 );
1900 T2_Pip1Pim1Pim2 = OrbitalTensors( Pim1Pim2Pip1, Pip1Pim1, Pim2, rRes, 2 );
1901 T2_Pip2Pim1Pim2 = OrbitalTensors( Pim1Pim2Pip2, Pip2Pim1, Pim2, rRes, 2 );
1902 T2_Pip1Pim2Pim1 = OrbitalTensors( Pim1Pim2Pip1, Pip1Pim2, Pim1, rRes, 2 );
1903 T2_Pip2Pim2Pim1 = OrbitalTensors( Pim1Pim2Pip2, Pip2Pim2, Pim1, rRes, 2 );
1904
1905 // D->XX Orbital
1906 vector<double> T1_2z11;
1907 T1_2z11.clear();
1908 vector<double> T1_2z12;
1909 T1_2z12.clear();
1910 vector<double> T1_2z21;
1911 T1_2z21.clear();
1912 vector<double> T1_2z22;
1913 T1_2z22.clear();
1914
1915 T1_2z11 = OrbitalTensors( D0, Pip1Pim1, Pip2Pim2, rD, 1 );
1916 T1_2z12 = OrbitalTensors( D0, Pip2Pim2, Pip1Pim1, rD, 1 );
1917 T1_2z21 = OrbitalTensors( D0, Pip1Pim2, Pip2Pim1, rD, 1 );
1918 T1_2z22 = OrbitalTensors( D0, Pip2Pim1, Pip1Pim2, rD, 1 );
1919
1920 vector<double> T2_2z11;
1921 T2_2z11.clear();
1922 vector<double> T2_2z12;
1923 T2_2z12.clear();
1924 vector<double> T2_2z21;
1925 T2_2z21.clear();
1926 vector<double> T2_2z22;
1927 T2_2z22.clear();
1928
1929 T2_2z11 = OrbitalTensors( D0, Pip1Pim1, Pip2Pim2, rD, 2 );
1930 T2_2z12 = OrbitalTensors( D0, Pip2Pim2, Pip1Pim1, rD, 2 );
1931 T2_2z21 = OrbitalTensors( D0, Pip1Pim2, Pip2Pim1, rD, 2 );
1932 T2_2z22 = OrbitalTensors( D0, Pip2Pim1, Pip1Pim2, rD, 2 );
1933
1934 // D->XP Orbital
1935 vector<double> T1_3p1;
1936 T1_3p1.clear();
1937 vector<double> T1_3p2;
1938 T1_3p2.clear();
1939 vector<double> T1_3m1;
1940 T1_3m1.clear();
1941 vector<double> T1_3m2;
1942 T1_3m2.clear();
1943
1944 T1_3p1 = OrbitalTensors( D0, Pip1Pip2Pim1, Pim2, rD, 1 );
1945 T1_3p2 = OrbitalTensors( D0, Pip1Pip2Pim2, Pim1, rD, 1 );
1946 T1_3m1 = OrbitalTensors( D0, Pim1Pim2Pip1, Pip2, rD, 1 );
1947 T1_3m2 = OrbitalTensors( D0, Pim1Pim2Pip2, Pip1, rD, 1 );
1948
1949 vector<double> T2_3p1;
1950 T2_3p1.clear();
1951 vector<double> T2_3p2;
1952 T2_3p2.clear();
1953 vector<double> T2_3m1;
1954 T2_3m1.clear();
1955 vector<double> T2_3m2;
1956 T2_3m2.clear();
1957
1958 T2_3p1 = OrbitalTensors( D0, Pip1Pip2Pim1, Pim2, rD, 2 );
1959 T2_3p2 = OrbitalTensors( D0, Pip1Pip2Pim2, Pim1, rD, 2 );
1960 T2_3m1 = OrbitalTensors( D0, Pim1Pim2Pip1, Pip2, rD, 2 );
1961 T2_3m2 = OrbitalTensors( D0, Pim1Pim2Pip2, Pip1, rD, 2 );
1962
1963 complex<double> amplitude( 0, 0 );
1964 vector<complex<double>> g_fitpara;
1965 g_fitpara.clear();
1966 g_fitpara.push_back( complex<double>( 100.0 * cos( 0.0 ), 100.0 * sin( 0.0 ) ) );
1967 g_fitpara.push_back(
1968 complex<double>( 7.95507 * cos( -0.0687407 ), 7.95507 * sin( -0.0687407 ) ) );
1969 g_fitpara.push_back(
1970 complex<double>( 37.5559 * cos( -1.74946 ), 37.5559 * sin( -1.74946 ) ) );
1971 g_fitpara.push_back( complex<double>( 61.2172 * cos( 2.98079 ), 61.2172 * sin( 2.98079 ) ) );
1972 g_fitpara.push_back( complex<double>( 187.79 * cos( 2.64471 ), 187.79 * sin( 2.64471 ) ) );
1973 g_fitpara.push_back(
1974 complex<double>( 385.474 * cos( -0.137107 ), 385.474 * sin( -0.137107 ) ) );
1975 g_fitpara.push_back(
1976 complex<double>( 0.330788 * cos( 0.268133 ), 0.330788 * sin( 0.268133 ) ) );
1977 g_fitpara.push_back(
1978 complex<double>( 127.158 * cos( -2.47773 ), 127.158 * sin( -2.47773 ) ) );
1979 g_fitpara.push_back( complex<double>( 339.914 * cos( 2.22856 ), 339.914 * sin( 2.22856 ) ) );
1980 g_fitpara.push_back(
1981 complex<double>( 0.320888 * cos( -2.6194 ), 0.320888 * sin( -2.6194 ) ) );
1982 g_fitpara.push_back(
1983 complex<double>( 0.366283 * cos( -0.26867 ), 0.366283 * sin( -0.26867 ) ) );
1984 g_fitpara.push_back(
1985 complex<double>( 86.0865 * cos( -2.49649 ), 86.0865 * sin( -2.49649 ) ) );
1986 g_fitpara.push_back( complex<double>( 6.1541 * cos( -1.18299 ), 6.1541 * sin( -1.18299 ) ) );
1987 g_fitpara.push_back(
1988 complex<double>( 56.6067 * cos( 0.142977 ), 56.6067 * sin( 0.142977 ) ) );
1989 g_fitpara.push_back(
1990 complex<double>( 92.3073 * cos( -2.15881 ), 92.3073 * sin( -2.15881 ) ) );
1991 g_fitpara.push_back(
1992 complex<double>( 10.5885 * cos( -3.03166 ), 10.5885 * sin( -3.03166 ) ) );
1993 g_fitpara.push_back( complex<double>( 8.36765 * cos( 1.8417 ), 8.36765 * sin( 1.8417 ) ) );
1994 g_fitpara.push_back(
1995 complex<double>( 6.56437 * cos( -2.93087 ), 6.56437 * sin( -2.93087 ) ) );
1996 g_fitpara.push_back( complex<double>( 15.7197 * cos( 0.96925 ), 15.7197 * sin( 0.96925 ) ) );
1997 g_fitpara.push_back(
1998 complex<double>( 21.4195 * cos( -1.23701 ), 21.4195 * sin( -1.23701 ) ) );
1999 g_fitpara.push_back(
2000 complex<double>( 56.8867 * cos( -0.385837 ), 56.8867 * sin( -0.385837 ) ) );
2001 g_fitpara.push_back( complex<double>( 231.626 * cos( 2.14842 ), 231.626 * sin( 2.14842 ) ) );
2002 g_fitpara.push_back(
2003 complex<double>( 2938.45 * cos( -0.693491 ), 2938.45 * sin( -0.693491 ) ) );
2004 g_fitpara.push_back( complex<double>( 7252.7 * cos( 2.23659 ), 7252.7 * sin( 2.23659 ) ) );
2005 g_fitpara.push_back(
2006 complex<double>( 5165.87 * cos( 0.913557 ), 5165.87 * sin( 0.913557 ) ) );
2007 g_fitpara.push_back(
2008 complex<double>( 11508.6 * cos( -1.07187 ), 11508.6 * sin( -1.07187 ) ) );
2009 g_fitpara.push_back( complex<double>( 2461.86 * cos( 1.8709 ), 2461.86 * sin( 1.8709 ) ) );
2010 g_fitpara.push_back( complex<double>( 8757.75 * cos( 2.40756 ), 8757.75 * sin( 2.40756 ) ) );
2011 g_fitpara.push_back( complex<double>( 19.7413 * cos( -1.0753 ), 19.7413 * sin( -1.0753 ) ) );
2012 g_fitpara.push_back( complex<double>( 66.3826 * cos( 2.34666 ), 66.3826 * sin( 2.34666 ) ) );
2013
2014 // D0 -> a1(1260)+ {rho(770)0 pi+[S]} pi-
2015 double SF_Ap_S_VP_1 = contract_11_0( contract_21_1( Proj1_3p1, T1_Pip2Pim1 ), T1_3p1 );
2016 double SF_Ap_S_VP_2 = contract_11_0( contract_21_1( Proj1_3p1, T1_Pip1Pim1 ), T1_3p1 );
2017 double SF_Ap_S_VP_3 = contract_11_0( contract_21_1( Proj1_3p2, T1_Pip2Pim2 ), T1_3p2 );
2018 double SF_Ap_S_VP_4 = contract_11_0( contract_21_1( Proj1_3p2, T1_Pip1Pim2 ), T1_3p2 );
2019 amplitude += g_fitpara[0] * ( SF_Ap_S_VP_1 * RBW_a11260p_1 * GS_rho770_21 +
2020 SF_Ap_S_VP_2 * RBW_a11260p_1 * GS_rho770_11 +
2021 SF_Ap_S_VP_3 * RBW_a11260p_2 * GS_rho770_22 +
2022 SF_Ap_S_VP_4 * RBW_a11260p_2 * GS_rho770_12 );
2023
2024 // D0 -> a1(1260)+ {rho(770)0 pi+[D]} pi-
2025 double SF_Ap_D_VP_1 = contract_11_0( contract_21_1( T2_Pip2Pim1Pip1, T1_Pip2Pim1 ), T1_3p1 );
2026 double SF_Ap_D_VP_2 = contract_11_0( contract_21_1( T2_Pip1Pim1Pip2, T1_Pip1Pim1 ), T1_3p1 );
2027 double SF_Ap_D_VP_3 = contract_11_0( contract_21_1( T2_Pip2Pim2Pip1, T1_Pip2Pim2 ), T1_3p2 );
2028 double SF_Ap_D_VP_4 = contract_11_0( contract_21_1( T2_Pip1Pim2Pip2, T1_Pip1Pim2 ), T1_3p2 );
2029
2030 amplitude += g_fitpara[1] * ( SF_Ap_D_VP_1 * RBW_a11260p_1 * GS_rho770_21 +
2031 SF_Ap_D_VP_2 * RBW_a11260p_1 * GS_rho770_11 +
2032 SF_Ap_D_VP_3 * RBW_a11260p_2 * GS_rho770_22 +
2033 SF_Ap_D_VP_4 * RBW_a11260p_2 * GS_rho770_12 );
2034
2035 // D0 -> a1(1260)+ {f2(1270) pi+ [P]} pi-
2036 double SF_Ap_P_TP_1 = contract_11_0(
2037 contract_21_1( contract_42_2( Proj2_3p1, T2_Pip2Pim1 ), T1_Pip2Pim1Pip1 ), T1_3p1 );
2038 double SF_Ap_P_TP_2 = contract_11_0(
2039 contract_21_1( contract_42_2( Proj2_3p1, T2_Pip1Pim1 ), T1_Pip1Pim1Pip2 ), T1_3p1 );
2040 double SF_Ap_P_TP_3 = contract_11_0(
2041 contract_21_1( contract_42_2( Proj2_3p2, T2_Pip2Pim2 ), T1_Pip2Pim2Pip1 ), T1_3p2 );
2042 double SF_Ap_P_TP_4 = contract_11_0(
2043 contract_21_1( contract_42_2( Proj2_3p2, T2_Pip1Pim2 ), T1_Pip1Pim2Pip2 ), T1_3p2 );
2044
2045 amplitude += g_fitpara[2] * ( SF_Ap_P_TP_1 * RBW_a11260p_1 * RBW_f21270_21 +
2046 SF_Ap_P_TP_2 * RBW_a11260p_1 * RBW_f21270_11 +
2047 SF_Ap_P_TP_3 * RBW_a11260p_2 * RBW_f21270_22 +
2048 SF_Ap_P_TP_4 * RBW_a11260p_2 * RBW_f21270_12 );
2049
2050 // D0 -> a1(1260)+ {f0 pi+ [P]} pi-
2051 double SF_Ap_P_SP_1 = contract_11_0( T1_3p1, T1_Pip2Pim1Pip1 );
2052 double SF_Ap_P_SP_2 = contract_11_0( T1_3p1, T1_Pip1Pim1Pip2 );
2053 double SF_Ap_P_SP_3 = contract_11_0( T1_3p2, T1_Pip2Pim2Pip1 );
2054 double SF_Ap_P_SP_4 = contract_11_0( T1_3p2, T1_Pip1Pim2Pip2 );
2055
2056 amplitude += g_fitpara[3] * ( SF_Ap_P_SP_1 * RBW_a11260p_1 * PiPiS_21_0 +
2057 SF_Ap_P_SP_2 * RBW_a11260p_1 * PiPiS_11_0 +
2058 SF_Ap_P_SP_3 * RBW_a11260p_2 * PiPiS_22_0 +
2059 SF_Ap_P_SP_4 * RBW_a11260p_2 * PiPiS_12_0 );
2060 amplitude += g_fitpara[4] * ( SF_Ap_P_SP_1 * RBW_a11260p_1 * PiPiS_21_1 +
2061 SF_Ap_P_SP_2 * RBW_a11260p_1 * PiPiS_11_1 +
2062 SF_Ap_P_SP_3 * RBW_a11260p_2 * PiPiS_22_1 +
2063 SF_Ap_P_SP_4 * RBW_a11260p_2 * PiPiS_12_1 );
2064 amplitude += g_fitpara[5] * ( SF_Ap_P_SP_1 * RBW_a11260p_1 * PiPiS_21_5 +
2065 SF_Ap_P_SP_2 * RBW_a11260p_1 * PiPiS_11_5 +
2066 SF_Ap_P_SP_3 * RBW_a11260p_2 * PiPiS_22_5 +
2067 SF_Ap_P_SP_4 * RBW_a11260p_2 * PiPiS_12_5 );
2068
2069 // D0->a1(1260)- {rho(770)0 pi- [S]} pi+
2070 double SF_Am_S_VP_1 = contract_11_0( contract_21_1( Proj1_3m1, T1_Pim2Pip1 ), T1_3m1 );
2071 double SF_Am_S_VP_2 = contract_11_0( contract_21_1( Proj1_3m1, T1_Pim1Pip1 ), T1_3m1 );
2072 double SF_Am_S_VP_3 = contract_11_0( contract_21_1( Proj1_3m2, T1_Pim2Pip2 ), T1_3m2 );
2073 double SF_Am_S_VP_4 = contract_11_0( contract_21_1( Proj1_3m2, T1_Pim1Pip2 ), T1_3m2 );
2074
2075 amplitude += g_fitpara[0] * g_fitpara[6] *
2076 ( SF_Am_S_VP_1 * RBW_a11260m_1 * GS_rho770_12 +
2077 SF_Am_S_VP_2 * RBW_a11260m_1 * GS_rho770_11 +
2078 SF_Am_S_VP_3 * RBW_a11260m_2 * GS_rho770_22 +
2079 SF_Am_S_VP_4 * RBW_a11260m_2 * GS_rho770_21 );
2080
2081 // D0 -> a1(1260)- {rho(770)0 pi-[D]} pi+
2082 double SF_Am_D_VP_1 = contract_11_0( contract_21_1( T2_Pip1Pim2Pim1, T1_Pim2Pip1 ), T1_3m1 );
2083 double SF_Am_D_VP_2 = contract_11_0( contract_21_1( T2_Pip1Pim1Pim2, T1_Pim1Pip1 ), T1_3m1 );
2084 double SF_Am_D_VP_3 = contract_11_0( contract_21_1( T2_Pip2Pim2Pim1, T1_Pim2Pip2 ), T1_3m2 );
2085 double SF_Am_D_VP_4 = contract_11_0( contract_21_1( T2_Pip2Pim1Pim2, T1_Pim1Pip2 ), T1_3m2 );
2086
2087 amplitude += g_fitpara[1] * g_fitpara[6] *
2088 ( SF_Am_D_VP_1 * RBW_a11260m_1 * GS_rho770_12 +
2089 SF_Am_D_VP_2 * RBW_a11260m_1 * GS_rho770_11 +
2090 SF_Am_D_VP_3 * RBW_a11260m_2 * GS_rho770_22 +
2091 SF_Am_D_VP_4 * RBW_a11260m_2 * GS_rho770_21 );
2092
2093 // D0 -> a1(1260)- {f2(1270) pi- [P]} pi+
2094 double SF_Am_P_TP_1 = contract_11_0(
2095 contract_21_1( contract_42_2( Proj2_3m1, T2_Pip1Pim2 ), T1_Pip1Pim2Pim1 ), T1_3m1 );
2096 double SF_Am_P_TP_2 = contract_11_0(
2097 contract_21_1( contract_42_2( Proj2_3m1, T2_Pip1Pim1 ), T1_Pip1Pim1Pim2 ), T1_3m1 );
2098 double SF_Am_P_TP_3 = contract_11_0(
2099 contract_21_1( contract_42_2( Proj2_3m2, T2_Pip2Pim2 ), T1_Pip2Pim2Pim1 ), T1_3m2 );
2100 double SF_Am_P_TP_4 = contract_11_0(
2101 contract_21_1( contract_42_2( Proj2_3m2, T2_Pip2Pim1 ), T1_Pip2Pim1Pim2 ), T1_3m2 );
2102
2103 amplitude += g_fitpara[2] * g_fitpara[6] *
2104 ( SF_Am_P_TP_1 * RBW_a11260m_1 * RBW_f21270_12 +
2105 SF_Am_P_TP_2 * RBW_a11260m_1 * RBW_f21270_11 +
2106 SF_Am_P_TP_3 * RBW_a11260m_2 * RBW_f21270_22 +
2107 SF_Am_P_TP_4 * RBW_a11260m_2 * RBW_f21270_21 );
2108
2109 // D0 -> a1(1260)- {f0 pi- [P]} pi+
2110 double SF_Am_P_SP_1 = contract_11_0( T1_3m1, T1_Pip1Pim2Pim1 );
2111 double SF_Am_P_SP_2 = contract_11_0( T1_3m1, T1_Pip1Pim1Pim2 );
2112 double SF_Am_P_SP_3 = contract_11_0( T1_3m2, T1_Pip2Pim2Pim1 );
2113 double SF_Am_P_SP_4 = contract_11_0( T1_3m2, T1_Pip2Pim1Pim2 );
2114
2115 amplitude +=
2116 g_fitpara[3] * g_fitpara[6] *
2117 ( SF_Am_P_SP_1 * RBW_a11260m_1 * PiPiS_12_0 + SF_Am_P_SP_2 * RBW_a11260m_1 * PiPiS_11_0 +
2118 SF_Am_P_SP_3 * RBW_a11260m_2 * PiPiS_22_0 +
2119 SF_Am_P_SP_4 * RBW_a11260m_2 * PiPiS_21_0 );
2120 amplitude +=
2121 g_fitpara[4] * g_fitpara[6] *
2122 ( SF_Am_P_SP_1 * RBW_a11260m_1 * PiPiS_12_1 + SF_Am_P_SP_2 * RBW_a11260m_1 * PiPiS_11_1 +
2123 SF_Am_P_SP_3 * RBW_a11260m_2 * PiPiS_22_1 +
2124 SF_Am_P_SP_4 * RBW_a11260m_2 * PiPiS_21_1 );
2125 amplitude +=
2126 g_fitpara[5] * g_fitpara[6] *
2127 ( SF_Am_P_SP_1 * RBW_a11260m_1 * PiPiS_12_5 + SF_Am_P_SP_2 * RBW_a11260m_1 * PiPiS_11_5 +
2128 SF_Am_P_SP_3 * RBW_a11260m_2 * PiPiS_22_5 +
2129 SF_Am_P_SP_4 * RBW_a11260m_2 * PiPiS_21_5 );
2130
2131 // D0 -> a1(1420)+ {rho(770)0 pi+[S]} pi-
2132 // amplitude += g_fitpara[7]*(SF_Ap_S_VP_1*RBW_a11420p_1*GS_rho770_21 +
2133 // SF_Ap_S_VP_2*RBW_a11420p_1*GS_rho770_11 + SF_Ap_S_VP_3*RBW_a11420p_2*GS_rho770_22 +
2134 // SF_Ap_S_VP_4*RBW_a11420p_2*GS_rho770_12);
2135
2136 // D0 -> a1(1420)+ {f0 pi+[P]} pi-
2137 amplitude += g_fitpara[7] * ( SF_Ap_P_SP_1 * RBW_a11420p_1 * PiPiS_21_5 +
2138 SF_Ap_P_SP_2 * RBW_a11420p_1 * PiPiS_11_5 +
2139 SF_Ap_P_SP_3 * RBW_a11420p_2 * PiPiS_22_5 +
2140 SF_Ap_P_SP_4 * RBW_a11420p_2 * PiPiS_12_5 );
2141 amplitude += g_fitpara[8] * ( SF_Ap_P_SP_1 * RBW_a11420p_1 * PiPiS_21_6 +
2142 SF_Ap_P_SP_2 * RBW_a11420p_1 * PiPiS_11_6 +
2143 SF_Ap_P_SP_3 * RBW_a11420p_2 * PiPiS_22_6 +
2144 SF_Ap_P_SP_4 * RBW_a11420p_2 * PiPiS_12_6 );
2145
2146 vector<double> m_epsilon_uvmn;
2147 m_epsilon_uvmn.clear();
2148 for ( int i = 0; i < 4; i++ )
2149 {
2150 for ( int j = 0; j < 4; j++ )
2151 {
2152 for ( int k = 0; k < 4; k++ )
2153 {
2154 for ( int l = 0; l < 4; l++ )
2155 {
2156 if ( i == j || i == k || i == l || j == k || j == l || k == l )
2157 { m_epsilon_uvmn.push_back( 0.0 ); }
2158 else
2159 {
2160 if ( i == 0 && j == 1 && k == 2 && l == 3 ) m_epsilon_uvmn.push_back( 1.0 );
2161 if ( i == 0 && j == 1 && k == 3 && l == 2 ) m_epsilon_uvmn.push_back( -1.0 );
2162 if ( i == 0 && j == 2 && k == 1 && l == 3 ) m_epsilon_uvmn.push_back( -1.0 );
2163 if ( i == 0 && j == 2 && k == 3 && l == 1 ) m_epsilon_uvmn.push_back( 1.0 );
2164 if ( i == 0 && j == 3 && k == 1 && l == 2 ) m_epsilon_uvmn.push_back( 1.0 );
2165 if ( i == 0 && j == 3 && k == 2 && l == 1 ) m_epsilon_uvmn.push_back( -1.0 );
2166
2167 if ( i == 1 && j == 0 && k == 2 && l == 3 ) m_epsilon_uvmn.push_back( -1.0 );
2168 if ( i == 1 && j == 0 && k == 3 && l == 2 ) m_epsilon_uvmn.push_back( 1.0 );
2169 if ( i == 1 && j == 2 && k == 0 && l == 3 ) m_epsilon_uvmn.push_back( 1.0 );
2170 if ( i == 1 && j == 2 && k == 3 && l == 0 ) m_epsilon_uvmn.push_back( -1.0 );
2171 if ( i == 1 && j == 3 && k == 0 && l == 2 ) m_epsilon_uvmn.push_back( -1.0 );
2172 if ( i == 1 && j == 3 && k == 2 && l == 0 ) m_epsilon_uvmn.push_back( 1.0 );
2173
2174 if ( i == 2 && j == 0 && k == 1 && l == 3 ) m_epsilon_uvmn.push_back( 1.0 );
2175 if ( i == 2 && j == 0 && k == 3 && l == 1 ) m_epsilon_uvmn.push_back( -1.0 );
2176 if ( i == 2 && j == 1 && k == 0 && l == 3 ) m_epsilon_uvmn.push_back( -1.0 );
2177 if ( i == 2 && j == 1 && k == 3 && l == 0 ) m_epsilon_uvmn.push_back( 1.0 );
2178 if ( i == 2 && j == 3 && k == 0 && l == 1 ) m_epsilon_uvmn.push_back( 1.0 );
2179 if ( i == 2 && j == 3 && k == 1 && l == 0 ) m_epsilon_uvmn.push_back( -1.0 );
2180
2181 if ( i == 3 && j == 0 && k == 1 && l == 2 ) m_epsilon_uvmn.push_back( -1.0 );
2182 if ( i == 3 && j == 0 && k == 2 && l == 1 ) m_epsilon_uvmn.push_back( 1.0 );
2183 if ( i == 3 && j == 1 && k == 0 && l == 2 ) m_epsilon_uvmn.push_back( 1.0 );
2184 if ( i == 3 && j == 1 && k == 2 && l == 0 ) m_epsilon_uvmn.push_back( -1.0 );
2185 if ( i == 3 && j == 2 && k == 0 && l == 1 ) m_epsilon_uvmn.push_back( -1.0 );
2186 if ( i == 3 && j == 2 && k == 1 && l == 0 ) m_epsilon_uvmn.push_back( 1.0 );
2187 }
2188 }
2189 }
2190 }
2191 }
2192 // D0 -> a2(1320)+ {rho(770)0 pi+[D]} pi-
2193 double SF_Tp_D_VP_1 = contract_22_0(
2194 contract_22_2( contract_31_2( contract_41_3( m_epsilon_uvmn,
2195 contract_21_1( Proj1_3p1, T1_Pip2Pim1 ) ),
2196 Pip1Pip2Pim1 ),
2197 contract_42_2( Proj2_3p1, T2_3p1 ) ),
2198 T2_Pip2Pim1Pip1 );
2199 double SF_Tp_D_VP_2 = contract_22_0(
2200 contract_22_2( contract_31_2( contract_41_3( m_epsilon_uvmn,
2201 contract_21_1( Proj1_3p1, T1_Pip1Pim1 ) ),
2202 Pip1Pip2Pim1 ),
2203 contract_42_2( Proj2_3p1, T2_3p1 ) ),
2204 T2_Pip1Pim1Pip2 );
2205 double SF_Tp_D_VP_3 = contract_22_0(
2206 contract_22_2( contract_31_2( contract_41_3( m_epsilon_uvmn,
2207 contract_21_1( Proj1_3p2, T1_Pip2Pim2 ) ),
2208 Pip1Pip2Pim2 ),
2209 contract_42_2( Proj2_3p2, T2_3p2 ) ),
2210 T2_Pip2Pim2Pip1 );
2211 double SF_Tp_D_VP_4 = contract_22_0(
2212 contract_22_2( contract_31_2( contract_41_3( m_epsilon_uvmn,
2213 contract_21_1( Proj1_3p2, T1_Pip1Pim2 ) ),
2214 Pip1Pip2Pim2 ),
2215 contract_42_2( Proj2_3p2, T2_3p2 ) ),
2216 T2_Pip1Pim2Pip2 );
2217
2218 amplitude += g_fitpara[9] * ( SF_Tp_D_VP_1 * RBW_a21320p_1 * GS_rho770_21 +
2219 SF_Tp_D_VP_2 * RBW_a21320p_1 * GS_rho770_11 +
2220 SF_Tp_D_VP_3 * RBW_a21320p_2 * GS_rho770_22 +
2221 SF_Tp_D_VP_4 * RBW_a21320p_2 * GS_rho770_12 );
2222
2223 // D0 -> a2(1320)- {rho(770)0 pi-[D]} pi+
2224 double SF_Tm_D_VP_1 = contract_22_0(
2225 contract_22_2( contract_31_2( contract_41_3( m_epsilon_uvmn,
2226 contract_21_1( Proj1_3m1, T1_Pim2Pip1 ) ),
2227 Pim1Pim2Pip1 ),
2228 contract_42_2( Proj2_3m1, T2_3m1 ) ),
2229 T2_Pip1Pim2Pim1 );
2230 double SF_Tm_D_VP_2 = contract_22_0(
2231 contract_22_2( contract_31_2( contract_41_3( m_epsilon_uvmn,
2232 contract_21_1( Proj1_3m1, T1_Pim1Pip1 ) ),
2233 Pim1Pim2Pip1 ),
2234 contract_42_2( Proj2_3m1, T2_3m1 ) ),
2235 T2_Pip1Pim1Pim2 );
2236 double SF_Tm_D_VP_3 = contract_22_0(
2237 contract_22_2( contract_31_2( contract_41_3( m_epsilon_uvmn,
2238 contract_21_1( Proj1_3m2, T1_Pim2Pip2 ) ),
2239 Pim1Pim2Pip2 ),
2240 contract_42_2( Proj2_3m2, T2_3m2 ) ),
2241 T2_Pip2Pim2Pim1 );
2242 double SF_Tm_D_VP_4 = contract_22_0(
2243 contract_22_2( contract_31_2( contract_41_3( m_epsilon_uvmn,
2244 contract_21_1( Proj1_3m2, T1_Pim1Pip2 ) ),
2245 Pim1Pim2Pip2 ),
2246 contract_42_2( Proj2_3m2, T2_3m2 ) ),
2247 T2_Pip2Pim1Pim2 );
2248
2249 amplitude += g_fitpara[10] * ( SF_Tm_D_VP_1 * RBW_a21320m_1 * GS_rho770_12 +
2250 SF_Tm_D_VP_2 * RBW_a21320m_1 * GS_rho770_11 +
2251 SF_Tm_D_VP_3 * RBW_a21320m_2 * GS_rho770_22 +
2252 SF_Tm_D_VP_4 * RBW_a21320m_2 * GS_rho770_21 );
2253
2254 // D0 -> pi(1300)- {rho(770)0 pi-} pi+
2255 double SF_Pm_P_VP_1 = contract_11_0( T1_Pim2Pip1, T1_Pip1Pim2Pim1 );
2256 double SF_Pm_P_VP_2 = contract_11_0( T1_Pim1Pip1, T1_Pip1Pim1Pim2 );
2257 double SF_Pm_P_VP_3 = contract_11_0( T1_Pim2Pip2, T1_Pip2Pim2Pim1 );
2258 double SF_Pm_P_VP_4 = contract_11_0( T1_Pim1Pip2, T1_Pip2Pim1Pim2 );
2259
2260 amplitude += g_fitpara[11] * ( SF_Pm_P_VP_1 * GS_rho770_12 * RBW_pi1300m_1 +
2261 SF_Pm_P_VP_2 * GS_rho770_11 * RBW_pi1300m_1 +
2262 SF_Pm_P_VP_3 * GS_rho770_22 * RBW_pi1300m_2 +
2263 SF_Pm_P_VP_4 * GS_rho770_21 * RBW_pi1300m_2 );
2264 /*
2265 // D0 -> pi(1300)- {f2(1270)0 pi-} pi+
2266 double SF_Pm_D_TP_1 = contract_22_0(T2_Pip1Pim2,T2_Pip1Pim2Pim1);
2267 double SF_Pm_D_TP_2 = contract_22_0(T2_Pip1Pim1,T2_Pip1Pim1Pim2);
2268 double SF_Pm_D_TP_3 = contract_22_0(T2_Pip2Pim2,T2_Pip2Pim2Pim1);
2269 double SF_Pm_D_TP_4 = contract_22_0(T2_Pip2Pim1,T2_Pip2Pim1Pim2);
2270
2271 amplitude += g_fitpara[12]*g_fitpara[11]*(SF_Pm_D_TP_1*RBW_f21270_12*RBW_pi1300m_1 +
2272 SF_Pm_D_TP_2*RBW_f21270_11*RBW_pi1300m_1 + SF_Pm_D_TP_3*RBW_f21270_22*RBW_pi1300m_2 +
2273 SF_Pm_D_TP_4*RBW_f21270_21*RBW_pi1300m_2);
2274 */
2275 // D0 -> pi(1300)- {f0 pi-} pi+
2276 amplitude += g_fitpara[12] * g_fitpara[11] *
2277 ( RBW_pi1300m_1 * PiPiS_12_0 + RBW_pi1300m_1 * PiPiS_11_0 +
2278 RBW_pi1300m_2 * PiPiS_22_0 + RBW_pi1300m_2 * PiPiS_21_0 );
2279 // amplitude += g_fitpara[13]*g_fitpara[11]*(RBW_pi1300m_1*PiPiS_12_5 +
2280 // RBW_pi1300m_1*PiPiS_11_5 + RBW_pi1300m_2*PiPiS_22_5 + RBW_pi1300m_2*PiPiS_21_5);
2281
2282 amplitude += g_fitpara[13] * g_fitpara[11] *
2283 ( RBW_pi1300m_1 * PiPiS_12_6 + RBW_pi1300m_1 * PiPiS_11_6 +
2284 RBW_pi1300m_2 * PiPiS_22_6 + RBW_pi1300m_2 * PiPiS_21_6 );
2285
2286 // D0 -> pi(1300)+ {rho(770)0 pi+} pi-
2287 double SF_Pp_P_VP_1 = contract_11_0( T1_Pip2Pim1, T1_Pip2Pim1Pip1 );
2288 double SF_Pp_P_VP_2 = contract_11_0( T1_Pip1Pim1, T1_Pip1Pim1Pip2 );
2289 double SF_Pp_P_VP_3 = contract_11_0( T1_Pip2Pim2, T1_Pip2Pim2Pip1 );
2290 double SF_Pp_P_VP_4 = contract_11_0( T1_Pip1Pim2, T1_Pip1Pim2Pip2 );
2291
2292 amplitude += g_fitpara[14] * ( SF_Pp_P_VP_1 * GS_rho770_21 * RBW_pi1300p_1 +
2293 SF_Pp_P_VP_2 * GS_rho770_11 * RBW_pi1300p_1 +
2294 SF_Pp_P_VP_3 * GS_rho770_22 * RBW_pi1300p_2 +
2295 SF_Pp_P_VP_4 * GS_rho770_12 * RBW_pi1300p_2 );
2296 /*
2297 // D0 -> pi(1300)+ {f2(1270)0 pi+} pi-
2298 double SF_Pp_D_TP_1 = contract_22_0(T2_Pip2Pim1,T2_Pip2Pim1Pip1);
2299 double SF_Pp_D_TP_2 = contract_22_0(T2_Pip1Pim1,T2_Pip1Pim1Pip2);
2300 double SF_Pp_D_TP_3 = contract_22_0(T2_Pip2Pim2,T2_Pip2Pim2Pip1);
2301 double SF_Pp_D_TP_4 = contract_22_0(T2_Pip1Pim2,T2_Pip1Pim2Pip2);
2302
2303 amplitude += g_fitpara[12]*g_fitpara[15]*(SF_Pp_D_TP_1*RBW_f21270_21*RBW_pi1300p_1 +
2304 SF_Pp_D_TP_2*RBW_f21270_11*RBW_pi1300p_1 + SF_Pp_D_TP_3*RBW_f21270_22*RBW_pi1300p_2 +
2305 SF_Pp_D_TP_4*RBW_f21270_12*RBW_pi1300p_2);
2306 */
2307 // D0 -> pi(1300)+ {f0 pi+} pi-
2308 amplitude += g_fitpara[12] * g_fitpara[14] *
2309 ( RBW_pi1300p_1 * PiPiS_21_0 + RBW_pi1300p_1 * PiPiS_11_0 +
2310 RBW_pi1300p_2 * PiPiS_22_0 + RBW_pi1300p_2 * PiPiS_12_0 );
2311 // amplitude += g_fitpara[13]*g_fitpara[15]*(RBW_pi1300p_1*PiPiS_21_5 +
2312 // RBW_pi1300p_1*PiPiS_11_5 + RBW_pi1300p_2*PiPiS_22_5 + RBW_pi1300p_2*PiPiS_12_5);
2313
2314 amplitude += g_fitpara[13] * g_fitpara[14] *
2315 ( RBW_pi1300p_1 * PiPiS_21_6 + RBW_pi1300p_1 * PiPiS_11_6 +
2316 RBW_pi1300p_2 * PiPiS_22_6 + RBW_pi1300p_2 * PiPiS_12_6 );
2317
2318 // D0 -> rho rho [S]
2319 double SF_VV_S_1 = contract_11_0( T1_Pip1Pim1, T1_Pip2Pim2 );
2320 double SF_VV_S_3 = contract_11_0( T1_Pip1Pim2, T1_Pip2Pim1 );
2321
2322 amplitude += g_fitpara[15] * ( SF_VV_S_1 * GS_rho770_11 * GS_rho770_22 +
2323 SF_VV_S_3 * GS_rho770_12 * GS_rho770_21 );
2324
2325 // D0 -> rho rho [P]
2326 double SF_VV_P_1 = contract_11_0(
2327 contract_21_1(
2328 contract_31_2( contract_41_3( m_epsilon_uvmn, T1_Pip1Pim1 ), T1_Pip2Pim2 ),
2329 T1_2z11 ),
2330 D0 );
2331 double SF_VV_P_2 = contract_11_0(
2332 contract_21_1(
2333 contract_31_2( contract_41_3( m_epsilon_uvmn, T1_Pip2Pim2 ), T1_Pip1Pim1 ),
2334 T1_2z12 ),
2335 D0 );
2336 double SF_VV_P_3 = contract_11_0(
2337 contract_21_1(
2338 contract_31_2( contract_41_3( m_epsilon_uvmn, T1_Pip1Pim2 ), T1_Pip2Pim1 ),
2339 T1_2z21 ),
2340 D0 );
2341 double SF_VV_P_4 = contract_11_0(
2342 contract_21_1(
2343 contract_31_2( contract_41_3( m_epsilon_uvmn, T1_Pip2Pim1 ), T1_Pip1Pim2 ),
2344 T1_2z22 ),
2345 D0 );
2346
2347 amplitude += g_fitpara[16] * ( SF_VV_P_1 * GS_rho770_11 * GS_rho770_22 +
2348 SF_VV_P_3 * GS_rho770_12 * GS_rho770_21 );
2349
2350 // D0 -> rho rho [D]
2351 double SF_VV_D_1 = contract_11_0( contract_21_1( T2_2z11, T1_Pip2Pim2 ), T1_Pip1Pim1 );
2352 double SF_VV_D_3 = contract_11_0( contract_21_1( T2_2z21, T1_Pip2Pim1 ), T1_Pip1Pim2 );
2353
2354 amplitude += g_fitpara[17] * ( SF_VV_D_1 * GS_rho770_11 * GS_rho770_22 +
2355 SF_VV_D_3 * GS_rho770_12 * GS_rho770_21 );
2356
2357 // D0 -> rho rho' [P]
2358 amplitude +=
2359 g_fitpara[18] *
2360 ( SF_VV_P_1 * GS_rho770_11 * GS_rho1450_22 + SF_VV_P_2 * GS_rho770_22 * GS_rho1450_11 +
2361 SF_VV_P_3 * GS_rho770_12 * GS_rho1450_21 + SF_VV_P_3 * GS_rho770_21 * GS_rho1450_12 );
2362
2363 // D0 ->rho f0
2364 double SF_VS_P_1 = contract_11_0( T1_Pip1Pim1, T1_2z11 );
2365 double SF_VS_P_2 = contract_11_0( T1_Pip2Pim2, T1_2z12 );
2366 double SF_VS_P_3 = contract_11_0( T1_Pip1Pim2, T1_2z21 );
2367 double SF_VS_P_4 = contract_11_0( T1_Pip2Pim1, T1_2z22 );
2368
2369 amplitude +=
2370 g_fitpara[19] *
2371 ( SF_VS_P_1 * GS_rho770_11 * PiPiS_22_0 + SF_VS_P_2 * GS_rho770_22 * PiPiS_11_0 +
2372 SF_VS_P_3 * GS_rho770_12 * PiPiS_21_0 + SF_VS_P_4 * GS_rho770_21 * PiPiS_12_0 );
2373 amplitude +=
2374 g_fitpara[20] *
2375 ( SF_VS_P_1 * GS_rho770_11 * PiPiS_22_5 + SF_VS_P_2 * GS_rho770_22 * PiPiS_11_5 +
2376 SF_VS_P_3 * GS_rho770_12 * PiPiS_21_5 + SF_VS_P_4 * GS_rho770_21 * PiPiS_12_5 );
2377 amplitude +=
2378 g_fitpara[21] *
2379 ( SF_VS_P_1 * GS_rho770_11 * PiPiS_22_6 + SF_VS_P_2 * GS_rho770_22 * PiPiS_11_6 +
2380 SF_VS_P_3 * GS_rho770_12 * PiPiS_21_6 + SF_VS_P_4 * GS_rho770_21 * PiPiS_12_6 );
2381
2382 // D0 -> f0f0
2383 // 00
2384 amplitude += g_fitpara[22] * ( PiPiS_11_0 * PiPiS_22_0 + PiPiS_12_0 * PiPiS_21_0 +
2385 PiPiS_22_0 * PiPiS_11_0 + PiPiS_21_0 * PiPiS_12_0 );
2386 amplitude += g_fitpara[23] * ( PiPiS_11_0 * PiPiS_22_1 + PiPiS_12_0 * PiPiS_21_1 +
2387 PiPiS_22_0 * PiPiS_11_1 + PiPiS_21_0 * PiPiS_12_1 );
2388 amplitude += g_fitpara[24] * ( PiPiS_11_1 * PiPiS_22_1 + PiPiS_12_1 * PiPiS_21_1 +
2389 PiPiS_22_1 * PiPiS_11_1 + PiPiS_21_1 * PiPiS_12_1 );
2390 amplitude += g_fitpara[25] * ( PiPiS_11_1 * PiPiS_22_5 + PiPiS_12_1 * PiPiS_21_5 +
2391 PiPiS_22_1 * PiPiS_11_5 + PiPiS_21_1 * PiPiS_12_5 );
2392 amplitude += g_fitpara[26] * ( PiPiS_11_5 * PiPiS_22_5 + PiPiS_12_5 * PiPiS_21_5 +
2393 PiPiS_22_5 * PiPiS_11_5 + PiPiS_21_5 * PiPiS_12_5 );
2394 amplitude += g_fitpara[27] * ( PiPiS_11_5 * PiPiS_22_6 + PiPiS_12_5 * PiPiS_21_6 +
2395 PiPiS_22_5 * PiPiS_11_6 + PiPiS_21_5 * PiPiS_12_6 );
2396
2397 // D0 -> f2 f0
2398 double SF_TS_D_1 = contract_22_0( T2_Pip1Pim1, T2_2z11 );
2399 double SF_TS_D_2 = contract_22_0( T2_Pip2Pim2, T2_2z12 );
2400 double SF_TS_D_3 = contract_22_0( T2_Pip1Pim2, T2_2z21 );
2401 double SF_TS_D_4 = contract_22_0( T2_Pip2Pim1, T2_2z22 );
2402
2403 amplitude +=
2404 g_fitpara[28] *
2405 ( SF_TS_D_1 * RBW_f21270_11 * PiPiS_22_5 + SF_TS_D_2 * RBW_f21270_22 * PiPiS_11_5 +
2406 SF_TS_D_3 * RBW_f21270_12 * PiPiS_21_5 + SF_TS_D_4 * RBW_f21270_21 * PiPiS_12_5 );
2407 amplitude +=
2408 g_fitpara[29] *
2409 ( SF_TS_D_1 * RBW_f21270_11 * PiPiS_22_6 + SF_TS_D_2 * RBW_f21270_22 * PiPiS_11_6 +
2410 SF_TS_D_3 * RBW_f21270_12 * PiPiS_21_6 + SF_TS_D_4 * RBW_f21270_21 * PiPiS_12_6 );
2411
2412 return amplitude;
2413}
sin
double sin(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:185
cos
double cos(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:187

Referenced by QCMCFilter::findD0Decay().

◆ init()

void D0To2pip2pim::init ( )

Definition at line 13 of file D0To2pip2pim.cxx.

13 {
14
15 // std::cout << "D0To2pip2pim ==> Initialization !" << std::endl;
16
17 g_mag[0] = 100.0;
18 g_pha[0] = 0.0;
19 g_mag[1] = 7.95507;
20 g_pha[1] = -0.0687407;
21 g_mag[2] = 37.5559;
22 g_pha[2] = -1.74946;
23 g_mag[3] = 61.2172;
24 g_pha[3] = 2.98079;
25 g_mag[4] = 187.79;
26 g_pha[4] = 2.64471;
27 g_mag[5] = 385.474;
28 g_pha[5] = -0.137107;
29 g_mag[6] = 0.330788;
30 g_pha[6] = 0.268133;
31 g_mag[7] = 127.158;
32 g_pha[7] = -2.47773;
33 g_mag[8] = 339.914;
34 g_pha[8] = 2.22856;
35 g_mag[9] = 0.320888;
36 g_pha[9] = -2.6194;
37 g_mag[10] = 0.366283;
38 g_pha[10] = -0.26867;
39 g_mag[11] = 86.0865;
40 g_pha[11] = -2.49649;
41 g_mag[12] = 6.1541;
42 g_pha[12] = -1.18299;
43 g_mag[13] = 56.6067;
44 g_pha[13] = 0.142977;
45 g_mag[14] = 92.3073;
46 g_pha[14] = -2.15881;
47 g_mag[15] = 10.5885;
48 g_pha[15] = -3.03166;
49 g_mag[16] = 8.36765;
50 g_pha[16] = 1.8417;
51 g_mag[17] = 6.56437;
52 g_pha[17] = -2.93087;
53 g_mag[18] = 15.7197;
54 g_pha[18] = 0.96925;
55 g_mag[19] = 21.4195;
56 g_pha[19] = -1.23701;
57 g_mag[20] = 56.8867;
58 g_pha[20] = -0.385837;
59 g_mag[21] = 231.626;
60 g_pha[21] = 2.14842;
61 g_mag[22] = 2938.45;
62 g_pha[22] = -0.693491;
63 g_mag[23] = 7252.7;
64 g_pha[23] = 2.23659;
65 g_mag[24] = 5165.87;
66 g_pha[24] = 0.913557;
67 g_mag[25] = 11508.6;
68 g_pha[25] = -1.07187;
69 g_mag[26] = 2461.86;
70 g_pha[26] = 1.8709;
71 g_mag[27] = 8757.75;
72 g_pha[27] = 2.40756;
73 g_mag[28] = 19.7413;
74 g_pha[28] = -1.0753;
75 g_mag[29] = 66.3826;
76 g_pha[29] = 2.34666;
77 fitpara.clear();
78 for ( int i = 0; i < 30; i++ )
79 {
80 complex<double> ctemp( g_mag[i] * cos( g_pha[i] ), g_mag[i] * sin( g_pha[i] ) );
81 fitpara.push_back( ctemp );
82 }
83
84 // g_uv.clear();
85 // for(int i=0; i<4; i++){
86 // for(int j=0; j<4; j++){
87 // if(i!=j){
88 // g_uv.push_back(0.0);
89 // }else if(i<3){
90 // g_uv.push_back(-1.0);
91 // }else if(i==3){
92 // g_uv.push_back(1.0);
93 // }
94 // }
95 // }
96
97 epsilon_uvmn.clear();
98 for ( int i = 0; i < 4; i++ )
99 {
100 for ( int j = 0; j < 4; j++ )
101 {
102 for ( int k = 0; k < 4; k++ )
103 {
104 for ( int l = 0; l < 4; l++ )
105 {
106 if ( i == j || i == k || i == l || j == k || j == l || k == l )
107 { epsilon_uvmn.push_back( 0.0 ); }
108 else
109 {
110 if ( i == 0 && j == 1 && k == 2 && l == 3 ) epsilon_uvmn.push_back( 1.0 );
111 if ( i == 0 && j == 1 && k == 3 && l == 2 ) epsilon_uvmn.push_back( -1.0 );
112 if ( i == 0 && j == 2 && k == 1 && l == 3 ) epsilon_uvmn.push_back( -1.0 );
113 if ( i == 0 && j == 2 && k == 3 && l == 1 ) epsilon_uvmn.push_back( 1.0 );
114 if ( i == 0 && j == 3 && k == 1 && l == 2 ) epsilon_uvmn.push_back( 1.0 );
115 if ( i == 0 && j == 3 && k == 2 && l == 1 ) epsilon_uvmn.push_back( -1.0 );
116
117 if ( i == 1 && j == 0 && k == 2 && l == 3 ) epsilon_uvmn.push_back( -1.0 );
118 if ( i == 1 && j == 0 && k == 3 && l == 2 ) epsilon_uvmn.push_back( 1.0 );
119 if ( i == 1 && j == 2 && k == 0 && l == 3 ) epsilon_uvmn.push_back( 1.0 );
120 if ( i == 1 && j == 2 && k == 3 && l == 0 ) epsilon_uvmn.push_back( -1.0 );
121 if ( i == 1 && j == 3 && k == 0 && l == 2 ) epsilon_uvmn.push_back( -1.0 );
122 if ( i == 1 && j == 3 && k == 2 && l == 0 ) epsilon_uvmn.push_back( 1.0 );
123
124 if ( i == 2 && j == 0 && k == 1 && l == 3 ) epsilon_uvmn.push_back( 1.0 );
125 if ( i == 2 && j == 0 && k == 3 && l == 1 ) epsilon_uvmn.push_back( -1.0 );
126 if ( i == 2 && j == 1 && k == 0 && l == 3 ) epsilon_uvmn.push_back( -1.0 );
127 if ( i == 2 && j == 1 && k == 3 && l == 0 ) epsilon_uvmn.push_back( 1.0 );
128 if ( i == 2 && j == 3 && k == 0 && l == 1 ) epsilon_uvmn.push_back( 1.0 );
129 if ( i == 2 && j == 3 && k == 1 && l == 0 ) epsilon_uvmn.push_back( -1.0 );
130
131 if ( i == 3 && j == 0 && k == 1 && l == 2 ) epsilon_uvmn.push_back( -1.0 );
132 if ( i == 3 && j == 0 && k == 2 && l == 1 ) epsilon_uvmn.push_back( 1.0 );
133 if ( i == 3 && j == 1 && k == 0 && l == 2 ) epsilon_uvmn.push_back( 1.0 );
134 if ( i == 3 && j == 1 && k == 2 && l == 0 ) epsilon_uvmn.push_back( -1.0 );
135 if ( i == 3 && j == 2 && k == 0 && l == 1 ) epsilon_uvmn.push_back( -1.0 );
136 if ( i == 3 && j == 2 && k == 1 && l == 0 ) epsilon_uvmn.push_back( 1.0 );
137 }
138 }
139 }
140 }
141 }
142
143 _nd = 4;
144 math_pi = 3.1415926;
145 mass_Pion = 0.13957;
146
147 rRes = 3.0 * 0.197321;
148 rD = 5.0 * 0.197321;
149 m_Pi = 0.139570;
150 m2_Pi = m_Pi * m_Pi;
151
152 m0_rho770 = 0.77526;
153 w0_rho770 = 0.1478;
154
155 m0_rho1450 = 1.465;
156 w0_rho1450 = 0.400;
157
158 m0_f21270 = 1.2755;
159 w0_f21270 = 0.1867;
160
161 m0_a11260 = 1.1337;
162 g1_a11260 = 0.00335;
163 g2_a11260 = 0.0;
164
165 m0_pi1300 = 1.498;
166 w0_pi1300 = 0.590;
167
168 m0_a11420 = 1.411;
169 w0_a11420 = 0.161;
170
171 m0_a11640 = 1.655;
172 w0_a11640 = 0.254;
173
174 m0_a21320 = 1.3186;
175 w0_a21320 = 0.105;
176
177 m0_pi11400 = 1.354;
178 w0_pi11400 = 0.330;
179
180 s0_prod = -5.0;
181
182 return;
183}

Referenced by QCMCFilter::findD0Decay().

The documentation for this class was generated from the following files: