BOSS 7.1.3
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | List of all members
D0ToKSpipi2018 Class Reference

#include <D0ToKSpipi2018.h>

Public Member Functions

 D0ToKSpipi2018 ()
 
virtual ~D0ToKSpipi2018 ()
 
void init ()
 
complex< double > Amp_PFT (vector< double > k0l, vector< double > pip, vector< double > pim)
 

Detailed Description

Definition at line 9 of file D0ToKSpipi2018.h.

Constructor & Destructor Documentation

◆ D0ToKSpipi2018()

D0ToKSpipi2018::D0ToKSpipi2018 ( )
inline

Definition at line 13 of file D0ToKSpipi2018.h.

13{}

◆ ~D0ToKSpipi2018()

D0ToKSpipi2018::~D0ToKSpipi2018 ( )
virtual

Definition at line 23 of file D0ToKSpipi2018.cxx.

23{}

Member Function Documentation

◆ Amp_PFT()

complex< double > D0ToKSpipi2018::Amp_PFT ( vector< double > k0l,
vector< double > pip,
vector< double > pim )

Definition at line 94 of file D0ToKSpipi2018.cxx.

94 {
95 // Breit-Wigner lineshapes
96 vector<double> pD;pD.clear();
97 if(k0l.size()!=4||pip.size()!=4||pim.size()!=4)cout<<"ERROR in KSPIPI daughter 4 momentum"<<endl;
98 for(int i=0;i<k0l.size();i++){
99 pD.push_back(k0l[i] + pip[i] + pim[i]);
100 }
101
102 complex<double> DK2piRes0 = Resonance2(pD, pip, pim, ar[0], phir[0], width_R[0], mass_R[0], spin_R[0]); //ar, phir, width, mass, spin Rho770
103 complex<double> DK2piRes1 = Resonance2(pD, pip, pim, ar[1], phir[1], width_R[1], mass_R[1], spin_R[1]); //ar, phir, width, mass, spin Omega782
104 complex<double> DK2piRes2 = Resonance2(pD, pip, pim, ar[2], phir[2], width_R[2], mass_R[2], spin_R[2]); //ar, phir, width, mass, spin ftwo1270
105 complex<double> DK2piRes3 = Resonance2(pD, pip, pim, ar[3], phir[3], width_R[3], mass_R[3], spin_R[3]); //ar, phir, width, mass, spin Rho1450
106 complex<double> DK2piRes4 = Resonance2(pD, k0l, pim, ar[4], phir[4], width_R[4], mass_R[4], spin_R[4]); //ar, phir, width, mass, spin Kstar892-
107 complex<double> DK2piRes5 = Resonance2(pD, k0l, pim, ar[5], phir[5], width_R[5], mass_R[5], spin_R[5]); //ar, phir, width, mass, spin K2star1430-
108 complex<double> DK2piRes6 = Resonance2(pD, k0l, pim, ar[6], phir[6], width_R[6], mass_R[6], spin_R[6]); //ar, phir, width, mass, spin Kstar1680-
109 complex<double> DK2piRes7 = Resonance2(pD, k0l, pim, ar[7], phir[7], width_R[7], mass_R[7], spin_R[7]); //ar, phir, width, mass, spin Kstar1410-
110 complex<double> DK2piRes8 = Resonance2(pD, k0l, pip, ar[8], phir[8], width_R[8], mass_R[8], spin_R[8]); //ar, phir, width, mass, spin Kstar892+
111 complex<double> DK2piRes9 = Resonance2(pD, k0l, pip, ar[9], phir[9], width_R[9], mass_R[9], spin_R[9]); //ar, phir, width, mass, spin K2star1430+
112 complex<double> DK2piRes10 = Resonance2(pD, k0l, pip, ar[10], phir[10], width_R[10], mass_R[10], spin_R[10]); //ar, phir, width, mass, spin Kstar1410+
113 // K-matrix for pipi S-wave
114 complex<double> pipi_s_wave = K_matrix(pip, pim);
115 if(pipi_s_wave == complex<double>(9999., 9999.)) return 1e-20;
116 // LASS parametrization for Kpi S-wave
117 complex<double> kpi_s_wave = amplitude_LASS(k0l, pip, pim, "k0spim", ar[11], phir[11]*pi180inv);
118 complex<double> dcs_kpi_s_wave = amplitude_LASS(k0l, pip, pim, "k0spip", ar[12], phir[12]*pi180inv);
119
120 complex<double> _tmpAmp = DK2piRes0 + DK2piRes1 + DK2piRes2 + DK2piRes3 + pipi_s_wave;
121 //complex<double> TOT_PFT_AMP = DK2piRes0+ DK2piRes1+ DK2piRes2+ DK2piRes3+ DK2piRes4+ DK2piRes5+ DK2piRes6+ DK2piRes7+ DK2piRes8+ DK2piRes9+ DK2piRes10+ pipi_s_wave + kpi_s_wave+ dcs_kpi_s_wave ;
122 complex<double> TOT_PFT_AMP = _tmpAmp + DK2piRes4+ DK2piRes5+ DK2piRes6+ DK2piRes7+ DK2piRes8+ DK2piRes9+ DK2piRes10 + kpi_s_wave+ dcs_kpi_s_wave ;
123 // Coherent sum for pure-flavor-tagged amplitudes (PFT)
124 return TOT_PFT_AMP;
125}

◆ init()

void D0ToKSpipi2018::init ( )

Definition at line 25 of file D0ToKSpipi2018.cxx.

25 {
26 //std::cout << "D0ToKSpipi2018 ==> Initialization !" << std::endl;
27
28 _nd = 3;
29 tan2thetaC = (0.22650*0.22650)/(1.-(0.22650*0.22650)) ; //sin(theta_C) = 0.22650 +/- 0.00048
30 pi180inv = 1.0*3.1415926/180;
31 double Pi = 3.1415926;
32
33 mass_R[0]= 0.77155; width_R[0]= 0.13469; spin_R[0]= 1; ar[0]= 1; phir[0]= 0;
34 mass_R[1]= 0.78265; width_R[1]= 0.00849; spin_R[1]= 1; ar[1]= 0.038791; phir[1]= (180./Pi)*2.1073;
35 mass_R[2]= 1.27510; width_R[2]= 0.18420; spin_R[2]= 2; ar[2]= 1.42887; phir[2]= (180./Pi)*-0.633296;
36 mass_R[3]= 1.46500; width_R[3]= 0.40000; spin_R[3]= 1; ar[3]= 2.85131; phir[3]= (180./Pi)*1.7820801;
37 mass_R[4]= 0.89371; width_R[4]= 0.04719; spin_R[4]= 1; ar[4]= 1.72044; phir[4]= (180./Pi)*2.38835877;
38 mass_R[5]= 1.42560; width_R[5]= 0.09850; spin_R[5]= 2; ar[5]= 1.27268; phir[5]= (180./Pi)*-0.769095;
39 mass_R[6]= 1.71700; width_R[6]= 0.3220; spin_R[6]= 1; ar[6]= 3.307642; phir[6]= (180./Pi)*-2.062227;
40 mass_R[7]= 1.41400; width_R[7]= 0.2320; spin_R[7]= 1; ar[7]= 0.286927; phir[7]= (180./Pi)*1.7346186;
41 mass_R[8]= 0.89371; width_R[8]= 0.04719; spin_R[8]= 1; ar[8]= 0.1641792;phir[8]= (180./Pi)*-0.735903;
42 mass_R[9]= 1.42560; width_R[9]= 0.0985; spin_R[9]= 2; ar[9]= 0.1025736;phir[9]= (180./Pi)*-1.56397;
43 mass_R[10]= 1.41400; width_R[10]= 0.2320; spin_R[10]= 1; ar[10]= 0.2090326;phir[10]= (180./Pi)*2.6208986;
44 ////////////////////////////////////////
45 mass_R[11]= 1.42500; width_R[11]= 0.2700; spin_R[11]= 1; ar[11]= 2.36; phir[11]= 99.4;//not found
46 mass_R[12]= 1.42500; width_R[12]= 0.2700; spin_R[12]= 1; ar[12]= 0.11267; phir[12]= -162.3;//rad
47
48 beta[0] = complex<double>( 8.521486*cos( 1.195641 ), 8.521486*sin( 1.195641));//
49 beta[1] = complex<double>( 12.1895 *cos( 0.41802), 12.1895 *sin( 0.41802));
50 beta[2] = complex<double>(29.14616 *cos(-0.0018386 ), 29.14616 *sin(-0.0018386 ));
51 beta[3] = complex<double>(10.745735 *cos(-0.9057014 ), 10.745735 *sin(-0.9057014 ));
52 beta[4] = complex<double>(0., 0.);
53
54 fprod[0] = complex<double>(8.04427*cos(-2.19847), 8.04427*sin(-2.19847));
55 fprod[1] = complex<double>(26.2986*cos(-2.65853), 26.2986*sin(-2.65853));
56 fprod[2] = complex<double>(33.0349*cos(-1.62714), 33.0349*sin(-1.62714));
57 fprod[3] = complex<double>(26.1741*cos(-2.11891), 26.1741*sin(-2.11891));
58 fprod[4] = complex<double>(0., 0.);
59 //beta.push_back( complex<double>( 0.255303*cos( 47.8861 *pi180inv), 0.255303*sin( 47.8861 *pi180inv)) );
60 //beta.push_back( complex<double>(13.4446 *cos( -5.11127*pi180inv), 13.4446 *sin( -5.11127*pi180inv)) );
61 //beta.push_back( complex<double>(38.8496 *cos(-30.06 *pi180inv), 38.8496 *sin(-30.06 *pi180inv)) );
62 //beta.push_back( complex<double>(13.1086 *cos(-81.4148 *pi180inv), 13.1086 *sin(-81.4148 *pi180inv)) );
63 //beta.push_back( complex<double>(0., 0.) );
64
65 //fprod.push_back( complex<double>(5.08049*cos(-182.312*pi180inv), 5.08049*sin(-182.312*pi180inv)));
66 //fprod.push_back( complex<double>(17.2388*cos(-219.209*pi180inv), 17.2388*sin(-219.209*pi180inv)));
67 //fprod.push_back( complex<double>(19.0145*cos(-76.9884*pi180inv), 19.0145*sin(-76.9884*pi180inv)));
68 //fprod.push_back( complex<double>(11.9875*cos(-190.502*pi180inv), 11.9875*sin(-190.502*pi180inv)));
69 //fprod.push_back( complex<double>(0., 0.));
70
71 ma[0]= 0.651; g[0][0]= 0.22889; g[0][1]= -0.55377; g[0][2]= 0; g[0][3]= -0.39899; g[0][4]= -0.34639;
72 ma[1]= 1.20360; g[1][0]= 0.94128; g[1][1]= 0.55095; g[1][2]= 0; g[1][3]= 0.39065; g[1][4]= 0.31503;
73 ma[2]= 1.55817; g[2][0]= 0.36856; g[2][1]= 0.23888; g[2][2]= 0.55639; g[2][3]= 0.18340; g[2][4]= 0.18681;
74 ma[3]= 1.21000; g[3][0]= 0.33650; g[3][1]= 0.40907; g[3][2]= 0.85679; g[3][3]= 0.19906; g[3][4]= -0.00984;
75 ma[4]= 1.82206; g[4][0]= 0.18171; g[4][1]= -0.17558; g[4][2]= -0.79658; g[4][3]= -0.00355; g[4][4]= 0.22358;
76
77 // Hadronic parameters for tag modes: 0=no-specified, 1=Kpi, 2=Kpipi0, 3=K3pi
78 rd[0] = 0.0;
79 rd[1] = 0.0586;
80 rd[2] = 0.0440;
81 rd[3] = 0.0546;
82 deltad[0] = 0.0;
83 deltad[1] = 194.7*pi180inv;
84 deltad[2] = 196.0*pi180inv;
85 deltad[3] = 167.0*pi180inv;
86 Rf[0] = 0.0;
87 Rf[1] = 1.0;
88 Rf[2] = 0.78;
89 Rf[3] = 0.52;
90
91 return;
92}
sin
double sin(const BesAngle a)
Definition BesAngle.h:210
cos
double cos(const BesAngle a)
Definition BesAngle.h:213
The documentation for this class was generated from the following files: