CGEM BOSS 6.6.5.f
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtBTo4piCP Class Reference

#include <EvtBTo4piCP.hh>

+ Inheritance diagram for EvtBTo4piCP:

Public Member Functions

 EvtBTo4piCP ()
 
virtual ~EvtBTo4piCP ()
 
void getName (std::string &name)
 
EvtDecayBaseclone ()
 
void init ()
 
void decay (EvtParticle *p)
 
- Public Member Functions inherited from EvtDecayAmp
void makeDecay (EvtParticle *p)
 
void setWeight (double weight)
 
void vertex (const EvtComplex &amp)
 
void vertex (int i1, const EvtComplex &amp)
 
void vertex (int i1, int i2, const EvtComplex &amp)
 
void vertex (int i1, int i2, int i3, const EvtComplex &amp)
 
void vertex (int *i1, const EvtComplex &amp)
 
virtual ~EvtDecayAmp ()
 
- Public Member Functions inherited from EvtDecayBase
virtual void getName (std::string &name)=0
 
virtual void decay (EvtParticle *p)=0
 
virtual void makeDecay (EvtParticle *p)=0
 
virtual EvtDecayBaseclone ()=0
 
virtual void init ()
 
virtual void initProbMax ()
 
virtual std::string commandName ()
 
virtual void command (std::string cmd)
 
double getProbMax (double prob)
 
double resetProbMax (double prob)
 
 EvtDecayBase ()
 
virtual ~EvtDecayBase ()
 
virtual bool matchingDecay (const EvtDecayBase &other) const
 
EvtId getParentId ()
 
double getBranchingFraction ()
 
void disableCheckQ ()
 
void checkQ ()
 
int getNDaug ()
 
EvtIdgetDaugs ()
 
EvtId getDaug (int i)
 
int getNArg ()
 
int getPHOTOS ()
 
void setPHOTOS ()
 
void setVerbose ()
 
void setSummary ()
 
double * getArgs ()
 
std::string * getArgsStr ()
 
double getArg (int j)
 
std::string getArgStr (int j)
 
std::string getModelName ()
 
int getDSum ()
 
int summary ()
 
int verbose ()
 
void saveDecayInfo (EvtId ipar, int ndaug, EvtId *daug, int narg, std::vector< std::string > &args, std::string name, double brfr)
 
void printSummary ()
 
void setProbMax (double prbmx)
 
void noProbMax ()
 
void checkNArg (int a1, int a2=-1, int a3=-1, int a4=-1)
 
void checkNDaug (int d1, int d2=-1)
 
void checkSpinParent (EvtSpinType::spintype sp)
 
void checkSpinDaughter (int d1, EvtSpinType::spintype sp)
 
virtual int nRealDaughters ()
 

Additional Inherited Members

- Static Public Member Functions inherited from EvtDecayBase
static void findMasses (EvtParticle *p, int ndaugs, EvtId daugs[10], double masses[10])
 
static void findMass (EvtParticle *p)
 
static double findMaxMass (EvtParticle *p)
 
- Protected Member Functions inherited from EvtDecayBase
bool daugsDecayedByParentModel ()
 
- Protected Attributes inherited from EvtDecayAmp
EvtAmp _amp2
 
- Protected Attributes inherited from EvtDecayBase
bool _daugsDecayedByParentModel
 

Detailed Description

Definition at line 28 of file EvtBTo4piCP.hh.

Constructor & Destructor Documentation

◆ EvtBTo4piCP()

EvtBTo4piCP::EvtBTo4piCP ( )
inline

Definition at line 32 of file EvtBTo4piCP.hh.

32{}

Referenced by clone().

◆ ~EvtBTo4piCP()

EvtBTo4piCP::~EvtBTo4piCP ( )
virtual

Definition at line 33 of file EvtBTo4piCP.cc.

33{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtBTo4piCP::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 156 of file EvtBTo4piCP.cc.

156 {
157
158 return new EvtBTo4piCP;
159
160}

◆ decay()

void EvtBTo4piCP::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 176 of file EvtBTo4piCP.cc.

176 {
177
178 //added by Lange Jan4,2000
179 static EvtId B0=EvtPDL::getId("B0");
180 static EvtId B0B=EvtPDL::getId("anti-B0");
181
182
183 double t;
184 EvtId other_b;
185
186 EvtCPUtil::OtherB(p,t,other_b);
187
189 EvtVector4R mom1 = p->getDaug(0)->getP4();
190 EvtVector4R mom2 = p->getDaug(1)->getP4();
191 EvtVector4R mom3 = p->getDaug(2)->getP4();
192 EvtVector4R mom4 = p->getDaug(3)->getP4();
193
194 // double alpha=getArg(0);
195 //double dm=getArg(1);
196
197 EvtComplex amp;
198
199 EvtComplex A,Abar;
200
201
202 EvtComplex A_a1p,Abar_a1p,A_a2p,Abar_a2p;
203 EvtComplex A_a1m,Abar_a1m,A_a2m,Abar_a2m;
204
205 A_a1p=EvtComplex(getArg(2)*cos(getArg(3)),getArg(2)*sin(getArg(3)));
206 Abar_a1p=EvtComplex(getArg(4)*cos(getArg(5)),getArg(4)*sin(getArg(5)));
207
208 A_a2p=EvtComplex(getArg(6)*cos(getArg(7)),getArg(6)*sin(getArg(7)));
209 Abar_a2p=EvtComplex(getArg(8)*cos(getArg(9)),getArg(8)*sin(getArg(9)));
210
211 A_a1m=EvtComplex(getArg(10)*cos(getArg(11)),getArg(10)*sin(getArg(11)));
212 Abar_a1m=EvtComplex(getArg(12)*cos(getArg(13)),getArg(12)*sin(getArg(13)));
213
214 A_a2m=EvtComplex(getArg(14)*cos(getArg(15)),getArg(14)*sin(getArg(15)));
215 Abar_a2m=EvtComplex(getArg(16)*cos(getArg(17)),getArg(16)*sin(getArg(17)));
216
217 EvtComplex a2p_amp=EvtAmpA2(mom1,mom2,mom3,mom4)+
218 EvtAmpA2(mom1,mom4,mom3,mom2)+
219 EvtAmpA2(mom3,mom2,mom1,mom4)+
220 EvtAmpA2(mom3,mom4,mom1,mom2);
221
222 EvtComplex a2m_amp=EvtAmpA2(mom2,mom3,mom4,mom1)+
223 EvtAmpA2(mom2,mom1,mom4,mom3)+
224 EvtAmpA2(mom4,mom3,mom2,mom1)+
225 EvtAmpA2(mom4,mom1,mom2,mom3);
226
227 EvtComplex a1p_amp=EvtAmpA1(mom1,mom2,mom3,mom4)+
228 EvtAmpA1(mom1,mom4,mom3,mom2)+
229 EvtAmpA1(mom3,mom2,mom1,mom4)+
230 EvtAmpA1(mom3,mom4,mom1,mom2);
231
232 EvtComplex a1m_amp=EvtAmpA1(mom2,mom3,mom4,mom1)+
233 EvtAmpA1(mom2,mom1,mom4,mom3)+
234 EvtAmpA1(mom4,mom3,mom2,mom1)+
235 EvtAmpA1(mom4,mom1,mom2,mom3);
236
237
238 A=A_a2p*a2p_amp+A_a1p*a1p_amp+
239 A_a2m*a2m_amp+A_a1m*a1m_amp;
240 Abar=Abar_a2p*a2p_amp+Abar_a1p*a1p_amp+
241 Abar_a2m*a2m_amp+Abar_a1m*a1m_amp;
242
243
244 if (other_b==B0B){
245 amp=A*cos(getArg(1)*t/(2*EvtConst::c))+
246 EvtComplex(cos(-2.0*getArg(0)),sin(-2.0*getArg(0)))*
247 getArg(2)*EvtComplex(0.0,1.0)*Abar*sin(getArg(1)*t/(2*EvtConst::c));
248 }
249 if (other_b==B0){
250 amp=A*EvtComplex(cos(2.0*getArg(0)),sin(2.0*getArg(0)))*
251 EvtComplex(0.0,1.0)*sin(getArg(1)*t/(2*EvtConst::c))+
252 getArg(2)*Abar*cos(getArg(1)*t/(2*EvtConst::c));
253 }
254
255 vertex(amp);
256
257 return ;
258}
EvtComplex EvtAmpA2(const EvtVector4R &p4pi1, const EvtVector4R &p4pi2, const EvtVector4R &p4pi3, const EvtVector4R &p4pi4)
Definition: EvtBTo4piCP.cc:36
EvtComplex EvtAmpA1(const EvtVector4R &p4pi1, const EvtVector4R &p4pi2, const EvtVector4R &p4pi3, const EvtVector4R &p4pi4)
Definition: EvtBTo4piCP.cc:95
double sin(const BesAngle a)
double cos(const BesAngle a)
static void OtherB(EvtParticle *p, double &t, EvtId &otherb)
Definition: EvtCPUtil.cc:229
static const double c
Definition: EvtConst.hh:32
void vertex(const EvtComplex &amp)
Definition: EvtDecayAmp.hh:37
double getArg(int j)
EvtId * getDaugs()
Definition: EvtDecayBase.hh:65
Definition: EvtId.hh:27
static EvtId getId(const std::string &name)
Definition: EvtPDL.cc:287
const EvtVector4R & getP4() const
Definition: EvtParticle.cc:121
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:85
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
int t()
Definition: t.c:1

◆ getName()

void EvtBTo4piCP::getName ( std::string &  name)
virtual

Implements EvtDecayBase.

Definition at line 149 of file EvtBTo4piCP.cc.

149 {
150
151 model_name="BTO4PI_CP";
152
153}

◆ init()

void EvtBTo4piCP::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 162 of file EvtBTo4piCP.cc.

162 {
163
164 // check that there are 18 arguments
165 checkNArg(18);
166 checkNDaug(4);
167
169
174}
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
void checkSpinParent(EvtSpinType::spintype sp)
void checkNDaug(int d1, int d2=-1)
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)

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