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

#include <EvtBto2piCPiso.hh>

+ Inheritance diagram for EvtBto2piCPiso:

Public Member Functions

 EvtBto2piCPiso ()
 
virtual ~EvtBto2piCPiso ()
 
void getName (std::string &name)
 
EvtDecayBase * clone ()
 
void init ()
 
void initProbMax ()
 
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 27 of file EvtBto2piCPiso.hh.

Constructor & Destructor Documentation

◆ EvtBto2piCPiso()

EvtBto2piCPiso::EvtBto2piCPiso ( )
inline

Definition at line 31 of file EvtBto2piCPiso.hh.

31{}

Referenced by clone().

◆ ~EvtBto2piCPiso()

EvtBto2piCPiso::~EvtBto2piCPiso ( )
virtual

Definition at line 34 of file EvtBto2piCPiso.cc.

34{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtBto2piCPiso::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 43 of file EvtBto2piCPiso.cc.

43 {
44
45 return new EvtBto2piCPiso;
46
47}

◆ decay()

void EvtBto2piCPiso::decay ( EvtParticle *  p)
virtual

Implements EvtDecayBase.

Definition at line 91 of file EvtBto2piCPiso.cc.

91 {
92
93 p->initializePhaseSpace(getNDaug(),getDaugs());
94
95 //added by Lange Jan4,2000
96 static EvtId B0=EvtPDL::getId("B0");
97 static EvtId B0B=EvtPDL::getId("anti-B0");
98 static EvtId PI0=EvtPDL::getId("pi0");
99 static EvtId PIP=EvtPDL::getId("pi+");
100 static EvtId PIM=EvtPDL::getId("pi-");
101
102 double t;
103 EvtId other_b;
104 int charged=0;
105
106//randomly generate the tag (B0 or B0B)
107
108 double tag = EvtRandom::Flat(0.0,1.0);
109 if (tag < 0.5) {
110
111 EvtCPUtil::OtherB(p,t,other_b,1.0);
112 other_b = B0;
113 }
114 else {
115
116 EvtCPUtil::OtherB(p,t,other_b,0.0);
117 other_b = B0B;
118 }
119
120 EvtComplex amp;
121
122 EvtComplex A,Abar;
123 EvtComplex A2, A2_bar, A0, A0_bar;
124
125 A2 = EvtComplex(getArg(2)*cos(getArg(3)),getArg(2)*sin(getArg(3)));
126 A2_bar = EvtComplex(getArg(4)*cos(getArg(5)),getArg(4)*sin(getArg(5)));
127
128 A0 = EvtComplex(getArg(6)*cos(getArg(7)),getArg(6)*sin(getArg(7)));
129 A0_bar = EvtComplex(getArg(8)*cos(getArg(9)),getArg(8)*sin(getArg(9)));
130
131 //depending on what combination of pi pi we have, there will be different
132 //A and Abar
133
134 if (((getDaugs()[0]==PIP) && (getDaugs()[1]==PI0)) || ((getDaugs()[0]==PI0) && (getDaugs()[1]==PIP))) {
135
136 //pi+ pi0, so just A_2
137
138 charged = 1;
139 A = 3.0*A2;
140
141 }
142
143 if (((getDaugs()[0]==PI0) && (getDaugs()[1]==PIM)) || ((getDaugs()[0]==PIM) && (getDaugs()[1]==PI0))) {
144
145 //pi- pi0, so just A2_bar
146
147 charged = 1;
148 A = 3.0*A2_bar;
149
150 }
151
152 if (((getDaugs()[0]==PIP) && (getDaugs()[1]==PIM)) || ((getDaugs()[0]==PIM) && (getDaugs()[1]==PIP))) {
153
154//pi+ pi-, so A_2 - A_0
155
156 charged = 0;
157 A=sqrt(2.0)*(A2-A0);
158 Abar=sqrt(2.0)*(A2_bar-A0_bar);
159 }
160
161 if ((getDaugs()[0]==PI0) && (getDaugs()[1]==PI0)) {
162
163 //pi0 pi0, so 2*A_2 + A_0
164
165 charged = 0;
166 A=2.0*A2 + A0;
167 Abar=2.0*A2_bar + A0_bar;
168 }
169
170 if(charged == 0) {
171
172 if (other_b==B0B){
173 amp=A*cos(getArg(1)*t/(2*EvtConst::c))+
174 EvtComplex(cos(-2.0*getArg(0)),sin(-2.0*getArg(0)))*
175 EvtComplex(0.0,1.0)*Abar*sin(getArg(1)*t/(2*EvtConst::c));
176 }
177 if (other_b==B0){
178 amp=A*EvtComplex(cos(2.0*getArg(0)),sin(2.0*getArg(0)))*
179 EvtComplex(0.0,1.0)*sin(getArg(1)*t/(2*EvtConst::c))+
180 Abar*cos(getArg(1)*t/(2*EvtConst::c));
181 }
182 }
183 else amp = A;
184
185 vertex(amp);
186
187 return ;
188}
sin
double sin(const BesAngle a)
Definition: BesAngle.h:210
cos
double cos(const BesAngle a)
Definition: BesAngle.h:213
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::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition: EvtParticle.cc:1069
EvtRandom::Flat
static double Flat()
Definition: EvtRandom.cc:73
t
int t()
Definition: t.c:1

◆ getName()

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

Implements EvtDecayBase.

Definition at line 36 of file EvtBto2piCPiso.cc.

36 {
37
38 model_name="BTO2PI_CP_ISO";
39
40}

◆ init()

void EvtBto2piCPiso::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 49 of file EvtBto2piCPiso.cc.

49 {
50
51 // check that there are 11 arguments
52
53 checkNArg(11);
54 checkNDaug(2);
55
56 checkSpinParent(EvtSpinType::SCALAR);
57
58 checkSpinDaughter(0,EvtSpinType::SCALAR);
59 checkSpinDaughter(1,EvtSpinType::SCALAR);
60
61}
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

◆ initProbMax()

void EvtBto2piCPiso::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 64 of file EvtBto2piCPiso.cc.

64 {
65
66 //added by Lange Jan4,2000
67 static EvtId PI0=EvtPDL::getId("pi0");
68 static EvtId PIP=EvtPDL::getId("pi+");
69 static EvtId PIM=EvtPDL::getId("pi-");
70
71//this may need to be revised
72
73if (((getDaugs()[0]==PIP) && (getDaugs()[1]==PIM)) || ((getDaugs()[0]==PIM) && (getDaugs()[1]==PIP))) {
74 setProbMax(4.0*(getArg(2)*getArg(2)+getArg(4)*getArg(4)));
75}
76
77if ((getDaugs()[0]==PI0) && (getDaugs()[1]==PI0)) {
78 setProbMax(2.0*(4.0*getArg(2)*getArg(2)+getArg(4)*getArg(4)));
79}
80
81if (((getDaugs()[0]==PIP) && (getDaugs()[1]==PI0)) || ((getDaugs()[0]==PI0) && (getDaugs()[1]==PIP))) {
82 setProbMax(6.0*getArg(2)*getArg(2));
83}
84
85if (((getDaugs()[0]==PI0) && (getDaugs()[1]==PIM)) || ((getDaugs()[0]==PIM) && (getDaugs()[1]==PI0))) {
86 setProbMax(6.0*getArg(4)*getArg(4));
87}
88
89}
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition: EvtDecayBase.cc:297
The documentation for this class was generated from the following files: