BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtChi1BB2 Class Reference

#include <EvtChi1BB2.hh>

Inheritance diagram for EvtChi1BB2:

Public Member Functions

 EvtChi1BB2 ()
virtual ~EvtChi1BB2 ()
void getName (std::string &name)
EvtDecayBase * clone ()
void decay (EvtParticle *p)
void init ()
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 ()
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 EvtChi1BB2.hh.

Constructor & Destructor Documentation

◆ EvtChi1BB2()

EvtChi1BB2::EvtChi1BB2 ( )
inline

Definition at line 31 of file EvtChi1BB2.hh.

31{}

Referenced by clone().

◆ ~EvtChi1BB2()

EvtChi1BB2::~EvtChi1BB2 ( )
virtual

Definition at line 36 of file EvtChi1BB2.cc.

36{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtChi1BB2::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 40 of file EvtChi1BB2.cc.

40{ return new EvtChi1BB2; }

◆ decay()

void EvtChi1BB2::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 60 of file EvtChi1BB2.cc.

60 {
61
62 p->initializePhaseSpace( getNDaug(), getDaugs() );
63
64 EvtParticle *v, *s1;
65 EvtVector4R pv, ps, ppr;
66
67 v = p->getDaug( 0 );
68 s1 = p->getDaug( 1 );
69 pv = v->getP4();
70 ps = s1->getP4();
71 ppr = p->getP4();
72
73 // Put phase space results into the daughters.
74 EvtHelSys angles( ppr, pv ); // using helicity sys. angles
75 double theta = angles.getHelAng( 1 );
76 double phi = angles.getHelAng( 2 );
77 double gamma = 0;
78
79 /* double alpha=getArg(0);
80 double F01=sqrt((1+alpha)/2);
81 double F00=sqrt((1-alpha)/4);
82 */
83 double F00 = getArg( 0 );
84 double F01 = getArg( 1 );
85 double F03 = getArg( 2 );
86 double F11 = getArg( 3 );
87 double F33 = -F00;
88 double F32 = -F01;
89 double F30 = -F03;
90 double F10 = F01;
91 double F23 = -F10;
92 double F22 = -F11;
93 // cout<<"********F00,F01,F03,F10,F11,F22,F23,F30,F32,F33="<<F00<<F01<<F03<<F10<<F11<<F22<<F23<<F30<<F32<<F33<<endl;
94 /*cout<<"********Djmn(1, 1, 0,phi,theta,gamma)="<<Djmn(1, 1, 0,phi,theta,gamma);
95 cout<<"********Djmn(1, 1,-1,phi,theta,gamma)="<<Djmn(1, 1,-1,phi,theta,gamma);
96 cout<<"********Djmn(1, 1, 1,phi,theta,gamma)="<<Djmn(1, 1, 1,phi,theta,gamma);
97 */
98 int j1, j2, j3;
99 for ( j1 = 0; j1 <= 2; )
100 {
101 for ( j2 = 0; j2 <= 3; )
102 {
103 for ( j3 = 0; j3 <= 3; )
104 {
105 vertex( j1, j2, j3, 0.0 );
106 j3++;
107 }
108 j2++;
109 }
110 j1++;
111 }
112
113 vertex( 0, 0, 0, Djmn( 1, 1, 0, phi, theta, gamma ) * F00 ); // Chi1 helicity =1
114 // corresponding index=0
115 vertex( 0, 0, 1, Djmn( 1, 1, -1, phi, theta, gamma ) * F01 );
116 vertex( 0, 0, 3, Djmn( 1, 1, 1, phi, theta, gamma ) * F03 );
117 vertex( 0, 1, 0, Djmn( 1, 1, 1, phi, theta, gamma ) * F10 );
118 vertex( 0, 1, 1, Djmn( 1, 1, 0, phi, theta, gamma ) * F11 );
119 vertex( 0, 2, 2, Djmn( 1, 1, 0, phi, theta, gamma ) * F22 );
120 vertex( 0, 2, 3, Djmn( 1, 1, -1, phi, theta, gamma ) * F23 );
121 vertex( 0, 3, 0, Djmn( 1, 1, -1, phi, theta, gamma ) * F30 );
122 vertex( 0, 3, 2, Djmn( 1, 1, 1, phi, theta, gamma ) * F32 );
123 vertex( 0, 3, 3, Djmn( 1, 1, 0, phi, theta, gamma ) * F33 );
124
125 vertex( 1, 0, 0, Djmn( 1, -1, 0, phi, theta, gamma ) * F00 ); // Chi1 helicity =-1
126 // corresponding index=1
127 vertex( 1, 0, 1, Djmn( 1, -1, -1, phi, theta, gamma ) * F01 );
128 vertex( 1, 0, 3, Djmn( 1, -1, 1, phi, theta, gamma ) * F03 );
129 vertex( 1, 1, 0, Djmn( 1, -1, 1, phi, theta, gamma ) * F10 );
130 vertex( 1, 1, 1, Djmn( 1, -1, 0, phi, theta, gamma ) * F11 );
131 vertex( 1, 2, 2, Djmn( 1, -1, 0, phi, theta, gamma ) * F22 );
132 vertex( 1, 2, 3, Djmn( 1, -1, -1, phi, theta, gamma ) * F23 );
133 vertex( 1, 3, 0, Djmn( 1, -1, -1, phi, theta, gamma ) * F30 );
134 vertex( 1, 3, 2, Djmn( 1, -1, 1, phi, theta, gamma ) * F32 );
135 vertex( 1, 3, 3, Djmn( 1, -1, 0, phi, theta, gamma ) * F33 );
136
137 vertex( 2, 0, 0, Djmn( 1, 0, 0, phi, theta, gamma ) * F00 ); // Chi1 helicity =0
138 // corresponding index=2
139 vertex( 2, 0, 1, Djmn( 1, 0, -1, phi, theta, gamma ) * F01 );
140 vertex( 2, 0, 3, Djmn( 1, 0, 1, phi, theta, gamma ) * F03 );
141 vertex( 2, 1, 0, Djmn( 1, 0, 1, phi, theta, gamma ) * F10 );
142 vertex( 2, 1, 1, Djmn( 1, 0, 0, phi, theta, gamma ) * F11 );
143 vertex( 2, 2, 2, Djmn( 1, 0, 0, phi, theta, gamma ) * F22 );
144 vertex( 2, 2, 3, Djmn( 1, 0, -1, phi, theta, gamma ) * F23 );
145 vertex( 2, 3, 0, Djmn( 1, 0, -1, phi, theta, gamma ) * F30 );
146 vertex( 2, 3, 2, Djmn( 1, 0, 1, phi, theta, gamma ) * F32 );
147 vertex( 2, 3, 3, Djmn( 1, 0, 0, phi, theta, gamma ) * F33 );
148
149 return;
150}
Djmn
EvtComplex Djmn(int j, int m, int n, double phi, double theta, double gamma)
Definition EvtHelSys.cc:165
v
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum v
Definition KarLud.h:35
EvtDecayAmp::vertex
void vertex(const EvtComplex &amp)
Definition EvtDecayAmp.hh:36
EvtDecayBase::getArg
double getArg(int j)
Definition EvtDecayBase.cc:520
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:68
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:108
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:80
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1042

◆ getName()

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

Implements EvtDecayBase.

Definition at line 38 of file EvtChi1BB2.cc.

38{ model_name = "Chi1BB2"; }

◆ init()

void EvtChi1BB2::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 42 of file EvtChi1BB2.cc.

42 {
43 // check that there are 1 arguments
44 checkNArg( 4 );
45 checkNDaug( 2 );
46 checkSpinParent( EvtSpinType::VECTOR );
47 checkSpinDaughter( 0, EvtSpinType::RARITASCHWINGER );
48 checkSpinDaughter( 1, EvtSpinType::RARITASCHWINGER );
49}
EvtDecayBase::checkSpinDaughter
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:492
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:479
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:465
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:449
The documentation for this class was generated from the following files: