Belle II Software
development
GenTauTag.cc
1
/**************************************************************************
2
* basf2 (Belle II Analysis Software Framework) *
3
* Author: The Belle II Collaboration *
4
* *
5
* See git log for contributors and copyright holders. *
6
* This file is licensed under LGPL-3.0, see LICENSE.md. *
7
**************************************************************************/
8
9
#include <analysis/utility/GenTauTag.h>
10
#include <analysis/utility/CheckDecayUtils.h>
11
12
using namespace
Belle2
;
13
14
GenTauTag::GenTauTag
() =
default
;
15
16
bool
GenTauTag::PcheckDecay
(std::vector<int> gp,
int
da1,
int
da2)
17
{
18
int
ndau = gp.size();
19
if
(ndau == 2) {
m_nPhotos
= 0;
return
CheckDecay::CheckDecay(gp, da1, da2);}
20
else
if
(ndau == 3) {
m_nPhotos
= 1;
return
CheckDecay::CheckDecay(gp, da1, da2, 22);}
21
else
if
(ndau == 4) {
m_nPhotos
= 2;
return
CheckDecay::CheckDecay(gp, da1, da2, 22, 22);}
22
else
if
(ndau == 5) {
m_nPhotos
= 3;
return
CheckDecay::CheckDecay(gp, da1, da2, 22, 22, 22);}
23
else
if
(ndau == 6) {
m_nPhotos
= 4;
return
CheckDecay::CheckDecay(gp, da1, da2, 22, 22, 22, 22);}
24
else
return
false
;
25
}
26
27
bool
GenTauTag::PcheckDecay
(std::vector<int> gp,
int
da1,
int
da2,
int
da3)
28
{
29
int
ndau = gp.size();
30
if
(ndau == 3) {
m_nPhotos
= 0;
return
CheckDecay::CheckDecay(gp, da1, da2, da3);}
31
else
if
(ndau == 4) {
m_nPhotos
= 1;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, 22);}
32
else
if
(ndau == 5) {
m_nPhotos
= 2;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, 22, 22);}
33
else
if
(ndau == 6) {
m_nPhotos
= 3;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, 22, 22, 22);}
34
else
if
(ndau == 7) {
m_nPhotos
= 4;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, 22, 22, 22, 22);}
35
else
return
false
;
36
}
37
38
bool
GenTauTag::PcheckDecay
(std::vector<int> gp,
int
da1,
int
da2,
int
da3,
int
da4)
39
{
40
int
ndau = gp.size();
41
if
(ndau == 4) {
m_nPhotos
= 0;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4);}
42
else
if
(ndau == 5) {
m_nPhotos
= 1;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, 22);}
43
else
if
(ndau == 6) {
m_nPhotos
= 2;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, 22, 22);}
44
else
if
(ndau == 7) {
m_nPhotos
= 3;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, 22, 22, 22);}
45
else
if
(ndau == 8) {
m_nPhotos
= 4;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, 22, 22, 22, 22);}
46
else
return
false
;
47
}
48
49
bool
GenTauTag::PcheckDecay
(std::vector<int> gp,
int
da1,
int
da2,
int
da3,
int
da4,
int
da5)
50
{
51
int
ndau = gp.size();
52
if
(ndau == 5) {
m_nPhotos
= 0;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5);}
53
else
if
(ndau == 6) {
m_nPhotos
= 1;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, 22);}
54
else
if
(ndau == 7) {
m_nPhotos
= 2;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, 22, 22);}
55
else
if
(ndau == 8) {
m_nPhotos
= 3;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, 22, 22, 22);}
56
else
if
(ndau == 9) {
m_nPhotos
= 4;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, 22, 22, 22, 22);}
57
else
return
false
;
58
}
59
60
bool
GenTauTag::PcheckDecay
(std::vector<int> gp,
int
da1,
int
da2,
int
da3,
int
da4,
int
da5,
int
da6)
61
{
62
int
ndau = gp.size();
63
if
(ndau == 6) {
m_nPhotos
= 0;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6);}
64
else
if
(ndau == 7) {
m_nPhotos
= 1;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, 22);}
65
else
if
(ndau == 8) {
m_nPhotos
= 2;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, 22, 22);}
66
else
if
(ndau == 9) {
m_nPhotos
= 3;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, 22, 22, 22);}
67
else
if
(ndau == 10) {
m_nPhotos
= 4;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, 22, 22, 22, 22);}
68
else
return
false
;
69
}
70
71
bool
GenTauTag::PcheckDecay
(std::vector<int> gp,
int
da1,
int
da2,
int
da3,
int
da4,
int
da5,
int
da6,
int
da7)
72
{
73
int
ndau = gp.size();
74
if
(ndau == 7) {
m_nPhotos
= 0;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, da7);}
75
else
if
(ndau == 8) {
m_nPhotos
= 1;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, da7, 22);}
76
else
if
(ndau == 9) {
m_nPhotos
= 2;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, da7, 22, 22);}
77
else
if
(ndau == 10) {
m_nPhotos
= 3;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, da7, 22, 22, 22);}
78
else
if
(ndau == 11) {
m_nPhotos
= 4;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, da7, 22, 22, 22, 22);}
79
else
return
false
;
80
}
81
82
bool
GenTauTag::PcheckDecay
(std::vector<int> gp,
int
da1,
int
da2,
int
da3,
int
da4,
int
da5,
int
da6,
int
da7,
int
da8)
83
{
84
int
ndau = gp.size();
85
if
(ndau == 8) {
m_nPhotos
= 0;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, da7, da8);}
86
else
if
(ndau == 9) {
m_nPhotos
= 1;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, da7, da8, 22);}
87
else
if
(ndau == 10) {
m_nPhotos
= 2;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, da7, da8, 22, 22);}
88
else
if
(ndau == 11) {
m_nPhotos
= 3;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, da7, da8, 22, 22, 22);}
89
else
if
(ndau == 12) {
m_nPhotos
= 4;
return
CheckDecay::CheckDecay(gp, da1, da2, da3, da4, da5, da6, da7, da8, 22, 22, 22, 22);}
90
else
return
false
;
91
}
92
93
//All the tag number for various decay mode is defined here
94
95
int
GenTauTag::Mode_tau_plus
(std::vector<int> genpart)
96
{
97
if
(
GenTauTag::PcheckDecay
(genpart, -11, 12, -16)) {
98
return
1000 *
m_nPhotos
+ 1;
99
}
//tau+ decays to e+ nu_e anti-nu_tau
100
if
(
GenTauTag::PcheckDecay
(genpart, -13, 14, -16)) {
101
return
1000 *
m_nPhotos
+ 2;
102
}
//tau+ decays to mu+ nu_mu anti-nu_tau
103
if
(
GenTauTag::PcheckDecay
(genpart, 211, -16)) {
104
return
1000 *
m_nPhotos
+ 3;
105
}
//tau+ decays to pi+ anti-nu_tau
106
if
(
GenTauTag::PcheckDecay
(genpart, 213, -16)) {
107
return
1000 *
m_nPhotos
+ 4;
108
}
//tau+ decays to rho+ anti-nu_tau
109
if
(
GenTauTag::PcheckDecay
(genpart, 20213, -16)) {
110
return
1000 *
m_nPhotos
+ 5;
111
}
//tau+ decays to a_1+ anti-nu_tau
112
if
(
GenTauTag::PcheckDecay
(genpart, 321, -16)) {
113
return
1000 *
m_nPhotos
+ 6;
114
}
//tau+ decays to K+ anti-nu_tau
115
if
(
GenTauTag::PcheckDecay
(genpart, 323, -16)) {
116
return
1000 *
m_nPhotos
+ 7;
117
}
//tau- decays to K*+ anti-nu_tau
118
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, -211, 211, 111)) {
119
return
1000 *
m_nPhotos
+ 8;
120
}
//tau- decays to anti-nu_tau pi+ pi- pi+ pi0
121
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 111, 111, 111)) {
122
return
1000 *
m_nPhotos
+ 9;
123
}
//tau- decays to anti-nu_tau pi+ pi0 pi0 pi0
124
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 211, -211, 111, 111)) {
125
return
1000 *
m_nPhotos
+ 10;
126
}
//tau- decays to anti-nu_tau pi+ pi+ pi- pi0 pi0
127
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 211, 211, -211, -211)) {
128
return
1000 *
m_nPhotos
+ 11;
129
}
//tau- decays to anti-nu_tau pi+ pi+ pi+ pi- pi-
130
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 211, 211, -211, -211, 111)) {
131
return
1000 *
m_nPhotos
+ 12;
132
}
//tau- decays to anti-nu_tau pi+ pi+ pi+ pi- pi- pi0
133
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 211, -211, 111, 111, 111)) {
134
return
1000 *
m_nPhotos
+ 13;
135
}
//tau- decays to anti-nu_tau pi+ pi+ pi- pi0 pi0 pi0
136
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, 211, -321)) {
137
return
1000 *
m_nPhotos
+ 14;
138
}
//tau- decays to anti-nu_tau K+ pi+ K-
139
if
(
GenTauTag::PcheckDecay
(genpart, -16, -10313, +211, 10313)) {
140
return
1000 *
m_nPhotos
+ 15;
141
}
//tau- decays to anti-nu_tau anti-K_10 pi+ K_10
142
if
(
GenTauTag::PcheckDecay
(genpart, -16, 311, 211, -311)) {
143
return
1000 *
m_nPhotos
+ 15;
144
}
//tau- decays to anti-nu_tau K0 pi+ anti-K0
145
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, -10313, 111)) {
146
return
1000 *
m_nPhotos
+ 16;
147
}
//tau- decays to anti-nu_tau K+ anti-K_10 pi0
148
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, 311, 111)) {
149
return
1000 *
m_nPhotos
+ 16;
150
}
//tau- decays to anti-nu_tau K+ K0 pi0
151
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, -311, 111)) {
152
return
1000 *
m_nPhotos
+ 16;
153
}
//tau- decays to anti-nu_tau K+ anti-K0 pi0
154
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, 111, 111)) {
155
return
1000 *
m_nPhotos
+ 17;
156
}
//tau- decays to anti-nu_tau K+ pi0 pi0
157
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, -211, 211)) {
158
return
1000 *
m_nPhotos
+ 18;
159
}
//tau- decays to anti-nu_tau K+ pi- pi+
160
if
(
GenTauTag::PcheckDecay
(genpart, -16, 10313, 211, 111)) {
161
return
1000 *
m_nPhotos
+ 19;
162
}
//tau- decays to anti-nu_tau K_10 pi+ pi0
163
if
(
GenTauTag::PcheckDecay
(genpart, -16, 311, 211, 111)) {
164
return
1000 *
m_nPhotos
+ 19;
165
}
//tau- decays to anti-nu_tau K0 pi+ pi0
166
if
(
GenTauTag::PcheckDecay
(genpart, -16, -311, 211, 111)) {
167
return
1000 *
m_nPhotos
+ 19;
168
}
//tau- decays to anti-nu_tau anti-K0 pi+ pi0
169
if
(
GenTauTag::PcheckDecay
(genpart, -16, 221, 211, 111)) {
170
return
1000 *
m_nPhotos
+ 20;
171
}
//tau- decays to anti-nu_tau eta pi+ pi0
172
if
(
GenTauTag::PcheckDecay
(genpart, 211, 111, -16)) {
173
return
1000 *
m_nPhotos
+ 21;
174
}
//tau+ decays to pi+ pi0 anti-nu_tau
175
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, -10313)) {
176
return
1000 *
m_nPhotos
+ 22;
177
}
//tau- decays to anti-nu_tau K+ anti-K_10
178
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, 311)) {
179
return
1000 *
m_nPhotos
+ 22;
180
}
//tau- decays to anti-nu_tau K+ K0
181
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, -311)) {
182
return
1000 *
m_nPhotos
+ 22;
183
}
//tau- decays to anti-nu_tau K+ anti-K0
184
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 111, 111, 111, 111)) {
185
return
1000 *
m_nPhotos
+ 23;
186
}
//tau- decays to anti-nu_tau pi+ pi0 pi0 pi0 pi0
187
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 223, 111)) {
188
return
1000 *
m_nPhotos
+ 24;
189
}
//tau- decays to anti-nu_tau pi+ omega pi0
190
if
(
GenTauTag::PcheckDecay
(genpart, 221, 211, -211, 211, -16)) {
191
return
1000 *
m_nPhotos
+ 25;
192
}
//tau- decays to eta pi+ pi- pi+ anti-nu_tau
193
if
(
GenTauTag::PcheckDecay
(genpart, 221, 211, 111, 111, -16)) {
194
return
1000 *
m_nPhotos
+ 26;
195
}
//tau- decays to eta pi+ pi0 pi0 anti-nu_tau
196
if
(
GenTauTag::PcheckDecay
(genpart, 221, 321, -16)) {
197
return
1000 *
m_nPhotos
+ 27;
198
}
//tau- decays to eta K+ anti-nu_tau
199
if
(
GenTauTag::PcheckDecay
(genpart, 221, 323, -16)) {
200
return
1000 *
m_nPhotos
+ 28;
201
}
//tau- decays to eta K*+ anti-nu_tau
202
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, -211, 211, 111)) {
203
return
1000 *
m_nPhotos
+ 29;
204
}
//tau- decays to anti-nu_tau K+ pi- pi+ pi0
205
if
(
GenTauTag::PcheckDecay
(genpart, 321, 111, 111, 111, -16)) {
206
return
1000 *
m_nPhotos
+ 30;
207
}
//tau- decays to K+ pi0 pi0 pi0 anti-nu_tau
208
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, -211, 211, -10313)) {
209
return
1000 *
m_nPhotos
+ 31;
210
}
//tau- decays to anti-nu_tau pi+ pi- pi+ anti-K_10
211
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, -211, 211, 311)) {
212
return
1000 *
m_nPhotos
+ 31;
213
}
//tau- decays to anti-nu_tau pi+ pi- pi+ K0
214
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, -211, 211, -311)) {
215
return
1000 *
m_nPhotos
+ 31;
216
}
//tau- decays to anti-nu_tau pi+ pi- pi+ anti-K0
217
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 111, 111, 10313)) {
218
return
1000 *
m_nPhotos
+ 32;
219
}
//tau- decays to anti-nu_tau pi+ pi0 pi0 K_10
220
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 111, 111, 311)) {
221
return
1000 *
m_nPhotos
+ 32;
222
}
//tau- decays to anti-nu_tau pi+ pi0 pi0 K0
223
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 111, 111, -311)) {
224
return
1000 *
m_nPhotos
+ 32;
225
}
//tau- decays to anti-nu_tau pi+ pi0 pi0 anti-K0
226
if
(
GenTauTag::PcheckDecay
(genpart, 321, -321, 211, 111, -16)) {
227
return
1000 *
m_nPhotos
+ 33;
228
}
//tau- decays to K+ K- pi+ pi0 anti-nu_tau
229
if
(
GenTauTag::PcheckDecay
(genpart, -16, -10313, 211, 111, 10313)) {
230
return
1000 *
m_nPhotos
+ 34;
231
}
//tau- decays to anti-nu_tau anti-K_10 pi+ pi0 K_10
232
if
(
GenTauTag::PcheckDecay
(genpart, -16, 311, 211, 111, -311)) {
233
return
1000 *
m_nPhotos
+ 34;
234
}
//tau- decays to anti-nu_tau K0 pi+ pi0 anti-K0
235
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 223, -211, 211)) {
236
return
1000 *
m_nPhotos
+ 35;
237
}
//tau- decays to anti-nu_tau pi+ omega pi- pi+
238
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 223, 111, 111)) {
239
return
1000 *
m_nPhotos
+ 36;
240
}
//tau- decays to anti-nu_tau pi+ omega pi0 pi0
241
if
(
GenTauTag::PcheckDecay
(genpart, -11, -11, 11, 12, -16)) {
242
return
1000 *
m_nPhotos
+ 37;
243
}
//tau- decays to e+ e+ e- nu_e anti-nu_tau
244
if
(
GenTauTag::PcheckDecay
(genpart, 20223, 211, -16)) {
245
return
1000 *
m_nPhotos
+ 38;
246
}
//tau- decays to f_1 pi+ anti-nu_tau
247
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, 223)) {
248
return
1000 *
m_nPhotos
+ 39;
249
}
//tau- decays to anti-nu_tau K+ omega
250
if
(
GenTauTag::PcheckDecay
(genpart, -16, -10313, 211, -211, 321)) {
251
return
1000 *
m_nPhotos
+ 40;
252
}
//tau- decays to anti-nu_tau anti-K_10 pi+ pi- K+
253
if
(
GenTauTag::PcheckDecay
(genpart, -16, 311, 211, -211, 321)) {
254
return
1000 *
m_nPhotos
+ 40;
255
}
//tau- decays to anti-nu_tau K0 pi+ pi- K+
256
if
(
GenTauTag::PcheckDecay
(genpart, -16, -311, 211, -211, 321)) {
257
return
1000 *
m_nPhotos
+ 40;
258
}
//tau- decays to anti-nu_tau anti-K0 pi+ pi- K+
259
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, 111, 111, -10313)) {
260
return
1000 *
m_nPhotos
+ 41;
261
}
//tau- decays to anti-nu_tau K+ pi0 pi0 anti-K_10
262
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, 111, 111, 311)) {
263
return
1000 *
m_nPhotos
+ 41;
264
}
//tau- decays to anti-nu_tau K+ pi0 pi0 K0
265
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, 111, 111, -311)) {
266
return
1000 *
m_nPhotos
+ 41;
267
}
//tau- decays to anti-nu_tau K+ pi0 pi0 anti-K0
268
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 211, -321, 10313)) {
269
return
1000 *
m_nPhotos
+ 42;
270
}
//tau- decays to anti-nu_tau pi+ pi+ K- K_10
271
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 211, -321, 311)) {
272
return
1000 *
m_nPhotos
+ 42;
273
}
//tau- decays to anti-nu_tau pi+ pi+ K- K0
274
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 211, -321, -311)) {
275
return
1000 *
m_nPhotos
+ 42;
276
}
//tau- decays to anti-nu_tau pi+ pi+ K- anti-K0
277
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 311)) {
278
return
1000 *
m_nPhotos
+ 43;
279
}
//tau- decays to anti-nu_tau pi+ K0
280
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, -311)) {
281
return
1000 *
m_nPhotos
+ 43;
282
}
//tau- decays to anti-nu_tau pi+ anti-K0
283
if
(
GenTauTag::PcheckDecay
(genpart, -16, 22, 211, 111)) {
284
return
1000 *
m_nPhotos
+ 44;
285
}
//tau- decays to anti-nu_tau gamma pi+ pi0
286
if
(
GenTauTag::PcheckDecay
(genpart, 111, 111, 211, -16)) {
287
return
1000 *
m_nPhotos
+ 45;
288
}
//tau- decays to pi0 pi0 pi+ anti-nu_tau
289
if
(
GenTauTag::PcheckDecay
(genpart, 211, 211, -211, -16)) {
290
return
1000 *
m_nPhotos
+ 46;
291
}
//tau- decays to pi+ pi+ pi- anti-nu_tau
292
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, -321, 321)) {
293
return
1000 *
m_nPhotos
+ 47;
294
}
//tau- decays to anti-nu_tau K+ K- K+
295
if
(
GenTauTag::PcheckDecay
(genpart, 313, 211, -16)) {
296
return
1000 *
m_nPhotos
+ 48;
297
}
//tau- decays to K*0 pi+ anti-nu_tau
298
if
(
GenTauTag::PcheckDecay
(genpart, -313, 321, -16)) {
299
return
1000 *
m_nPhotos
+ 49;
300
}
//tau- decays to anti-K*0 K+ anti-nu_tau
301
if
(
GenTauTag::PcheckDecay
(genpart, 333, 211, -16)) {
302
return
1000 *
m_nPhotos
+ 50;
303
}
//tau- decays to phi pi+ anti-nu_tau
304
if
(
GenTauTag::PcheckDecay
(genpart, 333, 321, -16)) {
305
return
1000 *
m_nPhotos
+ 51;
306
}
//tau- decays to phi K+ anti-nu_tau
307
if
(
GenTauTag::PcheckDecay
(genpart, -13, 14, -16, 22)) {
308
return
1000 *
m_nPhotos
+ 52;
309
}
//tau- decays to mu+ nu_mu anti-nu_tau gamma
310
if
(
GenTauTag::PcheckDecay
(genpart, -11, 12, -16, 22)) {
311
return
1000 *
m_nPhotos
+ 53;
312
}
//tau- decays to e+ nu_e anti-nu_tau gamma
313
if
(
GenTauTag::PcheckDecay
(genpart, 321, 111, -16)) {
314
return
1000 *
m_nPhotos
+ 54;
315
}
//tau- decays to K+ pi0 anti-nu_tau
316
if
(
GenTauTag::PcheckDecay
(genpart, 311, 213, -16)) {
317
return
1000 *
m_nPhotos
+ 55;
318
}
//tau- decays to K0 rho+ anti-nu_tau
319
if
(
GenTauTag::PcheckDecay
(genpart, -311, 213, -16)) {
320
return
1000 *
m_nPhotos
+ 55;
321
}
//tau- decays to anti-K0 rho+ anti-nu_tau
322
if
(
GenTauTag::PcheckDecay
(genpart, 10323, -16)) {
323
return
1000 *
m_nPhotos
+ 56;
324
}
//tau- decays to K_1+ anti-nu_tau
325
if
(
GenTauTag::PcheckDecay
(genpart, 20323, -16)) {
326
return
1000 *
m_nPhotos
+ 57;
327
}
//tau- decays to K'_1+ anti-nu_tau
328
if
(
GenTauTag::PcheckDecay
(genpart, 100323, -16)) {
329
return
1000 *
m_nPhotos
+ 58;
330
}
//tau- decays to K'*+ anti-nu_tau
331
if
(
GenTauTag::PcheckDecay
(genpart, 221, 321, 111, -16)) {
332
return
1000 *
m_nPhotos
+ 59;
333
}
//tau- decays to eta K+ pi0 anti-nu_tau
334
if
(
GenTauTag::PcheckDecay
(genpart, 221, 311, 211, -16)) {
335
return
1000 *
m_nPhotos
+ 60;
336
}
//tau- decays to eta K0 pi+ anti-nu_tau
337
if
(
GenTauTag::PcheckDecay
(genpart, 221, -311, 211, -16)) {
338
return
1000 *
m_nPhotos
+ 60;
339
}
//tau- decays to eta anti-K0 pi+ anti-nu_tau
340
if
(
GenTauTag::PcheckDecay
(genpart, -16, 310, 310, 211)) {
341
return
1000 *
m_nPhotos
+ 20015;
342
}
//tau- decays to anti-nu_tau K_S0 K_S0 pi+
343
if
(
GenTauTag::PcheckDecay
(genpart, -16, 310, 130, 211)) {
344
return
1000 *
m_nPhotos
+ 110015;
345
}
//tau- decays to anti-nu_tau K_S0 K_L0 pi+
346
if
(
GenTauTag::PcheckDecay
(genpart, -16, 130, 130, 211)) {
347
return
1000 *
m_nPhotos
+ 200015;
348
}
//tau- decays to anti-nu_tau K_L0 K_L0 pi+
349
if
(
GenTauTag::PcheckDecay
(genpart, -16, 310, 321, 111)) {
350
return
1000 *
m_nPhotos
+ 10016;
351
}
//tau- decays to anti-nu_tau K_S0 K+ pi0
352
if
(
GenTauTag::PcheckDecay
(genpart, -16, 130, 321, 111)) {
353
return
1000 *
m_nPhotos
+ 100016;
354
}
//tau- decays to anti-nu_tau K_L0 K+ pi0
355
if
(
GenTauTag::PcheckDecay
(genpart, -16, 310, 211, 111)) {
356
return
1000 *
m_nPhotos
+ 10019;
357
}
//tau- decays to anti-nu_tau K_S0 pi+ pi0
358
if
(
GenTauTag::PcheckDecay
(genpart, -16, 130, 211, 111)) {
359
return
1000 *
m_nPhotos
+ 100019;
360
}
//tau- decays to anti-nu_tau K_L0 pi+ pi0
361
if
(
GenTauTag::PcheckDecay
(genpart, -16, 310, 321)) {
362
return
1000 *
m_nPhotos
+ 10022;
363
}
//tau- decays to anti-nu_tau K_S0 K+
364
if
(
GenTauTag::PcheckDecay
(genpart, -16, 130, 321)) {
365
return
1000 *
m_nPhotos
+ 100022;
366
}
//tau- decays to anti-nu_tau K_L0 K+
367
if
(
GenTauTag::PcheckDecay
(genpart, -16, 310, 211, -211, 211)) {
368
return
1000 *
m_nPhotos
+ 10031;
369
}
//tau- decays to anti-nu_tau K_S0 pi+ pi- pi+
370
if
(
GenTauTag::PcheckDecay
(genpart, -16, 130, 211, -211, 211)) {
371
return
1000 *
m_nPhotos
+ 100031;
372
}
//tau- decays to anti-nu_tau K_L0 pi+ pi- pi+
373
if
(
GenTauTag::PcheckDecay
(genpart, -16, 310, 211, 111, 111)) {
374
return
1000 *
m_nPhotos
+ 10032;
375
}
//tau- decays to anti-nu_tau K_S0 pi+ pi0 pi0
376
if
(
GenTauTag::PcheckDecay
(genpart, -16, 130, 211, 111, 111)) {
377
return
1000 *
m_nPhotos
+ 100032;
378
}
//tau- decays to anti-nu_tau K_L0 pi+ pi0 pi0
379
if
(
GenTauTag::PcheckDecay
(genpart, -16, 310, 310, 211, 111)) {
380
return
1000 *
m_nPhotos
+ 20034;
381
}
//tau- decays to anti-nu_tau K_S0 K_S0 pi+ pi0
382
if
(
GenTauTag::PcheckDecay
(genpart, -16, 310, 130, 211, 111)) {
383
return
1000 *
m_nPhotos
+ 110034;
384
}
//tau- decays to anti-nu_tau K_S0 K_L0 pi+ pi0
385
if
(
GenTauTag::PcheckDecay
(genpart, -16, 130, 130, 211, 111)) {
386
return
1000 *
m_nPhotos
+ 200034;
387
}
//tau- decays to anti-nu_tau K_L0 K_L0 pi+ pi0
388
if
(
GenTauTag::PcheckDecay
(genpart, -16, 310, 211, -211, 321)) {
389
return
1000 *
m_nPhotos
+ 10040;
390
}
//tau- decays to anti-nu_tau K_S0 pi+ pi- K+
391
if
(
GenTauTag::PcheckDecay
(genpart, -16, 130, 211, -211, 321)) {
392
return
1000 *
m_nPhotos
+ 100040;
393
}
//tau- decays to anti-nu_tau K_L0 pi+ pi- K+
394
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, 111, 111, 310)) {
395
return
1000 *
m_nPhotos
+ 10041;
396
}
//tau- decays to anti-nu_tau K+ pi0 pi0 K_S0
397
if
(
GenTauTag::PcheckDecay
(genpart, -16, 321, 111, 111, 130)) {
398
return
1000 *
m_nPhotos
+ 100041;
399
}
//tau- decays to anti-nu_tau K+ pi0 pi0 K_L0
400
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 211, -321, 311)) {
401
return
1000 *
m_nPhotos
+ 10042;
402
}
//tau- decays to anti-nu_tau pi+ pi+ K- K_S0
403
if
(
GenTauTag::PcheckDecay
(genpart, -16, 211, 211, -321, 130)) {
404
return
1000 *
m_nPhotos
+ 100042;
405
}
//tau- decays to anti-nu_tau pi+ pi+ K- K_L0
406
return
0;
407
}
// not yet defined tau+ decay
408
409
410
int
GenTauTag::Mode_tau_minus
(std::vector<int> genpart)
411
{
412
if
(
GenTauTag::PcheckDecay
(genpart, 11, -12, 16)) {
413
return
1000 *
m_nPhotos
+ 1;
414
}
//tau- decays to e- anti-nu_e nu_tau
415
if
(
GenTauTag::PcheckDecay
(genpart, 13, -14, 16)) {
416
return
1000 *
m_nPhotos
+ 2;
417
}
//tau- decays to mu- anti-nu_mu nu_tau
418
if
(
GenTauTag::PcheckDecay
(genpart, -211, 16)) {
419
return
1000 *
m_nPhotos
+ 3;
420
}
//tau- decays to pi- nu_tau
421
if
(
GenTauTag::PcheckDecay
(genpart, -213, 16)) {
422
return
1000 *
m_nPhotos
+ 4;
423
}
//tau- decays to rho- nu_tau
424
if
(
GenTauTag::PcheckDecay
(genpart, -20213, 16)) {
425
return
1000 *
m_nPhotos
+ 5;
426
}
//tau- decays to a_1- nu_tau
427
if
(
GenTauTag::PcheckDecay
(genpart, -321, 16)) {
428
return
1000 *
m_nPhotos
+ 6;
429
}
//tau- decays to K- nu_tau
430
if
(
GenTauTag::PcheckDecay
(genpart, -323, 16)) {
431
return
1000 *
m_nPhotos
+ 7;
432
}
//tau- decays to K*- nu_tau
433
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 211, -211, 111)) {
434
return
1000 *
m_nPhotos
+ 8;
435
}
//tau- decays to nu_tau pi- pi+ pi- pi0
436
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 111, 111, 111)) {
437
return
1000 *
m_nPhotos
+ 9;
438
}
//tau- decays to nu_tau pi- pi0 pi0 pi0
439
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, -211, 211, 111, 111)) {
440
return
1000 *
m_nPhotos
+ 10;
441
}
//tau- decays to nu_tau pi- pi- pi+ pi0 pi0
442
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, -211, -211, 211, 211)) {
443
return
1000 *
m_nPhotos
+ 11;
444
}
//tau- decays to nu_tau pi- pi- pi- pi+ pi+
445
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, -211, -211, 211, 211, 111)) {
446
return
1000 *
m_nPhotos
+ 12;
447
}
//tau- decays to nu_tau pi- pi- pi- pi+ pi+ pi0
448
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, -211, 211, 111, 111, 111)) {
449
return
1000 *
m_nPhotos
+ 13;
450
}
//tau- decays to nu_tau pi- pi- pi+ pi0 pi0 pi0
451
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, -211, 321)) {
452
return
1000 *
m_nPhotos
+ 14;
453
}
//tau- decays to nu_tau K- pi- K+
454
if
(
GenTauTag::PcheckDecay
(genpart, 16, -10313, -211, 10313)) {
455
return
1000 *
m_nPhotos
+ 15;
456
}
//tau- decays to nu_tau anti-K_10 pi- K_10
457
if
(
GenTauTag::PcheckDecay
(genpart, 16, 311, -211, -311)) {
458
return
1000 *
m_nPhotos
+ 15;
459
}
//tau- decays to nu_tau K0 pi- anti-K0
460
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 10313, 111)) {
461
return
1000 *
m_nPhotos
+ 16;
462
}
//tau- decays to nu_tau K- K_10 pi0
463
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 311, 111)) {
464
return
1000 *
m_nPhotos
+ 16;
465
}
//tau- decays to nu_tau K- K0 pi0
466
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, -311, 111)) {
467
return
1000 *
m_nPhotos
+ 16;
468
}
//tau- decays to nu_tau K- anti-K0 pi0
469
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 111, 111)) {
470
return
1000 *
m_nPhotos
+ 17;
471
}
//tau- decays to nu_tau K- pi0 pi0
472
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 211, -211)) {
473
return
1000 *
m_nPhotos
+ 18;
474
}
//tau- decays to nu_tau K- pi+ pi-
475
if
(
GenTauTag::PcheckDecay
(genpart, 16, -10313, -211, 111)) {
476
return
1000 *
m_nPhotos
+ 19;
477
}
//tau- decays to nu_tau anti-K_10 pi- pi0
478
if
(
GenTauTag::PcheckDecay
(genpart, 16, 311, -211, 111)) {
479
return
1000 *
m_nPhotos
+ 19;
480
}
//tau- decays to nu_tau K0 pi- pi0
481
if
(
GenTauTag::PcheckDecay
(genpart, 16, -311, -211, 111)) {
482
return
1000 *
m_nPhotos
+ 19;
483
}
//tau- decays to nu_tau anti-K0 pi- pi0
484
if
(
GenTauTag::PcheckDecay
(genpart, 16, 221, -211, 111)) {
485
return
1000 *
m_nPhotos
+ 20;
486
}
//tau- decays to nu_tau eta pi- pi0
487
if
(
GenTauTag::PcheckDecay
(genpart, -211, 111, 16)) {
488
return
1000 *
m_nPhotos
+ 21;
489
}
//tau+ decays to pi- pi0 nu_tau
490
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 10313)) {
491
return
1000 *
m_nPhotos
+ 22;
492
}
//tau- decays to nu_tau K- K_10
493
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 311)) {
494
return
1000 *
m_nPhotos
+ 22;
495
}
//tau- decays to nu_tau K- K0
496
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, -311)) {
497
return
1000 *
m_nPhotos
+ 22;
498
}
//tau- decays to nu_tau K- anti-K0
499
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 111, 111, 111, 111)) {
500
return
1000 *
m_nPhotos
+ 23;
501
}
//tau- decays to nu_tau pi- pi0 pi0 pi0 pi0
502
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 223, 111)) {
503
return
1000 *
m_nPhotos
+ 24;
504
}
//tau- decays to nu_tau pi- omega pi0
505
if
(
GenTauTag::PcheckDecay
(genpart, 221, -211, 211, -211, 16)) {
506
return
1000 *
m_nPhotos
+ 25;
507
}
//tau- decays to eta pi- pi+ pi- nu_tau
508
if
(
GenTauTag::PcheckDecay
(genpart, 221, -211, 111, 111, 16)) {
509
return
1000 *
m_nPhotos
+ 26;
510
}
//tau- decays to eta pi- pi0 pi0 nu_tau
511
if
(
GenTauTag::PcheckDecay
(genpart, 221, -321, 16)) {
512
return
1000 *
m_nPhotos
+ 27;
513
}
//tau- decays to eta K- nu_tau
514
if
(
GenTauTag::PcheckDecay
(genpart, 221, -323, 16)) {
515
return
1000 *
m_nPhotos
+ 28;
516
}
//tau- decays to eta K*- nu_tau
517
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 211, -211, 111)) {
518
return
1000 *
m_nPhotos
+ 29;
519
}
//tau- decays to nu_tau K- pi+ pi- pi0
520
if
(
GenTauTag::PcheckDecay
(genpart, -321, 111, 111, 111, 16)) {
521
return
1000 *
m_nPhotos
+ 30;
522
}
//tau- decays to K- pi0 pi0 pi0 nu_tau
523
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 211, -211, 10313)) {
524
return
1000 *
m_nPhotos
+ 31;
525
}
//tau- decays to nu_tau pi- pi+ pi- K_10
526
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 211, -211, 311)) {
527
return
1000 *
m_nPhotos
+ 31;
528
}
//tau- decays to nu_tau pi- pi+ pi- K0
529
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 211, -211, -311)) {
530
return
1000 *
m_nPhotos
+ 31;
531
}
//tau- decays to nu_tau pi- pi+ pi- anti-K0
532
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 111, 111, -10313)) {
533
return
1000 *
m_nPhotos
+ 32;
534
}
//tau- decays to nu_tau pi- pi0 pi0 anti-K_10
535
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 111, 111, 311)) {
536
return
1000 *
m_nPhotos
+ 32;
537
}
//tau- decays to nu_tau pi- pi0 pi0 K0
538
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 111, 111, -311)) {
539
return
1000 *
m_nPhotos
+ 32;
540
}
//tau- decays to nu_tau pi- pi0 pi0 anti-K0
541
if
(
GenTauTag::PcheckDecay
(genpart, -321, 321, -211, 111, 16)) {
542
return
1000 *
m_nPhotos
+ 33;
543
}
//tau- decays to K- K+ pi- pi0 nu_tau
544
if
(
GenTauTag::PcheckDecay
(genpart, 16, -10313, -211, 111, 10313)) {
545
return
1000 *
m_nPhotos
+ 34;
546
}
//tau- decays to nu_tau anti-K_10 pi- pi0 K_10
547
if
(
GenTauTag::PcheckDecay
(genpart, 16, 311, -211, 111, -311)) {
548
return
1000 *
m_nPhotos
+ 34;
549
}
//tau- decays to nu_tau K0 pi- pi0 anti-K0
550
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 223, 211, -211)) {
551
return
1000 *
m_nPhotos
+ 35;
552
}
//tau- decays to nu_tau pi- omega pi+ pi-
553
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 223, 111, 111)) {
554
return
1000 *
m_nPhotos
+ 36;
555
}
//tau- decays to nu_tau pi- omega pi0 pi0
556
if
(
GenTauTag::PcheckDecay
(genpart, 11, 11, -11, -12, 16)) {
557
return
1000 *
m_nPhotos
+ 37;
558
}
//tau- decays to e- e- e+ anti-nu_e nu_tau
559
if
(
GenTauTag::PcheckDecay
(genpart, 20223, -211, 16)) {
560
return
1000 *
m_nPhotos
+ 38;
561
}
//tau- decays to f_1 pi- nu_tau
562
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 223)) {
563
return
1000 *
m_nPhotos
+ 39;
564
}
//tau- decays to nu_tau K- omega
565
if
(
GenTauTag::PcheckDecay
(genpart, 16, 10313, -211, 211, -321)) {
566
return
1000 *
m_nPhotos
+ 40;
567
}
//tau- decays to nu_tau K_10 pi- pi+ K-
568
if
(
GenTauTag::PcheckDecay
(genpart, 16, 311, -211, 211, -321)) {
569
return
1000 *
m_nPhotos
+ 40;
570
}
//tau- decays to nu_tau K0 pi- pi+ K-
571
if
(
GenTauTag::PcheckDecay
(genpart, 16, -311, -211, 211, -321)) {
572
return
1000 *
m_nPhotos
+ 40;
573
}
//tau- decays to nu_tau anti-K0 pi- pi+ K-
574
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 111, 111, 10313)) {
575
return
1000 *
m_nPhotos
+ 41;
576
}
//tau- decays to nu_tau K- pi0 pi0 K_10
577
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 111, 111, 311)) {
578
return
1000 *
m_nPhotos
+ 41;
579
}
//tau- decays to nu_tau K- pi0 pi0 K0
580
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 111, 111, -311)) {
581
return
1000 *
m_nPhotos
+ 41;
582
}
//tau- decays to nu_tau K- pi0 pi0 anti-K0
583
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, -211, 321, -10313)) {
584
return
1000 *
m_nPhotos
+ 42;
585
}
//tau- decays to nu_tau pi- pi- K+ anti-K_10
586
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, -211, 321, 311)) {
587
return
1000 *
m_nPhotos
+ 42;
588
}
//tau- decays to nu_tau pi- pi- K+ K0
589
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, -211, 321, -311)) {
590
return
1000 *
m_nPhotos
+ 42;
591
}
//tau- decays to nu_tau pi- pi- K+ anti-K0
592
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, 311)) {
593
return
1000 *
m_nPhotos
+ 43;
594
}
//tau- decays to nu_tau pi- K0
595
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, -311)) {
596
return
1000 *
m_nPhotos
+ 43;
597
}
//tau- decays to nu_tau pi- anti-K0
598
if
(
GenTauTag::PcheckDecay
(genpart, 16, 22, -211, 111)) {
599
return
1000 *
m_nPhotos
+ 44;
600
}
//tau- decays to nu_tau gamma pi- pi0
601
if
(
GenTauTag::PcheckDecay
(genpart, 111, 111, -211, 16)) {
602
return
1000 *
m_nPhotos
+ 45;
603
}
//tau- decays to pi0 pi0 pi- nu_tau
604
if
(
GenTauTag::PcheckDecay
(genpart, -211, -211, 211, 16)) {
605
return
1000 *
m_nPhotos
+ 46;
606
}
//tau- decays to pi- pi- pi+ nu_tau
607
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 321, -321)) {
608
return
1000 *
m_nPhotos
+ 47;
609
}
//tau- decays to nu_tau K- K+ K-
610
if
(
GenTauTag::PcheckDecay
(genpart, -313, -211, 16)) {
611
return
1000 *
m_nPhotos
+ 48;
612
}
//tau- decays to anti-K*0 pi- nu_tau
613
if
(
GenTauTag::PcheckDecay
(genpart, 313, -321, 16)) {
614
return
1000 *
m_nPhotos
+ 49;
615
}
//tau- decays to K*0 K- nu_tau
616
if
(
GenTauTag::PcheckDecay
(genpart, 333, -211, 16)) {
617
return
1000 *
m_nPhotos
+ 50;
618
}
//tau- decays to phi pi- nu_tau
619
if
(
GenTauTag::PcheckDecay
(genpart, 333, -321, 16)) {
620
return
1000 *
m_nPhotos
+ 51;
621
}
//tau- decays to phi K- nu_tau
622
if
(
GenTauTag::PcheckDecay
(genpart, 13, -14, 16, 22)) {
623
return
1000 *
m_nPhotos
+ 52;
624
}
//tau- decays to mu- anti-nu_mu nu_tau gamma
625
if
(
GenTauTag::PcheckDecay
(genpart, 11, -12, 16, 22)) {
626
return
1000 *
m_nPhotos
+ 53;
627
}
//tau- decays to e- anti-nu_e nu_tau gamma
628
if
(
GenTauTag::PcheckDecay
(genpart, -321, 111, 16)) {
629
return
1000 *
m_nPhotos
+ 54;
630
}
//tau- decays to K- pi0 nu_tau
631
if
(
GenTauTag::PcheckDecay
(genpart, 311, -213, 16)) {
632
return
1000 *
m_nPhotos
+ 55;
633
}
//tau- decays to K0 rho- nu_tau
634
if
(
GenTauTag::PcheckDecay
(genpart, -311, -213, 16)) {
635
return
1000 *
m_nPhotos
+ 55;
636
}
//tau- decays to anti-K0 rho- nu_tau
637
if
(
GenTauTag::PcheckDecay
(genpart, -10323, 16)) {
638
return
1000 *
m_nPhotos
+ 56;
639
}
//tau- decays to K_1- nu_tau
640
if
(
GenTauTag::PcheckDecay
(genpart, -20323, 16)) {
641
return
1000 *
m_nPhotos
+ 57;
642
}
//tau- decays to K'_1- nu_tau
643
if
(
GenTauTag::PcheckDecay
(genpart, -100323, 16)) {
644
return
1000 *
m_nPhotos
+ 58;
645
}
//tau- decays to K'*- nu_tau
646
if
(
GenTauTag::PcheckDecay
(genpart, 221, -321, 111, 16)) {
647
return
1000 *
m_nPhotos
+ 59;
648
}
//tau- decays to eta K- pi0 nu_tau
649
if
(
GenTauTag::PcheckDecay
(genpart, 221, 311, -211, 16)) {
650
return
1000 *
m_nPhotos
+ 60;
651
}
//tau- decays to eta K0 pi- nu_tau
652
if
(
GenTauTag::PcheckDecay
(genpart, 221, -311, -211, 16)) {
653
return
1000 *
m_nPhotos
+ 60;
654
}
//tau- decays to eta anti-K0 pi- nu_tau
655
if
(
GenTauTag::PcheckDecay
(genpart, 16, 310, 310, -211)) {
656
return
1000 *
m_nPhotos
+ 20015;
657
}
//tau- decays to nu_tau K_S0 K_S0 pi-
658
if
(
GenTauTag::PcheckDecay
(genpart, 16, 310, 130, -211)) {
659
return
1000 *
m_nPhotos
+ 110015;
660
}
//tau- decays to nu_tau K_S0 K_L0 pi-
661
if
(
GenTauTag::PcheckDecay
(genpart, 16, 130, 130, -211)) {
662
return
1000 *
m_nPhotos
+ 200015;
663
}
//tau- decays to nu_tau K_L0 K_L0 pi-
664
if
(
GenTauTag::PcheckDecay
(genpart, 16, 310, -321, 111)) {
665
return
1000 *
m_nPhotos
+ 10016;
666
}
//tau- decays to nu_tau K_S0 K- pi0
667
if
(
GenTauTag::PcheckDecay
(genpart, 16, 130, -321, 111)) {
668
return
1000 *
m_nPhotos
+ 100016;
669
}
//tau- decays to nu_tau K_L0 K- pi0
670
if
(
GenTauTag::PcheckDecay
(genpart, 16, 310, -211, 111)) {
671
return
1000 *
m_nPhotos
+ 10019;
672
}
//tau- decays to nu_tau K_S0 pi- pi0
673
if
(
GenTauTag::PcheckDecay
(genpart, 16, 130, -211, 111)) {
674
return
1000 *
m_nPhotos
+ 100019;
675
}
//tau- decays to nu_tau K_L0 pi- pi0
676
if
(
GenTauTag::PcheckDecay
(genpart, 16, 310, -321)) {
677
return
1000 *
m_nPhotos
+ 10022;
678
}
//tau- decays to nu_tau K_S0 K-
679
if
(
GenTauTag::PcheckDecay
(genpart, 16, 130, -321)) {
680
return
1000 *
m_nPhotos
+ 100022;
681
}
//tau- decays to nu_tau K_L0 K-
682
if
(
GenTauTag::PcheckDecay
(genpart, 16, 310, -211, 211, -211)) {
683
return
1000 *
m_nPhotos
+ 10031;
684
}
//tau- decays to nu_tau K_S0 pi- pi+ pi-
685
if
(
GenTauTag::PcheckDecay
(genpart, 16, 130, -211, 211, -211)) {
686
return
1000 *
m_nPhotos
+ 100031;
687
}
//tau- decays to nu_tau K_L0 pi- pi+ pi-
688
if
(
GenTauTag::PcheckDecay
(genpart, 16, 310, -211, 111, 111)) {
689
return
1000 *
m_nPhotos
+ 10032;
690
}
//tau- decays to nu_tau K_S0 pi- pi0 pi0
691
if
(
GenTauTag::PcheckDecay
(genpart, 16, 130, -211, 111, 111)) {
692
return
1000 *
m_nPhotos
+ 100032;
693
}
//tau- decays to nu_tau K_L0 pi- pi0 pi0
694
if
(
GenTauTag::PcheckDecay
(genpart, 16, 310, 310, -211, 111)) {
695
return
1000 *
m_nPhotos
+ 20034;
696
}
//tau- decays to nu_tau K_S0 K_S0 pi- pi0
697
if
(
GenTauTag::PcheckDecay
(genpart, 16, 310, 130, -211, 111)) {
698
return
1000 *
m_nPhotos
+ 110034;
699
}
//tau- decays to nu_tau K_S0 K_L0 pi- pi0
700
if
(
GenTauTag::PcheckDecay
(genpart, 16, 130, 130, -211, 111)) {
701
return
1000 *
m_nPhotos
+ 200034;
702
}
//tau- decays to nu_tau K_L0 K_L0 pi- pi0
703
if
(
GenTauTag::PcheckDecay
(genpart, 16, 310, -211, 211, -321)) {
704
return
1000 *
m_nPhotos
+ 10040;
705
}
//tau- decays to nu_tau K_S0 pi- pi+ K-
706
if
(
GenTauTag::PcheckDecay
(genpart, 16, 130, -211, 211, -321)) {
707
return
1000 *
m_nPhotos
+ 100040;
708
}
//tau- decays to nu_tau K_L0 pi- pi+ K-
709
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 111, 111, 310)) {
710
return
1000 *
m_nPhotos
+ 10041;
711
}
//tau- decays to nu_tau K- pi0 pi0 K_S0
712
if
(
GenTauTag::PcheckDecay
(genpart, 16, -321, 111, 111, 130)) {
713
return
1000 *
m_nPhotos
+ 100041;
714
}
//tau- decays to nu_tau K- pi0 pi0 K_L0
715
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, -211, 321, 310)) {
716
return
1000 *
m_nPhotos
+ 10042;
717
}
//tau- decays to nu_tau pi- pi- K+ K_S0
718
if
(
GenTauTag::PcheckDecay
(genpart, 16, -211, -211, 321, 130)) {
719
return
1000 *
m_nPhotos
+ 100042;
720
}
//tau- decays to nu_tau pi- pi- K+ K_L0
721
return
0;
722
}
// not yet defined tau- decay
Belle2::GenTauTag::Mode_tau_minus
int Mode_tau_minus(std::vector< int > genDAU)
returns tau- mode identifier
Definition
GenTauTag.cc:410
Belle2::GenTauTag::GenTauTag
GenTauTag()
Constructor.
Belle2::GenTauTag::m_nPhotos
int m_nPhotos
number of photons
Definition
GenTauTag.h:35
Belle2::GenTauTag::PcheckDecay
bool PcheckDecay(std::vector< int > gp, int da1, int da2)
check decay with two daughters
Definition
GenTauTag.cc:16
Belle2::GenTauTag::Mode_tau_plus
int Mode_tau_plus(std::vector< int > genDAU)
returns tau+ mode identifier
Definition
GenTauTag.cc:95
Belle2
Abstract base class for different kinds of events.
Definition
MillepedeAlgorithm.h:17
analysis
utility
src
GenTauTag.cc
Generated on Mon Sep 1 2025 02:46:21 for Belle II Software by
1.13.2