#include <BesMdcRes.hh>
Public Member Functions | |
| BesMdcRes (std::string &path) | |
| BesMdcRes (std::string &path) | |
| double | getD_dD (int ii, int jj) |
| double | getD_dD (int ii, int jj) |
| double | getD_effi (int ii, int jj) |
| double | getD_effi (int ii, int jj) |
| double | getD_iEntr (int ii, int jj, int kk) |
| double | getD_iEntr (int ii, int jj, int kk) |
| double | getD_lay (int ii, int jj) |
| double | getD_lay (int ii, int jj) |
| double | getD_theta (int ii, int jj) |
| double | getD_theta (int ii, int jj) |
| double | getM_dD (int ii, int jj) |
| double | getM_dD (int ii, int jj) |
| double | getM_iEntr (int ii, int jj, int kk) |
| double | getM_iEntr (int ii, int jj, int kk) |
| double | getM_lay (int ii, int jj) |
| double | getM_lay (int ii, int jj) |
| double | getM_theta (int ii, int jj) |
| double | getM_theta (int ii, int jj) |
| ~BesMdcRes () | |
| ~BesMdcRes () | |
Private Attributes | |
| double | par_lay1 [43][10] |
| double | par_lay11 [43][20] |
| double | par_lay2 [43][5] |
| double | par_lay3 [43][10] |
| double | par_lay33 [43][20] |
| double | par_lay4 [43][5] |
| double | par_lay5 [43][90] |
| double | par_lay6 [43][2][18] |
| double | par_lay7 [43][2][18] |
|
|
00007 {
00008 int layerId, bin,layId ;
00009 double sigma;
00010 double cons1,sigma1,cons2,sigma2,Sig,sigmaE;
00011 double effi,m_nhit_Raw,m_nhit_Exp;
00012
00013
00014 /*
00015 //read lay efficiency parameter
00016 ifstream infpar5;
00017 string logeff_raw= path+"/data/logeffi_0_42.dat";
00018 infpar5.open(logeff_raw.c_str());
00019 if( !infpar5.good()){
00020 std::cout<<"Error,mdc parameter file4 Data-effi-lay not exist "<<std::endl;
00021 exit (0);
00022 }else{
00023 std::cout<<" open mdc parameter file4 Data-effi-lay "<<std::endl;
00024 for(int lay =0; lay<43;lay++){
00025 infpar5 >> lay;
00026 for(int dD = 0; dD<91;dD++ ){
00027 infpar5 >> dD >> effi >> m_nhit_Raw >> m_nhit_Exp ;
00028 par_lay5[lay][dD] = effi;
00029 //cout<<"lay : "<<lay<<" dD "<<dD <<" =" <<effi<<endl;
00030 }
00031 }// end of lay loop
00032 }
00033
00034 */
00035 //read resolution Vs enter angle and doca
00036 ifstream infpar6;
00037 string logres_iEntr = path+"/data/logrVetr.dat";
00038 infpar6.open(logres_iEntr.c_str());
00039 if(!infpar6.good()){
00040 std::cout<<" Error , mdc parameters file6 Data_res_iEnter not exist "<<std::endl;
00041 exit (0);
00042 }else{
00043 std::cout<<" open mdc parameter file6 Data_res_iEnter "<<std::endl;
00044 for(int lay =0;lay<43;lay++){
00045 infpar6 >> lay;
00046 for(int iEntr =0;iEntr<2;iEntr++){
00047 infpar6 >> iEntr;
00048 for(int dD =0;dD<18;dD++){
00049 infpar6 >> dD >> sigmaE;
00050 par_lay6[lay][iEntr][dD] = sigmaE;
00051 // cout<<" lay : "<<lay<<" iEntr : "<<iEntr<<" dD : "<<dD<<" = "<<sigmaE<<endl;
00052 }
00053 }
00054 }// end of lay loop
00055 }
00056
00057 //read resolution Vs enter angle and doca
00058 ifstream infpar7;
00059 string logres_M_iEntr = path + "/mc/logrVetr.dat";
00060 infpar7.open(logres_M_iEntr.c_str());
00061 if(!infpar7.good()){
00062 std::cout<<"Error,mdc parameters file7 Mc_res_iEntr not exist"<<std::endl;
00063 exit (0);
00064 }else{
00065 std::cout<<"open mdc parameters file7 Mc_res_iEntr "<<std::endl;
00066 for(int lay =0;lay<43;lay++){
00067 infpar7 >> lay;
00068 for(int iEntr =0;iEntr<2;iEntr++){
00069 infpar7 >> iEntr;
00070 for(int dD =0;dD<18;dD++){
00071 infpar7 >> dD >> sigmaE;
00072 par_lay7[lay][iEntr][dD] = sigmaE;
00073 // std::cout<<"par_lay7["<<lay<<"]["<<iEntr<<"]["<<dD<<"] : "<<sigmaE<<std::endl;
00074 }
00075 }
00076 }//end of lay loop
00077 }
00078
00079 }
|
|
|
00008 {};
|
|
|
|
|
|
00008 {};
|
|
||||||||||||
|
00010 {return par_lay1[ii][jj];}
|
|
||||||||||||
|
00010 {return par_lay1[ii][jj];}
|
|
||||||||||||
|
00016 {return par_lay5[ii][jj];}
|
|
||||||||||||
|
00016 {return par_lay5[ii][jj];}
|
|
||||||||||||||||
|
00017 {return par_lay6[ii][jj][kk];}
|
|
||||||||||||||||
|
00017 {return par_lay6[ii][jj][kk];}
|
|
||||||||||||
|
00011 {return par_lay2[ii][jj];}
|
|
||||||||||||
|
00011 {return par_lay2[ii][jj];}
|
|
||||||||||||
|
00012 {return par_lay11[ii][jj];}
|
|
||||||||||||
|
00012 {return par_lay11[ii][jj];}
|
|
||||||||||||
|
00013 {return par_lay3[ii][jj];}
|
|
||||||||||||
|
00013 {return par_lay3[ii][jj];}
|
|
||||||||||||||||
|
00018 {return par_lay7[ii][jj][kk];}
|
|
||||||||||||||||
|
00018 {return par_lay7[ii][jj][kk];}
|
|
||||||||||||
|
00014 {return par_lay4[ii][jj];}
|
|
||||||||||||
|
00014 {return par_lay4[ii][jj];}
|
|
||||||||||||
|
00015 {return par_lay33[ii][jj];}
|
|
||||||||||||
|
00015 {return par_lay33[ii][jj];}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.3.9.1