BOSS 8.0.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 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 35 of file EvtD0ToKKpipi.cc.

35{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0ToKKpipi::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 39 of file EvtD0ToKKpipi.cc.

39{ return new EvtD0ToKKpipi; }

◆ decay()

void EvtD0ToKKpipi::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 54 of file EvtD0ToKKpipi.cc.

54 {
55 /*
56 double maxprob=0;
57 for(int ir=0;ir<=60000000;ir++){
58 p->initializePhaseSpace(getNDaug(),getDaugs());
59 EvtVector4R Kp = p->getDaug(0)->getP4();
60 EvtVector4R Km = p->getDaug(1)->getP4();
61 EvtVector4R pip = p->getDaug(2)->getP4();
62 EvtVector4R pim = p->getDaug(3)->getP4();
63
64 int id = EvtPDL::getStdHep(p->getId());
65 int charm;
66 if(id > 0) charm = 1;
67 else charm = -1;
68
69 double X[16];
70 std::complex<double> Total_amp;
71 if (charm==1) {
72 X[0] = Kp.get(1); X[4] = Km.get(1); X[8] = pip.get(1); X[12] = pim.get(1);
73 X[1] = Kp.get(2); X[5] = Km.get(2); X[9] = pip.get(2); X[13] = pim.get(2);
74 X[2] = Kp.get(3); X[6] = Km.get(3); X[10] = pip.get(3); X[14] = pim.get(3);
75 X[3] = Kp.get(0); X[7] = Km.get(0); X[11] = pip.get(0); X[15] = pim.get(0);
76 Total_amp = AMP(X, 1);
77 } else {
78 X[0] = -Km.get(1); X[4] = -Kp.get(1); X[8] = -pim.get(1); X[12] = -pip.get(1);
79 X[1] = -Km.get(2); X[5] = -Kp.get(2); X[9] = -pim.get(2); X[13] = -pip.get(2);
80 X[2] = -Km.get(3); X[6] = -Kp.get(3); X[10] = -pim.get(3); X[14] = -pip.get(3);
81 X[3] = Km.get(0); X[7] = Kp.get(0); X[11] = pim.get(0); X[15] = pip.get(0);
82 Total_amp = AMP(X,1);
83 }
84 double value1 = std::norm(Total_amp);
85 if(value1>maxprob) {
86 maxprob=value1;
87 cout << ir << " max= " << maxprob << endl;
88 }
89 }
90 cout << "Maxprob= " << maxprob << endl;
91 */
93 EvtVector4R Kp = p->getDaug( 0 )->getP4();
94 EvtVector4R Km = p->getDaug( 1 )->getP4();
95 EvtVector4R pip = p->getDaug( 2 )->getP4();
96 EvtVector4R pim = p->getDaug( 3 )->getP4();
97
98 int id = EvtPDL::getStdHep( p->getId() );
99 int charm;
100 if ( id > 0 ) charm = 1;
101 else charm = -1;
102
103 double X[16];
104 std::complex<double> Total_amp;
105 if ( charm == 1 )
106 {
107 X[0] = Kp.get( 1 );
108 X[4] = Km.get( 1 );
109 X[8] = pip.get( 1 );
110 X[12] = pim.get( 1 );
111 X[1] = Kp.get( 2 );
112 X[5] = Km.get( 2 );
113 X[9] = pip.get( 2 );
114 X[13] = pim.get( 2 );
115 X[2] = Kp.get( 3 );
116 X[6] = Km.get( 3 );
117 X[10] = pip.get( 3 );
118 X[14] = pim.get( 3 );
119 X[3] = Kp.get( 0 );
120 X[7] = Km.get( 0 );
121 X[11] = pip.get( 0 );
122 X[15] = pim.get( 0 );
123 Total_amp = AMP( X, 1 );
124 }
125 else
126 {
127 X[0] = -Km.get( 1 );
128 X[4] = -Kp.get( 1 );
129 X[8] = -pim.get( 1 );
130 X[12] = -pip.get( 1 );
131 X[1] = -Km.get( 2 );
132 X[5] = -Kp.get( 2 );
133 X[9] = -pim.get( 2 );
134 X[13] = -pip.get( 2 );
135 X[2] = -Km.get( 3 );
136 X[6] = -Kp.get( 3 );
137 X[10] = -pim.get( 3 );
138 X[14] = -pip.get( 3 );
139 X[3] = Km.get( 0 );
140 X[7] = Kp.get( 0 );
141 X[11] = pim.get( 0 );
142 X[15] = pip.get( 0 );
143 Total_amp = AMP( X, 1 );
144 }
145 double value = std::norm( Total_amp );
146 setProb( value );
147 return;
148}
EvtId * getDaugs()
void setProb(double prob)
static int getStdHep(EvtId id)
Definition EvtPDL.hh:61
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 get(int i) const

◆ getName()

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

Implements EvtDecayBase.

Definition at line 37 of file EvtD0ToKKpipi.cc.

37{ model_name = "D0ToKKpipi"; }

◆ init()

void EvtD0ToKKpipi::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 41 of file EvtD0ToKKpipi.cc.

41 {
42 checkNArg( 0 );
43 checkNDaug( 4 );
49 // cout << "Initializing D0ToKKpipi" << endl;
50}
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 52 of file EvtD0ToKKpipi.cc.

52{ setProbMax( 33.0 ); }
void setProbMax(double prbmx)

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