/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Tof/tofcalgsec/tofcalgsec-00-02-21/src/cut.cxx

Go to the documentation of this file.
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     // 090308-090413 psip
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     // 090528-090603 continue
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     // 090612-090727 jpsi
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     // 100116-100713 psipp
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     // 110504-1106 psi(4040)
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     // 111221-120331 psip
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     // 120407-120617 jpsi
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     // 121215 - 130606 Y4260/4360
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     // 131209 - 130529  R value scan
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     // 141231 - 150202 R value scan
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     // 150202 - R value scan
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     // 160106 - 4170
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     //    if( m_path<130. || m_path>170.) return false;
00147     if( m_zrhit<50. || m_zrhit>83.) return false;
00148 
00149     // 090308-090413 psip
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     // 090528-090603 continue
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     // 090612-090727 jpsi
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     // 100116-10713 psipp
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     // 110504-11060 psi(4040)
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     // 111221-120331 psip
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     // 120407-120617 jpsi
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     // 121215 - 130606      pisp4260
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     // 131209-130529  R value scan
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     // 141231-150203  R value scan
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     // 150202-  R value scan
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     // 150202 - R value scan
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     // 160106 - 160108  4170
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     // 160108  4170
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     // 160114  4170
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     // 160114  4170
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     // 160116  4170
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 }

Generated on Tue Nov 29 23:14:36 2016 for BOSS_7.0.2 by  doxygen 1.4.7