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

#include <EvtBTo4piCP.hh>

+ Inheritance diagram for EvtBTo4piCP:

Public Member Functions

 EvtBTo4piCP ()
 
virtual ~EvtBTo4piCP ()
 
void getName (std::string &name)
 
EvtDecayBase * clone ()
 
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 EvtDecayBase * clone ()=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 ()
 
EvtId * getDaugs ()
 
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
188 p->initializePhaseSpace(getNDaug(),getDaugs());
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}
sin
double sin(const BesAngle a)
Definition: BesAngle.h:210
cos
double cos(const BesAngle a)
Definition: BesAngle.h:213
EvtAmpA2
EvtComplex EvtAmpA2(const EvtVector4R &p4pi1, const EvtVector4R &p4pi2, const EvtVector4R &p4pi3, const EvtVector4R &p4pi4)
Definition: EvtBTo4piCP.cc:36
EvtAmpA1
EvtComplex EvtAmpA1(const EvtVector4R &p4pi1, const EvtVector4R &p4pi2, const EvtVector4R &p4pi3, const EvtVector4R &p4pi4)
Definition: EvtBTo4piCP.cc:95
t
TTree * t
Definition: binning.cxx:23
EvtCPUtil::OtherB
static void OtherB(EvtParticle *p, double &t, EvtId &otherb)
Definition: EvtCPUtil.cc:229
EvtConst::c
static const double c
Definition: EvtConst.hh:32
EvtDecayAmp::vertex
void vertex(const EvtComplex &amp)
Definition: EvtDecayAmp.hh:37
EvtDecayBase::getArg
double getArg(int j)
Definition: EvtDecayBase.cc:564
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition: EvtDecayBase.hh:65
EvtId
Definition: EvtId.hh:27
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition: EvtPDL.cc:287
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition: EvtParticle.cc:121
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:85
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition: EvtParticle.cc:1071

◆ 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
168 checkSpinParent(EvtSpinType::SCALAR);
169
170 checkSpinDaughter(0,EvtSpinType::SCALAR);
171 checkSpinDaughter(1,EvtSpinType::SCALAR);
172 checkSpinDaughter(2,EvtSpinType::SCALAR);
173 checkSpinDaughter(3,EvtSpinType::SCALAR);
174}
EvtDecayBase::checkSpinDaughter
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
Definition: EvtDecayBase.cc:533
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition: EvtDecayBase.cc:520
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition: EvtDecayBase.cc:504
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition: EvtDecayBase.cc:482
The documentation for this class was generated from the following files: