StRoot: StEmcPool/StEmcOfflineCalibrationMaker/macros/CalibrationHelperFunctions.cxx Source File

00001 #include "CalibrationHelperFunctions.h"
00002 #include <cmath>
00003 #include <iostream>
00004 #include <fstream>
00005 using namespace std;
00006 
00007 CalibrationHelperFunctions::CalibrationHelperFunctions(){       
00008         ifstream read_mapping;
00009         
00010         //try to read mapping from local disk
00011         read_mapping.open("btow_mapping.20030401.0.txt");
00012         if(!read_mapping.good()){
00013                 cout<<"CalibrationHelperFunctions - FATAL ERROR: could not read mapping!.  Answers will be nonsense"<<endl;
00014         }
00015         else{
00016                 while (1) {
00017                         int id;
00018                         double eta, phi, daqId, crate, crate_seq, module, etaDiv, subDiv, tdc, pmt, position, patchId;
00019                         read_mapping >> id >> eta >> phi >> daqId >> crate >> crate_seq >> module >> etaDiv >> subDiv >> tdc >> pmt >> position >> patchId;
00020                         if (!read_mapping.good() || id > ntowers) break;
00021                         tower_eta[id-1] = eta;
00022                         tower_phi[id-1] = phi;
00023                         tower_theta[id-1] = 2*atan(exp(-1*eta));
00024                 }
00025         }
00026 }
00027 
00028 CalibrationHelperFunctions::~CalibrationHelperFunctions(){
00029         
00030 }
00031 
00032 bool CalibrationHelperFunctions::isGoodTower2004(int towerid){
00033         switch(towerid){
00034                 default:
00035                         return true;
00036         }
00037 }
00038 
00039 bool CalibrationHelperFunctions::isGoodTower2005(int towerid){
00040         if(isBadIso2005(towerid) || isBadPoverE2005(towerid)) return false;
00041         
00042         switch(towerid){
00043                 //bad towers
00044                 case(36): case(45): case(46): case(47): case(48): case(50): case(59): case(63): case(139): case(220):
00045                 case(297): case(303): case(389): case(411): case(426): case(446): case(448): case(485): case(492):
00046                 case(528): case(559): case(565): case(638): case(671): case(691): case(738): case(744): case(761):
00047                 case(846): case(855): case(875): case(897): case(916): case(1028): case(1080): case(1100): case(1104):
00048                 case(1160): case(1176): case(1220): case(1280): case(1317): case(1397): case(1398): case(1400): 
00049                 case(1417): case(1419): case(1420): case(1440): case(1471): case(1505): case(1612): case(1668):
00050                 case(1676): case(1679): case(1720): case(1764): case(1856): case(1866): case(1880): case(1909): 
00051                 case(2069): case(2074): case(2075): case(2079): case(2092): case(2161): case(2168): case(2241):
00052                 case(2257): case(2285): case(2378): case(2394): case(2403): case(2409): case(2458): case(2460):
00053                 case(2470): case(2504): case(2529): case(2592): case(2610): case(2658): case(2794): case(2834):
00054                 case(2835): case(2863): case(2865): case(2897): case(2961): case(2969): case(2972): case(3069):
00055                 case(3086): case(3093): case(3097): case(3232): case(3255): case(3283): case(3289): case(3309): 
00056                 case(3372): case(3407): case(3515): case(4171): case(4217): case(4232): case(4240): case(4312): 
00057                 case(4357): case(4377): case(4506): case(4507): case(4508): case(4514): case(4519): case(4543): 
00058                 case(4558): case(4560): case(4580): case(4585): case(4639): case(4660): case(4671): case(4678): 
00059                 case(4766): case(4768):
00060                         return false;
00061                         
00062                 default:
00063                         return true;
00064         }
00065 }
00066 
00067 bool CalibrationHelperFunctions::isBadIso2005(int id){
00068         switch(id){
00069                 //isolation cut failures
00070                 case(34): case(35): case(266): case(267): case(286): case(287): case(561): case(562):
00071                 case(615): case(616): case(633): case(653): case(637): case(657): case(649): case(650):
00072                 case(673): case(674): case(813): case(814): case(837): case(857): case(953): case(954): 
00073                 case(1026): case(1046): case(1353): case(1354): case(1574): case(1575): case(1753): case(1773):
00074                 case(1765): case(1766): case(1897): case(1898): case(2073): case(2093): case(2077): case(2097): 
00075                 case(2440): case(2460): case(2589): case(2590): case(3070): case(3071): case(3494): case(3495):
00076                 case(4677): case(4678):
00077                         return true;
00078                         
00079                 default:
00080                         return false;
00081         }
00082 }
00083 
00084 bool CalibrationHelperFunctions::isBadPoverE2005(int id){
00085         switch(id){
00086                 //poverE failures
00087                 case(30): case(95): case(108): case(162): case(308): case(533): case(555): case(750):
00088                 case(762): case(779): case(873): case(882): case(899): case(1024): case(1130): case(1132): 
00089                 case(1197): case(1204): case(1217): case(1237): case(1257): case(1294): case(1306): case(1375):
00090                 case(1434): case(1487): case(1537): case(1709): case(1984): case(2043): case(2162): case(2339):
00091                 case(2392):
00092                         return true;
00093                         
00094                 default:
00095                         return false;
00096         }
00097 }
00098 
00099 bool CalibrationHelperFunctions::isGoodTower2006(int id){
00100         switch(id){
00101                 //copy bad towers from 2005, remove a bunch of good ones
00102                 case(50): case(139): case(220):
00103                 case(389): case(411): case(426): case(446): case(485): case(492):
00104                 case(638): case(671):
00105                 case(846): case(855): case(875): case(916): case(1028): case(1080): case(1100):
00106                 case(1160): case(1176): case(1220): case(1280): case(1397):  
00107                 case(1612): case(1668):
00108                 case(1720): case(1856): case(1880):  
00109                 case(2074): case(2075): case(2079): case(2092): case(2168):
00110                 case(2257): case(2394): case(2403): case(2409): case(2458):
00111                 case(2470): case(2504): case(2529): case(2610): case(2658): case(2794): case(2834):
00112                 case(2865): case(2897): case(2961): case(2969): case(2972):
00113                 case(3097): case(3255): case(3289):
00114                 case(3372): case(3515): case(4171): case(4217): case(4240):
00115                 case(4357): case(4377): case(4506): case(4507): case(4508): case(4514): case(4543): 
00116                 case(4560): case(4580): case(4585): case(4671):
00117                 case(4766): case(4768):
00118                         return false;
00119                         
00120                 //copy bad iso from 2005, dropped 1897-8 (should this have been 1877-8?)
00121                 case(34): case(35): case(266): case(267): case(286): case(287): case(561): case(562):
00122                 case(615): case(616): case(633): case(653): case(637): case(657): case(649): case(650):
00123                 case(673): case(674): case(813): case(814): case(837): case(857): case(953): case(954): 
00124                 case(1026): case(1046): case(1353): case(1354): case(1574): case(1575): case(1753): case(1773):
00125                 case(1765): case(1766): case(2073): case(2093): case(2077): case(2097): 
00126                 case(2440): case(2460): case(2589): case(2590): case(3070): case(3071): case(3494): case(3495):
00127                 case(4677): case(4678):
00128                         return false;
00129                 
00130                 //new bad towers
00131                 case(240): case(390): case(391): case(392): case(409): case(410): case(412): case(504): case(541): case(594):
00132                 case(629): case(639): case(647): case(681): case(749): case(760): case(839): case(840): case(844): 
00133                 case(859): case(873): case(880): case(933): case(1018): case(1125): case(1142): case(1143): 
00134                 case(1158): case(1159): case(1161): case(1162): case(1163): case(1171): case(1180): case(1198):
00135                 case(1217): case(1224): case(1225): case(1237): case(1240): case(1244): case(1250):
00136                 case(1301): case(1319): case(1321): case(1341): case(1342): case(1348): case(1375): case(1381): case(1401):
00137                 case(1422): case(1507): case(1588): case(1608): case(1654): case(1732): case(1779): case(1838): 
00138                 case(1892): case(1893): case(1949): case(1985): case(2005): case(2021): case(2025): case(2070): case(2085):
00139                 case(2094): case(2095): case(2101): case(2105): case(2108): case(2116): case(2177): case(2188): case(2196):
00140                 case(2222): case(2262): case(2301): case(2305): case(2337): case(2453): case(2580): case(2633): case(2652):
00141                 case(2727): case(3007): case(3017): case(3154): case(3171): case(3220): case(3231): case(3287): case(3290):
00142                 case(3296): case(3328): case(3333): case(3493): case(3508): case(3544): case(3557): case(3588): case(3604):
00143                 case(3611): case(3653): case(3678): case(3679): case(3709): case(3715): case(3727): case(3745): case(3746):
00144                 case(3761): case(3769): case(3795): case(3803): case(3986): case(4014): case(4105): case(4016): case(4107): 
00145                 case(4018): case(4019): case(4020): case(4053): case(4054): case(4055): case(4056): case(4057): case(4075): 
00146                 case(4077): case(4078): case(4080): case(4100): case(4119): case(4120): case(4130): case(4174): case(4326):
00147                 case(4388): case(4459): case(4464): case(4505): case(4549): case(4569): case(4596): case(4672): case(4684):
00148                 case(4765):
00149                         return false;
00150                         
00151                 //a few more I missed the first time around
00152                 case(341):  case(900):  case(1044): case(1063): case(1078): case(1850): case(1877): case(1878): case(3738):
00153                 case(4015): case(4017): case(1048): case(3690): case(3718):
00154                         return false;
00155                         
00156                 default:
00157                         return true;
00158         }
00159 }