BOSS 7.0.3
BESIII Offline Software System
Loading...
Searching...
No Matches
cut.cxx
Go to the documentation of this file.
1#include "tofcalgsec/TofDataSet.h"
2#include <cmath>
3bool Record::cutBarrel() {
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<80000 ) {
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 }
129 else {
130 if( m_texp<2.7 || m_texp>4.9 ) return false;
131 if( m_tleft<15. || m_tleft>33.5 ) return false;
132 if( m_tright<15. || m_tright>33.5 ) return false;
133 if( m_qleft<200. || m_qleft>20000. ) return false;
134 if( m_qright<200.|| m_qright>20000.) return false;
135 if( fabs(m_zrhit)>115.0 ) return false;
136 }
137
138 return true;
139}
140
141
142bool Record::cutEndcap() {
143 if( m_run>0 ) {
144 if( m_texp<4.6 || m_texp>5.4 ) return false;
145 if( m_qleft<50.|| m_qleft>20000.) return false;
146 // if( m_path<130. || m_path>170.) return false;
147 if( m_zrhit<50. || m_zrhit>83.) return false;
148
149 // 090308-090413 psip
150 if( m_run>=8107 && m_run<9671 ) {
151 if( m_tleft<24.0 || m_tleft>32.0 ) return false;
152 if( fabs(m_tleft-21.4-0.1*m_zrhit)>4.0 ) return false;
153 }
154
155 // 090528-090603 continue
156 else if( m_run>=9671 && m_run<9947 ) {
157 if( m_tleft<20.0 || m_tleft>28.0 ) return false;
158 if( fabs(m_tleft-17.3-0.1*m_zrhit)>4.0 ) return false;
159 }
160
161 // 090612-090727 jpsi
162 else if( m_run>=9947 && m_run<=10878 ) {
163 if( m_tleft<20.5 || m_tleft>28.5 ) return false;
164 if( fabs(m_tleft-17.6-0.1*m_zrhit)>4.0 ) return false;
165 }
166
167 // 100116-10713 psipp
168 else if( m_run>=11397 && m_run<12002 ) {
169 if( m_tleft<20.5 || m_tleft>28.5 ) return false;
170 if( fabs(m_tleft-17.8-0.1*m_zrhit)>4.0 ) return false;
171 }
172
173 // 110504-11060 psi(4040)
174 else if( m_run>=23463 && m_run<=24896 ) {
175 if( m_tleft<20.5 || m_tleft>28.5 ) return false;
176 if( fabs(m_tleft-18.0-0.1*m_zrhit)>4.0 ) return false;
177 }
178
179 // 111221-120331 psip
180 else if( m_run>=24897 && m_run<=27101 ) {
181 if( m_tleft<21.0 || m_tleft>30.0 ) return false;
182 if( fabs(m_tleft-18.60-0.1*m_zrhit)>4.0 ) return false;
183 }
184
185 // 120407-120617 jpsi
186 else if( m_run>=27102 && m_run<=29675 ) {
187 if( m_tleft<22.0 || m_tleft>30.0 ) return false;
188 if( fabs(m_tleft-20.17-0.1*m_zrhit)>4.0 ) return false;
189 }
190
191 // 121215 - 130606 pisp4260
192 else if( m_run>=29676 && m_run<33996 ) {
193 if( m_tleft<18.0 || m_tleft>27.0 ) return false;
194 if( fabs(m_tleft-15.63-0.1*m_zrhit)>4.0 ) return false;
195 }
196
197 // 131209-130529 R value scan
198 else if( m_run>=33997 && m_run<39355 ) {
199 if( m_tleft<17.5 || m_tleft>26.5 ) return false;
200 if( fabs(m_tleft-15.44-0.1*m_zrhit)>4.0 ) return false;
201 }
202
203 // 141231-150203 R value scan
204 else if( m_run>=39355 && m_run<40208 ) {
205 if( m_tleft<18.0 || m_tleft>27.0 ) return false;
206 if( fabs(m_tleft-15.8-0.1*m_zrhit)>3.0 ) return false;
207 }
208
209 // 150202- R value scan
210 else if( m_run>=40208 && m_run<80000 ) {
211 if( m_tleft<18.5 || m_tleft>27.5 ) return false;
212 if( fabs(m_tleft-16.5-0.1*m_zrhit)>3.0 ) return false;
213 }
214
215 }
216 else {
217 if( m_tleft<16.|| m_tleft>22. ) return false;
218 if( m_texp<4.5 || m_texp>5.4 ) return false;
219 if( m_qleft<200.|| m_qleft>20000.) return false;
220 if( m_path<120. || m_path>180.) return false;
221 if( m_zrhit<50. || m_zrhit>83.) return false;
222 }
223 return true;
224}
225
226
227bool Record::cutEtf() {
228 if( m_run>0 ) {
229 if( m_texp<4.5 || m_texp>5.6 ) return false;
230 if( fabs(m_zrhit)>7.8 ) return false;
231 // 150202 - R value scan
232 if( m_run>=40208 && m_run<43646 ) {
233 if( m_qleft<5.|| m_qleft>60. ) return false;
234 if( m_qright<5.|| m_qright>60.) return false;
235 if( m_tleft<21.5 || m_tleft>27.5 ) return false;
236 if( m_tright<21.5 || m_tright>27.5 ) return false;
237 }
238 // 160106 - 160108 4170
239 else if( m_run>=43646 && m_run<43682 ) {
240 if( m_qleft<5.|| m_qleft>60. ) return false;
241 if( m_qright<5.|| m_qright>60.) return false;
242 if( m_tleft<9.5 || m_tleft>15.5 ) return false;
243 if( m_tright<9.5 || m_tright>15.5 ) return false;
244 }
245 // 160108 4170
246 else if( m_run>=43682 && m_run<43809 ) {
247 if( m_qleft<5.|| m_qleft>60. ) return false;
248 if( m_qright<5.|| m_qright>60.) return false;
249 if( m_tleft<-254.5 || m_tleft>-248.5 ) return false;
250 if( m_tright<-254.5 || m_tright>-248.5 ) return false;
251 }
252 // 160114 4170
253 else if( m_run>=43809 && m_run<43810 ) {
254 if( m_qleft<5.|| m_qleft>60. ) return false;
255 if( m_qright<5.|| m_qright>60.) return false;
256 if( m_tleft<-14.5 || m_tleft>-8.5 ) return false;
257 if( m_tright<-14.5 || m_tright>-8.5 ) return false;
258 }
259 // 160114 4170
260 else if( m_run>=43810 && m_run<43879 ) {
261 if( m_qleft<5.|| m_qleft>60. ) return false;
262 if( m_qright<5.|| m_qright>60.) return false;
263 if( m_tleft<33.5 || m_tleft>39.5 ) return false;
264 if( m_tright<33.5 || m_tright>39.5 ) return false;
265 }
266 // 160116 4170
267 else if( m_run>=43879 && m_run<80000 ) {
268 if( m_qleft<5.|| m_qleft>60. ) return false;
269 if( m_qright<5.|| m_qright>60.) return false;
270 if( m_tleft<9.5 || m_tleft>15.5 ) return false;
271 if( m_tright<9.5 || m_tright>15.5 ) return false;
272 }
273
274 }
275 else {
276 if( m_texp<4.5 || m_texp>5.6 ) return false;
277 if( fabs(m_zrhit)>7.8 ) return false;
278 if( m_qleft<5.|| m_qleft>20. ) return false;
279 if( m_qright<5.|| m_qright>20.) return false;
280 if( m_tleft<4.5 || m_tleft>7.5 ) return false;
281 if( m_tright<4.5 || m_tright>7.5 ) return false;
282 }
283
284 return true;
285}
Record::cutEtf
bool cutEtf()
Definition: cut.cxx:227
Record::cutBarrel
bool cutBarrel()
Definition: cut.cxx:3
Record::cutEndcap
bool cutEndcap()
Definition: cut.cxx:142