BOSS 7.1.2
BESIII Offline Software System
Loading...
Searching...
No Matches
cut.cxx
Go to the documentation of this file.
2#include <cmath>
4 if( m_run>0 ) {
5 if( m_texp<2.7 || m_texp>4.9 ) return false;
6 if( fabs(m_zrhit)>115.0 ) return false;
7
8 // 090308-090413 psip
9 if( m_run>=8107 && m_run<9671 ) {
10 if( ( m_tofid != 141 ) && ( m_qleft<200.|| m_qleft>20000. ) ) return false;
11 if( m_qright<200.|| m_qright>20000. ) return false;
12 if( ( m_tofid != 141 ) && ( m_tleft<17.5 || m_tleft>39.0 ) ) return false;
13 if( m_tright<17.5 || m_tright>39.0 ) return false;
14 if( ( m_tofid != 141 ) && ( fabs(m_tleft-27.4+0.07*m_zrhit ) > 4.0 ) ) return false;
15 if( fabs(m_tright-27.6-0.07*m_zrhit) > 4.0 ) return false;
16 }
17
18 // 090528-090603 continue
19 else if( m_run>=9671 && m_run<9947 ) {
20 if( ( m_tofid != 141 ) && ( m_qleft<200.|| m_qleft>20000. ) ) return false;
21 if( m_qright<200.|| m_qright>20000. ) return false;
22 if( ( m_tofid != 141 ) && ( m_tleft<13.5 || m_tleft>35.0 ) ) return false;
23 if( m_tright<13.5 || m_tright>35.0 ) return false;
24 if( ( m_tofid != 141 ) && ( fabs(m_tleft-23.7+0.07*m_zrhit ) > 4.0 ) ) return false;
25 if( fabs(m_tright-23.7-0.07*m_zrhit) > 4.0 ) return false;
26 }
27
28 // 090612-090727 jpsi
29 else if( m_run>=9947 && m_run<=10878 ) {
30 if( ( m_tofid != 141 ) && ( m_qleft<200.|| m_qleft>20000. ) ) return false;
31 if( m_qright<200.|| m_qright>20000. ) return false;
32 if( ( m_tofid != 141 ) && ( m_tleft<14.0 || m_tleft>35.0 ) ) return false;
33 if( m_tright<14.0 || m_tright>35.0 ) return false;
34 if( ( m_tofid != 141 ) && ( fabs(m_tleft-23.6+0.07*m_zrhit ) > 4.0 ) ) return false;
35 if( fabs(m_tright-23.6-0.07*m_zrhit) > 4.0 ) return false;
36 }
37
38 // 100116-100713 psipp
39 else if( m_run>=11397 && m_run<12002 ) {
40 if( m_qleft<200.|| m_qleft>20000. ) return false;
41 if( m_qright<200.|| m_qright>20000. ) return false;
42 if( m_tleft<13.5 || m_tleft>34.5 ) return false;
43 if( m_tright<13.5 || m_tright>34.5 ) return false;
44 if( fabs(m_tleft-23.5+0.07*m_zrhit ) > 4.0 ) return false;
45 if( fabs(m_tright-23.5-0.07*m_zrhit) > 4.0 ) return false;
46 }
47
48 // 110504-1106 psi(4040)
49 else if( m_run>=23463 && m_run<=24896 ) {
50 if( m_qleft<200.|| m_qleft>20000. ) return false;
51 if( m_qright<200.|| m_qright>20000. ) return false;
52 if( m_tleft<14.0 || m_tleft>35.0 ) return false;
53 if( m_tright<14.0 || m_tright>35.0 ) return false;
54 if( fabs(m_tleft-23.75+0.07*m_zrhit ) > 4.0 ) return false;
55 if( fabs(m_tright-23.75-0.07*m_zrhit) > 4.0 ) return false;
56 }
57
58 // 111221-120331 psip
59 else if( m_run>=24897 && m_run<=27101 ) {
60 if( m_qleft<200.|| m_qleft>20000. ) return false;
61 if((m_tofid != 137 ) && ( m_qright<200.|| m_qright>20000.) ) return false;
62 if( m_tleft<14.5 || m_tleft>35.5 ) return false;
63 if((m_tofid != 137 ) && ( m_tright<14.5 || m_tright>35.5 )) return false;
64 if( fabs(m_tleft-24.35+0.07*m_zrhit ) > 4.0 ) return false;
65 if(( m_tofid != 137 ) && (fabs(m_tright-24.35-0.07*m_zrhit) > 4.0 )) return false;
66 }
67
68 // 120407-120617 jpsi
69 else if( m_run>=27102 && m_run<=29675 ) {
70 if( m_qleft<200.|| m_qleft>20000. ) return false;
71 if((m_tofid != 137 ) && ( m_qright<200.|| m_qright>20000.) ) return false;
72 if( m_tleft<14.5 || m_tleft>35.5 ) return false;
73 if((m_tofid != 137 ) && ( m_tright<14.5 || m_tright>35.5 )) return false;
74 if( fabs(m_tleft-25.40+0.07*m_zrhit ) > 4.0 ) return false;
75 if(( m_tofid != 137 ) && (fabs(m_tright-25.40-0.07*m_zrhit) > 4.0 )) return false;
76 }
77
78 // 121215 - 130606 Y4260/4360
79 else if( m_run>=29676 && m_run<33996 ) {
80 if( m_qleft<200.|| m_qleft>20000. ) return false;
81 if((m_tofid != 137 ) && ( m_qright<200.|| m_qright>20000.) ) return false;
82 if( m_tleft<11 || m_tleft>32.5 ) return false;
83 if((m_tofid != 137 ) && ( m_tright<11 || m_tright>32.5 )) return false;
84 if( fabs(m_tleft-21.00+0.07*m_zrhit ) > 4.0 ) return false;
85 if(( m_tofid != 137 ) && (fabs(m_tright-21.01-0.07*m_zrhit) > 4.0 )) return false;
86 }
87
88 // 131209 - 130529 R value scan
89 else if( m_run>=33997 && m_run<39355 ) {
90 if( m_qleft<200.|| m_qleft>20000. ) return false;
91 if((m_tofid != 8 && m_tofid != 137 ) && ( m_qright<200.|| m_qright>20000.) ) return false;
92 if( m_tleft<11 || m_tleft>32.5 ) return false;
93 if((m_tofid !=8 && m_tofid != 137 ) && ( m_tright<11 || m_tright>32.5 )) return false;
94 if( fabs(m_tleft-21.00+0.07*m_zrhit ) > 4.0 ) return false;
95 if((m_tofid !=8 && m_tofid != 137 ) && (fabs(m_tright-21.00-0.07*m_zrhit) > 4.0 )) return false;
96 }
97
98 // 141231 - 150202 R value scan
99 else if( m_run>=39355 && m_run<40208 ) {
100 if( m_qleft<200.|| m_qleft>20000. ) return false;
101 if( m_qright<200.|| m_qright>20000.) return false;
102 if( m_tleft<11 || m_tleft>32.5 ) return false;
103 if( m_tright<11 || m_tright>32.5 ) return false;
104 if( fabs(m_tleft-21.00+0.07*m_zrhit ) > 3.0 ) return false;
105 if( fabs(m_tright-21.00-0.07*m_zrhit) > 3.0 ) return false;
106 }
107
108 // 150202 - R value scan
109 else if( m_run>=40208 && m_run<43646 ) {
110 if( m_qleft<200.|| m_qleft>20000. ) return false;
111 if( m_qright<200.|| m_qright>20000.) return false;
112 if( m_tleft<11.5 || m_tleft>33.0 ) return false;
113 if( m_tright<11.5 || m_tright>33.0 ) return false;
114 if( fabs(m_tleft-21.50+0.07*m_zrhit ) > 3.0 ) return false;
115 if( fabs(m_tright-21.50-0.07*m_zrhit) > 3.0 ) return false;
116 }
117
118 // 160106 - 4170
119 else if( m_run>=43646 && m_run<47432 ) {
120 if( m_qleft<200.|| m_qleft>20000. ) return false;
121 if( m_qright<200.|| m_qright>20000.) return false;
122 if( m_tleft<9.0 || m_tleft>30.0 ) return false;
123 if( m_tright<9.0 || m_tright>30.0 ) return false;
124 if( fabs(m_tleft-18.90+0.07*m_zrhit ) > 3.5 ) return false;
125 if( fabs(m_tright-18.90-0.07*m_zrhit) > 3.5 ) return false;
126 }
127
128 // 161211 - XYZ
129 else if( m_run>=47432 && m_run<52940 ) {
130 if( m_qleft<200.|| m_qleft>20000. ) return false;
131 if( m_qright<200.|| m_qright>20000.) return false;
132 if( m_tleft<8.5 || m_tleft>29.5 ) return false;
133 if( m_tright<8.5 || m_tright>29.5 ) return false;
134 if( fabs(m_tleft-18.50+0.07*m_zrhit ) > 3.5 ) return false;
135 if( fabs(m_tright-18.50-0.07*m_zrhit) > 3.5 ) return false;
136 }
137 // 171208 - Jpsi
138 else if( m_run>=52940 && m_run<53500 ) {
139 if( m_qleft<200.|| m_qleft>20000. ) return false;
140 if( m_qright<200.|| m_qright>20000.) return false;
141 if( m_tleft<8.5 || m_tleft>29.5 ) return false;
142 if( m_tright<8.5 || m_tright>29.5 ) return false;
143 if( fabs(m_tleft-18.50+0.07*m_zrhit ) > 3.5 ) return false;
144 if( fabs(m_tright-18.50-0.07*m_zrhit) > 3.5 ) return false;
145 }
146 // 180202 - Jpsi
147 else if( m_run>=53500 && m_run<53586 ) {
148 if( m_qleft<200.|| m_qleft>20000. ) return false;
149 if( m_qright<200.|| m_qright>20000.) return false;
150 if( ( m_tofid>=9 && m_tofid<=12 ) || ( m_tofid>=97 && m_tofid<=100 ) ) {
151 if( m_tleft<-12.0 || m_tleft>9.0 ) return false;
152 if( m_tright<-12.0 || m_tright>9.0 ) return false;
153 if( fabs(m_tleft+2.7+0.07*m_zrhit ) > 3.5 ) return false;
154 if( fabs(m_tright+2.7-0.07*m_zrhit) > 3.5 ) return false;
155 }
156 else {
157 if( m_tleft<9.5 || m_tleft>31.0 ) return false;
158 if( m_tright<9.5 || m_tright>31.0 ) return false;
159 if( fabs(m_tleft-18.50+0.07*m_zrhit ) > 3.5 ) return false;
160 if( fabs(m_tright-18.50-0.07*m_zrhit) > 3.5 ) return false;
161 }
162 }
163 // 180207 - 180405 Jpsi
164 else if( m_run>=53586 && m_run<54780 ) {
165 if( m_qleft<200.|| m_qleft>20000. ) return false;
166 if( m_qright<200.|| m_qright>20000.) return false;
167 if( m_tleft<8.5 || m_tleft>31.0 ) return false;
168 if( m_tright<8.5 || m_tright>31.0 ) return false;
169 if( (m_tleft+0.07*m_zrhit)<15.0 || (m_tleft+0.07*m_zrhit)>24.0 ) return false;
170 if( (m_tright-0.07*m_zrhit)<15.0 || (m_tright-0.07*m_zrhit)>24.0 ) return false;
171 }
172 // 180405 Jpsi MDC t0 changed
173 else if( m_run>=54780 && m_run<55115 ) {
174 if( m_qleft<200.|| m_qleft>20000. ) return false;
175 if( m_qright<200.|| m_qright>20000.) return false;
176 if( m_tleft<9.4 || m_tleft>32.0 ) return false;
177 if( m_tright<9.4 || m_tright>32.0 ) return false;
178 if( (m_tleft+0.07*m_zrhit)<16.0 || (m_tleft+0.07*m_zrhit)>25.0 ) return false;
179 if( (m_tright-0.07*m_zrhit)<16.0 || (m_tright-0.07*m_zrhit)>25.0 ) return false;
180 }
181 // 180420 - Tau and psip scan
182 else if( m_run>=55115 && m_run<55861 ) {
183 if( m_qleft<200.|| m_qleft>20000. ) return false;
184 if( m_qright<200.|| m_qright>20000. ) return false;
185 if( m_tleft<9.0 || m_tleft>30.5 ) return false;
186 if( m_tright<9.0 || m_tright>30.5 ) return false;
187 if( fabs(m_tleft+0.07*m_zrhit-19.0 ) > 3.5 ) return false;
188 if( fabs(m_tright-0.07*m_zrhit-19.0) > 3.5 ) return false;
189 }
190 // 180526 - Jpsi
191 else if( m_run>=55861 && m_run<56769 ) {
192 if( m_qleft<200.|| m_qleft>20000. ) return false;
193 if( m_qright<200.|| m_qright>20000. ) return false;
194 if( m_tleft<10.0 || m_tleft>31.0 ) return false;
195 if( m_tright<10.0 || m_tright>31.0 ) return false;
196 if( fabs(m_tleft+0.07*m_zrhit-20.0 ) > 3.5 ) return false;
197 if( fabs(m_tright-0.07*m_zrhit-20.0) > 3.5 ) return false;
198 }
199 // 181118 - Jpsi
200 else if( m_run>=56769 && m_run<1000000 ) {
201 if( m_qleft<200.|| m_qleft>20000. ) return false;
202 if( m_qright<200.|| m_qright>20000. ) return false;
203 if( m_tleft<10.5 || m_tleft>31.5 ) return false;
204 if( m_tright<10.5 || m_tright>31.5 ) return false;
205 if( fabs(m_tleft+0.07*m_zrhit-20.0 ) > 3.5 ) return false;
206 if( fabs(m_tright-0.07*m_zrhit-20.0) > 3.5 ) return false;
207 }
208
209 }
210 else {
211 if( m_texp<2.7 || m_texp>4.9 ) return false;
212 if( m_tleft<15. || m_tleft>33.5 ) return false;
213 if( m_tright<15. || m_tright>33.5 ) return false;
214 if( m_qleft<200. || m_qleft>20000. ) return false;
215 if( m_qright<200.|| m_qright>20000.) return false;
216 if( fabs(m_zrhit)>115.0 ) return false;
217 }
218
219 return true;
220}
221
222
224 if( m_run>0 ) {
225 if( m_texp<4.6 || m_texp>5.4 ) return false;
226 if( m_qleft<50.|| m_qleft>20000.) return false;
227 // if( m_path<130. || m_path>170.) return false;
228 if( m_zrhit<50. || m_zrhit>83.) return false;
229
230 // 090308-090413 psip
231 if( m_run>=8107 && m_run<9671 ) {
232 if( m_tleft<24.0 || m_tleft>32.0 ) return false;
233 if( fabs(m_tleft-21.4-0.1*m_zrhit)>4.0 ) return false;
234 }
235
236 // 090528-090603 continue
237 else if( m_run>=9671 && m_run<9947 ) {
238 if( m_tleft<20.0 || m_tleft>28.0 ) return false;
239 if( fabs(m_tleft-17.3-0.1*m_zrhit)>4.0 ) return false;
240 }
241
242 // 090612-090727 jpsi
243 else if( m_run>=9947 && m_run<=10878 ) {
244 if( m_tleft<20.5 || m_tleft>28.5 ) return false;
245 if( fabs(m_tleft-17.6-0.1*m_zrhit)>4.0 ) return false;
246 }
247
248 // 100116-10713 psipp
249 else if( m_run>=11397 && m_run<12002 ) {
250 if( m_tleft<20.5 || m_tleft>28.5 ) return false;
251 if( fabs(m_tleft-17.8-0.1*m_zrhit)>4.0 ) return false;
252 }
253
254 // 110504-11060 psi(4040)
255 else if( m_run>=23463 && m_run<=24896 ) {
256 if( m_tleft<20.5 || m_tleft>28.5 ) return false;
257 if( fabs(m_tleft-18.0-0.1*m_zrhit)>4.0 ) return false;
258 }
259
260 // 111221-120331 psip
261 else if( m_run>=24897 && m_run<=27101 ) {
262 if( m_tleft<21.0 || m_tleft>30.0 ) return false;
263 if( fabs(m_tleft-18.60-0.1*m_zrhit)>4.0 ) return false;
264 }
265
266 // 120407-120617 jpsi
267 else if( m_run>=27102 && m_run<=29675 ) {
268 if( m_tleft<22.0 || m_tleft>30.0 ) return false;
269 if( fabs(m_tleft-20.17-0.1*m_zrhit)>4.0 ) return false;
270 }
271
272 // 121215 - 130606 pisp4260
273 else if( m_run>=29676 && m_run<33996 ) {
274 if( m_tleft<18.0 || m_tleft>27.0 ) return false;
275 if( fabs(m_tleft-15.63-0.1*m_zrhit)>4.0 ) return false;
276 }
277
278 // 131209-130529 R value scan
279 else if( m_run>=33997 && m_run<39355 ) {
280 if( m_tleft<17.5 || m_tleft>26.5 ) return false;
281 if( fabs(m_tleft-15.44-0.1*m_zrhit)>4.0 ) return false;
282 }
283
284 // 141231-150203 R value scan
285 else if( m_run>=39355 && m_run<40208 ) {
286 if( m_tleft<18.0 || m_tleft>27.0 ) return false;
287 if( fabs(m_tleft-15.8-0.1*m_zrhit)>3.0 ) return false;
288 }
289
290 // 150202- R value scan
291 else if( m_run>=40208 && m_run<1000000 ) {
292 if( m_tleft<18.5 || m_tleft>27.5 ) return false;
293 if( fabs(m_tleft-16.5-0.1*m_zrhit)>3.0 ) return false;
294 }
295
296 }
297 else {
298 if( m_tleft<16.|| m_tleft>22. ) return false;
299 if( m_texp<4.5 || m_texp>5.4 ) return false;
300 if( m_qleft<200.|| m_qleft>20000.) return false;
301 if( m_path<120. || m_path>180.) return false;
302 if( m_zrhit<50. || m_zrhit>83.) return false;
303 }
304 return true;
305}
306
307
309 if( m_run>0 ) {
310 if( m_texp<4.5 || m_texp>5.6 ) return false;
311 if( fabs(m_zrhit)>7.8 ) return false;
312 // dead channels
313 bool rightDead = ( ( m_tofid==3 && ( m_strip==1 || m_strip==3 || ( m_strip>=5 && m_strip<=11 ) ) ) || ( m_tofid==13 && m_strip==7 ) );
314 bool leftDead = ( m_tofid==34 && ( m_strip>=0 && m_strip<=5 ) );
315 bool leftSpecial = ( m_tofid==13 && m_strip==5 );
316
317 // 150202 - R value scan
318 if( m_run>=40208 && m_run<43646 ) {
319 if( m_qleft<5.|| m_qleft>60. ) return false;
320 if( m_qright<5.|| m_qright>60.) return false;
321 if( m_tleft<21.5 || m_tleft>27.5 ) return false;
322 if( m_tright<21.5 || m_tright>27.5 ) return false;
323 }
324 // 160106 - 160108 4170
325 else if( m_run>=43646 && m_run<43682 ) {
326 if( m_qleft<5.|| m_qleft>60. ) return false;
327 if( m_qright<5.|| m_qright>60.) return false;
328 if( m_tleft<9.5 || m_tleft>15.5 ) return false;
329 if( m_tright<9.5 || m_tright>15.5 ) return false;
330 }
331 // 160108 4170
332 else if( m_run>=43682 && m_run<43809 ) {
333 if( m_qleft<5.|| m_qleft>60. ) return false;
334 if( m_qright<5.|| m_qright>60.) return false;
335 if( m_tleft<-254.5 || m_tleft>-248.5 ) return false;
336 if( m_tright<-254.5 || m_tright>-248.5 ) return false;
337 }
338 // 160114 4170
339 else if( m_run>=43809 && m_run<43810 ) {
340 if( m_qleft<5.|| m_qleft>60. ) return false;
341 if( m_qright<5.|| m_qright>60.) return false;
342 if( m_tleft<-14.5 || m_tleft>-8.5 ) return false;
343 if( m_tright<-14.5 || m_tright>-8.5 ) return false;
344 }
345 // 160114 4170
346 else if( m_run>=43810 && m_run<43879 ) {
347 if( m_qleft<5.|| m_qleft>60. ) return false;
348 if( m_qright<5.|| m_qright>60.) return false;
349 if( m_tleft<33.5 || m_tleft>39.5 ) return false;
350 if( m_tright<33.5 || m_tright>39.5 ) return false;
351 }
352 // 160116 4170
353 else if( m_run>=43879 && m_run<47432 ) {
354 if( m_qleft<5.|| m_qleft>60. ) return false;
355 if( m_qright<5.|| m_qright>60.) return false;
356 if( m_tleft<9.5 || m_tleft>15.5 ) return false;
357 if( m_tright<9.5 || m_tright>15.5 ) return false;
358 }
359 // 161211 XYZ
360 else if( m_run>=47432 && m_run<52940 ) {
361 if( m_qleft<5.|| m_qleft>60. ) return false;
362 if( m_qright<5.|| m_qright>60.) return false;
363 if( m_tleft<9.5 || m_tleft>15.5 ) return false;
364 if( m_tright<9.5 || m_tright>15.5 ) return false;
365 }
366 // 171208 - 180405 Jpsi also add dead channel
367 else if( m_run>=52940 && m_run<54780 ) {
368 if( !leftDead && ( m_qleft<5.|| m_qleft>60. ) ) return false;
369 if( !rightDead && ( m_qright<5.|| m_qright>60.) ) return false;
370 if( !leftDead && ( m_tleft<9.5 || m_tleft>15.5 ) ) return false;
371 if( !rightDead && ( m_tright<9.5 || m_tright>15.5 ) ) return false;
372 }
373 // 180405 Jpsi MDC t0 changed
374 else if( m_run>=54780 && m_run<55115 ) {
375 if( !leftDead && ( m_qleft<5.|| m_qleft>60. ) ) return false;
376 if( !rightDead && ( m_qright<5.|| m_qright>60.) ) return false;
377 if( !leftDead && ( m_tleft<9.5 || m_tleft>15.5 ) ) return false;
378 if( !rightDead && ( m_tright<9.5 || m_tright>15.5 ) ) return false;
379 }
380 // 180420 - Tau and psip scan
381 else if( m_run>=55115 && m_run<55861 ) {
382 if( !leftDead && ( m_qleft<5.|| m_qleft>60. ) ) return false;
383 if( !rightDead && ( m_qright<5.|| m_qright>60.) ) return false;
384 if( !leftDead && ( m_tleft<10.0 || m_tleft>16.0 ) ) return false;
385 if( !rightDead && ( m_tright<10.0 || m_tright>16.0 ) ) return false;
386 }
387 // 180526 - Tau and psip scan
388 else if( m_run>=55861 && m_run<56769 ) {
389 if( !leftDead && ( m_qleft<5.|| m_qleft>60. ) ) return false;
390 if( !rightDead && ( m_qright<5.|| m_qright>60.) ) return false;
391 if( !leftDead && ( m_tleft<11.0 || m_tleft>17.0 ) ) return false;
392 if( !rightDead && ( m_tright<11.0 || m_tright>17.0 ) ) return false;
393 }
394 // 181118 - 190108 Jpsi
395 else if( m_run>=56769 && m_run<58239 ) {
396 if( !leftDead && ( m_qleft<5.|| m_qleft>60. ) ) return false;
397 if( !rightDead && ( m_qright<5.|| m_qright>60.) ) return false;
398 if( !leftDead && ( m_tleft<11.0 || m_tleft>18.0 ) ) return false;
399 if( !rightDead && ( m_tright<11.0 || m_tright>18.0 ) ) return false;
400 }
401 // 190109 Jpsi
402 else if( m_run>=58239 && m_run<1000000 ) {
403 if( !leftDead && ( m_qleft<5.|| m_qleft>60. ) ) return false;
404 if( !rightDead && ( m_qright<5.|| m_qright>60.) ) return false;
405 if( !leftDead && ( m_tleft<11.0 || m_tleft>17.5 ) ) return false;
406 if( !rightDead && ( m_tright<11.0 || m_tright>17.5 ) ) return false;
407 }
408
409 }
410 else {
411 if( m_texp<4.5 || m_texp>5.6 ) return false;
412 if( fabs(m_zrhit)>7.8 ) return false;
413 if( m_qleft<5.|| m_qleft>20. ) return false;
414 if( m_qright<5.|| m_qright>20.) return false;
415 if( m_tleft<4.5 || m_tleft>7.5 ) return false;
416 if( m_tright<4.5 || m_tright>7.5 ) return false;
417 }
418
419 return true;
420}
bool cutEtf()
Definition cut.cxx:308
bool cutBarrel()
Definition cut.cxx:3
bool cutEndcap()
Definition cut.cxx:223