CGEM BOSS 6.6.5.g
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtBtoKpiCPiso Class Reference

#include <EvtBtoKpiCPiso.hh>

+ Inheritance diagram for EvtBtoKpiCPiso:

Public Member Functions

 EvtBtoKpiCPiso ()
 
virtual ~EvtBtoKpiCPiso ()
 
void getName (std::string &name)
 
EvtDecayBaseclone ()
 
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 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 EvtBtoKpiCPiso.hh.

Constructor & Destructor Documentation

◆ EvtBtoKpiCPiso()

EvtBtoKpiCPiso::EvtBtoKpiCPiso ( )
inline

Definition at line 32 of file EvtBtoKpiCPiso.hh.

32{}

Referenced by clone().

◆ ~EvtBtoKpiCPiso()

EvtBtoKpiCPiso::~EvtBtoKpiCPiso ( )
virtual

Definition at line 30 of file EvtBtoKpiCPiso.cc.

30{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtBtoKpiCPiso::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 39 of file EvtBtoKpiCPiso.cc.

39 {
40
41 return new EvtBtoKpiCPiso;
42
43}

◆ decay()

void EvtBtoKpiCPiso::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 105 of file EvtBtoKpiCPiso.cc.

105 {
106
107
109 //added by Lange Jan4,2000
110 static EvtId PI0=EvtPDL::getId("pi0");
111 static EvtId PIP=EvtPDL::getId("pi+");
112 static EvtId PIM=EvtPDL::getId("pi+");
113 static EvtId K0=EvtPDL::getId("K0");
114 static EvtId KB=EvtPDL::getId("anti-K0");
115 static EvtId KP=EvtPDL::getId("K+");
116 static EvtId KM=EvtPDL::getId("K-");
117
118
120 EvtComplex U, Ubar, V, Vbar, W, Wbar;
121
122 U = EvtComplex(getArg(2)*cos(getArg(3)),getArg(2)*sin(getArg(3)));
123 Ubar = EvtComplex(getArg(4)*cos(getArg(5)),getArg(4)*sin(getArg(5)));
124 V = EvtComplex(getArg(6)*cos(getArg(7)),getArg(6)*sin(getArg(7)));
125 Vbar = EvtComplex(getArg(8)*cos(getArg(9)),getArg(8)*sin(getArg(9)));
126 W = EvtComplex(getArg(10)*cos(getArg(11)),getArg(10)*sin(getArg(11)));
127 Wbar = EvtComplex(getArg(12)*cos(getArg(13)),getArg(12)*sin(getArg(13)));
128
129//depending on what combination of K pi we have, there will be different
130//A and Abar (only A given in comments!)
131
132if (((getDaug(0)==PI0) && (getDaug(1)==KP)) || ((getDaug(0)==KP) && (getDaug(1)==PI0))) {
133
134//pi0 K+, so U - W
135
136 A = U - W;
137 }
138
139if (((getDaug(0)==PI0) && (getDaug(1)==KM)) || ((getDaug(0)==KM) && (getDaug(1)==PI0))) {
140
141//pi0 K-, so Ubar - Wbar
142
143 A = Ubar - Wbar;
144 }
145
146 if (((getDaug(0)==PIP) && (getDaug(1)==K0)) || ((getDaug(0)==K0) && (getDaug(1)==PIP))) {
147
148//pi+ K0, so V + W
149
150 A = sqrt(2.0)*(V + W);
151 }
152
153 if (((getDaug(0)==PIM) && (getDaug(1)==KB)) || ((getDaug(0)==KB) && (getDaug(1)==PIM))) {
154
155//pi- K0bar, so Vbar + Wbar
156 A = sqrt(2.0)*(Vbar + Wbar);
157 }
158
159 if (((getDaug(0)==PI0) && (getDaug(1)==K0)) || ((getDaug(0)==K0) && (getDaug(1)==PI0))) {
160
161//pi0 K0, so U + W
162
163 A= U + W;
164 }
165
166 if (((getDaug(0)==PI0) && (getDaug(1)==KB)) || ((getDaug(0)==KB) && (getDaug(1)==PI0))) {
167
168 A = Ubar + Wbar;
169 }
170
171 if (((getDaug(0)==PIM) && (getDaug(1)==KP)) || ((getDaug(0)==KP) && (getDaug(1)==PIM))) {
172
173//pi- K+, so V - W
174
175 A = sqrt(2.0)*(V - W);
176 }
177
178 if (((getDaug(0)==PIP) && (getDaug(1)==KM)) || ((getDaug(0)==KM) && (getDaug(1)==PIP))) {
179
180 A = sqrt(2.0)*(Vbar - Wbar);
181 }
182
183 vertex(A);
184
185 return ;
186}
double sin(const BesAngle a)
Definition: BesAngle.h:210
double cos(const BesAngle a)
Definition: BesAngle.h:213
void vertex(const EvtComplex &amp)
Definition: EvtDecayAmp.hh:37
double getArg(int j)
EvtId * getDaugs()
Definition: EvtDecayBase.hh:65
EvtId getDaug(int i)
Definition: EvtDecayBase.hh:66
Definition: EvtId.hh:27
static EvtId getId(const std::string &name)
Definition: EvtPDL.cc:287
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)

◆ getName()

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

Implements EvtDecayBase.

Definition at line 32 of file EvtBtoKpiCPiso.cc.

32 {
33
34 model_name="BTOKPI_CP_ISO";
35
36}

◆ init()

void EvtBtoKpiCPiso::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 46 of file EvtBtoKpiCPiso.cc.

46 {
47
48 // check that there are 15 arguments
49 checkNArg(15);
50 checkNDaug(2);
51
53
56}
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)

◆ initProbMax()

void EvtBtoKpiCPiso::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 59 of file EvtBtoKpiCPiso.cc.

59 {
60
61//this might need to be revised
62
63 //added by Lange Jan4,2000
64 static EvtId PI0=EvtPDL::getId("pi0");
65 static EvtId PIP=EvtPDL::getId("pi+");
66 static EvtId PIM=EvtPDL::getId("pi+");
67 static EvtId K0=EvtPDL::getId("K0");
68 static EvtId KB=EvtPDL::getId("anti-K0");
69 static EvtId KP=EvtPDL::getId("K+");
70 static EvtId KM=EvtPDL::getId("K-");
71
72if (((getDaug(0)==PI0) && (getDaug(1)==KP)) || ((getDaug(0)==KP) && (getDaug(1)==PI0))) {
73 setProbMax(2.0*(getArg(2)*getArg(2) + getArg(10)*getArg(10)));
74 }
75
76if (((getDaug(0)==PI0) && (getDaug(1)==KM)) || ((getDaug(0)==KM) && (getDaug(1)==PI0))) {
77 setProbMax( 2.0*(getArg(4)*getArg(4) + getArg(12)*getArg(12)));
78 }
79
80if (((getDaug(0)==PIP) && (getDaug(1)==K0)) || ((getDaug(0)==K0) && (getDaug(1)==PIP))) {
81 setProbMax( 4.0*(getArg(6)*getArg(6) + getArg(10)*getArg(10)));
82 }
83
84if (((getDaug(0)==PIM) && (getDaug(1)==KB)) || ((getDaug(0)==KB) && (getDaug(1)==PIM))) {
85 setProbMax( 4.0*(getArg(8)*getArg(8) + getArg(12)*getArg(12)));
86}
87
88if (((getDaug(0)==PI0) && (getDaug(1)==K0)) || ((getDaug(0)==K0) && (getDaug(1)==PI0))) {
89 setProbMax( 2.0*(getArg(2)*getArg(2) + getArg(10)*getArg(10)));
90}
91
92if (((getDaug(0)==PI0) && (getDaug(1)==KB)) || ((getDaug(0)==KB) && (getDaug(1)==PI0))) {
93 setProbMax( 2.0*(getArg(4)*getArg(4) + getArg(12)*getArg(12)));
94}
95
96if (((getDaug(0)==PIM) && (getDaug(1)==KP)) || ((getDaug(0)==KP) && (getDaug(1)==PIM))) {
97 setProbMax( 4.0*(getArg(6)*getArg(6) + getArg(10)*getArg(10)));
98}
99
100if (((getDaug(0)==PIP) && (getDaug(1)==KM)) || ((getDaug(0)==KM) && (getDaug(1)==PIP))) {
101 setProbMax( 4.0*(getArg(8)*getArg(8) + getArg(12)*getArg(12)));
102}
103}
void setProbMax(double prbmx)

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