00001 #include "tofcalgsec/TofDataSet.h"
00002 #include <cmath>
00003 bool Record::cutBarrel() {
00004 if( m_run>0 ) {
00005 if( m_texp<2.7 || m_texp>4.9 ) return false;
00006 if( fabs(m_zrhit)>115.0 ) return false;
00007
00008
00009 if( m_run>=8107 && m_run<9671 ) {
00010 if( ( m_tofid != 141 ) && ( m_qleft<200.|| m_qleft>20000. ) ) return false;
00011 if( m_qright<200.|| m_qright>20000. ) return false;
00012 if( ( m_tofid != 141 ) && ( m_tleft<17.5 || m_tleft>39.0 ) ) return false;
00013 if( m_tright<17.5 || m_tright>39.0 ) return false;
00014 if( ( m_tofid != 141 ) && ( fabs(m_tleft-27.4+0.07*m_zrhit ) > 4.0 ) ) return false;
00015 if( fabs(m_tright-27.6-0.07*m_zrhit) > 4.0 ) return false;
00016 }
00017
00018
00019 else if( m_run>=9671 && m_run<9947 ) {
00020 if( ( m_tofid != 141 ) && ( m_qleft<200.|| m_qleft>20000. ) ) return false;
00021 if( m_qright<200.|| m_qright>20000. ) return false;
00022 if( ( m_tofid != 141 ) && ( m_tleft<13.5 || m_tleft>35.0 ) ) return false;
00023 if( m_tright<13.5 || m_tright>35.0 ) return false;
00024 if( ( m_tofid != 141 ) && ( fabs(m_tleft-23.7+0.07*m_zrhit ) > 4.0 ) ) return false;
00025 if( fabs(m_tright-23.7-0.07*m_zrhit) > 4.0 ) return false;
00026 }
00027
00028
00029 else if( m_run>=9947 && m_run<=10878 ) {
00030 if( ( m_tofid != 141 ) && ( m_qleft<200.|| m_qleft>20000. ) ) return false;
00031 if( m_qright<200.|| m_qright>20000. ) return false;
00032 if( ( m_tofid != 141 ) && ( m_tleft<14.0 || m_tleft>35.0 ) ) return false;
00033 if( m_tright<14.0 || m_tright>35.0 ) return false;
00034 if( ( m_tofid != 141 ) && ( fabs(m_tleft-23.6+0.07*m_zrhit ) > 4.0 ) ) return false;
00035 if( fabs(m_tright-23.6-0.07*m_zrhit) > 4.0 ) return false;
00036 }
00037
00038
00039 else if( m_run>=11397 && m_run<12002 ) {
00040 if( m_qleft<200.|| m_qleft>20000. ) return false;
00041 if( m_qright<200.|| m_qright>20000. ) return false;
00042 if( m_tleft<13.5 || m_tleft>34.5 ) return false;
00043 if( m_tright<13.5 || m_tright>34.5 ) return false;
00044 if( fabs(m_tleft-23.5+0.07*m_zrhit ) > 4.0 ) return false;
00045 if( fabs(m_tright-23.5-0.07*m_zrhit) > 4.0 ) return false;
00046 }
00047
00048
00049 else if( m_run>=23463 && m_run<=24896 ) {
00050 if( m_qleft<200.|| m_qleft>20000. ) return false;
00051 if( m_qright<200.|| m_qright>20000. ) return false;
00052 if( m_tleft<14.0 || m_tleft>35.0 ) return false;
00053 if( m_tright<14.0 || m_tright>35.0 ) return false;
00054 if( fabs(m_tleft-23.75+0.07*m_zrhit ) > 4.0 ) return false;
00055 if( fabs(m_tright-23.75-0.07*m_zrhit) > 4.0 ) return false;
00056 }
00057
00058
00059 else if( m_run>=24897 && m_run<=27101 ) {
00060 if( m_qleft<200.|| m_qleft>20000. ) return false;
00061 if((m_tofid != 137 ) && ( m_qright<200.|| m_qright>20000.) ) return false;
00062 if( m_tleft<14.5 || m_tleft>35.5 ) return false;
00063 if((m_tofid != 137 ) && ( m_tright<14.5 || m_tright>35.5 )) return false;
00064 if( fabs(m_tleft-24.35+0.07*m_zrhit ) > 4.0 ) return false;
00065 if(( m_tofid != 137 ) && (fabs(m_tright-24.35-0.07*m_zrhit) > 4.0 )) return false;
00066 }
00067
00068
00069 else if( m_run>=27102 && m_run<=29675 ) {
00070 if( m_qleft<200.|| m_qleft>20000. ) return false;
00071 if((m_tofid != 137 ) && ( m_qright<200.|| m_qright>20000.) ) return false;
00072 if( m_tleft<14.5 || m_tleft>35.5 ) return false;
00073 if((m_tofid != 137 ) && ( m_tright<14.5 || m_tright>35.5 )) return false;
00074 if( fabs(m_tleft-25.40+0.07*m_zrhit ) > 4.0 ) return false;
00075 if(( m_tofid != 137 ) && (fabs(m_tright-25.40-0.07*m_zrhit) > 4.0 )) return false;
00076 }
00077
00078
00079 else if( m_run>=29676 && m_run<33996 ) {
00080 if( m_qleft<200.|| m_qleft>20000. ) return false;
00081 if((m_tofid != 137 ) && ( m_qright<200.|| m_qright>20000.) ) return false;
00082 if( m_tleft<11 || m_tleft>32.5 ) return false;
00083 if((m_tofid != 137 ) && ( m_tright<11 || m_tright>32.5 )) return false;
00084 if( fabs(m_tleft-21.00+0.07*m_zrhit ) > 4.0 ) return false;
00085 if(( m_tofid != 137 ) && (fabs(m_tright-21.01-0.07*m_zrhit) > 4.0 )) return false;
00086 }
00087
00088
00089 else if( m_run>=33997 && m_run<39355 ) {
00090 if( m_qleft<200.|| m_qleft>20000. ) return false;
00091 if((m_tofid != 8 && m_tofid != 137 ) && ( m_qright<200.|| m_qright>20000.) ) return false;
00092 if( m_tleft<11 || m_tleft>32.5 ) return false;
00093 if((m_tofid !=8 && m_tofid != 137 ) && ( m_tright<11 || m_tright>32.5 )) return false;
00094 if( fabs(m_tleft-21.00+0.07*m_zrhit ) > 4.0 ) return false;
00095 if((m_tofid !=8 && m_tofid != 137 ) && (fabs(m_tright-21.00-0.07*m_zrhit) > 4.0 )) return false;
00096 }
00097
00098
00099 else if( m_run>=39355 && m_run<40208 ) {
00100 if( m_qleft<200.|| m_qleft>20000. ) return false;
00101 if( m_qright<200.|| m_qright>20000.) return false;
00102 if( m_tleft<11 || m_tleft>32.5 ) return false;
00103 if( m_tright<11 || m_tright>32.5 ) return false;
00104 if( fabs(m_tleft-21.00+0.07*m_zrhit ) > 3.0 ) return false;
00105 if( fabs(m_tright-21.00-0.07*m_zrhit) > 3.0 ) return false;
00106 }
00107
00108
00109 else if( m_run>=40208 && m_run<43646 ) {
00110 if( m_qleft<200.|| m_qleft>20000. ) return false;
00111 if( m_qright<200.|| m_qright>20000.) return false;
00112 if( m_tleft<11.5 || m_tleft>33.0 ) return false;
00113 if( m_tright<11.5 || m_tright>33.0 ) return false;
00114 if( fabs(m_tleft-21.50+0.07*m_zrhit ) > 3.0 ) return false;
00115 if( fabs(m_tright-21.50-0.07*m_zrhit) > 3.0 ) return false;
00116 }
00117
00118
00119 else if( m_run>=43646 && m_run<80000 ) {
00120 if( m_qleft<200.|| m_qleft>20000. ) return false;
00121 if( m_qright<200.|| m_qright>20000.) return false;
00122 if( m_tleft<9.0 || m_tleft>30.0 ) return false;
00123 if( m_tright<9.0 || m_tright>30.0 ) return false;
00124 if( fabs(m_tleft-18.90+0.07*m_zrhit ) > 3.5 ) return false;
00125 if( fabs(m_tright-18.90-0.07*m_zrhit) > 3.5 ) return false;
00126 }
00127
00128 }
00129 else {
00130 if( m_texp<2.7 || m_texp>4.9 ) return false;
00131 if( m_tleft<15. || m_tleft>33.5 ) return false;
00132 if( m_tright<15. || m_tright>33.5 ) return false;
00133 if( m_qleft<200. || m_qleft>20000. ) return false;
00134 if( m_qright<200.|| m_qright>20000.) return false;
00135 if( fabs(m_zrhit)>115.0 ) return false;
00136 }
00137
00138 return true;
00139 }
00140
00141
00142 bool Record::cutEndcap() {
00143 if( m_run>0 ) {
00144 if( m_texp<4.6 || m_texp>5.4 ) return false;
00145 if( m_qleft<50.|| m_qleft>20000.) return false;
00146
00147 if( m_zrhit<50. || m_zrhit>83.) return false;
00148
00149
00150 if( m_run>=8107 && m_run<9671 ) {
00151 if( m_tleft<24.0 || m_tleft>32.0 ) return false;
00152 if( fabs(m_tleft-21.4-0.1*m_zrhit)>4.0 ) return false;
00153 }
00154
00155
00156 else if( m_run>=9671 && m_run<9947 ) {
00157 if( m_tleft<20.0 || m_tleft>28.0 ) return false;
00158 if( fabs(m_tleft-17.3-0.1*m_zrhit)>4.0 ) return false;
00159 }
00160
00161
00162 else if( m_run>=9947 && m_run<=10878 ) {
00163 if( m_tleft<20.5 || m_tleft>28.5 ) return false;
00164 if( fabs(m_tleft-17.6-0.1*m_zrhit)>4.0 ) return false;
00165 }
00166
00167
00168 else if( m_run>=11397 && m_run<12002 ) {
00169 if( m_tleft<20.5 || m_tleft>28.5 ) return false;
00170 if( fabs(m_tleft-17.8-0.1*m_zrhit)>4.0 ) return false;
00171 }
00172
00173
00174 else if( m_run>=23463 && m_run<=24896 ) {
00175 if( m_tleft<20.5 || m_tleft>28.5 ) return false;
00176 if( fabs(m_tleft-18.0-0.1*m_zrhit)>4.0 ) return false;
00177 }
00178
00179
00180 else if( m_run>=24897 && m_run<=27101 ) {
00181 if( m_tleft<21.0 || m_tleft>30.0 ) return false;
00182 if( fabs(m_tleft-18.60-0.1*m_zrhit)>4.0 ) return false;
00183 }
00184
00185
00186 else if( m_run>=27102 && m_run<=29675 ) {
00187 if( m_tleft<22.0 || m_tleft>30.0 ) return false;
00188 if( fabs(m_tleft-20.17-0.1*m_zrhit)>4.0 ) return false;
00189 }
00190
00191
00192 else if( m_run>=29676 && m_run<33996 ) {
00193 if( m_tleft<18.0 || m_tleft>27.0 ) return false;
00194 if( fabs(m_tleft-15.63-0.1*m_zrhit)>4.0 ) return false;
00195 }
00196
00197
00198 else if( m_run>=33997 && m_run<39355 ) {
00199 if( m_tleft<17.5 || m_tleft>26.5 ) return false;
00200 if( fabs(m_tleft-15.44-0.1*m_zrhit)>4.0 ) return false;
00201 }
00202
00203
00204 else if( m_run>=39355 && m_run<40208 ) {
00205 if( m_tleft<18.0 || m_tleft>27.0 ) return false;
00206 if( fabs(m_tleft-15.8-0.1*m_zrhit)>3.0 ) return false;
00207 }
00208
00209
00210 else if( m_run>=40208 && m_run<80000 ) {
00211 if( m_tleft<18.5 || m_tleft>27.5 ) return false;
00212 if( fabs(m_tleft-16.5-0.1*m_zrhit)>3.0 ) return false;
00213 }
00214
00215 }
00216 else {
00217 if( m_tleft<16.|| m_tleft>22. ) return false;
00218 if( m_texp<4.5 || m_texp>5.4 ) return false;
00219 if( m_qleft<200.|| m_qleft>20000.) return false;
00220 if( m_path<120. || m_path>180.) return false;
00221 if( m_zrhit<50. || m_zrhit>83.) return false;
00222 }
00223 return true;
00224 }
00225
00226
00227 bool Record::cutEtf() {
00228 if( m_run>0 ) {
00229 if( m_texp<4.5 || m_texp>5.6 ) return false;
00230 if( fabs(m_zrhit)>7.8 ) return false;
00231
00232 if( m_run>=40208 && m_run<43646 ) {
00233 if( m_qleft<5.|| m_qleft>60. ) return false;
00234 if( m_qright<5.|| m_qright>60.) return false;
00235 if( m_tleft<21.5 || m_tleft>27.5 ) return false;
00236 if( m_tright<21.5 || m_tright>27.5 ) return false;
00237 }
00238
00239 else if( m_run>=43646 && m_run<43682 ) {
00240 if( m_qleft<5.|| m_qleft>60. ) return false;
00241 if( m_qright<5.|| m_qright>60.) return false;
00242 if( m_tleft<9.5 || m_tleft>15.5 ) return false;
00243 if( m_tright<9.5 || m_tright>15.5 ) return false;
00244 }
00245
00246 else if( m_run>=43682 && m_run<43809 ) {
00247 if( m_qleft<5.|| m_qleft>60. ) return false;
00248 if( m_qright<5.|| m_qright>60.) return false;
00249 if( m_tleft<-254.5 || m_tleft>-248.5 ) return false;
00250 if( m_tright<-254.5 || m_tright>-248.5 ) return false;
00251 }
00252
00253 else if( m_run>=43809 && m_run<43810 ) {
00254 if( m_qleft<5.|| m_qleft>60. ) return false;
00255 if( m_qright<5.|| m_qright>60.) return false;
00256 if( m_tleft<-14.5 || m_tleft>-8.5 ) return false;
00257 if( m_tright<-14.5 || m_tright>-8.5 ) return false;
00258 }
00259
00260 else if( m_run>=43810 && m_run<43879 ) {
00261 if( m_qleft<5.|| m_qleft>60. ) return false;
00262 if( m_qright<5.|| m_qright>60.) return false;
00263 if( m_tleft<33.5 || m_tleft>39.5 ) return false;
00264 if( m_tright<33.5 || m_tright>39.5 ) return false;
00265 }
00266
00267 else if( m_run>=43879 && m_run<80000 ) {
00268 if( m_qleft<5.|| m_qleft>60. ) return false;
00269 if( m_qright<5.|| m_qright>60.) return false;
00270 if( m_tleft<9.5 || m_tleft>15.5 ) return false;
00271 if( m_tright<9.5 || m_tright>15.5 ) return false;
00272 }
00273
00274 }
00275 else {
00276 if( m_texp<4.5 || m_texp>5.6 ) return false;
00277 if( fabs(m_zrhit)>7.8 ) return false;
00278 if( m_qleft<5.|| m_qleft>20. ) return false;
00279 if( m_qright<5.|| m_qright>20.) return false;
00280 if( m_tleft<4.5 || m_tleft>7.5 ) return false;
00281 if( m_tright<4.5 || m_tright>7.5 ) return false;
00282 }
00283
00284 return true;
00285 }