CGEM BOSS 6.6.5.g
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | List of all members
EvtSemiLeptonicTensorAmp Class Reference

#include <EvtSemiLeptonicTensorAmp.hh>

+ Inheritance diagram for EvtSemiLeptonicTensorAmp:

Public Member Functions

void CalcAmp (EvtParticle *parent, EvtAmp &amp, EvtSemiLeptonicFF *FormFactors)
 
- Public Member Functions inherited from EvtSemiLeptonicAmp
virtual void CalcAmp (EvtParticle *parent, EvtAmp &amp, EvtSemiLeptonicFF *FormFactors)
 
double CalcMaxProb (EvtId parent, EvtId meson, EvtId lepton, EvtId nudaug, EvtSemiLeptonicFF *FormFactors)
 

Detailed Description

Definition at line 30 of file EvtSemiLeptonicTensorAmp.hh.

Member Function Documentation

◆ CalcAmp()

void EvtSemiLeptonicTensorAmp::CalcAmp ( EvtParticle *  parent,
EvtAmp &  amp,
EvtSemiLeptonicFF *  FormFactors 
)
virtual

Reimplemented from EvtSemiLeptonicAmp.

Definition at line 36 of file EvtSemiLeptonicTensorAmp.cc.

38 {
39 static EvtId EM=EvtPDL::getId("e-");
40 static EvtId MUM=EvtPDL::getId("mu-");
41 static EvtId TAUM=EvtPDL::getId("tau-");
42 static EvtId EP=EvtPDL::getId("e+");
43 static EvtId MUP=EvtPDL::getId("mu+");
44 static EvtId TAUP=EvtPDL::getId("tau+");
45
46 static EvtId D0=EvtPDL::getId("D0");
47 static EvtId D0B=EvtPDL::getId("anti-D0");
48 static EvtId DP=EvtPDL::getId("D+");
49 static EvtId DM=EvtPDL::getId("D-");
50 static EvtId DSM=EvtPDL::getId("D_s-");
51 static EvtId DSP=EvtPDL::getId("D_s+");
52
53 //Add the lepton and neutrino 4 momenta to find q2
54
55 EvtVector4R q = parent->getDaug(1)->getP4()
56 + parent->getDaug(2)->getP4();
57 double q2 = (q.mass2());
58
59 double hf,kf,bpf,bmf;
60
61 FormFactors->gettensorff(parent->getId(),
62 parent->getDaug(0)->getId(),
63 q2,
64 parent->getDaug(0)->mass(),
65 &hf,
66 &kf,
67 &bpf,
68 &bmf);
69
70
71 double costhl_flag = 1.0;
72
73 if(parent->getId()==D0||parent->getId()==D0B||
74 parent->getId()==DP||parent->getId()==DM) {
75 costhl_flag = -1.0;
76 }
77 if(parent->getId()==DSP||parent->getId()==DSM) {
78 costhl_flag = -1.0;
79 }
80 hf = hf * costhl_flag;
81
82 EvtVector4R p4b;
83 p4b.set(parent->mass(),0.0,0.0,0.0);
84
85 EvtVector4R p4meson = parent->getDaug(0)->getP4();
86
87 EvtVector4C l1,l2;
88
89 EvtId l_num = parent->getDaug(1)->getId();
90
91 EvtVector4C ep_meson_b[5];
92
93 ep_meson_b[0] = ((parent->getDaug(0)->epsTensorParent(0)).cont2(p4b)).conj();
94 ep_meson_b[1] = ((parent->getDaug(0)->epsTensorParent(1)).cont2(p4b)).conj();
95 ep_meson_b[2] = ((parent->getDaug(0)->epsTensorParent(2)).cont2(p4b)).conj();
96 ep_meson_b[3] = ((parent->getDaug(0)->epsTensorParent(3)).cont2(p4b)).conj();
97 ep_meson_b[4] = ((parent->getDaug(0)->epsTensorParent(4)).cont2(p4b)).conj();
98
99 EvtVector4R pp,pm;
100
101 pp=p4b+p4meson;
102 pm=p4b-p4meson;
103
104 //lange - October 31,2002 - try to lessen the mass dependence of probmax
105 double q2max = p4b.mass2() + p4meson.mass2() - 2.0*p4b.mass()*p4meson.mass();
106 double q2maxin=1.0/q2max;
107
108 EvtComplex ep_meson_bb[5];
109
110 ep_meson_bb[0]=ep_meson_b[0]*(p4b);
111 ep_meson_bb[1]=ep_meson_b[1]*(p4b);
112 ep_meson_bb[2]=ep_meson_b[2]*(p4b);
113 ep_meson_bb[3]=ep_meson_b[3]*(p4b);
114 ep_meson_bb[4]=ep_meson_b[4]*(p4b);
115
116
117 EvtVector4C tds0,tds1,tds2,tds3,tds4;
118
119 EvtTensor4C tds;
120 if (l_num==EM||l_num==MUM||l_num==TAUM){
121 EvtTensor4C tdual=EvtComplex(0.0,hf)*dual(directProd(pp,pm));
122 tds0=tdual.cont2(ep_meson_b[0])
123 -kf*ep_meson_b[0]
124 -bpf*ep_meson_bb[0]*pp-bmf*ep_meson_bb[0]*pm;
125 tds0*=q2maxin;
126
127 tds1=tdual.cont2(ep_meson_b[1])
128 -kf*ep_meson_b[1]
129 -bpf*ep_meson_bb[1]*pp-bmf*ep_meson_bb[1]*pm;
130 tds1*=q2maxin;
131
132 tds2=tdual.cont2(ep_meson_b[2])
133 -kf*ep_meson_b[2]
134 -bpf*ep_meson_bb[2]*pp-bmf*ep_meson_bb[2]*pm;
135 tds2*=q2maxin;
136
137 tds3=tdual.cont2(ep_meson_b[3])
138 -kf*ep_meson_b[3]
139 -bpf*ep_meson_bb[3]*pp-bmf*ep_meson_bb[3]*pm;
140 tds3*=q2maxin;
141
142 tds4=tdual.cont2(ep_meson_b[4])
143 -kf*ep_meson_b[4]
144 -bpf*ep_meson_bb[4]*pp-bmf*ep_meson_bb[4]*pm;
145 tds4*=q2maxin;
146
147
148 l1=EvtLeptonVACurrent(parent->getDaug(1)->spParent(0),
149 parent->getDaug(2)->spParentNeutrino());
150 l2=EvtLeptonVACurrent(parent->getDaug(1)->spParent(1),
151 parent->getDaug(2)->spParentNeutrino());
152 }
153 else{
154 if (l_num==EP||l_num==MUP||l_num==TAUP){
155 EvtTensor4C tdual=EvtComplex(0.0,-hf)*dual(directProd(pp,pm));
156 tds0=tdual.cont2(ep_meson_b[0])
157 -kf*ep_meson_b[0]
158 -bpf*ep_meson_bb[0]*pp-bmf*ep_meson_bb[0]*pm;
159 tds0*=q2maxin;
160
161 tds1=tdual.cont2(ep_meson_b[1])
162 -kf*ep_meson_b[1]
163 -bpf*ep_meson_bb[1]*pp-bmf*ep_meson_bb[1]*pm;
164 tds1*=q2maxin;
165
166 tds2=tdual.cont2(ep_meson_b[2])
167 -kf*ep_meson_b[2]
168 -bpf*ep_meson_bb[2]*pp-bmf*ep_meson_bb[2]*pm;
169 tds2*=q2maxin;
170
171 tds3=tdual.cont2(ep_meson_b[3])
172 -kf*ep_meson_b[3]
173 -bpf*ep_meson_bb[3]*pp-bmf*ep_meson_bb[3]*pm;
174 tds3*=q2maxin;
175
176 tds4=tdual.cont2(ep_meson_b[4])
177 -kf*ep_meson_b[4]
178 -bpf*ep_meson_bb[4]*pp-bmf*ep_meson_bb[4]*pm;
179 tds4*=q2maxin;
180
181 l1=EvtLeptonVACurrent(parent->getDaug(2)->spParentNeutrino(),
182 parent->getDaug(1)->spParent(0));
183 l2=EvtLeptonVACurrent(parent->getDaug(2)->spParentNeutrino(),
184 parent->getDaug(1)->spParent(1));
185 }
186 else{
187 report(ERROR,"EvtGen") << "dfnb89agngri wrong lepton number\n";
188 }
189 }
190
191 amp.vertex(0,0,l1*tds0);
192 amp.vertex(0,1,l2*tds0);
193
194 amp.vertex(1,0,l1*tds1);
195 amp.vertex(1,1,l2*tds1);
196
197 amp.vertex(2,0,l1*tds2);
198 amp.vertex(2,1,l2*tds2);
199
200 amp.vertex(3,0,l1*tds3);
201 amp.vertex(3,1,l2*tds3);
202
203 amp.vertex(4,0,l1*tds4);
204 amp.vertex(4,1,l2*tds4);
205
206 return;
207
208}
conj
Evt3Rank3C conj(const Evt3Rank3C &t2)
Definition: Evt3Rank3C.cc:175
directProd
Evt3Rank3C directProd(const EvtVector3C &c1, const EvtVector3C &c2, const EvtVector3C &c3)
EvtLeptonVACurrent
EvtVector4C EvtLeptonVACurrent(const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
Definition: EvtDiracSpinor.cc:187
report
ostream & report(Severity severity, const char *facility)
Definition: EvtReport.cc:36
ERROR
@ ERROR
Definition: EvtReport.hh:49
dual
EvtTensor4C dual(const EvtTensor4C &t2)
Definition: EvtTensor4C.cc:379
q
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_Chfin DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_sinw2 DOUBLE PRECISION m_GFermi DOUBLE PRECISION m_MfinMin DOUBLE PRECISION m_ta2 INTEGER m_out INTEGER m_KeyFSR INTEGER m_KeyQCD *COMMON c_Semalib $ !copy of input $ !CMS energy $ !beam mass $ !final mass $ !beam charge $ !final charge $ !smallest final mass $ !Z mass $ !Z width $ !EW mixing angle $ !Gmu Fermi $ alphaQED at q
Definition: KKsem.h:33
EvtAmp::vertex
void vertex(const EvtComplex &amp)
Definition: EvtAmp.cc:477
EvtId
Definition: EvtId.hh:27
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition: EvtPDL.cc:287
EvtParticle::epsTensorParent
virtual EvtTensor4C epsTensorParent(int i) const
Definition: EvtParticle.cc:655
EvtParticle::getId
EvtId getId() const
Definition: EvtParticle.cc:113
EvtParticle::spParentNeutrino
virtual EvtDiracSpinor spParentNeutrino() const
Definition: EvtParticle.cc:633
EvtParticle::spParent
virtual EvtDiracSpinor spParent(int) const
Definition: EvtParticle.cc:607
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition: EvtParticle.cc:121
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:85
EvtParticle::mass
double mass() const
Definition: EvtParticle.cc:127
EvtSemiLeptonicFF::gettensorff
virtual void gettensorff(EvtId parent, EvtId daught, double t, double mass, double *a1f, double *a2f, double *vf, double *a0f)
Definition: EvtSemiLeptonicFF.hh:37
EvtTensor4C::cont2
EvtVector4C cont2(const EvtVector4C &v4) const
Definition: EvtTensor4C.cc:474
EvtVector4R::mass
double mass() const
Definition: EvtVector4R.cc:39
EvtVector4R::mass2
double mass2() const
Definition: EvtVector4R.hh:116
EvtVector4R::set
void set(int i, double d)
Definition: EvtVector4R.hh:183
The documentation for this class was generated from the following files: