CGEM BOSS 6.6.5.h
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtSTSCP Class Reference

#include <EvtSTSCP.hh>

+ Inheritance diagram for EvtSTSCP:

Public Member Functions

 EvtSTSCP ()
 
virtual ~EvtSTSCP ()
 
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 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 EvtSTSCP.hh.

Constructor & Destructor Documentation

◆ EvtSTSCP()

EvtSTSCP::EvtSTSCP ( )
inline

Definition at line 32 of file EvtSTSCP.hh.

32{}

Referenced by clone().

◆ ~EvtSTSCP()

EvtSTSCP::~EvtSTSCP ( )
virtual

Definition at line 36 of file EvtSTSCP.cc.

36{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtSTSCP::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 53 of file EvtSTSCP.cc.

53 {
54
55 return new EvtSTSCP;
56
57}

◆ decay()

void EvtSTSCP::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 73 of file EvtSTSCP.cc.

73 {
74
75 //added by Lange Jan4,2000
76 static EvtId B0=EvtPDL::getId("B0");
77 static EvtId B0B=EvtPDL::getId("anti-B0");
78
79 double t;
80 EvtId other_b;
81
82 EvtParticle *t1;
84 t1 = p->getDaug(0);
85 EvtVector4R momt = t1->getP4();
86 EvtVector4R moms = p->getDaug(1)->getP4();
87 double masst = t1->mass();
88
89 EvtCPUtil::OtherB(p,t,other_b);
90
91 double m_parent = p->mass();
92
93 EvtComplex amp;
94
95 EvtComplex A,Abar;
96
98 Abar=EvtComplex(getArg(5)*cos(getArg(6)),getArg(5)*sin(getArg(6)));
99
100 if (other_b==B0B){
101 amp=A*cos(getArg(1)*t/(2*EvtConst::c))+
102 EvtComplex(cos(-2.0*getArg(0)),sin(-2.0*getArg(0)))*
103 getArg(2)*EvtComplex(0.0,1.0)*Abar*sin(getArg(1)*t/(2*EvtConst::c));
104 }
105 if (other_b==B0){
106 amp=A*EvtComplex(cos(2.0*getArg(0)),sin(2.0*getArg(0)))*
107 EvtComplex(0.0,1.0)*sin(getArg(1)*t/(2*EvtConst::c))+
108 getArg(2)*Abar*cos(getArg(1)*t/(2*EvtConst::c));
109 }
110
111 EvtVector4R p4_parent;
112
113 p4_parent=momt+moms;
114
115 double norm=masst*masst/(m_parent*momt.d3mag()*momt.d3mag());
116
117
118 vertex(0,amp*norm*t1->epsTensorParent(0).cont1(p4_parent)*p4_parent);
119 vertex(1,amp*norm*t1->epsTensorParent(1).cont1(p4_parent)*p4_parent);
120 vertex(2,amp*norm*t1->epsTensorParent(2).cont1(p4_parent)*p4_parent);
121 vertex(3,amp*norm*t1->epsTensorParent(3).cont1(p4_parent)*p4_parent);
122 vertex(4,amp*norm*t1->epsTensorParent(4).cont1(p4_parent)*p4_parent);
123
124 return ;
125}
double sin(const BesAngle a)
Definition BesAngle.h:210
double cos(const BesAngle a)
Definition BesAngle.h:213
static void OtherB(EvtParticle *p, double &t, EvtId &otherb)
Definition EvtCPUtil.cc:229
static const double c
Definition EvtConst.hh:32
void vertex(const EvtComplex &amp)
double getArg(int j)
EvtId * getDaugs()
Definition EvtId.hh:27
static EvtId getId(const std::string &name)
Definition EvtPDL.cc:287
virtual EvtTensor4C epsTensorParent(int i) const
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double mass() const
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
EvtVector4C cont1(const EvtVector4C &v4) const
double d3mag() const
int t()
Definition t.c:1

◆ getName()

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

Implements EvtDecayBase.

Definition at line 38 of file EvtSTSCP.cc.

38 {
39
40 model_name="STS_CP";
41
42}

◆ init()

void EvtSTSCP::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 59 of file EvtSTSCP.cc.

59 {
60
61 // check that there are 7 arguments
62 checkNArg(7);
63 checkNDaug(2);
64
66
69
70}
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 EvtSTSCP::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 44 of file EvtSTSCP.cc.

44 {
45
46 //This is probably not quite right, but it should do as a start...
47 //Anders
48
49 setProbMax(20*(getArg(3)*getArg(3)+getArg(5)*getArg(5)));
50
51}
void setProbMax(double prbmx)

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