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

#include <EvtSVPCP.hh>

+ Inheritance diagram for EvtSVPCP:

Public Member Functions

 EvtSVPCP ()
 
virtual ~EvtSVPCP ()
 
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 ()
 

Static Public Member Functions

static void SVPHel (EvtParticle *parent, EvtAmp &amp, EvtId n_v1, EvtId n_v2, const EvtComplex &hp, const EvtComplex &hm)
 
- 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)
 

Additional Inherited Members

- 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 36 of file EvtSVPCP.hh.

Constructor & Destructor Documentation

◆ EvtSVPCP()

EvtSVPCP::EvtSVPCP ( )
inline

Definition at line 40 of file EvtSVPCP.hh.

40{}

Referenced by clone().

◆ ~EvtSVPCP()

EvtSVPCP::~EvtSVPCP ( )
virtual

Definition at line 38 of file EvtSVPCP.cc.

38{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtSVPCP::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 47 of file EvtSVPCP.cc.

47 {
48
49 return new EvtSVPCP;
50
51}
EvtSVPCP::EvtSVPCP
EvtSVPCP()
Definition: EvtSVPCP.hh:40

◆ decay()

void EvtSVPCP::decay ( EvtParticle *  p)
virtual

Implements EvtDecayBase.

Definition at line 73 of file EvtSVPCP.cc.

73 {
74
75 static EvtId B0=EvtPDL::getId("B0");
76 static EvtId B0B=EvtPDL::getId("anti-B0");
77
78 double t;
79 EvtId other_b;
80
81 EvtCPUtil::OtherB(p,t,other_b);
82
83 EvtComplex G1P,G1M, G1_T_even, G1_T_odd;
84
85 double norm = getArg(3)*getArg(3)+getArg(5)*getArg(5);
86
87 G1P=EvtComplex(getArg(3)*cos(getArg(4))/norm,getArg(3)*sin(getArg(4))/norm);
88 G1M=EvtComplex(getArg(5)*cos(getArg(6))/norm,getArg(5)*sin(getArg(6))/norm);
89
90 G1_T_even = (G1P+G1M)/sqrt(2.0);
91 G1_T_odd = (G1P-G1M)/sqrt(2.0);
92
93 EvtComplex lambda_km =EvtComplex(cos(-2*getArg(0)),sin(-2*getArg(0)));
94
95 double cdmt=cos(getArg(1)*t/(2*EvtConst::c));
96 double sdmt=sin(getArg(1)*t/(2*EvtConst::c));
97
98 EvtComplex cG1_T_even,cG1_T_odd;
99
100 if (other_b==B0B){
101 cG1_T_even = G1_T_even*(cdmt+lambda_km*EvtComplex(0.0,getArg(2)*sdmt));
102 cG1_T_odd = G1_T_odd*(cdmt-lambda_km*EvtComplex(0.0,getArg(2)*sdmt));
103 }
104 if (other_b==B0){
105 cG1_T_even = G1_T_even*(cdmt+(1.0/lambda_km)*EvtComplex(0.0,getArg(2)*sdmt));
106 cG1_T_odd =-G1_T_odd*(cdmt-(1.0/lambda_km)*EvtComplex(0.0,getArg(2)*sdmt));
107 }
108
109 EvtComplex hp, hm, h0;
110
111 // This part is adopted from EvtSVVHel and since there is
112 // a photon that can not have helicity 0 this is put in by
113 // setting the h0 amplitude to 0.
114 hm=(cG1_T_even-cG1_T_odd)/sqrt(2.0);
115 hp=(cG1_T_even+cG1_T_odd)/sqrt(2.0);
116 h0=EvtComplex(0.0,0.0);
117
118 EvtParticle *v1,*ph;
119
120 p->initializePhaseSpace(getNDaug(),getDaugs());
121 v1 = p->getDaug(0);
122 ph = p->getDaug(1);
123 EvtVector4R momv1 = v1->getP4();
124 EvtVector4R momph = ph->getP4();
125
126 EvtTensor4C d,g;
127
128 g.setdiag(1.0,-1.0,-1.0,-1.0);
129
130 EvtVector4R v,vp;
131
132 v=momv1/momv1.d3mag();
133 vp=(momv1+momph)/(momv1+momph).mass();
134
135 d=((1.0/sqrt(3.0))*(h0-(hp+hm))*(-1.0/sqrt(3.0)))*g+
136 ((1.0/sqrt(2.0))*(hp-hm)*EvtComplex(0.0,1.0)*(sqrt(1.0/2.0)))*dual(directProd(v,vp))+
137 (sqrt(2.0/3.0)*(h0+0.5*(hp+hm))*sqrt(3.0/2.0))*(directProd(v,v)+(1.0/3.0)*g);
138
139 EvtVector4C ep0,ep1,ep2;
140
141 ep0=d.cont1(v1->eps(0).conj());
142 ep1=d.cont1(v1->eps(1).conj());
143 ep2=d.cont1(v1->eps(2).conj());
144
145 EvtVector4C ep20,ep21,ep22;
146
147 ep20=ph->epsParentPhoton(0).conj();
148 ep21=ph->epsParentPhoton(1).conj();
149
150 vertex(0,0,ep0*ep20);
151 vertex(0,1,ep0*ep21);
152
153 vertex(1,0,ep1*ep20);
154 vertex(1,1,ep1*ep21);
155
156 vertex(2,0,ep2*ep20);
157 vertex(2,1,ep2*ep21);
158
159
160 return ;
161
162}
sin
double sin(const BesAngle a)
Definition: BesAngle.h:210
cos
double cos(const BesAngle a)
Definition: BesAngle.h:213
directProd
Evt3Rank3C directProd(const EvtVector3C &c1, const EvtVector3C &c2, const EvtVector3C &c3)
dual
EvtTensor4C dual(const EvtTensor4C &t2)
Definition: EvtTensor4C.cc:379
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
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::epsParentPhoton
virtual EvtVector4C epsParentPhoton(int i)
Definition: EvtParticle.cc:585
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition: EvtParticle.cc:121
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:85
EvtParticle::eps
virtual EvtVector4C eps(int i) const
Definition: EvtParticle.cc:574
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition: EvtParticle.cc:1069
EvtTensor4C::setdiag
void setdiag(double t00, double t11, double t22, double t33)
Definition: EvtTensor4C.cc:218
EvtTensor4C::cont1
EvtVector4C cont1(const EvtVector4C &v4) const
Definition: EvtTensor4C.cc:461
EvtVector4C::conj
EvtVector4C conj() const
Definition: EvtVector4C.hh:206
EvtVector4R::mass
double mass() const
Definition: EvtVector4R.cc:39
EvtVector4R::d3mag
double d3mag() const
Definition: EvtVector4R.cc:186
t
int t()
Definition: t.c:1

◆ getName()

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

Implements EvtDecayBase.

Definition at line 40 of file EvtSVPCP.cc.

40 {
41
42 model_name="SVP_CP";
43
44}

◆ init()

void EvtSVPCP::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 60 of file EvtSVPCP.cc.

60 {
61
62 // check that there are 7 arguments
63 checkNArg(7);
64 checkNDaug(2);
65
66 checkSpinParent(EvtSpinType::SCALAR);
67
68 checkSpinDaughter(0,EvtSpinType::VECTOR);
69 checkSpinDaughter(1,EvtSpinType::PHOTON);
70
71}
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 EvtSVPCP::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 53 of file EvtSVPCP.cc.

53 {
54
55 setProbMax(2*(getArg(3)*getArg(3)+getArg(5)*getArg(5)));
56
57}
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition: EvtDecayBase.cc:297

◆ SVPHel()

static void EvtSVPCP::SVPHel ( EvtParticle *  parent,
EvtAmp &  amp,
EvtId  n_v1,
EvtId  n_v2,
const EvtComplex &  hp,
const EvtComplex &  hm 
)
static
The documentation for this class was generated from the following files: