BOSS 7.1.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtD0ToKKpipi Class Reference

#include <EvtD0ToKKpipi.hh>

+ Inheritance diagram for EvtD0ToKKpipi:

Public Member Functions

 EvtD0ToKKpipi ()
 
virtual ~EvtD0ToKKpipi ()
 
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 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 EvtDecayBase
bool _daugsDecayedByParentModel
 

Detailed Description

Definition at line 8 of file EvtD0ToKKpipi.hh.

Constructor & Destructor Documentation

◆ EvtD0ToKKpipi()

EvtD0ToKKpipi::EvtD0ToKKpipi ( )
inline

Definition at line 10 of file EvtD0ToKKpipi.hh.

10{}

Referenced by clone().

◆ ~EvtD0ToKKpipi()

EvtD0ToKKpipi::~EvtD0ToKKpipi ( )
virtual

Definition at line 34 of file EvtD0ToKKpipi.cc.

34{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0ToKKpipi::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 40 of file EvtD0ToKKpipi.cc.

40 {
41 return new EvtD0ToKKpipi;
42}

◆ decay()

void EvtD0ToKKpipi::decay ( EvtParticle p)
virtual

Implements EvtDecayBase.

Definition at line 59 of file EvtD0ToKKpipi.cc.

59 {
60/*
61 double maxprob=0;
62 for(int ir=0;ir<=60000000;ir++){
63 p->initializePhaseSpace(getNDaug(),getDaugs());
64 EvtVector4R Kp = p->getDaug(0)->getP4();
65 EvtVector4R Km = p->getDaug(1)->getP4();
66 EvtVector4R pip = p->getDaug(2)->getP4();
67 EvtVector4R pim = p->getDaug(3)->getP4();
68
69 int id = EvtPDL::getStdHep(p->getId());
70 int charm;
71 if(id > 0) charm = 1;
72 else charm = -1;
73
74 double X[16];
75 std::complex<double> Total_amp;
76 if (charm==1) {
77 X[0] = Kp.get(1); X[4] = Km.get(1); X[8] = pip.get(1); X[12] = pim.get(1);
78 X[1] = Kp.get(2); X[5] = Km.get(2); X[9] = pip.get(2); X[13] = pim.get(2);
79 X[2] = Kp.get(3); X[6] = Km.get(3); X[10] = pip.get(3); X[14] = pim.get(3);
80 X[3] = Kp.get(0); X[7] = Km.get(0); X[11] = pip.get(0); X[15] = pim.get(0);
81 Total_amp = AMP(X, 1);
82 } else {
83 X[0] = -Km.get(1); X[4] = -Kp.get(1); X[8] = -pim.get(1); X[12] = -pip.get(1);
84 X[1] = -Km.get(2); X[5] = -Kp.get(2); X[9] = -pim.get(2); X[13] = -pip.get(2);
85 X[2] = -Km.get(3); X[6] = -Kp.get(3); X[10] = -pim.get(3); X[14] = -pip.get(3);
86 X[3] = Km.get(0); X[7] = Kp.get(0); X[11] = pim.get(0); X[15] = pip.get(0);
87 Total_amp = AMP(X,1);
88 }
89 double value1 = std::norm(Total_amp);
90 if(value1>maxprob) {
91 maxprob=value1;
92 cout << ir << " max= " << maxprob << endl;
93 }
94 }
95 cout << "Maxprob= " << maxprob << endl;
96*/
98 EvtVector4R Kp = p->getDaug(0)->getP4();
99 EvtVector4R Km = p->getDaug(1)->getP4();
100 EvtVector4R pip = p->getDaug(2)->getP4();
101 EvtVector4R pim = p->getDaug(3)->getP4();
102
103 int id = EvtPDL::getStdHep(p->getId());
104 int charm;
105 if(id > 0) charm = 1;
106 else charm = -1;
107
108 double X[16];
109 std::complex<double> Total_amp;
110 if (charm==1) {
111 X[0] = Kp.get(1); X[4] = Km.get(1); X[8] = pip.get(1); X[12] = pim.get(1);
112 X[1] = Kp.get(2); X[5] = Km.get(2); X[9] = pip.get(2); X[13] = pim.get(2);
113 X[2] = Kp.get(3); X[6] = Km.get(3); X[10] = pip.get(3); X[14] = pim.get(3);
114 X[3] = Kp.get(0); X[7] = Km.get(0); X[11] = pip.get(0); X[15] = pim.get(0);
115 Total_amp = AMP(X, 1);
116 } else {
117 X[0] = -Km.get(1); X[4] = -Kp.get(1); X[8] = -pim.get(1); X[12] = -pip.get(1);
118 X[1] = -Km.get(2); X[5] = -Kp.get(2); X[9] = -pim.get(2); X[13] = -pip.get(2);
119 X[2] = -Km.get(3); X[6] = -Kp.get(3); X[10] = -pim.get(3); X[14] = -pip.get(3);
120 X[3] = Km.get(0); X[7] = Kp.get(0); X[11] = pim.get(0); X[15] = pip.get(0);
121 Total_amp = AMP(X,1);
122 }
123 double value = std::norm(Total_amp);
124 setProb(value);
125 return;
126}
EvtId * getDaugs()
Definition: EvtDecayBase.hh:65
void setProb(double prob)
Definition: EvtDecayProb.hh:34
static int getStdHep(EvtId id)
Definition: EvtPDL.hh:56
EvtId getId() const
Definition: EvtParticle.cc:112
const EvtVector4R & getP4() const
Definition: EvtParticle.cc:120
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:84
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
double get(int i) const
Definition: EvtVector4R.hh:179

◆ getName()

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

Implements EvtDecayBase.

Definition at line 36 of file EvtD0ToKKpipi.cc.

36 {
37 model_name="D0ToKKpipi";
38}

◆ init()

void EvtD0ToKKpipi::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 44 of file EvtD0ToKKpipi.cc.

44 {
45 checkNArg(0);
46 checkNDaug(4);
52 //cout << "Initializing D0ToKKpipi" << endl;
53}
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 EvtD0ToKKpipi::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 55 of file EvtD0ToKKpipi.cc.

55 {
56 setProbMax(33.0);
57}
void setProbMax(double prbmx)

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