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

#include <D0Topippim2pi0.h>

Public Member Functions

 D0Topippim2pi0 ()
virtual ~D0Topippim2pi0 ()
void init ()
complex< double > Amp (vector< double > Pip, vector< double > Pim, vector< double > Pi01, vector< double > Pi02)

Detailed Description

Definition at line 9 of file D0Topippim2pi0.h.

Constructor & Destructor Documentation

◆ D0Topippim2pi0()

D0Topippim2pi0::D0Topippim2pi0 ( )
inline

Definition at line 12 of file D0Topippim2pi0.h.

12{}

◆ ~D0Topippim2pi0()

D0Topippim2pi0::~D0Topippim2pi0 ( )
virtual

Definition at line 11 of file D0Topippim2pi0.cxx.

11{}

Member Function Documentation

◆ Amp()

complex< double > D0Topippim2pi0::Amp ( vector< double > Pip,
vector< double > Pim,
vector< double > Pi01,
vector< double > Pi02 )

Definition at line 1626 of file D0Topippim2pi0.cxx.

1627 {
1628
1629 vector<double> PipPim;
1630 PipPim.clear();
1631 vector<double> PipPi01;
1632 PipPi01.clear();
1633 vector<double> PipPi02;
1634 PipPi02.clear();
1635 vector<double> PimPi01;
1636 PimPi01.clear();
1637 vector<double> PimPi02;
1638 PimPi02.clear();
1639 vector<double> Pi01Pi02;
1640 Pi01Pi02.clear();
1641
1642 PipPim = sum_tensor( Pip, Pim );
1643 PipPi01 = sum_tensor( Pip, Pi01 );
1644 PipPi02 = sum_tensor( Pip, Pi02 );
1645 PimPi01 = sum_tensor( Pim, Pi01 );
1646 PimPi02 = sum_tensor( Pim, Pi02 );
1647 Pi01Pi02 = sum_tensor( Pi01, Pi02 );
1648
1649 vector<double> PipPimPi01;
1650 PipPimPi01.clear();
1651 vector<double> PipPimPi02;
1652 PipPimPi02.clear();
1653 vector<double> PipPi01Pi02;
1654 PipPi01Pi02.clear();
1655 vector<double> PimPi01Pi02;
1656 PimPi01Pi02.clear();
1657
1658 PipPimPi01 = sum_tensor( PipPim, Pi01 );
1659 PipPimPi02 = sum_tensor( PipPim, Pi02 );
1660 PipPi01Pi02 = sum_tensor( PipPi01, Pi02 );
1661 PimPi01Pi02 = sum_tensor( PimPi01, Pi02 );
1662
1663 vector<double> D0;
1664 D0.clear();
1665 D0 = sum_tensor( PipPimPi01, Pi02 );
1666
1667 double M2_PipPim = contract_11_0( PipPim, PipPim );
1668 double M2_PipPi01 = contract_11_0( PipPi01, PipPi01 );
1669 double M2_PipPi02 = contract_11_0( PipPi02, PipPi02 );
1670 double M2_PimPi01 = contract_11_0( PimPi01, PimPi01 );
1671 double M2_PimPi02 = contract_11_0( PimPi02, PimPi02 );
1672 double M2_Pi01Pi02 = contract_11_0( Pi01Pi02, Pi01Pi02 );
1673
1674 double M2_PipPimPi01 = contract_11_0( PipPimPi01, PipPimPi01 );
1675 double M2_PipPimPi02 = contract_11_0( PipPimPi02, PipPimPi02 );
1676 double M2_PipPi01Pi02 = contract_11_0( PipPi01Pi02, PipPi01Pi02 );
1677 double M2_PimPi01Pi02 = contract_11_0( PimPi01Pi02, PimPi01Pi02 );
1678 double M2_D0 = contract_11_0( D0, D0 );
1679
1680 complex<double> GS_rho770_pm =
1681 GS( M2_PipPim, m0_rho7700, w0_rho7700, m2_Pi, m2_Pi, rRes, 1 );
1682 complex<double> GS_rho770_p1 =
1683 GS( M2_PipPi01, m0_rho770p, w0_rho770p, m2_Pi, m2_Pi0, rRes, 1 );
1684 complex<double> GS_rho770_p2 =
1685 GS( M2_PipPi02, m0_rho770p, w0_rho770p, m2_Pi, m2_Pi0, rRes, 1 );
1686 complex<double> GS_rho770_m1 =
1687 GS( M2_PimPi01, m0_rho770p, w0_rho770p, m2_Pi, m2_Pi0, rRes, 1 );
1688 complex<double> GS_rho770_m2 =
1689 GS( M2_PimPi02, m0_rho770p, w0_rho770p, m2_Pi, m2_Pi0, rRes, 1 );
1690
1691 complex<double> RBW_f21270_pm =
1692 RBW( M2_PipPim, m0_f21270, w0_f21270, m2_Pi, m2_Pi, rRes, 2 );
1693 complex<double> RBW_f21270_00 =
1694 RBW( M2_Pi01Pi02, m0_f21270, w0_f21270, m2_Pi0, m2_Pi0, rRes, 2 );
1695
1696 complex<double> PiPiS_pm_0 = Fvector( M2_PipPim, s0_prod, 0 );
1697 complex<double> PiPiS_00_0 = Fvector( M2_Pi01Pi02, s0_prod, 0 );
1698
1699 complex<double> PiPiS_pm_1 = Fvector( M2_PipPim, s0_prod, 1 );
1700 complex<double> PiPiS_00_1 = Fvector( M2_Pi01Pi02, s0_prod, 1 );
1701
1702 complex<double> PiPiS_pm_5 = Fvector( M2_PipPim, s0_prod, 5 );
1703 complex<double> PiPiS_00_5 = Fvector( M2_Pi01Pi02, s0_prod, 5 );
1704
1705 complex<double> PiPiS_pm_6 = Fvector( M2_PipPim, s0_prod, 6 );
1706 complex<double> PiPiS_00_6 = Fvector( M2_Pi01Pi02, s0_prod, 6 );
1707
1708 complex<double> RBW_a11260_p = RBWa1260( M2_PipPi01Pi02, m0_a11260, g1_a11260, g2_a11260 );
1709 complex<double> RBW_a11260_m = RBWa1260( M2_PimPi01Pi02, m0_a11260, g1_a11260, g2_a11260 );
1710 complex<double> RBW_a11260_01 = RBWa1260( M2_PipPimPi01, m0_a11260, g1_a11260, g2_a11260 );
1711 complex<double> RBW_a11260_02 = RBWa1260( M2_PipPimPi02, m0_a11260, g1_a11260, g2_a11260 );
1712
1713 complex<double> RBW_a11420_p = RBW( M2_PipPi01Pi02, m0_a11420, w0_a11420, -1, -1, -1, -1 );
1714 complex<double> RBW_a11420_m = RBW( M2_PimPi01Pi02, m0_a11420, w0_a11420, -1, -1, -1, -1 );
1715 complex<double> RBW_a11420_01 = RBW( M2_PipPimPi01, m0_a11420, w0_a11420, -1, -1, -1, -1 );
1716 complex<double> RBW_a11420_02 = RBW( M2_PipPimPi02, m0_a11420, w0_a11420, -1, -1, -1, -1 );
1717
1718 complex<double> RBW_omega_01 = RBW( M2_PipPimPi01, m0_omega, w0_omega, -1, -1, -1, -1 );
1719 complex<double> RBW_omega_02 = RBW( M2_PipPimPi02, m0_omega, w0_omega, -1, -1, -1, -1 );
1720
1721 complex<double> RBW_phi_01 = RBW( M2_PipPimPi01, m0_phi, w0_phi, -1, -1, -1, -1 );
1722 complex<double> RBW_phi_02 = RBW( M2_PipPimPi02, m0_phi, w0_phi, -1, -1, -1, -1 );
1723
1724 complex<double> RBW_a21320_p = RBW( M2_PipPi01Pi02, m0_a21320, w0_a21320, -1, -1, -1, -1 );
1725 complex<double> RBW_a21320_m = RBW( M2_PimPi01Pi02, m0_a21320, w0_a21320, -1, -1, -1, -1 );
1726
1727 complex<double> RBW_pi1300_p = RBWpi1300( M2_PipPi01Pi02, m0_pi1300, w0_pi1300 );
1728 complex<double> RBW_pi1300_m = RBWpi1300( M2_PimPi01Pi02, m0_pi1300, w0_pi1300 );
1729 complex<double> RBW_pi1300_01 = RBWpi1300( M2_PipPimPi01, m0_pi1300, w0_pi1300 );
1730 complex<double> RBW_pi1300_02 = RBWpi1300( M2_PipPimPi02, m0_pi1300, w0_pi1300 );
1731
1732 complex<double> RBW_h11170_01 = RBW( M2_PipPimPi01, m0_h11170, w0_h11170, -1, -1, -1, -1 );
1733 complex<double> RBW_h11170_02 = RBW( M2_PipPimPi02, m0_h11170, w0_h11170, -1, -1, -1, -1 );
1734
1735 complex<double> RBW_pi21670_01 =
1736 RBW( M2_PipPimPi01, m0_pi21670, w0_pi21670, -1, -1, -1, -1 );
1737 complex<double> RBW_pi21670_02 =
1738 RBW( M2_PipPimPi02, m0_pi21670, w0_pi21670, -1, -1, -1, -1 );
1739
1740 // D->XX Projection
1741 vector<double> Proj1_3p;
1742 Proj1_3p.clear();
1743 vector<double> Proj1_3m;
1744 Proj1_3m.clear();
1745 vector<double> Proj1_3z1;
1746 Proj1_3z1.clear();
1747 vector<double> Proj1_3z2;
1748 Proj1_3z2.clear();
1749
1750 Proj1_3p = ProjectionTensors( PipPi01Pi02, 1 );
1751 Proj1_3m = ProjectionTensors( PimPi01Pi02, 1 );
1752 Proj1_3z1 = ProjectionTensors( PipPimPi01, 1 );
1753 Proj1_3z2 = ProjectionTensors( PipPimPi02, 1 );
1754
1755 vector<double> Proj2_3p;
1756 Proj2_3p.clear();
1757 vector<double> Proj2_3m;
1758 Proj2_3m.clear();
1759 vector<double> Proj2_3z1;
1760 Proj2_3z1.clear();
1761 vector<double> Proj2_3z2;
1762 Proj2_3z2.clear();
1763
1764 Proj2_3p = ProjectionTensors( PipPi01Pi02, 2 );
1765 Proj2_3m = ProjectionTensors( PimPi01Pi02, 2 );
1766 Proj2_3z1 = ProjectionTensors( PipPimPi01, 2 );
1767 Proj2_3z2 = ProjectionTensors( PipPimPi02, 2 );
1768
1769 // X->PP Orbital
1770 vector<double> T1_PipPim;
1771 T1_PipPim.clear();
1772 vector<double> T1_PipPi01;
1773 T1_PipPi01.clear();
1774 vector<double> T1_PipPi02;
1775 T1_PipPi02.clear();
1776 vector<double> T1_PimPi01;
1777 T1_PimPi01.clear();
1778 vector<double> T1_PimPi02;
1779 T1_PimPi02.clear();
1780 vector<double> T1_Pi01Pi02;
1781 T1_Pi01Pi02.clear();
1782
1783 T1_PipPim = OrbitalTensors( PipPim, Pip, Pim, rRes, 1 );
1784 T1_PipPi01 = OrbitalTensors( PipPi01, Pip, Pi01, rRes, 1 );
1785 T1_PipPi02 = OrbitalTensors( PipPi02, Pip, Pi02, rRes, 1 );
1786 T1_PimPi01 = OrbitalTensors( PimPi01, Pim, Pi01, rRes, 1 );
1787 T1_PimPi02 = OrbitalTensors( PimPi02, Pim, Pi02, rRes, 1 );
1788 T1_Pi01Pi02 = OrbitalTensors( Pi01Pi02, Pi01, Pi02, rRes, 1 );
1789
1790 vector<double> T2_PipPim;
1791 T2_PipPim.clear();
1792 vector<double> T2_Pi01Pi02;
1793 T2_Pi01Pi02.clear();
1794
1795 T2_PipPim = OrbitalTensors( PipPim, Pip, Pim, rRes, 2 );
1796 T2_Pi01Pi02 = OrbitalTensors( Pi01Pi02, Pi01, Pi02, rRes, 2 );
1797
1798 // X->YP Orbital
1799 vector<double> T1_PipPimPi01;
1800 T1_PipPimPi01.clear();
1801 vector<double> T1_PipPimPi02;
1802 T1_PipPimPi02.clear();
1803 vector<double> T1_PipPi01Pi02;
1804 T1_PipPi01Pi02.clear();
1805 vector<double> T1_PipPi02Pi01;
1806 T1_PipPi02Pi01.clear();
1807 vector<double> T1_PimPi01Pi02;
1808 T1_PimPi01Pi02.clear();
1809 vector<double> T1_PimPi02Pi01;
1810 T1_PimPi02Pi01.clear();
1811 vector<double> T1_PipPi01Pim;
1812 T1_PipPi01Pim.clear();
1813 vector<double> T1_PipPi02Pim;
1814 T1_PipPi02Pim.clear();
1815 vector<double> T1_PimPi01Pip;
1816 T1_PimPi01Pip.clear();
1817 vector<double> T1_PimPi02Pip;
1818 T1_PimPi02Pip.clear();
1819 vector<double> T1_Pi01Pi02Pip;
1820 T1_Pi01Pi02Pip.clear();
1821 vector<double> T1_Pi01Pi02Pim;
1822 T1_Pi01Pi02Pim.clear();
1823
1824 T1_PipPimPi01 = OrbitalTensors( PipPimPi01, PipPim, Pi01, rRes, 1 );
1825 T1_PipPimPi02 = OrbitalTensors( PipPimPi02, PipPim, Pi02, rRes, 1 );
1826 T1_PipPi01Pi02 = OrbitalTensors( PipPi01Pi02, PipPi01, Pi02, rRes, 1 );
1827 T1_PipPi02Pi01 = OrbitalTensors( PipPi01Pi02, PipPi02, Pi01, rRes, 1 );
1828 T1_PimPi01Pi02 = OrbitalTensors( PimPi01Pi02, PimPi01, Pi02, rRes, 1 );
1829 T1_PimPi02Pi01 = OrbitalTensors( PimPi01Pi02, PimPi02, Pi01, rRes, 1 );
1830 T1_PipPi01Pim = OrbitalTensors( PipPimPi01, PipPi01, Pim, rRes, 1 );
1831 T1_PipPi02Pim = OrbitalTensors( PipPimPi02, PipPi02, Pim, rRes, 1 );
1832 T1_PimPi01Pip = OrbitalTensors( PipPimPi01, PimPi01, Pip, rRes, 1 );
1833 T1_PimPi02Pip = OrbitalTensors( PipPimPi02, PimPi02, Pip, rRes, 1 );
1834 T1_Pi01Pi02Pip = OrbitalTensors( PipPi01Pi02, Pi01Pi02, Pip, rRes, 1 );
1835 T1_Pi01Pi02Pim = OrbitalTensors( PimPi01Pi02, Pi01Pi02, Pim, rRes, 1 );
1836
1837 vector<double> T2_PipPimPi01;
1838 T2_PipPimPi01.clear();
1839 vector<double> T2_PipPimPi02;
1840 T2_PipPimPi02.clear();
1841 vector<double> T2_PipPi01Pi02;
1842 T2_PipPi01Pi02.clear();
1843 vector<double> T2_PipPi02Pi01;
1844 T2_PipPi02Pi01.clear();
1845 vector<double> T2_PimPi01Pi02;
1846 T2_PimPi01Pi02.clear();
1847 vector<double> T2_PimPi02Pi01;
1848 T2_PimPi02Pi01.clear();
1849 vector<double> T2_PipPi01Pim;
1850 T2_PipPi01Pim.clear();
1851 vector<double> T2_PipPi02Pim;
1852 T2_PipPi02Pim.clear();
1853 vector<double> T2_PimPi01Pip;
1854 T2_PimPi01Pip.clear();
1855 vector<double> T2_PimPi02Pip;
1856 T2_PimPi02Pip.clear();
1857 vector<double> T2_Pi01Pi02Pip;
1858 T2_Pi01Pi02Pip.clear();
1859 vector<double> T2_Pi01Pi02Pim;
1860 T2_Pi01Pi02Pim.clear();
1861
1862 T2_PipPimPi01 = OrbitalTensors( PipPimPi01, PipPim, Pi01, rRes, 2 );
1863 T2_PipPimPi02 = OrbitalTensors( PipPimPi02, PipPim, Pi02, rRes, 2 );
1864 T2_PipPi01Pi02 = OrbitalTensors( PipPi01Pi02, PipPi01, Pi02, rRes, 2 );
1865 T2_PipPi02Pi01 = OrbitalTensors( PipPi01Pi02, PipPi02, Pi01, rRes, 2 );
1866 T2_PimPi01Pi02 = OrbitalTensors( PimPi01Pi02, PimPi01, Pi02, rRes, 2 );
1867 T2_PimPi02Pi01 = OrbitalTensors( PimPi01Pi02, PimPi02, Pi01, rRes, 2 );
1868 T2_PipPi01Pim = OrbitalTensors( PipPimPi01, PipPi01, Pim, rRes, 2 );
1869 T2_PipPi02Pim = OrbitalTensors( PipPimPi02, PipPi02, Pim, rRes, 2 );
1870 T2_PimPi01Pip = OrbitalTensors( PipPimPi01, PimPi01, Pip, rRes, 2 );
1871 T2_PimPi02Pip = OrbitalTensors( PipPimPi02, PimPi02, Pip, rRes, 2 );
1872 T2_Pi01Pi02Pip = OrbitalTensors( PipPi01Pi02, Pi01Pi02, Pip, rRes, 2 );
1873 T2_Pi01Pi02Pim = OrbitalTensors( PimPi01Pi02, Pi01Pi02, Pim, rRes, 2 );
1874
1875 // D->XX Orbital
1876 vector<double> T1_2pm12;
1877 T1_2pm12.clear();
1878 vector<double> T1_2p1m2;
1879 T1_2p1m2.clear();
1880 vector<double> T1_2p2m1;
1881 T1_2p2m1.clear();
1882
1883 T1_2pm12 = OrbitalTensors( D0, PipPim, Pi01Pi02, rD, 1 );
1884 T1_2p1m2 = OrbitalTensors( D0, PipPi01, PimPi02, rD, 1 );
1885 T1_2p2m1 = OrbitalTensors( D0, PipPi02, PimPi01, rD, 1 );
1886
1887 vector<double> T2_2pm12;
1888 T2_2pm12.clear();
1889 vector<double> T2_2p1m2;
1890 T2_2p1m2.clear();
1891 vector<double> T2_2p2m1;
1892 T2_2p2m1.clear();
1893
1894 T2_2pm12 = OrbitalTensors( D0, PipPim, Pi01Pi02, rD, 2 );
1895 T2_2p1m2 = OrbitalTensors( D0, PipPi01, PimPi02, rD, 2 );
1896 T2_2p2m1 = OrbitalTensors( D0, PipPi02, PimPi01, rD, 2 );
1897
1898 // D->XP Orbital
1899 vector<double> T1_3pm;
1900 T1_3pm.clear();
1901 vector<double> T1_3mp;
1902 T1_3mp.clear();
1903 vector<double> T1_3z12;
1904 T1_3z12.clear();
1905 vector<double> T1_3z21;
1906 T1_3z21.clear();
1907
1908 T1_3pm = OrbitalTensors( D0, PipPi01Pi02, Pim, rD, 1 );
1909 T1_3mp = OrbitalTensors( D0, PimPi01Pi02, Pip, rD, 1 );
1910 T1_3z12 = OrbitalTensors( D0, PipPimPi01, Pi02, rD, 1 );
1911 T1_3z21 = OrbitalTensors( D0, PipPimPi02, Pi01, rD, 1 );
1912
1913 vector<double> T2_3pm;
1914 T2_3pm.clear();
1915 vector<double> T2_3mp;
1916 T2_3mp.clear();
1917 vector<double> T2_3z12;
1918 T2_3z12.clear();
1919 vector<double> T2_3z21;
1920 T2_3z21.clear();
1921
1922 T2_3pm = OrbitalTensors( D0, PipPi01Pi02, Pim, rD, 2 );
1923 T2_3mp = OrbitalTensors( D0, PimPi01Pi02, Pip, rD, 2 );
1924 T2_3z12 = OrbitalTensors( D0, PipPimPi01, Pi02, rD, 2 );
1925 T2_3z21 = OrbitalTensors( D0, PipPimPi02, Pi01, rD, 2 );
1926
1927 complex<double> amplitude( 0, 0 );
1928
1929 // D0 -> a1(1260)+ {rho(770)+ pi0[S]} pi-
1930 double SF_Ap_S_Vp1P = contract_11_0( contract_21_1( Proj1_3p, T1_PipPi01 ), T1_3pm );
1931 double SF_Ap_S_Vp2P = contract_11_0( contract_21_1( Proj1_3p, T1_PipPi02 ), T1_3pm );
1932
1933 amplitude += fitpara[0] * ( SF_Ap_S_Vp1P * RBW_a11260_p * GS_rho770_p1 +
1934 SF_Ap_S_Vp2P * RBW_a11260_p * GS_rho770_p2 );
1935
1936 // D0 -> a1(1260)+ {rho(770)+ pi0[D]} pi-
1937 double SF_Ap_D_Vp1P = contract_11_0( contract_21_1( T2_PipPi01Pi02, T1_PipPi01 ), T1_3pm );
1938 double SF_Ap_D_Vp2P = contract_11_0( contract_21_1( T2_PipPi02Pi01, T1_PipPi02 ), T1_3pm );
1939
1940 // cout<<SF_Ap_D_VP_1<<","<<SF_Ap_D_VP_2<<","<<SF_Ap_D_VP_3<<","<<SF_Ap_D_VP_4<<","<<endl;
1941 // cout<<"-------"<<endl;
1942 amplitude += fitpara[1] * ( SF_Ap_D_Vp1P * RBW_a11260_p * GS_rho770_p1 +
1943 SF_Ap_D_Vp2P * RBW_a11260_p * GS_rho770_p2 );
1944
1945 // D0 -> a1(1260)+ {f2(1270) pi+ [P]} pi-
1946 double SF_Ap_P_TP = contract_11_0(
1947 contract_21_1( contract_42_2( Proj2_3p, T2_Pi01Pi02 ), T1_Pi01Pi02Pip ), T1_3pm );
1948
1949 amplitude += fitpara[2] * ( SF_Ap_P_TP * RBW_a11260_p * RBW_f21270_00 );
1950
1951 // D0 -> a1(1260)+ {f0 pi+ [P]} pi-
1952 double SF_Ap_P_SP = contract_11_0( T1_3pm, T1_Pi01Pi02Pip );
1953
1954 amplitude += fitpara[3] * ( SF_Ap_P_SP * RBW_a11260_p * PiPiS_00_0 );
1955 amplitude += fitpara[4] * ( SF_Ap_P_SP * RBW_a11260_p * PiPiS_00_1 );
1956 amplitude += fitpara[5] * ( SF_Ap_P_SP * RBW_a11260_p * PiPiS_00_5 );
1957
1958 // D0 -> a1(1260)- { rho(770)- pi0 [S]} pi+
1959 double SF_Am_S_Vm1P = contract_11_0( contract_21_1( Proj1_3m, T1_PimPi01 ), T1_3mp );
1960 double SF_Am_S_Vm2P = contract_11_0( contract_21_1( Proj1_3m, T1_PimPi02 ), T1_3mp );
1961
1962 amplitude += fitpara[6] * fitpara[0] *
1963 ( SF_Am_S_Vm1P * RBW_a11260_m * GS_rho770_m1 +
1964 SF_Am_S_Vm2P * RBW_a11260_m * GS_rho770_m2 );
1965
1966 // D0 -> a1(1260)- {rho(770)- pi0[D]} pi+
1967 double SF_Am_D_Vm1P = contract_11_0( contract_21_1( T2_PimPi01Pi02, T1_PimPi01 ), T1_3mp );
1968 double SF_Am_D_Vm2P = contract_11_0( contract_21_1( T2_PimPi02Pi01, T1_PimPi02 ), T1_3mp );
1969
1970 amplitude += fitpara[6] * fitpara[1] *
1971 ( SF_Am_D_Vm1P * RBW_a11260_m * GS_rho770_m1 +
1972 SF_Am_D_Vm2P * RBW_a11260_m * GS_rho770_m2 );
1973
1974 // D0 -> a1(1260)- {f2(1270) pi- [P]} pi+
1975 double SF_Am_P_TP = contract_11_0(
1976 contract_21_1( contract_42_2( Proj2_3m, T2_Pi01Pi02 ), T1_Pi01Pi02Pim ), T1_3mp );
1977
1978 amplitude += fitpara[6] * fitpara[2] * ( SF_Am_P_TP * RBW_a11260_m * RBW_f21270_00 );
1979
1980 // D0 -> a1(1260)- {f0 pi- [P]} pi+
1981 double SF_Am_P_SP = contract_11_0( T1_3mp, T1_Pi01Pi02Pim );
1982
1983 amplitude += fitpara[6] * fitpara[3] * ( SF_Am_P_SP * RBW_a11260_m * PiPiS_00_0 );
1984 amplitude += fitpara[6] * fitpara[4] * ( SF_Am_P_SP * RBW_a11260_m * PiPiS_00_1 );
1985 amplitude += fitpara[6] * fitpara[5] * ( SF_Am_P_SP * RBW_a11260_m * PiPiS_00_5 );
1986
1987 // D -> a1(1260)0 pi0
1988 double SF_A01_S_Vp1P = contract_11_0( contract_21_1( Proj1_3z1, T1_PipPi01 ), T1_3z12 );
1989 double SF_A02_S_Vp2P = contract_11_0( contract_21_1( Proj1_3z2, T1_PipPi02 ), T1_3z21 );
1990 double SF_A01_S_Vm1P = contract_11_0( contract_21_1( Proj1_3z1, T1_PimPi01 ), T1_3z12 );
1991 double SF_A02_S_Vm2P = contract_11_0( contract_21_1( Proj1_3z2, T1_PimPi02 ), T1_3z21 );
1992 double SF_A01_S_VzP = contract_11_0( contract_21_1( Proj1_3z1, T1_PipPim ), T1_3z12 );
1993 double SF_A02_S_VzP = contract_11_0( contract_21_1( Proj1_3z2, T1_PipPim ), T1_3z21 );
1994
1995 amplitude += fitpara[7] * fitpara[0] *
1996 ( SF_A01_S_Vp1P * RBW_a11260_01 * GS_rho770_p1 +
1997 SF_A02_S_Vp2P * RBW_a11260_02 * GS_rho770_p2 +
1998 SF_A01_S_Vm1P * RBW_a11260_01 * GS_rho770_m1 +
1999 SF_A02_S_Vm2P * RBW_a11260_02 * GS_rho770_m2 );
2000
2001 double SF_A01_D_Vp1P = contract_11_0( contract_21_1( T2_PipPi01Pim, T1_PipPi01 ), T1_3z12 );
2002 double SF_A02_D_Vp2P = contract_11_0( contract_21_1( T2_PipPi02Pim, T1_PipPi02 ), T1_3z21 );
2003 double SF_A01_D_Vm1P = contract_11_0( contract_21_1( T2_PimPi01Pip, T1_PimPi01 ), T1_3z12 );
2004 double SF_A02_D_Vm2P = contract_11_0( contract_21_1( T2_PimPi02Pip, T1_PimPi02 ), T1_3z21 );
2005
2006 amplitude += fitpara[7] * fitpara[1] *
2007 ( SF_A01_D_Vp1P * RBW_a11260_01 * GS_rho770_p1 +
2008 SF_A02_D_Vp2P * RBW_a11260_02 * GS_rho770_p2 +
2009 SF_A01_D_Vm1P * RBW_a11260_01 * GS_rho770_m1 +
2010 SF_A02_D_Vm2P * RBW_a11260_02 * GS_rho770_m2 );
2011
2012 double SF_A01_P_TP = contract_11_0(
2013 contract_21_1( contract_42_2( Proj2_3z1, T2_PipPim ), T1_PipPimPi01 ), T1_3z12 );
2014 double SF_A02_P_TP = contract_11_0(
2015 contract_21_1( contract_42_2( Proj2_3z2, T2_PipPim ), T1_PipPimPi02 ), T1_3z21 );
2016
2017 amplitude += fitpara[7] * fitpara[2] * ( -1.0 ) *
2018 ( SF_A01_P_TP * RBW_a11260_01 * RBW_f21270_pm +
2019 SF_A02_P_TP * RBW_a11260_02 * RBW_f21270_pm );
2020
2021 double SF_A01_P_SP = contract_11_0( T1_3z12, T1_PipPimPi01 );
2022 double SF_A02_P_SP = contract_11_0( T1_3z21, T1_PipPimPi02 );
2023
2024 amplitude +=
2025 fitpara[7] * fitpara[3] * ( -1.0 ) *
2026 ( SF_A01_P_SP * RBW_a11260_01 * PiPiS_pm_0 + SF_A02_P_SP * RBW_a11260_02 * PiPiS_pm_0 );
2027 amplitude +=
2028 fitpara[7] * fitpara[4] * ( -1.0 ) *
2029 ( SF_A01_P_SP * RBW_a11260_01 * PiPiS_pm_1 + SF_A02_P_SP * RBW_a11260_02 * PiPiS_pm_1 );
2030 amplitude +=
2031 fitpara[7] * fitpara[5] * ( -1.0 ) *
2032 ( SF_A01_P_SP * RBW_a11260_01 * PiPiS_pm_5 + SF_A02_P_SP * RBW_a11260_02 * PiPiS_pm_5 );
2033
2034 // D0 -> a1(1420)+ {rho(770)+ pi0[S]} pi-
2035 // amplitude += fitpara[8]*(SF_Ap_S_Vp1P*RBW_a11420_p*GS_rho770_p1 +
2036 // SF_Ap_S_Vp2P*RBW_a11420_p*GS_rho770_p2); amplitude +=
2037 // fitpara[9]*(SF_A01_S_Vp1P*RBW_a11420_01*GS_rho770_p1 +
2038 // SF_A02_S_Vp2P*RBW_a11420_02*GS_rho770_p2 + SF_A01_S_Vm1P*RBW_a11420_01*GS_rho770_m1 +
2039 // SF_A02_S_Vm2P*RBW_a11420_02*GS_rho770_m2);
2040
2041 // D0 -> a1(1420)+ {f0 pi0+[S]} pi-
2042 amplitude += fitpara[8] * ( SF_Ap_P_SP * RBW_a11420_p * PiPiS_00_5 );
2043 amplitude += fitpara[9] * ( SF_Ap_P_SP * RBW_a11420_p * PiPiS_00_6 );
2044
2045 // D0 -> a2(1320)+ {rho+ pi0} pi-
2046 double SF_Tp_D_Vp1P = contract_22_0(
2047 contract_22_2(
2048 contract_31_2( contract_41_3( epsilon_uvmn, contract_21_1( Proj1_3p, T1_PipPi01 ) ),
2049 PipPi01Pi02 ),
2050 contract_42_2( Proj2_3p, T2_3pm ) ),
2051 T2_PipPi01Pi02 );
2052 double SF_Tp_D_Vp2P = contract_22_0(
2053 contract_22_2(
2054 contract_31_2( contract_41_3( epsilon_uvmn, contract_21_1( Proj1_3p, T1_PipPi02 ) ),
2055 PipPi01Pi02 ),
2056 contract_42_2( Proj2_3p, T2_3pm ) ),
2057 T2_PipPi02Pi01 );
2058
2059 amplitude += fitpara[10] * ( SF_Tp_D_Vp1P * GS_rho770_p1 * RBW_a21320_p +
2060 SF_Tp_D_Vp2P * GS_rho770_p2 * RBW_a21320_p );
2061
2062 // D0 -> a2(1320)- {rho- pi0} pi+
2063 double SF_Tm_D_Vm1P = contract_22_0(
2064 contract_22_2(
2065 contract_31_2( contract_41_3( epsilon_uvmn, contract_21_1( Proj1_3m, T1_PimPi01 ) ),
2066 PimPi01Pi02 ),
2067 contract_42_2( Proj2_3m, T2_3mp ) ),
2068 T2_PimPi01Pi02 );
2069 double SF_Tm_D_Vm2P = contract_22_0(
2070 contract_22_2(
2071 contract_31_2( contract_41_3( epsilon_uvmn, contract_21_1( Proj1_3m, T1_PimPi02 ) ),
2072 PimPi01Pi02 ),
2073 contract_42_2( Proj2_3m, T2_3mp ) ),
2074 T2_PimPi02Pi01 );
2075 amplitude += fitpara[11] * ( SF_Tm_D_Vm1P * GS_rho770_m1 * RBW_a21320_m +
2076 SF_Tm_D_Vm2P * GS_rho770_m2 * RBW_a21320_m );
2077
2078 // D0 -> h1(1170)0 {rho pi0} pi0
2079 amplitude += fitpara[12] * ( SF_A01_S_Vp1P * RBW_h11170_01 * GS_rho770_p1 +
2080 SF_A02_S_Vp2P * RBW_h11170_02 * GS_rho770_p2 -
2081 SF_A01_S_Vm1P * RBW_h11170_01 * GS_rho770_m1 -
2082 SF_A02_S_Vm2P * RBW_h11170_02 * GS_rho770_m2 -
2083 SF_A01_S_VzP * RBW_h11170_01 * GS_rho770_pm -
2084 SF_A02_S_VzP * RBW_h11170_02 * GS_rho770_pm );
2085
2086 // D0 -> pi(1300)- {rho(770)- pi0} pi+
2087 double SF_Pm_P_Vm1P = contract_11_0( T1_PimPi01, T1_PimPi01Pi02 );
2088 double SF_Pm_P_Vm2P = contract_11_0( T1_PimPi02, T1_PimPi02Pi01 );
2089
2090 amplitude += fitpara[13] * ( SF_Pm_P_Vm1P * GS_rho770_m1 * RBW_pi1300_m +
2091 SF_Pm_P_Vm2P * GS_rho770_m2 * RBW_pi1300_m );
2092
2093 // D0 -> pi(1300)- {f2 pi-} pi+
2094 // double SF_Pm_D_TP = contract_22_0(T2_Pi01Pi02,T2_Pi01Pi02Pim);
2095 // amplitude += fitpara[14]*fitpara[13]*(SF_Pm_D_TP*RBW_f21270_00*RBW_pi1300_m);
2096
2097 // D0 -> pi(1300)- {f0 pi-} pi+
2098 amplitude += fitpara[14] * fitpara[13] * ( RBW_pi1300_m * PiPiS_00_0 );
2099 // amplitude += fitpara[15]*fitpara[13]*(RBW_pi1300_m*PiPiS_00_5);
2100 amplitude += fitpara[15] * fitpara[13] * ( RBW_pi1300_m * PiPiS_00_6 );
2101
2102 // D0 -> pi(1300)+ {rho(770)+ pi0} pi-
2103 double SF_Pp_P_Vp1P = contract_11_0( T1_PipPi01, T1_PipPi01Pi02 );
2104 double SF_Pp_P_Vp2P = contract_11_0( T1_PipPi02, T1_PipPi02Pi01 );
2105
2106 amplitude += fitpara[16] * ( SF_Pp_P_Vp1P * GS_rho770_p1 * RBW_pi1300_p +
2107 SF_Pp_P_Vp2P * GS_rho770_p2 * RBW_pi1300_p );
2108
2109 // D0 -> pi(1300)+ {f2 pi-} pi+
2110 // double SF_Pp_D_TP = contract_22_0(T2_Pi01Pi02,T2_Pi01Pi02Pip);
2111 // amplitude += fitpara[14]*fitpara[17]*(SF_Pp_D_TP*RBW_f21270_00*RBW_pi1300_p);
2112
2113 // D0 -> pi(1300)+ {f0 pi+} pi-
2114 amplitude += fitpara[14] * fitpara[16] * ( RBW_pi1300_p * PiPiS_00_0 );
2115 // amplitude += fitpara[15]*fitpara[17]*(RBW_pi1300_p*PiPiS_00_5);
2116 amplitude += fitpara[15] * fitpara[16] * ( RBW_pi1300_p * PiPiS_00_6 );
2117
2118 // D0 -> pi(1300)0 {rho pi} pi0
2119 double SF_P01_P_Vp1P = contract_11_0( T1_PipPi01, T1_PipPi01Pim );
2120 double SF_P02_P_Vp2P = contract_11_0( T1_PipPi02, T1_PipPi02Pim );
2121 double SF_P01_P_Vm1P = contract_11_0( T1_PimPi01, T1_PimPi01Pip );
2122 double SF_P02_P_Vm2P = contract_11_0( T1_PimPi02, T1_PimPi02Pip );
2123
2124 amplitude += fitpara[17] * ( SF_P01_P_Vp1P * RBW_pi1300_01 * GS_rho770_p1 +
2125 SF_P02_P_Vp2P * RBW_pi1300_02 * GS_rho770_p2 +
2126 SF_P01_P_Vm1P * RBW_pi1300_01 * GS_rho770_m1 +
2127 SF_P02_P_Vm2P * RBW_pi1300_02 * GS_rho770_m2 );
2128
2129 // double SF_P01_D_TP = contract_22_0(T2_PipPim,T2_PipPimPi01);
2130 // double SF_P02_D_TP = contract_22_0(T2_PipPim,T2_PipPimPi02);
2131 // amplitude += fitpara[14]*fitpara[18]*(-1.0)*(SF_P01_D_TP*RBW_f21270_pm*RBW_pi1300_01 +
2132 // SF_P02_D_TP*RBW_f21270_pm*RBW_pi1300_02);
2133
2134 amplitude += fitpara[14] * fitpara[17] * ( -1.0 ) *
2135 ( RBW_pi1300_01 * PiPiS_pm_0 + RBW_pi1300_02 * PiPiS_pm_0 );
2136 // amplitude += fitpara[15]*fitpara[18]*(-1.0)*(RBW_pi1300_01*PiPiS_pm_5 +
2137 // RBW_pi1300_02*PiPiS_pm_5);
2138 amplitude += fitpara[15] * fitpara[17] * ( -1.0 ) *
2139 ( RBW_pi1300_01 * PiPiS_pm_6 + RBW_pi1300_02 * PiPiS_pm_6 );
2140
2141 // D0 -> rho+ rho- [S]
2142 double SF_Vp1Vm2_S = contract_11_0( T1_PipPi01, T1_PimPi02 );
2143 double SF_Vp2Vm1_S = contract_11_0( T1_PipPi02, T1_PimPi01 );
2144
2145 amplitude += fitpara[18] * ( SF_Vp1Vm2_S * GS_rho770_p1 * GS_rho770_m2 +
2146 SF_Vp2Vm1_S * GS_rho770_p2 * GS_rho770_m1 );
2147
2148 // D0 -> rho+ rho- [P]
2149 double SF_Vp1Vm2_P = contract_11_0(
2150 contract_21_1( contract_31_2( contract_41_3( epsilon_uvmn, T1_PipPi01 ), T1_PimPi02 ),
2151 T1_2p1m2 ),
2152 D0 );
2153 double SF_Vp2Vm1_P = contract_11_0(
2154 contract_21_1( contract_31_2( contract_41_3( epsilon_uvmn, T1_PipPi02 ), T1_PimPi01 ),
2155 T1_2p2m1 ),
2156 D0 );
2157
2158 amplitude += fitpara[19] * ( SF_Vp1Vm2_P * GS_rho770_p1 * GS_rho770_m2 +
2159 SF_Vp2Vm1_P * GS_rho770_p2 * GS_rho770_m1 );
2160
2161 // D0 -> rho+ rho- [D]
2162 double SF_Vp1Vm2_D = contract_11_0( contract_21_1( T2_2p1m2, T1_PipPi01 ), T1_PimPi02 );
2163 double SF_Vp2Vm1_D = contract_11_0( contract_21_1( T2_2p2m1, T1_PipPi02 ), T1_PimPi01 );
2164 amplitude += fitpara[20] * ( SF_Vp1Vm2_D * GS_rho770_p1 * GS_rho770_m2 +
2165 SF_Vp2Vm1_D * GS_rho770_p2 * GS_rho770_m1 );
2166
2167 // D0 -> rho0 (Pi0 Pi0)S
2168 double SF_VpmS12_P = contract_11_0( T1_PipPim, T1_2pm12 );
2169
2170 amplitude += fitpara[21] * ( SF_VpmS12_P * GS_rho770_pm * PiPiS_00_0 );
2171 amplitude += fitpara[22] * ( SF_VpmS12_P * GS_rho770_pm * PiPiS_00_5 );
2172 amplitude += fitpara[23] * ( SF_VpmS12_P * GS_rho770_pm * PiPiS_00_6 );
2173
2174 // D0 -> f0f0
2175 // 00
2176 amplitude += fitpara[24] * ( PiPiS_pm_0 * PiPiS_00_0 + PiPiS_00_0 * PiPiS_pm_0 );
2177 amplitude += fitpara[25] * ( PiPiS_pm_0 * PiPiS_00_1 + PiPiS_00_0 * PiPiS_pm_1 );
2178 amplitude += fitpara[26] * ( PiPiS_pm_1 * PiPiS_00_1 + PiPiS_00_1 * PiPiS_pm_1 );
2179 amplitude += fitpara[27] * ( PiPiS_pm_1 * PiPiS_00_5 + PiPiS_00_1 * PiPiS_pm_5 );
2180 amplitude += fitpara[28] * ( PiPiS_pm_5 * PiPiS_00_5 + PiPiS_00_5 * PiPiS_pm_5 );
2181 amplitude += fitpara[29] * ( PiPiS_pm_5 * PiPiS_00_6 + PiPiS_00_5 * PiPiS_pm_6 );
2182
2183 // D0 -> f2(1270) f0
2184 double SF_TpmS00_D = contract_22_0( T2_PipPim, T2_2pm12 );
2185 double SF_T00Spm_D = contract_22_0( T2_Pi01Pi02, T2_2pm12 );
2186
2187 amplitude += fitpara[30] * ( SF_TpmS00_D * RBW_f21270_pm * PiPiS_00_5 +
2188 SF_T00Spm_D * RBW_f21270_00 * PiPiS_pm_5 );
2189 amplitude += fitpara[31] * ( SF_TpmS00_D * RBW_f21270_pm * PiPiS_00_6 +
2190 SF_T00Spm_D * RBW_f21270_00 * PiPiS_pm_6 );
2191
2192 // D0 -> pi2(1670)0[f2(1270) pi0] pi0
2193 double SF_PT01_S_TP = contract_22_0( contract_42_2( Proj2_3z1, T2_PipPim ), T2_3z12 );
2194 double SF_PT02_S_TP = contract_22_0( contract_42_2( Proj2_3z2, T2_PipPim ), T2_3z21 );
2195
2196 amplitude += fitpara[32] * ( -1.0 ) *
2197 ( SF_PT01_S_TP * RBW_f21270_pm * RBW_pi21670_01 +
2198 SF_PT02_S_TP * RBW_f21270_pm * RBW_pi21670_02 );
2199
2200 // D -> 1--(rho pi) pi0
2201 double SF_V1_Vz = contract_11_0(
2202 contract_21_1( contract_31_2( contract_41_3( epsilon_uvmn, PipPimPi01 ), T1_PipPimPi01 ),
2203 T1_PipPim ),
2204 contract_21_1( Proj1_3z1, T1_3z12 ) );
2205 double SF_V1_Vp1 = contract_11_0(
2206 contract_21_1( contract_31_2( contract_41_3( epsilon_uvmn, PipPimPi01 ), T1_PipPi01Pim ),
2207 T1_PipPi01 ),
2208 contract_21_1( Proj1_3z1, T1_3z12 ) );
2209 double SF_V1_Vm1 = contract_11_0(
2210 contract_21_1( contract_31_2( contract_41_3( epsilon_uvmn, PipPimPi01 ), T1_PimPi01Pip ),
2211 T1_PimPi01 ),
2212 contract_21_1( Proj1_3z1, T1_3z12 ) );
2213
2214 double SF_V2_Vz = contract_11_0(
2215 contract_21_1( contract_31_2( contract_41_3( epsilon_uvmn, PipPimPi02 ), T1_PipPimPi02 ),
2216 T1_PipPim ),
2217 contract_21_1( Proj1_3z2, T1_3z21 ) );
2218 double SF_V2_Vp2 = contract_11_0(
2219 contract_21_1( contract_31_2( contract_41_3( epsilon_uvmn, PipPimPi02 ), T1_PipPi02Pim ),
2220 T1_PipPi02 ),
2221 contract_21_1( Proj1_3z2, T1_3z21 ) );
2222 double SF_V1_Vm2 = contract_11_0(
2223 contract_21_1( contract_31_2( contract_41_3( epsilon_uvmn, PipPimPi02 ), T1_PimPi02Pip ),
2224 T1_PimPi02 ),
2225 contract_21_1( Proj1_3z2, T1_3z21 ) );
2226
2227 // D0 -> omega(rho pi) pi0
2228 amplitude +=
2229 ( -1.0 ) * fitpara[33] *
2230 ( SF_V1_Vp1 * RBW_omega_01 * GS_rho770_p1 - SF_V1_Vz * RBW_omega_01 * GS_rho770_pm -
2231 SF_V1_Vm1 * RBW_omega_01 * GS_rho770_m1 + SF_V2_Vp2 * RBW_omega_02 * GS_rho770_p2 -
2232 SF_V2_Vz * RBW_omega_02 * GS_rho770_pm - SF_V1_Vm2 * RBW_omega_02 * GS_rho770_m2 );
2233
2234 // D0 -> phi(rho pi) pi0
2235 amplitude +=
2236 ( -1.0 ) * fitpara[34] *
2237 ( SF_V1_Vp1 * RBW_phi_01 * GS_rho770_p1 - SF_V1_Vz * RBW_phi_01 * GS_rho770_pm -
2238 SF_V1_Vm1 * RBW_phi_01 * GS_rho770_m1 + SF_V2_Vp2 * RBW_phi_02 * GS_rho770_p2 -
2239 SF_V2_Vz * RBW_phi_02 * GS_rho770_pm - SF_V1_Vm2 * RBW_phi_02 * GS_rho770_m2 );
2240
2241 return amplitude;
2242}

Referenced by QCMCFilter::findD0Decay().

◆ init()

void D0Topippim2pi0::init ( )

Definition at line 13 of file D0Topippim2pi0.cxx.

13 {
14
15 // std::cout << "D0Topippim2pi0 ==> Initialization !" << std::endl;
16
17 double mag[35], pha[35];
18 mag[0] = 100.0;
19 pha[0] = 0.0;
20 mag[1] = 7.95507;
21 pha[1] = -0.0687407;
22 mag[2] = 37.5559;
23 pha[2] = -1.74946;
24 mag[3] = 61.2172;
25 pha[3] = 2.98079;
26 mag[4] = 187.79;
27 pha[4] = 2.64471;
28 mag[5] = 385.474;
29 pha[5] = -0.137107;
30 mag[6] = 0.330788;
31 pha[6] = 0.268133;
32 mag[7] = 0.584175;
33 pha[7] = -2.89693;
34 mag[8] = 127.158;
35 pha[8] = -2.47773;
36 mag[9] = 339.914;
37 pha[9] = 2.22856;
38 mag[10] = 0.320888;
39 pha[10] = -2.6194;
40 mag[11] = 0.366283;
41 pha[11] = -0.26867;
42 mag[12] = 14.1344;
43 pha[12] = -0.41164;
44 mag[13] = 86.0865;
45 pha[13] = -2.49649;
46 mag[14] = 6.1541;
47 pha[14] = -1.18299;
48 mag[15] = 56.6067;
49 pha[15] = 0.142977;
50 mag[16] = 92.3073;
51 pha[16] = -2.15881;
52 mag[17] = 80.9453;
53 pha[17] = 0.825815;
54 mag[18] = 16.9555;
55 pha[18] = -2.98994;
56 mag[19] = 9.72524;
57 pha[19] = -1.39929;
58 mag[20] = 5.71448;
59 pha[20] = 0.271902;
60 mag[21] = 21.4195;
61 pha[21] = -1.23701;
62 mag[22] = 56.8867;
63 pha[22] = -0.385837;
64 mag[23] = 231.626;
65 pha[23] = 2.14842;
66 mag[24] = 2938.45;
67 pha[24] = -0.693491;
68 mag[25] = 7252.7;
69 pha[25] = 2.23659;
70 mag[26] = 5165.87;
71 pha[26] = 0.913557;
72 mag[27] = 11508.6;
73 pha[27] = -1.07187;
74 mag[28] = 2461.86;
75 pha[28] = 1.8709;
76 mag[29] = 8757.75;
77 pha[29] = 2.40756;
78 mag[30] = 19.7413;
79 pha[30] = -1.0753;
80 mag[31] = 66.3826;
81 pha[31] = 2.34666;
82 mag[32] = 11.2904;
83 pha[32] = -0.822345;
84 mag[33] = 2.04576;
85 pha[33] = -0.281429;
86 mag[34] = 0.57927;
87 pha[34] = 2.7182;
88
89 fitpara.clear();
90 for ( int i = 0; i < 35; i++ )
91 {
92 complex<double> ctemp( mag[i] * cos( pha[i] ), mag[i] * sin( pha[i] ) );
93 fitpara.push_back( ctemp );
94 }
95
96 g_uv.clear();
97 for ( int i = 0; i < 4; i++ )
98 {
99 for ( int j = 0; j < 4; j++ )
100 {
101 if ( i != j ) { g_uv.push_back( 0.0 ); }
102 else if ( i < 3 ) { g_uv.push_back( -1.0 ); }
103 else if ( i == 3 ) { g_uv.push_back( 1.0 ); }
104 }
105 }
106
107 epsilon_uvmn.clear();
108 for ( int i = 0; i < 4; i++ )
109 {
110 for ( int j = 0; j < 4; j++ )
111 {
112 for ( int k = 0; k < 4; k++ )
113 {
114 for ( int l = 0; l < 4; l++ )
115 {
116 if ( i == j || i == k || i == l || j == k || j == l || k == l )
117 { epsilon_uvmn.push_back( 0.0 ); }
118 else
119 {
120 if ( i == 0 && j == 1 && k == 2 && l == 3 ) epsilon_uvmn.push_back( 1.0 );
121 if ( i == 0 && j == 1 && k == 3 && l == 2 ) epsilon_uvmn.push_back( -1.0 );
122 if ( i == 0 && j == 2 && k == 1 && l == 3 ) epsilon_uvmn.push_back( -1.0 );
123 if ( i == 0 && j == 2 && k == 3 && l == 1 ) epsilon_uvmn.push_back( 1.0 );
124 if ( i == 0 && j == 3 && k == 1 && l == 2 ) epsilon_uvmn.push_back( 1.0 );
125 if ( i == 0 && j == 3 && k == 2 && l == 1 ) epsilon_uvmn.push_back( -1.0 );
126
127 if ( i == 1 && j == 0 && k == 2 && l == 3 ) epsilon_uvmn.push_back( -1.0 );
128 if ( i == 1 && j == 0 && k == 3 && l == 2 ) epsilon_uvmn.push_back( 1.0 );
129 if ( i == 1 && j == 2 && k == 0 && l == 3 ) epsilon_uvmn.push_back( 1.0 );
130 if ( i == 1 && j == 2 && k == 3 && l == 0 ) epsilon_uvmn.push_back( -1.0 );
131 if ( i == 1 && j == 3 && k == 0 && l == 2 ) epsilon_uvmn.push_back( -1.0 );
132 if ( i == 1 && j == 3 && k == 2 && l == 0 ) epsilon_uvmn.push_back( 1.0 );
133
134 if ( i == 2 && j == 0 && k == 1 && l == 3 ) epsilon_uvmn.push_back( 1.0 );
135 if ( i == 2 && j == 0 && k == 3 && l == 1 ) epsilon_uvmn.push_back( -1.0 );
136 if ( i == 2 && j == 1 && k == 0 && l == 3 ) epsilon_uvmn.push_back( -1.0 );
137 if ( i == 2 && j == 1 && k == 3 && l == 0 ) epsilon_uvmn.push_back( 1.0 );
138 if ( i == 2 && j == 3 && k == 0 && l == 1 ) epsilon_uvmn.push_back( 1.0 );
139 if ( i == 2 && j == 3 && k == 1 && l == 0 ) epsilon_uvmn.push_back( -1.0 );
140
141 if ( i == 3 && j == 0 && k == 1 && l == 2 ) epsilon_uvmn.push_back( -1.0 );
142 if ( i == 3 && j == 0 && k == 2 && l == 1 ) epsilon_uvmn.push_back( 1.0 );
143 if ( i == 3 && j == 1 && k == 0 && l == 2 ) epsilon_uvmn.push_back( 1.0 );
144 if ( i == 3 && j == 1 && k == 2 && l == 0 ) epsilon_uvmn.push_back( -1.0 );
145 if ( i == 3 && j == 2 && k == 0 && l == 1 ) epsilon_uvmn.push_back( -1.0 );
146 if ( i == 3 && j == 2 && k == 1 && l == 0 ) epsilon_uvmn.push_back( 1.0 );
147 }
148 }
149 }
150 }
151 }
152
153 _nd = 4;
154 math_pi = 3.1415926f;
155 mass_Pion = 0.13957f;
156
157 rRes = 3.0 * 0.197321;
158 rD = 5.0 * 0.197321;
159 m_Pi = 0.139570;
160 m_Pi0 = 0.134977;
161 m2_Pi = m_Pi * m_Pi;
162 m2_Pi0 = m_Pi0 * m_Pi0;
163
164 m0_rho7700 = 0.77526;
165 w0_rho7700 = 0.1478;
166
167 m0_rho770p = 0.77511;
168 w0_rho770p = 0.1491;
169
170 m0_rho1450 = 1.465;
171 w0_rho1450 = 0.400;
172
173 m0_f21270 = 1.2755;
174 w0_f21270 = 0.1867;
175
176 m0_a11260 = 1.1337;
177 g1_a11260 = 0.00335;
178 g2_a11260 = 0.0;
179
180 m0_pi1300 = 1.498;
181 w0_pi1300 = 0.590;
182
183 m0_a11420 = 1.411;
184 w0_a11420 = 0.161;
185
186 m0_a11640 = 1.655;
187 w0_a11640 = 0.254;
188
189 m0_a21320 = 1.3186;
190 w0_a21320 = 0.105;
191
192 m0_pi21670 = 1.6706;
193 w0_pi21670 = 0.258;
194
195 m0_pi11400 = 1.354;
196 w0_pi11400 = 0.330;
197
198 m0_h11170 = 1.166;
199 w0_h11170 = 0.375;
200
201 m0_omega = 0.78265;
202 w0_omega = 0.00849;
203
204 m0_phi = 1.019461;
205 w0_phi = 0.004249;
206
207 s0_prod = -5.0;
208}
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().

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