BOSS 6.6.4.p03
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtPsi3Sdecay.hh
Go to the documentation of this file.
1//--------------------------------------------------------------------------
2//
3// Environment:
4// This software is part of models developed at BES collaboration
5// based on the EvtGen framework. If you use all or part
6// of it, please give an appropriate acknowledgement.
7//
8// Copyright Information: See EvtGen/BesCopyright
9// Copyright (A) 2006 Ping Rong-Gang, Pang Cai-Ying@IHEP
10//
11// Module: EvtPsi3Sdecay.hh
12//
13// Description: Routine to re-select the psi(4040) decay according the .
14// measured xsection at different energy point, see CLEOc measurement:
15// PRD 80, 072001
16// Modification history:
17//
18// Ping R.-G. December, 2010 Module created
19//
20//------------------------------------------------------------------------
21#ifndef EVTPsi3Sdecay_HH
22#define EVTPsi3Sdecay_HH
23
25#include "EvtGenBase/EvtId.hh"
26#include "EvtGenBase/EvtPDL.hh"
28#include <stdlib.h>
29#include <string>
30#include <vector>
31using std::endl;
32
34public:
35 // constructor
36 EvtPsi3Sdecay(double ecms, EvtParticle *parent) { //for 2-body decays
37 //initializer
38 Ecms = ecms;
39 theParent = parent;
40 Ndaugs=parent->getNDaug();
41 nsize = 24;
42
43 _excflag=0;
44 x.clear();
45 // open charm cross section, see PRD 80, 072001, xs in unit pb.
46 double xx[24]={3.72968, 3.97, 3.99, 4.01, 4.015, 4.03, 4.06, 4.12, 4.14, 4.16, 4.17, 4.18, 4.2, 4.26, 4.30, 4.34, 4.38, 4.42, 4.46, 4.50, 4.54, 4.58, 4.62, 4.66 }; // 14 energy points
47 double y0[24]={0., 86, 133, 76, 10, 334, 410, 303, 177, 167, 177, 179, 180, 86, 31, 49, 65, 196, 52, 87, 166, 14, 33, 49};// 0) D0D0bar cross section
48 double y1[24]={0., 137, 90, 135, 38, 196, 480, 310, 200, 200, 182, 197, 181, 94, 108, 96, 154, 165, 171, 106, 27, 144, 36, 22}; // 1) D+D-
49 double y2[24]={0., 1140, 1370, 1660, 1920, 1600, 1115, 700, 675, 626, 636, 605, 515, 540, 748, 880, 556, 657, 477, 494, 320, 616, 575, 373};// 2)D0D*0bar
50 double y3[24]={0., 1140, 1370, 1660, 1920, 1600, 1115, 700, 675, 626, 636, 605, 515, 540, 748, 880, 556, 657, 477, 494, 320, 616, 575, 373};// 3)D0bar D*0
51 double y4[24]={0., 0, 0, 0, 213, 2000, 2290, 2550, 2443, 2566, 2363, 2173, 1830, 269, 822, 1045, 1020, 820, 398, 588, 690, 459, 360, 398};// 4)D*0 D*0bar
52 double y5[24]={0., 1115, 1375, 1650, 1851, 1650, 1085, 780, 688, 688, 642, 648, 535, 511, 748, 880, 556, 657, 477, 494, 320, 616, 575, 373};// 5)D*+D-
53 double y6[24]={0., 1115, 1375, 1650, 1851, 1650, 1085, 780, 688, 688, 642, 648, 535, 511, 748, 880, 556, 657, 477, 494, 320, 616, 575, 373};// 6)D*-D+
54 double y7[24]={0., 0, 0, 0, 0, 1400, 2390, 2280, 2556, 2479, 2357, 2145, 1564, 237, 822, 1045, 1020, 820, 398, 588, 690, 459, 360, 398};// 7)D*+D*-
55 double y8[24]={0., 102, 133, 269, 250, 174, 51, 26, 25, 15, 34, 7, 15, 47, 106, 70, 36, 10, 2, 28, 60, 60, 48, 36}; // 8)Ds+ Ds-
56 double y9[24]={0., 0, 0, 0, 0, 0, 0, 478, 684, 905, 916, 889, 812, 34, 314, 368, 318, 357, 292, 171, 66, 103, 190, 272};// 9)Ds*+ Ds-
57 double y10[24]={0., 0, 0, 0, 0, 0, 0, 478, 684, 905, 916, 889, 812, 34, 314, 368, 318, 357, 292, 171, 66, 103, 190, 272};//10)Ds*- Ds+
58 double y11[24]={0., 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 440, 398, 428, 310, 131, 0, 45, 126, 98, 39, 0}; //11)Ds*+ Ds*-
59 double y12[24]={0., 0, 0, 0, 0, 0, 12, 3.8, 32.4, 32.4, 37, 48, 61.2, 53.6, 49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 };//12)D*+ D- pi0 //------ DD* pi----
60 double y13[24]={0., 0, 0, 0, 0, 0, 12, 3.8, 32.4, 32.4, 37, 48, 61.2, 53.6, 49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 }; //13)D*- D+ pi0
61 double y14[24]={0., 0, 0, 0, 0, 0, 24, 7.5, 68.6, 64.8, 73.3, 95.8, 122.5,106.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110}; //14)D*+ anti-D0 pi-
62 double y15[24]={0., 0, 0, 0, 0, 0, 24, 7.5, 68.6, 64.8, 73.3, 95.8, 122.5,106.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110}; //15)D*- D0 pi+
63 double y16[24]={0., 0, 0, 0, 0, 0, 24, 7.5, 68.6, 64.8, 73.3, 95.8, 122.5,106.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110}; //16)D+ anti-D*0 pi-
64 double y17[24]={0., 0, 0, 0, 0, 0, 24, 7.5, 68.6, 64.8, 73.3, 95.8, 122.5,106.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110}; //17)D- D*0 pi+
65 double y18[24]={0., 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 53.6,49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 }; //18) D*+ D*- pi0 //------D*D*pi, above 4.26Gev, assumed xs as D*D pi
66 double y19[24]={0., 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,107.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110 }; //19) anti-D*0 D*+ pi-
67 double y20[24]={0., 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,107.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110 }; //20) D*0 D*- pi+
68 double y21[24]={0., 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 53.6,49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 }; //21) D*0 D*0bar pi0
69 double y22[24]={0., 0, 0, 0, 0, 0, 12, 3.8, 32.4, 32.4, 37, 48, 61.2, 53.6,49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 }; //22)D*0 D0bar pi0 //------ DD* pi----
70 double y23[24]={0., 0, 0, 0, 0, 0, 12, 3.8, 32.4, 32.4, 37, 48, 61.2, 53.6,49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 }; //23)D*bar D0 pi0
71 d0d0bar.clear();
72 dpdm.clear();
73 d0dst0bar.clear();
74 dst0dst0bar.clear();
75 d0bardst0.clear();
76 dstpdm.clear();
77 dstmdp.clear();
78 dstpdstm.clear();
79 dspdsm.clear();
80 dsspdsm.clear();
81 dssmdsp.clear();
82 dsspdssm.clear();
83 xs12.clear();
84 xs13.clear();
85 xs14.clear();
86 xs15.clear();
87 xs16.clear();
88 xs17.clear();
89 xs18.clear();
90 xs19.clear();
91 xs20.clear();
92 xs21.clear();
93 xs22.clear();
94 xs23.clear();
95
96 for(int i=0;i<24;i++){
97 x.push_back(xx[i]);
98 d0d0bar.push_back(y0[i]);
99 dpdm.push_back(y1[i]);
100 d0dst0bar.push_back(y2[i]);
101 d0bardst0.push_back(y3[i]);
102 dst0dst0bar.push_back(y4[i]);
103 dstpdm.push_back( y5[i]);
104 dstmdp.push_back( y6[i]);
105 dstpdstm.push_back(y7[i]);
106 dspdsm.push_back( y8[i]);
107 dsspdsm.push_back( y9[i]);
108 dssmdsp.push_back( y10[i]);
109 dsspdssm.push_back( y11[i]);
110 xs12.push_back( y12[i] );
111 xs13.push_back( y13[i] );
112 xs14.push_back( y14[i] );
113 xs15.push_back( y15[i] );
114 xs16.push_back( y16[i] );
115 xs17.push_back( y17[i] );
116 xs18.push_back( y18[i] );
117 xs19.push_back( y19[i] );
118 xs20.push_back( y20[i] );
119 xs21.push_back( y21[i] );
120 xs22.push_back( y22[i] );
121 xs23.push_back( y23[i] );
122 }
123 }
124
125//---
126 EvtPsi3Sdecay() {//for 2-body and 3-body decays
127 //initializer
128 // Ecms = ecms;
129 nsize = 24;
130
131 x.clear();
132 // open charm cross section, see PRD 80, 072001, xross section in pb
133 double xx[24]={3.72968, 3.97, 3.99, 4.01, 4.015, 4.03, 4.06, 4.12, 4.14, 4.16, 4.17, 4.18, 4.2, 4.26, 4.30, 4.34, 4.38, 4.42, 4.46, 4.50, 4.54, 4.58, 4.62, 4.66 }; // 14 energy points
134 double y0[24]={0., 86, 133, 76, 10, 334, 410, 303, 177, 167, 177, 179, 180, 86, 31, 49, 65, 196, 52, 87, 166, 14, 33, 49};// 0) D0D0bar cross section
135 double y1[24]={0., 137, 90, 135, 38, 196, 480, 310, 200, 200, 182, 197, 181, 94, 108, 96, 154, 165, 171, 106, 27, 144, 36, 22}; // 1) D+D-
136 double y2[24]={0., 1140, 1370, 1660, 1920, 1600, 1115, 700, 675, 626, 636, 605, 515, 540, 748, 880, 556, 657, 477, 494, 320, 616, 575, 373};// 2)D0D*0bar
137 double y3[24]={0., 1140, 1370, 1660, 1920, 1600, 1115, 700, 675, 626, 636, 605, 515, 540, 748, 880, 556, 657, 477, 494, 320, 616, 575, 373};// 3)D0bar D*0
138 double y4[24]={0., 0, 0, 0, 213, 2000, 2290, 2550, 2443, 2566, 2363, 2173, 1830, 269, 822, 1045, 1020, 820, 398, 588, 690, 459, 360, 398};// 4)D*0 D*0bar
139 double y5[24]={0., 1115, 1375, 1650, 1851, 1650, 1085, 780, 688, 688, 642, 648, 535, 511, 748, 880, 556, 657, 477, 494, 320, 616, 575, 373};// 5)D*+D-
140 double y6[24]={0., 1115, 1375, 1650, 1851, 1650, 1085, 780, 688, 688, 642, 648, 535, 511, 748, 880, 556, 657, 477, 494, 320, 616, 575, 373};// 6)D*-D+
141 double y7[24]={0., 0, 0, 0, 0, 1400, 2390, 2280, 2556, 2479, 2357, 2145, 1564, 237, 822, 1045, 1020, 820, 398, 588, 690, 459, 360, 398};// 7)D*+D*-
142 double y8[24]={0., 102, 133, 269, 250, 174, 51, 26, 25, 15, 34, 7, 15, 47, 106, 70, 36, 10, 2, 28, 60, 60, 48, 36}; // 8)Ds+ Ds-
143 double y9[24]={0., 0, 0, 0, 0, 0, 0, 478, 684, 905, 916, 889, 812, 34, 314, 368, 318, 357, 292, 171, 66, 103, 190, 272};// 9)Ds*+ Ds-
144 double y10[24]={0., 0, 0, 0, 0, 0, 0, 478, 684, 905, 916, 889, 812, 34, 314, 368, 318, 357, 292, 171, 66, 103, 190, 272};//10)Ds*- Ds+
145 double y11[24]={0., 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 440, 398, 428, 310, 131, 0, 45, 126, 98, 39, 0}; //11)Ds*+ Ds*-
146 double y12[24]={0., 0, 0, 0, 0, 0, 12, 3.8, 32.4, 32.4, 37, 48, 61.2, 53.6, 49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 };//12)D*+ D- pi0 //------ DD* pi----
147 double y13[24]={0., 0, 0, 0, 0, 0, 12, 3.8, 32.4, 32.4, 37, 48, 61.2, 53.6, 49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 }; //13)D*- D+ pi0
148 double y14[24]={0., 0, 0, 0, 0, 0, 24, 7.5, 68.6, 64.8, 73.3, 95.8, 122.5,106.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110}; //14)D*+ anti-D0 pi-
149 double y15[24]={0., 0, 0, 0, 0, 0, 24, 7.5, 68.6, 64.8, 73.3, 95.8, 122.5,106.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110}; //15)D*- D0 pi+
150 double y16[24]={0., 0, 0, 0, 0, 0, 24, 7.5, 68.6, 64.8, 73.3, 95.8, 122.5,106.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110}; //16)D+ anti-D*0 pi-
151 double y17[24]={0., 0, 0, 0, 0, 0, 24, 7.5, 68.6, 64.8, 73.3, 95.8, 122.5,106.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110}; //17)D- D*0 pi+
152 double y18[24]={0., 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 53.6,49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 }; //18) D*+ D*- pi0 //------D*D*pi, above 4.26Gev, assumed xs as D*D pi
153 double y19[24]={0., 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,107.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110 }; //19) anti-D*0 D*+ pi-
154 double y20[24]={0., 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,107.3,98, 378, 650, 564, 167, 304, 359, 381, 497, 110 }; //20) D*0 D*- pi+
155 double y21[24]={0., 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 53.6,49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 }; //21) D*0 D*0bar pi0
156 double y22[24]={0., 0, 0, 0, 0, 0, 12, 3.8, 32.4, 32.4, 37, 48, 61.2, 53.6,49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 }; //22)D*0 D0bar pi0 //------ DD* pi----
157 double y23[24]={0., 0, 0, 0, 0, 0, 12, 3.8, 32.4, 32.4, 37, 48, 61.2, 53.6,49, 169, 325, 282, 83.5, 152, 179, 190, 248, 55 }; //23)D*bar D0 pi0
158
159 d0d0bar.clear();
160 dpdm.clear();
161 d0dst0bar.clear();
162 dst0dst0bar.clear();
163 d0bardst0.clear();
164 dstpdm.clear();
165 dstmdp.clear();
166 dstpdstm.clear();
167 dspdsm.clear();
168 dsspdsm.clear();
169 dssmdsp.clear();
170 dsspdssm.clear();
171 xs12.clear();
172 xs13.clear();
173 xs14.clear();
174 xs15.clear();
175 xs16.clear();
176 xs17.clear();
177 xs18.clear();
178 xs19.clear();
179 xs20.clear();
180 xs21.clear();
181 xs22.clear();
182 xs23.clear();
183 for(int i=0;i<24;i++){
184 x.push_back(xx[i]);
185 d0d0bar.push_back(y0[i]);
186 dpdm.push_back(y1[i]);
187 d0dst0bar.push_back(y2[i]);
188 d0bardst0.push_back(y3[i]);
189 dst0dst0bar.push_back(y4[i]);
190 dstpdm.push_back( y5[i]);
191 dstmdp.push_back( y6[i]);
192 dstpdstm.push_back(y7[i]);
193 dspdsm.push_back( y8[i]);
194 dsspdsm.push_back( y9[i]);
195 dssmdsp.push_back( y10[i]);
196 dsspdssm.push_back( y11[i]);
197 xs12.push_back( y12[i] );
198 xs13.push_back( y13[i] );
199 xs14.push_back( y14[i] );
200 xs15.push_back( y15[i] );
201 xs16.push_back( y16[i] );
202 xs17.push_back( y17[i] );
203 xs18.push_back( y18[i] );
204 xs19.push_back( y19[i] );
205 xs20.push_back( y20[i] );
206 xs21.push_back( y21[i] );
207 xs22.push_back( y22[i] );
208 xs23.push_back( y23[i] );
209 }
210
211//---- initilize Vmode
212
213 VmodeSon.clear();
214 //0: D0 anti-D0
215 Vson.clear();
216 Vson.push_back("D0"); Vson.push_back("anti-D0");
217 VmodeSon.push_back(Vson);
218
219 //1: D+ D-
220 Vson.clear();
221 Vson.push_back("D+"); Vson.push_back("D-");
222 VmodeSon.push_back(Vson);
223
224 //2: D0 anti-D*0
225 Vson.clear();
226 Vson.push_back("D0"); Vson.push_back("anti-D*0");
227 VmodeSon.push_back(Vson);
228
229 //3: anti-D0 D*0
230 Vson.clear();
231 Vson.push_back("anti-D0"); Vson.push_back("D*0");
232 VmodeSon.push_back(Vson);
233
234 //4: D*0 anti-D*0
235 Vson.clear();
236 Vson.push_back("D*0"); Vson.push_back("anti-D*0");
237 VmodeSon.push_back(Vson);
238
239 //5: D*+ D-
240 Vson.clear();
241 Vson.push_back("D*+"); Vson.push_back("D-");
242 VmodeSon.push_back(Vson);
243
244 //6: D*- D+
245 Vson.clear();
246 Vson.push_back("D*-"); Vson.push_back("D+");
247 VmodeSon.push_back(Vson);
248
249 //7: D*+ D*-
250 Vson.clear();
251 Vson.push_back("D*+"); Vson.push_back("D*-");
252 VmodeSon.push_back(Vson);
253
254 //8: D_s+ D_s-
255 Vson.clear();
256 Vson.push_back("D_s+"); Vson.push_back("D_s-");
257 VmodeSon.push_back(Vson);
258
259 //9: D_s*+ D_s-
260 Vson.clear();
261 Vson.push_back("D_s*+"); Vson.push_back("D_s-");
262 VmodeSon.push_back(Vson);
263
264 //10: D_s*- D_s+
265 Vson.clear();
266 Vson.push_back("D_s*-"); Vson.push_back("D_s+");
267 VmodeSon.push_back(Vson);
268
269 //11: D_s*+ D_s*-
270 Vson.clear();
271 Vson.push_back("D_s*+"); Vson.push_back("D_s*-");
272 VmodeSon.push_back(Vson);
273
274 //12: D*+ D- pi0
275 Vson.clear();
276 Vson.push_back("D*+"); Vson.push_back("D-");Vson.push_back("pi0");
277 VmodeSon.push_back(Vson);
278
279 //13: D*- D+ pi0
280 Vson.clear();
281 Vson.push_back("D*-"); Vson.push_back("D+");Vson.push_back("pi0");
282 VmodeSon.push_back(Vson);
283
284 //14: D*+ anti-D0 pi-
285 Vson.clear();
286 Vson.push_back("D*+"); Vson.push_back("anti-D0");Vson.push_back("pi-");
287 VmodeSon.push_back(Vson);
288
289 //15: D*- D0 pi+
290 Vson.clear();
291 Vson.push_back("D*-"); Vson.push_back("D0");Vson.push_back("pi+");
292 VmodeSon.push_back(Vson);
293
294 //16: D+ anti-D*0 pi-
295 Vson.clear();
296 Vson.push_back("D+"); Vson.push_back("anti-D*0");Vson.push_back("pi-");
297 VmodeSon.push_back(Vson);
298
299 //17: D- D*0 pi+
300 Vson.clear();
301 Vson.push_back("D-"); Vson.push_back("D*0");Vson.push_back("pi+");
302 VmodeSon.push_back(Vson);
303
304 //18: D*+ D*- pi0
305 Vson.clear();
306 Vson.push_back("D*+"); Vson.push_back("D*-");Vson.push_back("pi0");
307 VmodeSon.push_back(Vson);
308
309 //19: anti-D*0 D*+ pi-
310 Vson.clear();
311 Vson.push_back("anti-D*0"); Vson.push_back("D*+");Vson.push_back("pi-");
312 VmodeSon.push_back(Vson);
313
314 //20: D*0 D*- pi+
315 Vson.clear();
316 Vson.push_back("D*0"); Vson.push_back("D*-");Vson.push_back("pi+");
317 VmodeSon.push_back(Vson);
318
319 //21: D*0 D*0bar pi0
320 Vson.clear();
321 Vson.push_back("D*0"); Vson.push_back("anti-D*0");Vson.push_back("pi0");
322 VmodeSon.push_back(Vson);
323
324 //22: D0bar D*0 pi0
325 Vson.clear();
326 Vson.push_back("anti-D0"); Vson.push_back("D*0");Vson.push_back("pi0");
327 VmodeSon.push_back(Vson);
328
329 //23: D*0bar D0 pi0
330 Vson.clear();
331 Vson.push_back("anti-D*0"); Vson.push_back("D0");Vson.push_back("pi0");
332 VmodeSon.push_back(Vson);
333
334 }
335
336 // destructor
337 virtual ~EvtPsi3Sdecay(){}
338
339 int findPoints(); //return the lower size at the end of the segment containing Ecms points
340 double polint( std::vector <double> points);
341 bool choseDecay(); //determing accept or reject a generated decay
343 int getDecay(double ecms); //determing accept or reject a generated decay
344 double theProb(std::vector<double> myxs,int ich);
345 int findMode( );
346 int getMode(){return _themode;};
347 std::vector<EvtId> getVId(int mode);
348
349
350 void PHSPDecay(EvtParticle * par);
351 std::vector<EvtId> getDaugId(){return Vid;}
352 std::vector<EvtVector4R> getDaugP4(){return v_p4;}
353 bool AngSam(EvtVector4R parent_p4cm,EvtVector4R son_p4cm,double alpha);
354 void setMode(int m){
355 if(m<0 || m>24) {std::cout<<"EvtPsi3Decay::setMode: bad mode"<<std::endl;abort();}
356 _themode = m;_excflag=1;}
357private:
358 double Ecms;
359 EvtParticle *theParent;
360 EvtId *Daughters;
361 int Ndaugs;
362 std::vector <double> x;
363 std::vector <double> d0d0bar; // 0)
364 std::vector <double> dpdm; // 1)
365 std::vector <double> d0dst0bar; // 2)
366 std::vector <double> d0bardst0; // 3)
367 std::vector <double> dst0dst0bar; // 4)
368 std::vector <double> dstpdm; // 5)
369 std::vector <double> dstmdp; // 6)
370 std::vector <double> dstpdstm; // 7)
371 std::vector <double> dspdsm; // 8)
372 std::vector <double> dsspdsm; // 9)
373 std::vector <double> dssmdsp; //10)
374 std::vector <double> dsspdssm; //11)
375 std::vector <double> xs12; //12)
376 std::vector <double> xs13; //13)
377 std::vector <double> xs14; //14)
378 std::vector <double> xs15; //15)
379 std::vector <double> xs16; //16)
380 std::vector <double> xs17; //17)
381 std::vector <double> xs18; //18)
382 std::vector <double> xs19; //19)
383 std::vector <double> xs20; //20)
384 std::vector <double> xs21; //21)
385 std::vector <double> xs22; //22)
386 std::vector <double> xs23; //23)
387
388 int nsize; //the simze of the energy points
389 double theXsection[50];
390 int theLocation;
391 std::vector<std::string> Vson;
392 std::vector<EvtId> Vid;
393
394 std::vector< std::vector<std::string> > VmodeSon;
395 std::vector< std::vector<EvtId> > VmodeId;
396 int _themode,_excflag;
397 //
398 std::vector<EvtVector4R> v_p4;
399
400};
401
402
403#endif
const double alpha
Definition: EvtId.hh:27
int getNDaug() const
Definition: EvtParticle.cc:125
virtual ~EvtPsi3Sdecay()
int getDecay(double ecms)
double theProb(std::vector< double > myxs, int ich)
void PHSPDecay(EvtParticle *par)
bool AngSam(EvtVector4R parent_p4cm, EvtVector4R son_p4cm, double alpha)
std::vector< EvtId > getVId(int mode)
void setMode(int m)
double polint(std::vector< double > points)
EvtPsi3Sdecay(double ecms, EvtParticle *parent)
std::vector< EvtVector4R > getDaugP4()
std::vector< EvtId > getDaugId()
const double ecms
Definition: inclkstar.cxx:42