BOSS 7.1.2
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDsToKKenu Class Reference

#include <EvtDsToKKenu.hh>

+ Inheritance diagram for EvtDsToKKenu:

Public Member Functions

 EvtDsToKKenu ()
 
virtual ~EvtDsToKKenu ()
 
void getName (std::string &name)
 
EvtDecayBaseclone ()
 
void init ()
 
void initProbMax ()
 
void decay (EvtParticle *p)
 
- Public Member Functions inherited from EvtDecayProb
void makeDecay (EvtParticle *p)
 
void setProb (double prob)
 
double getProb ()
 
void setWeight (double weight)
 
virtual ~EvtDecayProb ()
 
- 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 EvtDecayBase
bool _daugsDecayedByParentModel
 

Detailed Description

Definition at line 8 of file EvtDsToKKenu.hh.

Constructor & Destructor Documentation

◆ EvtDsToKKenu()

EvtDsToKKenu::EvtDsToKKenu ( )
inline

Definition at line 12 of file EvtDsToKKenu.hh.

12{}

Referenced by clone().

◆ ~EvtDsToKKenu()

EvtDsToKKenu::~EvtDsToKKenu ( )
virtual

Definition at line 27 of file EvtDsToKKenu.cc.

27{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtDsToKKenu::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 33 of file EvtDsToKKenu.cc.

33 {
34 return new EvtDsToKKenu;
35}

◆ decay()

void EvtDsToKKenu::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 112 of file EvtDsToKKenu.cc.

112 {
113/*
114 double maxprob = 0.0;
115 for(int ir=0;ir<=30000000;ir++){
116 p->initializePhaseSpace(getNDaug(),getDaugs());
117 EvtVector4R _K = p->getDaug(0)->getP4();
118 EvtVector4R _pi = p->getDaug(1)->getP4();
119 EvtVector4R _e = p->getDaug(2)->getP4();
120 EvtVector4R _nu = p->getDaug(3)->getP4();
121 int pid = EvtPDL::getStdHep(p->getDaug(0)->getId());
122 int charm;
123 if(pid == -321) charm = 1;
124 else charm = -1;
125 double m2, q2, cosV, cosL, chi;
126 KinVGen(_K, _pi, _e, _nu, charm, m2, q2, cosV, cosL, chi);
127 double _prob = calPDF(m2, q2, cosV, cosL, chi);
128 if(_prob>maxprob) {
129 maxprob=_prob;
130 std::cout << "Max PDF = " << ir << " charm= " << charm << " prob= " << _prob << std::endl;
131 }
132 }
133 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
134*/
135
137 EvtVector4R K = p->getDaug(0)->getP4();
138 EvtVector4R pi = p->getDaug(1)->getP4();
139 EvtVector4R e = p->getDaug(2)->getP4();
140 EvtVector4R nu = p->getDaug(3)->getP4();
141
142 int pid = EvtPDL::getStdHep(p->getDaug(0)->getId());
143 int charm;
144 if(pid == -321) charm = 1;
145 else charm = -1;
146 double m2, q2, cosV, cosL, chi;
147 KinVGen(K, pi, e, nu, charm, m2, q2, cosV, cosL, chi);
148 double prob = calPDF(m2, q2, cosV, cosL, chi);
149 setProb(prob);
150 return;
151}
EvtId * getDaugs()
void setProb(double prob)
static int getStdHep(EvtId id)
Definition EvtPDL.hh:56
EvtId getId() const
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
double double * m2
Definition qcdloop1.h:75
const float pi
Definition vector3.h:133

◆ getName()

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

Implements EvtDecayBase.

Definition at line 29 of file EvtDsToKKenu.cc.

29 {
30 model_name="DsToKKenu";
31}

◆ init()

void EvtDsToKKenu::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 37 of file EvtDsToKKenu.cc.

37 {
38 checkNArg(0);
39 checkNDaug(4);
43
44 //std::cout << "EvtDsToKKenu ==> Initialization !" << std::endl;
45 nAmps = 2;
46
47 //rS = -11.57; // S-wave
48 //rS1 = 0.08;
49 //a_delta = 1.94;
50 //b_delta = -0.81;
51 //m0_1430_S = 1.425;
52 //width0_1430_S = 0.270;
53 //type[0] = 0;
54
55 mV = 2.1;
56 mA = 2.5;
57
58//for Ds to f0 e nu
59 mf0 = 0.965;
60 m2f0 = mf0*mf0;
61 g1 = 0.165*mf0;
62 g2 = 4.21*0.165*mf0;
63 rho_f0 = 7.1762;
64 phi_f0 = -3.3642;
65 //mDs = 1.96834;
66 type[0] = 0;
67
68 V_0 = 1.6456;
69 A1_0 = 1;
70 A2_0 = 0.83083;
71
72 m0 = 1.0195; // P-wave K*
73 width0 = 0.004249;
74 rBW = 3.0;
75 rho = 1;
76 phi = 0;
77 type[1] = 1;
78
79 m0_1410 = 1.414; // P-wave K*(1410)
80 width0_1410 = 0.232;
81 rho_1410 = 0.1;
82 phi_1410 = 0.;
83 type[2] = 2;
84
85 TV_0 = 1; // D-wave K*2(1430)
86 T1_0 = 1;
87 T2_0 = 1;
88 m0_1430 = 1.4324;
89 width0_1430 = 0.109;
90 rho_1430 = 15;
91 phi_1430 = 0;
92 type[3] = 3;
93
94 mPi = 0.13957;
95 mPi0 = 0.1349766;
96 mK0 = 0.497611;
97
98 mD = 1.96834;//Ds
99 mK1 = 0.49368;//k
100 mK2 = 0.49368;//k
101 Pi = atan2(0.0,-1.0);
102 root2 = sqrt(2.);
103 root2d3 = sqrt(2./3);
104 root1d2 = sqrt(0.5);
105 root3d2 = sqrt(1.5);
106}
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 EvtDsToKKenu::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 108 of file EvtDsToKKenu.cc.

108 {
109 setProbMax(450880.0);
110}
void setProbMax(double prbmx)

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