BOSS 7.1.2
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtD0ToKSpi0pi0pi0.hh
Go to the documentation of this file.
1#ifndef EVTD0TOKSPI0PI0PI0_HH
2#define EVTD0TOKSPI0PI0PI0_HH
3
4#include "EvtGenBase/EvtDecayProb.hh"
5class EvtParticle;
6
7class EvtD0ToKSpi0pi0pi0:public EvtDecayProb {
8
9 public:
10
11 EvtD0ToKSpi0pi0pi0() {}
12 virtual ~EvtD0ToKSpi0pi0pi0();
13
14 void getName(std::string& name);
15 EvtDecayBase* clone();
16
17 void init();
18 void initProbMax();
19 void decay(EvtParticle *p);
20
21 private:
22 double calPDF(double KS[], double Kp[], double Pim[], double Pi0[], double & Result);
23
24 void Com_Multi(double a1[2], double a2[2], double res[2]);
25 void Com_Divide(double a1[2], double a2[2], double res[2]);
26 double SCADot(double a1[4], double a2[4]);
27 double Barrier(double mass2, int l, double sa, double sb, double sc, double r2);
28 double barrier(int l, double sa, double sb, double sc, double r2, double Smass);
29 void calt1(double daug1[4], double daug2[4], double t1[4]);
30 void calt2(double daug1[4], double daug2[4], double t2[4][4]);
31 void propagator(double mass2, double mass, double width, double sx, double prop[2]);
32 double wid(double mass2, double mass, double sa, double sb, double sc, double r2, int l);
33 double widl1(double mass2, double mass, double sa, double sb, double sc, double r2);
34 void propagatorRBW(double mass2, double mass, double width, double sa, double sb, double sc, double r2, int l, double prop[2]);
35 void propagatorRBWl1(double mass2, double mass, double width, double sa, double sb, double sc, double r2, double prop[2]);
36 void propagatorRBW_a1(double mass2, double mass, double width, double sa, double sb, double sc, double r2, int l, double prop[2]);
37 void propagatorRBW_K1(double mass2, double mass, double width, double sa, double sb, double sc, double r2, int l, double prop[2]);
38 void propagatorGS(double mass2, double mass, double width, double sa, double sb, double sc, double r2, double prop[2]);
39 void rhoab(double sa, double sb, double sc, double res[2]);
40 void rho4Pi(double sa, double res[2]);
41 void propagatorsigma500(double sa, double sb, double sc, double prop[2]);
42 void Flatte_rhoab(double sa, double sb, double sc, double rho[2]);
43 void propagator980(double mass, double sx, double *sb, double *sc, double prop[2]);
44 void propagatorFlatte(double mass, double width, double sa, double prop[2]);
45 void KPiSLASS(double sa, double sb, double sc, double prop[2]);
46 //double a1_width[2000], K1270_width[2000], K1400_width[2000];
47 void PiPiSWAVE(double sa, double sb, double sc, double prop[2]);
48 //pipi S-wave K-matrix
49 double CalRho4pi(double_t s);
50 void rhoMTX(int i, int j, double s, double Rho[2]);
51 void KMTX(int i, int j, double s,double KM[2]);
52 void IMTX(int i, int j, double IMTX[2]);
53 void FMTX(double Kijx, double Kijy, double rhojjx, double rhojjy, int i, int j, double FM[2]);
54 void PVTR(int ID, double s, double PV[2],double sp0, double f1, double f2, double f3, double f4, double f5, double f6, double f7, double f8, double f9, double f10, double b1, double b2, double b3, double b4, double b5, double b6, double b7, double b8, double b9, double b10);
55 void FINVMTX(double s, double *FINVx, double *FINVy);
56 void Fvector(double sa, double s0, double Fv[2],double sp0, double f1, double f2, double f3, double f4, double f5, double f6, double f7, double f8, double f9, double f10, double b1, double b2, double b3, double b4, double b5, double b6, double b7, double b8, double b9, double b10);
57
58 //Decay process
59 void DtoVP_S(double Ks_e, double Ks_px, double Ks_py, double Ks_pz,
60 double Pi01_e, double Pi01_px, double Pi01_py, double Pi01_pz,
61 double Pi02_e, double Pi02_px, double Pi02_py, double Pi02_pz,
62 double Pi03_e, double Pi03_px, double Pi03_py, double Pi03_pz,
63 double mass1, double width1, int ang1, int ang2, int g0, double propagator[], double& tmp_PDF,
64 double sp0, double f1, double f2, double f3, double f4, double f5, double f6, double f7, double f8, double f9, double f10,
65 double b1, double b2, double b3, double b4, double b5, double b6, double b7, double b8, double b9, double b10 );
66
67 void DtoAP( double Ks_e, double Ks_px, double Ks_py, double Ks_pz,
68 double Pi01_e, double Pi01_px, double Pi01_py, double Pi01_pz,
69 double Pi02_e, double Pi02_px, double Pi02_py, double Pi02_pz,
70 double Pi03_e, double Pi03_px, double Pi03_py, double Pi03_pz,
71 double mass1, double mass2, double width1, double width2, int ang1, int ang2, int ang3, int f0, double propagator[], double& tmp_PDF,
72 double akeMag, double a32Mag, double akePhs, double a32Phs);
73
74 void DtoPP(double Ks_e, double Ks_px, double Ks_py, double Ks_pz,
75 double Pi01_e, double Pi01_px, double Pi01_py, double Pi01_pz,
76 double Pi02_e, double Pi02_px, double Pi02_py, double Pi02_pz,
77 double Pi03_e, double Pi03_px, double Pi03_py, double Pi03_pz,
78 double mass1, double mass2, double width1, double width2,int ang1, int ang2, int ang3, double propagator[], double& tmp_PDF,
79 double sp0, double f1, double f2, double f3, double f4, double f5, double f6, double f7, double f8, double f9, double f10,
80 double b1, double b2, double b3, double b4, double b5, double b6, double b7, double b8, double b9, double b10);
81
82 void DtoaP( double Ks_e, double Ks_px, double Ks_py, double Ks_pz,
83 double Pi01_e, double Pi01_px, double Pi01_py, double Pi01_pz,
84 double Pi02_e, double Pi02_px, double Pi02_py, double Pi02_pz,
85 double Pi03_e, double Pi03_px, double Pi03_py, double Pi03_pz,
86 double mass1, double width1, int ang1, int ang2, int g0, double propagator[], double& tmp_PDF,
87 double sp0, double f1, double f2, double f3, double f4, double f5, double f6, double f7, double f8, double f9, double f10,
88 double b1, double b2, double b3, double b4, double b5, double b6, double b7, double b8, double b9, double b10);
89
90 void kpiSwave(double m2,
91 double akeMag, double a32Mag, double akePhs, double a32Phs, double propagator[]);
92
93 void DtoKPi( double Ks_e, double Ks_px, double Ks_py, double Ks_pz,
94 double Pi01_e, double Pi01_px, double Pi01_py, double Pi01_pz,
95 double Pi02_e, double Pi02_px, double Pi02_py, double Pi02_pz,
96 double Pi03_e, double Pi03_px, double Pi03_py, double Pi03_pz,
97 int g0, int g1, double propagator[], double& tmp_PDF,
98 double sp0, double f1, double f2, double f3, double f4, double f5, double f6, double f7, double f8, double f9, double f10,
99 double b1, double b2, double b3, double b4, double b5, double b6, double b7, double b8, double b9, double b10,
100 double akeMag, double a32Mag, double akePhs, double a32Phs);
101
102
103 int mother_c;
104 double G[4][4],E[4][4][4][4];
105 //double rho[24];
106 //double phi[24];
107 int modetype[20];
108
109 double rho01, rho02, rho03, rho04, rho10, rho06, rho07, rho08, rho126;
110 double phi01, phi02, phi03, phi04, phi10, phi06, phi07, phi08, phi126;
111 double Kstf1, Kstf6, Kstb2, Kstb7;
112 double K1460f2, K1460f7, K1460b1, K1460b2, K1460b6, K1460b7;
113 double a1260f1, a1260f6, a1260b1, a1260b2, a1260b6, a1260b7;
114 double Phspf1, Phspf6, Phspb2, Phspb7;
115 double ake1, ake2, a321, a322;
116
117 double mrhop, Grhop, mrho0, Grho0, mKstr0, GKstr0, mKstrp, GKstrp, msigma, Gsigma, mOmega, GOmega, ma980, Ga980;
118 double ma1_1260, Ga1_1260, mK1_1270, GK1_1270, mK1_1400, GK1_1400, mK1460, GK1460, mK1650, GK1650, mpi_1300, Gpi_1300, mD, mPi, mKa;
119 double mf0_500, Gf0_500;
120 double math_pi;
121 double mass_Pion, mass_Eta, rD2, rRes2;
122};
123
124#endif
mass
double mass
Definition CosmicGenerator.cxx:138
f1
TFile * f1
Definition DataBase/tau_mode.c:4
g1
TF1 * g1
Definition DataBase/tau_mode.c:66
ID
int ID[no]
Definition DocaEangle/fitDocaeangle.h:2
s
XmlRpcServer s
Definition HelloServer.cpp:11
Rho
***************************************************************************************Pseudo Class RRes *****************************************************************************************Parameters and physical constants **Maarten sept ************************************************************************DOUBLE PRECISION xsmu **************************************************************************PARTICLE DATA ** Rho(770) and Omega(782) are taken from CMD-2 F_pi fit *(hep-ex/9904027)
EvtD0ToKSpi0pi0pi0::getName
void getName(std::string &name)
Definition EvtD0ToKSpi0pi0pi0.cc:35
EvtD0ToKSpi0pi0pi0::decay
void decay(EvtParticle *p)
Definition EvtD0ToKSpi0pi0pi0.cc:158
EvtD0ToKSpi0pi0pi0::clone
EvtDecayBase * clone()
Definition EvtD0ToKSpi0pi0pi0.cc:39
m2
double double * m2
Definition qcdloop1.h:75