BOSS 7.0.4
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtSVSNONCPEIGEN Class Reference

#include <EvtSVSNONCPEIGEN.hh>

+ Inheritance diagram for EvtSVSNONCPEIGEN:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtSVSNONCPEIGEN()

EvtSVSNONCPEIGEN::EvtSVSNONCPEIGEN ( )
inline

Definition at line 32 of file EvtSVSNONCPEIGEN.hh.

32{}

Referenced by clone().

◆ ~EvtSVSNONCPEIGEN()

EvtSVSNONCPEIGEN::~EvtSVSNONCPEIGEN ( )
virtual

Definition at line 35 of file EvtSVSNONCPEIGEN.cc.

35{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtSVSNONCPEIGEN::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 44 of file EvtSVSNONCPEIGEN.cc.

44 {
45
46 return new EvtSVSNONCPEIGEN;
47
48}

◆ decay()

void EvtSVSNONCPEIGEN::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 84 of file EvtSVSNONCPEIGEN.cc.

84 {
85
86 //added by Lange Jan4,2000
87 static EvtId B0=EvtPDL::getId("B0");
88 static EvtId B0B=EvtPDL::getId("anti-B0");
89
90 double t;
91 EvtId other_b;
92 EvtId daugs[2];
93
94 // MB: flip selects the final of the decay
95 int flip = ((p->getId() == B0) ? 0 : 1);
96 daugs[0]=getDaug(0);
97 daugs[1]=getDaug(1);
98 p->initializePhaseSpace(2, daugs);
99
100 EvtCPUtil::OtherB(p,t,other_b,0.5);
101
102 EvtComplex amp;
103 double dmt2 = (_dm * t) / (2 * EvtConst::c);
104 EvtComplex ePlusIPhi(cos(_phickm), sin(_phickm));
105 EvtComplex eMinusIPhi(cos(-_phickm), -sin(_phickm));
106
107 // flip == 0 : D-rho+
108 // flip == 1 : D+rho-
109
110 if (!flip) {
111 if (other_b==B0B){
112 // At t=0 we have a B0
113 amp = cos(dmt2)*_A_f + eMinusIPhi*EvtComplex(0.0,sin(dmt2))*_Abar_f;
114 }
115 if (other_b==B0){
116 // At t=0 we have a B0bar
117 amp = ePlusIPhi*EvtComplex(0.0,sin(dmt2))*_A_f + cos(dmt2)*_Abar_f;
118 }
119 }
120 else{
121 if (other_b==B0B){
122 // At t=0 we have a B0
123 amp = cos(dmt2)*_A_fbar + eMinusIPhi*EvtComplex(0.0,sin(dmt2))*_Abar_fbar;
124 }
125 if (other_b==B0){
126 // At t=0 we have a B0bar
127 amp = ePlusIPhi*EvtComplex(0.0,sin(dmt2))*_A_fbar + cos(dmt2)*_Abar_fbar;
128 }
129 }
130
131
132 EvtParticle *v;
133 v= p->getDaug(0);
134
135 EvtVector4R momv = p->getDaug(0)->getP4();
136 EvtVector4R moms = p->getDaug(1)->getP4();
137 EvtVector4R p4_parent=momv+moms;
138
139 double norm=momv.mass()/(momv.d3mag()*p->mass());
140
141 vertex(0,amp*norm*p4_parent*(v->epsParent(0)));
142 vertex(1,amp*norm*p4_parent*(v->epsParent(1)));
143 vertex(2,amp*norm*p4_parent*(v->epsParent(2)));
144
145 return ;
146}
double sin(const BesAngle a)
double cos(const BesAngle a)
**********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
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)
Definition: EvtDecayAmp.hh:37
EvtId getDaug(int i)
Definition: EvtDecayBase.hh:66
Definition: EvtId.hh:27
static EvtId getId(const std::string &name)
Definition: EvtPDL.cc:287
EvtId getId() const
Definition: EvtParticle.cc:113
const EvtVector4R & getP4() const
Definition: EvtParticle.cc:121
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:85
double mass() const
Definition: EvtParticle.cc:127
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
double mass() const
Definition: EvtVector4R.cc:39
double d3mag() const
Definition: EvtVector4R.cc:186
int t()
Definition: t.c:1

◆ getName()

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

Implements EvtDecayBase.

Definition at line 37 of file EvtSVSNONCPEIGEN.cc.

37 {
38
39 model_name="SVS_NONCPEIGEN";
40
41}

◆ init()

void EvtSVSNONCPEIGEN::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 50 of file EvtSVSNONCPEIGEN.cc.

50 {
51
52 // check that there are 11 arguments
53 checkNArg(11,7);
54 checkNDaug(2);
55
58
59 _dm=getArg(1);
60 _phickm=2*getArg(0)+getArg(2);
61
62 _A_f=EvtComplex(getArg(3)*cos(getArg(4)),getArg(3)*sin(getArg(4)));
63 _Abar_f=EvtComplex(getArg(5)*cos(getArg(6)),getArg(5)*sin(getArg(6)));
64
65 _A_fbar=_Abar_f;
66 _Abar_fbar=_A_f;
67
68 if (getNArg()==11){
69 _A_fbar=EvtComplex(getArg(7)*cos(getArg(8)),getArg(7)*sin(getArg(8)));
70 _Abar_fbar=EvtComplex(getArg(9)*cos(getArg(10)),getArg(9)*sin(getArg(10)));
71 }
72}
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
double getArg(int j)
void checkNDaug(int d1, int d2=-1)
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)

◆ initProbMax()

void EvtSVSNONCPEIGEN::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 74 of file EvtSVSNONCPEIGEN.cc.

74 {
75 double theProbMax =
76 abs(_A_f) * abs(_A_f) +
77 abs(_Abar_f) * abs(_Abar_f) +
78 abs(_A_fbar) * abs(_A_fbar) +
79 abs(_Abar_fbar) * abs(_Abar_fbar);
80
81 setProbMax(theProbMax);
82}
void setProbMax(double prbmx)

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