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

#include <D0ToKLpipi.h>

Public Member Functions

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

Detailed Description

Definition at line 9 of file D0ToKLpipi.h.

Constructor & Destructor Documentation

◆ D0ToKLpipi()

D0ToKLpipi::D0ToKLpipi ( )
inline

Definition at line 13 of file D0ToKLpipi.h.

13{}

◆ ~D0ToKLpipi()

D0ToKLpipi::~D0ToKLpipi ( )
virtual

Definition at line 24 of file D0ToKLpipi.cxx.

24{}

Member Function Documentation

◆ Amp_PFT()

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

Definition at line 100 of file D0ToKLpipi.cxx.

100 {
101
102 vector<double> pD;pD.clear();
103 if(k0l.size()!=4||pip.size()!=4||pim.size()!=4)cout<<"ERROR in KSPIPI daughter 4 momentum"<<endl;
104 for(int i=0;i<k0l.size();i++){
105 pD.push_back(k0l[i] + pip[i] + pim[i]);
106 }
107
108 //EvtResonance2 DK2piRes0(pD, pip, pim, ar[0], phir[0], width_R[0], mass_R[0], spin[0]); //ar, phir, width, mass, spin //Rho770
109 //EvtResonance2 DK2piRes1(pD, pip, pim, ar[1], phir[1], width_R[1], mass_R[1], spin[1]); //ar, phir, width, mass, spin //Omega782
110 //EvtResonance2 DK2piRes2(pD, pip, pim, ar[2], phir[2], width_R[2], mass_R[2], spin[2]); //ar, phir, width, mass, spin //ftwo1270
111 //EvtResonance2 DK2piRes3(pD, pip, pim, ar[3], phir[3], width_R[3], mass_R[3], spin[3]); //ar, phir, width, mass, spin //Rho1450
112 //EvtResonance2 DK2piRes4(pD, k0l, pim, ar[4], phir[4], width_R[4], mass_R[4], spin[4]); //ar, phir, width, mass, spin //Kstar892minus
113 //EvtResonance2 DK2piRes5(pD, k0l, pim, ar[5], phir[5], width_R[5], mass_R[5], spin[5]); //ar, phir, width, mass, spin //K2star1430minus
114 //EvtResonance2 DK2piRes6(pD, k0l, pim, ar[6], phir[6], width_R[6], mass_R[6], spin[6]); //ar, phir, width, mass, spin //Kstar1680minus
115 //EvtResonance2 DK2piRes7(pD, k0l, pim, ar[7], phir[7], width_R[7], mass_R[7], spin[7]); //ar, phir, width, mass, spin //Kstar1410minus
116 //EvtResonance2 DK2piRes8(pD, k0l, pip, ar[8], phir[8], width_R[8], mass_R[8], spin[8]); //ar, phir, width, mass, spin //Kstar892plus
117 //EvtResonance2 DK2piRes9(pD, k0l, pip, ar[9], phir[9], width_R[9], mass_R[9], spin[9]); //ar, phir, width, mass, spin //K2star1430plus
118 //EvtResonance2 DK2piRes10(pD, k0l, pip, ar[10], phir[10], width_R[10], mass_R[10], spin[10]); //ar, phir, width, mass, spin //Kstar1410plus
119 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
120 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
121 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
122 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
123 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-
124 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-
125 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-
126 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-
127 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+
128 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+
129 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+
130
131 complex<double> pipi_s_wave = K_matrix(pip, pim) ;
132 if(pipi_s_wave == complex<double>(9999., 9999.)) return 1e-20 ;
133
134 complex<double> kpi_s_wave = amplitude_LASS(k0l, pip, pim, "k0lpim", ar[11], phir[11]*pi180inv) ;
135 //complex<double> kpi_s_wave_dcs = amplitude_LASS(k0l, pip, pim, "k0spip", ar[12], phir[12]*pi180inv) ; should be there but not observed yet GUESS
136
137 complex<double> TOT_PFT_AMP = DK2piRes0 * CP_mult[0]
138 + DK2piRes1 * CP_mult[1]
139 + DK2piRes2 * CP_mult[2]
140 + DK2piRes3 * CP_mult[3]
141 + DK2piRes4
142 + DK2piRes5
143 + DK2piRes6
144 + DK2piRes7
145 + DK2piRes8 * (-1.)
146 + DK2piRes9 * (-1.)
147 + DK2piRes10* (-1.)
148 + pipi_s_wave * CP_mult[4]
149 + kpi_s_wave ;
150
151
152 return TOT_PFT_AMP;
153
154}

Referenced by QCMCFilter::findD0Decay().

◆ init()

void D0ToKLpipi::init ( )

Definition at line 26 of file D0ToKLpipi.cxx.

26 {
27 //std::cout << "D0ToKLpipi ==> Initialization !" << std::endl;
28
29 _nd = 3;
30 tan2thetaC = (0.22650*0.22650)/(1.-(0.22650*0.22650)) ; //sin(theta_C) = 0.22650 +/- 0.00048
31 pi180inv = 1.0*3.1415926/180;
32
33 //spin[0] = 1 ; //Rho770
34 //spin[1] = 1 ; //Omega782
35 //spin[2] = 2 ; //ftwo1270
36 //spin[3] = 1 ; //Rho1450
37 //spin[4] = 1 ; //Kstar892minus
38 //spin[5] = 2 ; //K2star1430minus
39 //spin[6] = 1 ; //Kstar1680minus
40 //spin[7] = 1 ; //Kstar1410minus
41 //spin[8] = 1 ; //Kstar892plus
42 //spin[9] = 2 ; //K2star1430plus
43 //spin[10]= 1 ; //Kstar1410plus
44 mass_R[0]= 0.77526; width_R[0]= 0.14740; spin_R[0]= 1; ar[0]= 1 ; phir[0]= 0 ;
45 mass_R[1]= 0.78266; width_R[1]= 0.00868; spin_R[1]= 1; ar[1]= 0.03928; phir[1]= 106.1 ;
46 mass_R[2]= 1.27550; width_R[2]= 0.18670; spin_R[2]= 2; ar[2]= 1.302 ; phir[2]= -39.97;
47 mass_R[3]= 1.46500; width_R[3]= 0.40000; spin_R[3]= 1; ar[3]= 1.656 ; phir[3]= 113 ;
48 mass_R[4]= 0.89167; width_R[4]= 0.0514; spin_R[4]= 1; ar[4]= 1.843 ; phir[4]= 138 ;
49 mass_R[5]= 1.42730; width_R[5]= 0.10000; spin_R[5]= 2; ar[5]= 1.537 ; phir[5]= -49.45;
50 mass_R[6]= 1.71800; width_R[6]= 0.3220; spin_R[6]= 1; ar[6]= 2.251 ; phir[6]= -196.3;
51 mass_R[7]= 1.41400; width_R[7]= 0.2320; spin_R[7]= 1; ar[7]= 0.4733 ; phir[7]= 1659 ;
52 mass_R[8]= 0.89167; width_R[8]= 0.0514; spin_R[8]= 1; ar[8]= 0.1617 ; phir[8]= -36.25;
53 mass_R[9]= 1.42730; width_R[9]= 0.1000; spin_R[9]= 2; ar[9]= 0.2267 ; phir[9]= -84.64;
54 mass_R[10]= 1.41400; width_R[10]= 0.2320; spin_R[10]= 1; ar[10]= 0.2185 ; phir[10]= 73.12 ;
55 mass_R[11]= 1.42500; width_R[11]= 0.2700; spin_R[11]= 1; ar[11]= 2.396 ; phir[11]= 98.14 ;
56 // beta[kb] = EvtComplex(mag*cos(phase*pi180inv), mag*sin(phase*pi180inv)) ;
57 // fprod[kf] = EvtComplex(mag*cos(phase*pi180inv), mag*sin(phase*pi180inv)) ;
58 beta[0] = complex<double>(8.5 *cos(68.5 *pi180inv), 8.5 *sin(68.5 *pi180inv));
59 beta[1] = complex<double>(12.2*cos(24 *pi180inv), 12.2 *sin(24 *pi180inv));
60 beta[2] = complex<double>(29.2*cos(-0.1 *pi180inv), 29.2 *sin(-0.1 *pi180inv));
61 beta[3] = complex<double>(10.8*cos(-51.9*pi180inv), 10.8 *sin(-51.9*pi180inv));
62 beta[4] = complex<double>(0., 0.);
63
64 fprod[0] = complex<double>(8 *cos(-126 *pi180inv), 8 *sin(-126 *pi180inv));
65 fprod[1] = complex<double>(26.3*cos(-152.3*pi180inv), 26.3*sin(-152.3*pi180inv));
66 fprod[2] = complex<double>(33 *cos(-93.2 *pi180inv), 33 *sin(-93.2 *pi180inv));
67 fprod[3] = complex<double>(26.2*cos(-121.4*pi180inv), 26.2*sin(-121.4*pi180inv));
68 fprod[4] = complex<double>(0., 0.);
69
70 //CP_mult[w] = (EvtComplex(1., 0.) - 2.*tan2thetaC*EvtComplex(r*cos(delta), r*sin(delta))) ;
71 CP_mult[0] = complex<double>(1.,0.)-2.*tan2thetaC*complex<double>(1.851 *cos(-94.07 *pi180inv),1.851 *sin(-94.07 *pi180inv));
72 CP_mult[1] = complex<double>(1.,0.)-2.*tan2thetaC*complex<double>(6.332 *cos(2.103 *pi180inv),6.332 *sin(2.103 *pi180inv));
73 CP_mult[2] = complex<double>(1.,0.)-2.*tan2thetaC*complex<double>(3.229 *cos(-60.05 *pi180inv),3.229 *sin(-60.05 *pi180inv));
74 CP_mult[3] = complex<double>(1.,0.)-2.*tan2thetaC*complex<double>(12.75 *cos(73.62 *pi180inv),12.75 *sin(73.62 *pi180inv));
75 CP_mult[4] = complex<double>(1.,0.)-2.*tan2thetaC*complex<double>(0.7116*cos(-177.149*pi180inv),0.7116*sin(-177.149*pi180inv));
76
77 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;
78 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;
79 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;
80 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;
81 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;
82
83 // Hadronic parameters for tag modes: 0=no-specified, 1=Kpi, 2=Kpipi0, 3=K3pi
84 rd[0] = 0.0;
85 rd[1] = 0.0586;
86 rd[2] = 0.0440;
87 rd[3] = 0.0546;
88 deltad[0] = 0.0;
89 deltad[1] = 194.7*pi180inv;
90 deltad[2] = 196.0*pi180inv;
91 deltad[3] = 167.0*pi180inv;
92 Rf[0] = 0.0;
93 Rf[1] = 1.0;
94 Rf[2] = 0.78;
95 Rf[3] = 0.52;
96
97 return;
98}
sin
double sin(const BesAngle a)
Definition BesAngle.h:210
cos
double cos(const BesAngle a)
Definition BesAngle.h:213

Referenced by QCMCFilter::findD0Decay().

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