1800 {
1801
1802 initPar();
1803
1804
1805
1806
1807
1808
1809
1810 createPropagator( ks0, pip, pim, pi0 );
1811 createSpinfactor( ks0, pip, pim, pi0 );
1812
1813
1814
1815
1816
1817 double f_isoConj = -1. * sqrt( 2.0 );
1818 double one = 1.0, minus_one = -1.0;
1819 addPartialWave( spinfactor["D2VV_ks0pim_pippi0_0"], propagator["kstarm_892"],
1820 propagator["rhop_770"], coefficient["<kstarm_892|rhop_770>_0_mag"],
1821 coefficient["<kstarm_892|rhop_770>_0_phase"] );
1822 addPartialWave( spinfactor["D2VV_ks0pim_pippi0_1"], propagator["kstarm_892"],
1823 propagator["rhop_770"], coefficient["<kstarm_892|rhop_770>_1_mag"],
1824 coefficient["<kstarm_892|rhop_770>_1_phase"] );
1825 addPartialWave( spinfactor["D2VV_ks0pim_pippi0_2"], propagator["kstarm_892"],
1826 propagator["rhop_770"], coefficient["<kstarm_892|rhop_770>_2_mag"],
1827 coefficient["<kstarm_892|rhop_770>_2_phase"] );
1828 addPartialWave( spinfactor["D2VV_ks0pi0_pippim_0"], propagator["kstar0_892"],
1829 propagator["rho0_770"], coefficient["<kstar0_892|rho0_770>_0_mag"],
1830 coefficient["<kstar0_892|rho0_770>_0_phase"] );
1831 addPartialWave( spinfactor["D2VV_ks0pi0_pippim_1"], propagator["kstar0_892"],
1832 propagator["rho0_770"], coefficient["<kstar0_892|rho0_770>_1_mag"],
1833 coefficient["<kstar0_892|rho0_770>_1_phase"] );
1834 addPartialWave( spinfactor["D2VV_ks0pi0_pippim_2"], propagator["kstar0_892"],
1835 propagator["rho0_770"], coefficient["<kstar0_892|rho0_770>_2_mag"],
1836 coefficient["<kstar0_892|rho0_770>_2_phase"] );
1837 complex<double> omega_ks0_propagator =
1838 propagator["rhop_770"] * spinfactor["D2VP_V2VP_pippimpi0_pippi0_1"] +
1839 minus_one * propagator["rho0_770"] * spinfactor["D2VP_V2VP_pippimpi0_pippim_1"] +
1840 propagator["rhom_770"] * spinfactor["D2VP_V2VP_pippimpi0_pimpi0_1"];
1841 addPartialWave( omega_ks0_propagator, propagator["omega3_782"],
1842 coefficient["<omega3_782|rho_770>_1_mag"],
1843 coefficient["<omega3_782|rho_770>_1_phase"] );
1844 complex<double> phi_ks0_propagator =
1845 propagator["rhop_770"] * spinfactor["D2VP_V2VP_pippimpi0_pippi0_1"] +
1846 minus_one * propagator["rho0_770"] * spinfactor["D2VP_V2VP_pippimpi0_pippim_1"] +
1847 propagator["rhom_770"] * spinfactor["D2VP_V2VP_pippimpi0_pimpi0_1"];
1848 addPartialWave( phi_ks0_propagator, propagator["phi_1020"],
1849 coefficient["<phi_1020|rho_770>_1_mag"],
1850 coefficient["<phi_1020|rho_770>_1_phase"] );
1851 addPartialWave( spinfactor["D2PP_P2VP_pippimpi0_pimpi0_1"], propagator["eta_547"],
1852 propagator["rhom_phsp"], coefficient["<eta_547|rhom_phsp>_1_mag"],
1853 coefficient["<eta_547|rhom_phsp>_1_phase"] );
1854 addPartialWave( spinfactor["D2PP_P2SP_pippimpi0_pimpi0_0"], propagator["eta_547"],
1855 propagator["rhom_phsp"], coefficient["<eta_547|rhom_phsp>_0_mag"],
1856 coefficient["<eta_547|rhom_phsp>_0_phase"] );
1857 addPartialWave( spinfactor["D2AP_A2VP_ks0pimpi0_pimpi0_0"], propagator["k1m_1270"],
1858 propagator["rhom_770"], coefficient["<k1m_1270|rhom_770>_0_mag"],
1859 coefficient["<k1m_1270|rhom_770>_0_phase"] );
1860 addPartialWave( spinfactor["D2AP_A2VP_ks0pimpi0_ks0pim_0"], propagator["k1m_1270"],
1861 propagator["kstarm_892"], spinfactor["D2AP_A2VP_ks0pimpi0_ks0pi0_0"],
1862 propagator["k1m_1270"], propagator["kstar0_892"], f_isoConj,
1863 coefficient["<k1m_1270|kstar_892>_0_mag"],
1864 coefficient["<k1m_1270|kstar_892>_0_phase"] );
1865 addPartialWave( spinfactor["D2AP_A2SP_ks0pimpi0_ks0pim_1"], propagator["k1m_1270"],
1866 propagator["k0starm_1430"], spinfactor["D2AP_A2SP_ks0pimpi0_ks0pi0_1"],
1867 propagator["k1m_1270"], propagator["k0star0_1430"], f_isoConj,
1868 coefficient["<k1m_1270|k0star_1430>_1_mag"],
1869 coefficient["<k1m_1270|k0star_1430>_1_phase"] );
1870 addPartialWave( spinfactor["D2AP_A2VP_ks0pippim_pippim_0"], propagator["k10_1270"],
1871 propagator["rho0_770"], coefficient["<k10_1270|rho0_770>_0_mag"],
1872 coefficient["<k10_1270|rho0_770>_0_phase"] );
1873 addPartialWave( spinfactor["D2AP_A2VP_ks0pippim_ks0pim_0"], propagator["k10_1270"],
1874 propagator["kstarm_892"], coefficient["<k10_1270|kstarm_892>_0_mag"],
1875 coefficient["<k10_1270|kstarm_892>_0_phase"] );
1876 addPartialWave( spinfactor["D2AP_A2SP_ks0pippim_ks0pim_1"], propagator["k10_1270"],
1877 propagator["k0starm_1430"], coefficient["<k10_1270|k0starm_1430>_1_mag"],
1878 coefficient["<k10_1270|k0starm_1430>_1_phase"] );
1879 addPartialWave( spinfactor["D2AP_A2VP_ks0pimpi0_ks0pim_0"], propagator["k1m_1400"],
1880 propagator["kstarm_892"], spinfactor["D2AP_A2VP_ks0pimpi0_ks0pi0_0"],
1881 propagator["k1m_1400"], propagator["kstar0_892"], f_isoConj,
1882 coefficient["<k1m_1400|kstar_892>_0_mag"],
1883 coefficient["<k1m_1400|kstar_892>_0_phase"] );
1884 addPartialWave( spinfactor["D2AP_A2VP_ks0pippim_ks0pim_0"], propagator["k10_1400"],
1885 propagator["kstarm_892"], coefficient["<k10_1400|kstarm_892>_0_mag"],
1886 coefficient["<k10_1400|kstarm_892>_0_phase"] );
1887 addPartialWave( spinfactor["D2VV_ks0pi0_pippim_1"], propagator["kstar02_1410"],
1888 propagator["rho0_770"], coefficient["<kstar02_1410|rho0_770>_1_mag"],
1889 coefficient["<kstar02_1410|rho0_770>_1_phase"] );
1890 addPartialWave( spinfactor["D2VV_ks0pim_pippi0_0"], propagator["kstarm2_1410"],
1891 propagator["rhop_770"], coefficient["<kstarm2_1410|rhop_770>_0_mag"],
1892 coefficient["<kstarm2_1410|rhop_770>_0_phase"] );
1893 addPartialWave( spinfactor["D2VV_ks0pim_pippi0_1"], propagator["kstarm2_1410"],
1894 propagator["rhop_770"], coefficient["<kstarm2_1410|rhop_770>_1_mag"],
1895 coefficient["<kstarm2_1410|rhop_770>_1_phase"] );
1896 addPartialWave( spinfactor["D2VV_ks0pim_pippi0_2"], propagator["kstarm2_1410"],
1897 propagator["rhop_770"], coefficient["<kstarm2_1410|rhop_770>_2_mag"],
1898 coefficient["<kstarm2_1410|rhop_770>_2_phase"] );
1899 addPartialWave( spinfactor["D2VP_V2VP_ks0pimpi0_ks0pim_1"], propagator["kstarm3_1410"],
1900 propagator["kstarm_892"], spinfactor["D2VP_V2VP_ks0pimpi0_ks0pi0_1"],
1901 propagator["kstarm3_1410"], propagator["kstar0_892"], f_isoConj,
1902 coefficient["<kstarm3_1410|kstar_892>_1_mag"],
1903 coefficient["<kstarm3_1410|kstar_892>_1_phase"] );
1904 addPartialWave( spinfactor["D2VP_V2VP_ks0pimpi0_pimpi0_1"], propagator["kstarm3_1410"],
1905 propagator["rhom_770"], coefficient["<kstarm3_1410|rhom_770>_1_mag"],
1906 coefficient["<kstarm3_1410|rhom_770>_1_phase"] );
1907 addPartialWave( spinfactor["D2PP_P2VP_ks0pimpi0_ks0pim_1"], propagator["km_1460"],
1908 propagator["kstarm_892"], spinfactor["D2PP_P2VP_ks0pimpi0_ks0pi0_1"],
1909 propagator["km_1460"], propagator["kstar0_892"], f_isoConj,
1910 coefficient["<km_1460|kstar_892>_1_mag"],
1911 coefficient["<km_1460|kstar_892>_1_phase"] );
1912 addPartialWave( spinfactor["D2PP_P2SP_ks0pimpi0_ks0pim_0"], propagator["km_1460"],
1913 propagator["k0starm_1430"], spinfactor["D2PP_P2SP_ks0pimpi0_ks0pi0_0"],
1914 propagator["km_1460"], propagator["k0star0_1430"], f_isoConj,
1915 coefficient["<km_1460|k0star_1430>_0_mag"],
1916 coefficient["<km_1460|k0star_1430>_0_phase"] );
1917 addPartialWave( spinfactor["D2PP_P2VP_ks0pippim_ks0pim_1"], propagator["k0_1460"],
1918 propagator["kstarm_892"], coefficient["<k0_1460|kstarm_892>_1_mag"],
1919 coefficient["<k0_1460|kstarm_892>_1_phase"] );
1920 addPartialWave( spinfactor["D2AP_A2VP_ks0pimpi0_ks0pim_0"], propagator["k1m_1650"],
1921 propagator["kstarm_892"], spinfactor["D2AP_A2VP_ks0pimpi0_ks0pi0_0"],
1922 propagator["k1m_1650"], propagator["kstar0_892"], f_isoConj,
1923 coefficient["<k1m_1650|kstar_892>_0_mag"],
1924 coefficient["<k1m_1650|kstar_892>_0_phase"] );
1925 addPartialWave( spinfactor["D2AP_A2VP_ks0pippim_ks0pim_0"], propagator["k10_1650"],
1926 propagator["kstarm_892"], coefficient["<k10_1650|kstarm_892>_0_mag"],
1927 coefficient["<k10_1650|kstarm_892>_0_phase"] );
1928 addPartialWave( spinfactor["D2AP_A2VP_ks0pimpi0_pimpi0_0"], propagator["k1m_1650"],
1929 propagator["rhom_770"], coefficient["<k1m_1650|rhom_770>_0_mag"],
1930 coefficient["<k1m_1650|rhom_770>_0_phase"] );
1931 addPartialWave( spinfactor["D2VP_V2VP_ks0pippim_ks0pim_1"], propagator["kstar03_1680"],
1932 propagator["kstarm_892"], coefficient["<kstar03_1680|kstarm_892>_1_mag"],
1933 coefficient["<kstar03_1680|kstarm_892>_1_phase"] );
1934 addPartialWave( spinfactor["D2VP_V2VP_ks0pimpi0_ks0pim_1"], propagator["kstarm3_1680"],
1935 propagator["kstarm_892"], spinfactor["D2VP_V2VP_ks0pimpi0_ks0pi0_1"],
1936 propagator["kstarm3_1680"], propagator["kstar0_892"], f_isoConj,
1937 coefficient["<kstarm3_1680|kstar_892>_1_mag"],
1938 coefficient["<kstarm3_1680|kstar_892>_1_phase"] );
1939 addPartialWave( spinfactor["D2SS_ks0pi0_pippim_0"], propagator["k0star0_1430"],
1940 propagator["rho0_phsp"], coefficient["<k0star0_1430|rho0_phsp>_0_mag"],
1941 coefficient["<k0star0_1430|rho0_phsp>_0_phase"] );
1942 addPartialWave( spinfactor["D2PP_P2VP_ks0pippi0_pippi0_1"], propagator["k1p_phsp"],
1943 propagator["rhop_770"], coefficient["<k1p_phsp|rhop_770>_1_mag"],
1944 coefficient["<k1p_phsp|rhop_770>_1_phase"] );
1945 addPartialWave( spinfactor["D2PP_P2VP_ks0pimpi0_pimpi0_1"], propagator["k1m_phsp"],
1946 propagator["rhom_770"], coefficient["<k1m_phsp|rhom_770>_1_mag"],
1947 coefficient["<k1m_phsp|rhom_770>_1_phase"] );
1948 addPartialWave( spinfactor["D2AP_A2VP_ks0pimpi0_pimpi0_2"], propagator["k1m_phsp"],
1949 propagator["rhom_770"], coefficient["<k1m_phsp|rhom_770>_2_mag"],
1950 coefficient["<k1m_phsp|rhom_770>_2_phase"] );
1951 addPartialWave( spinfactor["D2SS_ks0pim_pippi0_0"], propagator["k0starm_1430"],
1952 propagator["rhop_phsp"], coefficient["<k0starm_1430|rhop_phsp>_0_mag"],
1953 coefficient["<k0starm_1430|rhop_phsp>_0_phase"] );
1954 addPartialWave( spinfactor["D2PP_P2SP_ks0pippim_pippim_0"], propagator["k0_1460"],
1955 propagator["sigma_500"], coefficient["<k0_1460|sigma_500>_0_mag"],
1956 coefficient["<k0_1460|sigma_500>_0_phase"] );
1957 addPartialWave( spinfactor["D2AP_A2SP_ks0pippim_pippim_1"], propagator["k10_1270"],
1958 propagator["f0_980"], coefficient["<k10_1270|f0_980>_1_mag"],
1959 coefficient["<k10_1270|f0_980>_1_phase"] );
1960 addPartialWave( spinfactor["D2AP_A2SP_ks0pippim_pippim_1"], propagator["k10_1400"],
1961 propagator["f0_1370"], coefficient["<k10_1400|f0_1370>_1_mag"],
1962 coefficient["<k10_1400|f0_1370>_1_phase"] );
1963 addPartialWave( spinfactor["D2VS_pippim_ks0pi0_1"], propagator["rho0_770"],
1964 propagator["k0star0_1430"], coefficient["<rho0_770|k0star0_1430>_1_mag"],
1965 coefficient["<rho0_770|k0star0_1430>_1_phase"] );
1966 addPartialWave( spinfactor["D2TP_T2VP_ks0pimpi0_ks0pim_2"], propagator["k2starm3_1430"],
1967 propagator["kstarm_892"], spinfactor["D2TP_T2VP_ks0pimpi0_ks0pi0_2"],
1968 propagator["k2starm3_1430"], propagator["kstar0_892"], f_isoConj,
1969 coefficient["<k2starm3_1430|kstar_892>_2_mag"],
1970 coefficient["<k2starm3_1430|kstar_892>_2_phase"] );
1971 addPartialWave( spinfactor["D2TT_ks0pi0_pippim_1"], propagator["k2star02_1430"],
1972 propagator["f2_1270"], coefficient["<k2star02_1430|f2_1270>_1_mag"],
1973 coefficient["<k2star02_1430|f2_1270>_1_phase"] );
1974 addPartialWave( spinfactor["D2VS_ks0pi0_pippim_1"], propagator["kstar0_892"],
1975 propagator["f0_980"], coefficient["<kstar0_892|f0_980>_1_mag"],
1976 coefficient["<kstar0_892|f0_980>_1_phase"] );
1977 addPartialWave( spinfactor["D2SS_ks0pim_pippi0_0"], propagator["kstarm_phsp"],
1978 propagator["rhop_phsp"], coefficient["<kstarm_phsp|rhop_phsp>_0_mag"],
1979 coefficient["<kstarm_phsp|rhop_phsp>_0_phase"] );
1980 addPartialWave( spinfactor["D2VS_pippim_ks0pi0_1"], propagator["rho0_770"],
1981 propagator["kstar0_phsp"], coefficient["<rho0_770|kstar0_phsp>_1_mag"],
1982 coefficient["<rho0_770|kstar0_phsp>_1_phase"] );
1983
1984
1985
1986 return totalAmp;
1987}