40 {
44
45
46
47
48
49
50 std::cout << "EvtDToKSpipipi ==> Initialization !" << std::endl;
51
52 mK1270 = 1.272; mK1400 = 1.403;
53 GK1270 = 0.09; GK1400 = 0.174;
54 mKstr = 0.89166; mrho = 0.77549;
55 GKstr = 0.0462; Grho = 0.1491;
56 msigma = 0.472;
57 Gsigma = 0.542;
58 phi_omega = -0.02;
59 mK1650 = 1.65;
60 GK1650 = 0.158;
61 rho[0] = 1.0;
62 phi[0] = 0.0;
63
64 ma1 = 1.22;
65 Ga1 = 0.4282;
66 mK1460 = 1.4152;
67 GK1460 = 0.2485;
68 rho_omega = 0.00294;
69
70 phi[1] = -1.55;
71 rho[1] = 0.5843;
72
73 phi[2] = -1.8223;
74 rho[2] = 2.0974;
75
76 phi[3] = -2.6751;
77 rho[3] = 0.46642;
78
79 phi[4] = -2.2429;
80 rho[4] = 0.33334;
81
82 phi[5] = -0.55888;
83 rho[5] = 0.15549;
84
85 phi[6] = -1.8778;
86 rho[6] = 0.94452;
87
88 phi[7] = 2.7724;
89 rho[7] = 0.99411;
90
91 phi[8] = -2.6461;
92 rho[8] = 0.21231;
93
94 phi[9] = -0.95137;
95 rho[9] = 0.29895;
96
97 phi[10] = -3.0843;
98 rho[10] = 3.6361;
99
100 phi[11] = 2.0954;
101 rho[11] = 0.96472;
102
103 phi[12] = -2.4965;
104 rho[12] = 0.48470;
105
106 for (int i=0; i<13; i++) {
107 cout << i << "rho,phi = " << rho[i] << ", "<< phi[i] << endl;
108 }
109
110 mD = 1.86486;
111 rD = 5;
112 metap = 0.95778;
113 mkstr = 0.89594;
114 mk0 = 0.497614;
115 mass_Kaon = 0.49368;
116 mass_Pion = 0.13957;
117 math_pi = 3.1415926;
118
119 pi = 3.1415926;
120 mpi = 0.13957;
121 g1 = 0.5468;
122 g2 = 0.23;
123
124 int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
125 int EE[4][4][4][4] =
126 { { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
127 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
128 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
129 {{0,0,0,0}, {0,0,1,0}, {0,-1,0,0}, {0,0,0,0} } },
130 { {{0,0,0,0}, {0,0,0,0}, {0,0,0,-1}, {0,0,1,0} },
131 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
132 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
133 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
134 { {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
135 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
136 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
137 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
138 { {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
139 {{0,0,1,0}, {0,0,0,0}, {-1,0,0,0}, {0,0,0,0} },
140 {{0,-1,0,0}, {1,0,0,0}, {0,0,0,0}, {0,0,0,0} },
141 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} } } };
142 for (int i=0; i<4; i++) {
143 for (int j=0; j<4; j++) {
144 G[i][j] = GG[i][j];
145 for (int k=0; k<4; k++) {
146 for (int l=0; l<4; l++) {
147 E[i][j][k][l] = EE[i][j][k][l];
148 }
149 }
150 }
151 }
152}
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)