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

#include <EvtHypWK.hh>

Inheritance diagram for EvtHypWK:

Public Member Functions

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

Constructor & Destructor Documentation

◆ EvtHypWK()

EvtHypWK::EvtHypWK ( )
inline

Definition at line 31 of file EvtHypWK.hh.

31{}

Referenced by clone().

◆ ~EvtHypWK()

EvtHypWK::~EvtHypWK ( )
virtual

Definition at line 39 of file EvtHypWK.cc.

39{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtHypWK::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 43 of file EvtHypWK.cc.

43{ return new EvtHypWK; }

◆ decay()

void EvtHypWK::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 59 of file EvtHypWK.cc.

59 {
60
61 static EvtId PIP = EvtPDL::getId( "pi+" );
62 static EvtId PIM = EvtPDL::getId( "pi-" );
63 static EvtId PI0 = EvtPDL::getId( "pi0" );
64 static EvtId PP = EvtPDL::getId( "p+" );
65 static EvtId PM = EvtPDL::getId( "anti-p-" );
66 static EvtId NP = EvtPDL::getId( "n0" );
67 static EvtId NM = EvtPDL::getId( "anti-n0" );
68
69 static EvtId LP = EvtPDL::getId( "Lambda0" );
70 static EvtId LM = EvtPDL::getId( "anti-Lambda0" );
71 static EvtId SPP = EvtPDL::getId( "Sigma+" );
72 static EvtId SPM = EvtPDL::getId( "anti-Sigma-" );
73 static EvtId SMP = EvtPDL::getId( "Sigma-" );
74 static EvtId SMM = EvtPDL::getId( "anti-Sigma+" );
75
76 static EvtId XP = EvtPDL::getId( "Xi-" );
77 static EvtId XM = EvtPDL::getId( "anti-Xi+" );
78 static EvtId X0P = EvtPDL::getId( "Xi0" );
79 static EvtId X0M = EvtPDL::getId( "anti-Xi0" );
80
81 static int count = 0;
82
83 // if(count==0) findMaxPro(p); //using the default getProbMax() in
84 // EvtGenBase/EvtDecayProb.cc
85 count += 1;
86
87 p->initializePhaseSpace( getNDaug(), getDaugs() );
88
89 EvtParticle *l1, *l2;
90 l1 = p->getDaug( 0 );
91 l2 = p->getDaug( 1 );
92 EvtVector4R DRP4 = l1->getP4(); // Dirac Daughter
93 EvtVector4R MEP4 = l2->getP4(); // Meson Daughter
94 EvtVector4R PRP4 = p->getP4(); // Parent Momentum
95 EvtHelSys angles( PRP4, DRP4 ); // using helicity sys.angles
96 double tht = angles.getHelAng( 1 );
97 double ph = angles.getHelAng( 2 );
98 double gmm = 0;
99
100 EvtVector4R pDaug = l1->getP4();
101 double p3 = pDaug.d3mag();
102
103 double M = p->mass();
104 double M2 = M * M;
105 double m = l1->mass();
106 double m2 = m * m;
107 double theta, phi, alpha, beta, ga; // the following is the
108
109 EvtId p1, d1, d2;
110 p1 = p->getId();
111 d1 = l1->getId();
112 d2 = l2->getId();
113
114 if ( p1 == XP )
115 {
116 alpha = -0.46;
117 phi = -2.1;
118 }
119 else if ( p1 == XM )
120 {
121 alpha = 0.46;
122 phi = -2.1;
123 }
124 else if ( p1 == X0P )
125 {
126 alpha = -0.41;
127 phi = 21.0;
128 }
129 else if ( p1 == X0M )
130 {
131 alpha = 0.41;
132 phi = 21.0;
133 }
134 else if ( p1 == LP )
135 {
136 alpha = 0.64;
137 phi = -6.5;
138 }
139 else if ( p1 == LM )
140 {
141 alpha = -0.64;
142 phi = -6.5;
143 }
144 else if ( p1 == SMP )
145 {
146 alpha = -0.07;
147 phi = 10.0;
148 }
149 else if ( p1 == SMM )
150 {
151 alpha = 0.07;
152 phi = 10.0;
153 }
154 else if ( p1 == SPP )
155 {
156 if ( d1 == PP )
157 {
158 alpha = -0.98;
159 phi = 36.0;
160 }
161 else if ( d1 == NP )
162 {
163 alpha = 0.07;
164 phi = 167.0;
165 }
166 else { cout << "The decay is not discussed." << endl; }
167 }
168 else if ( p1 == SPM )
169 {
170 if ( d1 == PM )
171 {
172 alpha = 0.98;
173 phi = 36.0;
174 }
175 else if ( d1 == NM )
176 {
177 alpha = -0.07;
178 phi = 167.0;
179 }
180 else { cout << "The decay is not discussed." << endl; }
181 }
182 else { cout << "The decay is not discussed."; }
183 if ( getNArg() > 0 ) { alpha = getArg( 0 ); }
184
185 double FP = sqrt( ( 1 + alpha ) / 2. );
186 double FM = sqrt( ( 1 - alpha ) / 2. );
187
188 vertex( 0, 0, Djmn( 0.5, 0.5, 0.5, ph, tht, gmm ) * FP );
189 vertex( 0, 1, Djmn( 0.5, 0.5, -0.5, ph, tht, gmm ) * FM );
190 vertex( 1, 0, Djmn( 0.5, -0.5, 0.5, ph, tht, gmm ) * FP );
191 vertex( 1, 1, Djmn( 0.5, -0.5, -0.5, ph, tht, gmm ) * FM );
192
193 return;
194}
p1
double p1[4]
Definition BabayagaNLO.cxx:41
Djmn
EvtComplex Djmn(int j, int m, int n, double phi, double theta, double gamma)
Definition EvtHelSys.cc:165
count
DOUBLE_PRECISION count[3]
Definition EvtPhokharaDef.hh:63
alpha
double alpha
Definition EvtPhokharaDef.hh:17
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
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition EvtPDL.cc:272
EvtParticle::getId
EvtId getId() const
Definition EvtParticle.cc:100
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:108
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:80
EvtParticle::mass
double mass() const
Definition EvtParticle.cc:114
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1042
EvtVector4R::d3mag
double d3mag() const
Definition EvtVector4R.cc:168
m2
double double * m2
Definition qcdloop1.h:83

◆ findMaxPro()

void EvtHypWK::findMaxPro ( EvtParticle * p)

◆ getName()

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

Implements EvtDecayBase.

Definition at line 41 of file EvtHypWK.cc.

41{ model_name = "HypWK"; }

◆ init()

void EvtHypWK::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 45 of file EvtHypWK.cc.

45 {
46
47 // check that there are 0 arguments
48 // checkNArg(0);
49 checkNDaug( 2 );
50
51 // checkSpinParent(EvtSpinType::VECTOR);
52 checkSpinParent( EvtSpinType::DIRAC );
53
54 checkSpinDaughter( 0, EvtSpinType::DIRAC );
55 checkSpinDaughter( 1, EvtSpinType::SCALAR );
56 // checkSpinDaughter(1,EvtSpinType::DIRAC);
57}
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
The documentation for this class was generated from the following files: