BOSS 8.0.0
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 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 31 of file EvtBTo4piCP.hh.

31{}

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 146 of file EvtBTo4piCP.cc.

146{ return new EvtBTo4piCP; }

◆ decay()

void EvtBTo4piCP::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 162 of file EvtBTo4piCP.cc.

162 {
163
164 // added by Lange Jan4,2000
165 static EvtId B0 = EvtPDL::getId( "B0" );
166 static EvtId B0B = EvtPDL::getId( "anti-B0" );
167
168 double t;
169 EvtId other_b;
170
171 EvtCPUtil::OtherB( p, t, other_b );
172
174 EvtVector4R mom1 = p->getDaug( 0 )->getP4();
175 EvtVector4R mom2 = p->getDaug( 1 )->getP4();
176 EvtVector4R mom3 = p->getDaug( 2 )->getP4();
177 EvtVector4R mom4 = p->getDaug( 3 )->getP4();
178
179 // double alpha=getArg(0);
180 // double dm=getArg(1);
181
182 EvtComplex amp;
183
184 EvtComplex A, Abar;
185
186 EvtComplex A_a1p, Abar_a1p, A_a2p, Abar_a2p;
187 EvtComplex A_a1m, Abar_a1m, A_a2m, Abar_a2m;
188
189 A_a1p = EvtComplex( getArg( 2 ) * cos( getArg( 3 ) ), getArg( 2 ) * sin( getArg( 3 ) ) );
190 Abar_a1p = EvtComplex( getArg( 4 ) * cos( getArg( 5 ) ), getArg( 4 ) * sin( getArg( 5 ) ) );
191
192 A_a2p = EvtComplex( getArg( 6 ) * cos( getArg( 7 ) ), getArg( 6 ) * sin( getArg( 7 ) ) );
193 Abar_a2p = EvtComplex( getArg( 8 ) * cos( getArg( 9 ) ), getArg( 8 ) * sin( getArg( 9 ) ) );
194
195 A_a1m = EvtComplex( getArg( 10 ) * cos( getArg( 11 ) ), getArg( 10 ) * sin( getArg( 11 ) ) );
196 Abar_a1m =
197 EvtComplex( getArg( 12 ) * cos( getArg( 13 ) ), getArg( 12 ) * sin( getArg( 13 ) ) );
198
199 A_a2m = EvtComplex( getArg( 14 ) * cos( getArg( 15 ) ), getArg( 14 ) * sin( getArg( 15 ) ) );
200 Abar_a2m =
201 EvtComplex( getArg( 16 ) * cos( getArg( 17 ) ), getArg( 16 ) * sin( getArg( 17 ) ) );
202
203 EvtComplex a2p_amp = EvtAmpA2( mom1, mom2, mom3, mom4 ) +
204 EvtAmpA2( mom1, mom4, mom3, mom2 ) +
205 EvtAmpA2( mom3, mom2, mom1, mom4 ) + EvtAmpA2( mom3, mom4, mom1, mom2 );
206
207 EvtComplex a2m_amp = EvtAmpA2( mom2, mom3, mom4, mom1 ) +
208 EvtAmpA2( mom2, mom1, mom4, mom3 ) +
209 EvtAmpA2( mom4, mom3, mom2, mom1 ) + EvtAmpA2( mom4, mom1, mom2, mom3 );
210
211 EvtComplex a1p_amp = EvtAmpA1( mom1, mom2, mom3, mom4 ) +
212 EvtAmpA1( mom1, mom4, mom3, mom2 ) +
213 EvtAmpA1( mom3, mom2, mom1, mom4 ) + EvtAmpA1( mom3, mom4, mom1, mom2 );
214
215 EvtComplex a1m_amp = EvtAmpA1( mom2, mom3, mom4, mom1 ) +
216 EvtAmpA1( mom2, mom1, mom4, mom3 ) +
217 EvtAmpA1( mom4, mom3, mom2, mom1 ) + EvtAmpA1( mom4, mom1, mom2, mom3 );
218
219 A = A_a2p * a2p_amp + A_a1p * a1p_amp + A_a2m * a2m_amp + A_a1m * a1m_amp;
220 Abar = Abar_a2p * a2p_amp + Abar_a1p * a1p_amp + Abar_a2m * a2m_amp + Abar_a1m * a1m_amp;
221
222 if ( other_b == B0B )
223 {
224 amp = A * cos( getArg( 1 ) * t / ( 2 * EvtConst::c ) ) +
225 EvtComplex( cos( -2.0 * getArg( 0 ) ), sin( -2.0 * getArg( 0 ) ) ) * getArg( 2 ) *
226 EvtComplex( 0.0, 1.0 ) * Abar * sin( getArg( 1 ) * t / ( 2 * EvtConst::c ) );
227 }
228 if ( other_b == B0 )
229 {
230 amp = A * EvtComplex( cos( 2.0 * getArg( 0 ) ), sin( 2.0 * getArg( 0 ) ) ) *
231 EvtComplex( 0.0, 1.0 ) * sin( getArg( 1 ) * t / ( 2 * EvtConst::c ) ) +
232 getArg( 2 ) * Abar * cos( getArg( 1 ) * t / ( 2 * EvtConst::c ) );
233 }
234
235 vertex( amp );
236
237 return;
238}
EvtComplex EvtAmpA2(const EvtVector4R &p4pi1, const EvtVector4R &p4pi2, const EvtVector4R &p4pi3, const EvtVector4R &p4pi4)
EvtComplex EvtAmpA1(const EvtVector4R &p4pi1, const EvtVector4R &p4pi2, const EvtVector4R &p4pi3, const EvtVector4R &p4pi4)
static void OtherB(EvtParticle *p, double &t, EvtId &otherb)
Definition EvtCPUtil.cc:225
static const double c
Definition EvtConst.hh:31
void vertex(const EvtComplex &amp)
double getArg(int j)
EvtId * getDaugs()
static EvtId getId(const std::string &name)
Definition EvtPDL.cc:272
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
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 144 of file EvtBTo4piCP.cc.

144{ model_name = "BTO4PI_CP"; }

◆ init()

void EvtBTo4piCP::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 148 of file EvtBTo4piCP.cc.

148 {
149
150 // check that there are 18 arguments
151 checkNArg( 18 );
152 checkNDaug( 4 );
153
155
160}
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: