00088 {
00089
00090 MsgStream log(msgSvc(), name());
00091
00092
00093 std::string fullPath = "/Calib/DedxCal";
00094
00095 log << MSG::INFO<<"execute() fullPath = "<<fullPath<< endreq;
00096
00097 SmartDataPtr<CalibData::DedxCalibData> test(m_pCalibDataSvc, fullPath);
00098 int i;
00099 int j;
00100 double ddg;
00101 double ggs;
00102 double wireg;
00103 double gain;
00104 double zdep;
00105 double resol;
00106 double layerg;
00107 double m_id_doca;
00108 double m_iner_chi;
00109 double m_iner_gain;
00110 double m_iner_hits;
00111 double m_ip_eangle;
00112 double m_out_chi;
00113 double m_out_gain;
00114 double m_out_hits;
00115
00116 for(i=0;i<4;i++){
00117 for(j=0;j<43;j++){
00118 ddg = test->getddg(i,j);
00119 ggs = test->getggs(i,j);
00120 zdep = test->getzdep(i,j);
00121
00122
00123
00124
00125 }
00126 }
00127
00128 for(int i=0;i<1600;i++){
00129 m_id_doca = test->get_id_doca(i);
00130 m_iner_chi= test->get_iner_chi(i);
00131 m_iner_gain= test->get_iner_gain(i);
00132 m_iner_hits= test->get_iner_hits(i);
00133 m_ip_eangle= test->get_ip_eangle(i);
00134 m_out_chi= test->get_out_chi(i);
00135 m_out_gain= test->get_out_gain(i);
00136 m_out_hits= test->get_out_hits(i);
00137 std::cout<<"m_id_doca: "<<m_id_doca<<" m_iner_chi: "<<m_iner_chi<<" m_iner_gain: "<<m_iner_gain<<" m_iner_hits:"<<m_iner_hits<<" m_ip_eangle:"<<m_ip_eangle<<" m_out_chi="<<m_out_chi<<" m_out_gain:"<<m_out_gain<<" m_out_hits="<<m_out_hits<<std::endl;
00138 }
00139 gain = test->getgain();
00140
00141 double m_costheta;
00142 for(int i=0;i<80;i++){
00143 m_costheta = test->get_costheta(i);
00144 std::cout<<" m_costheta="<<m_costheta;
00145 if(i%5==0) std::cout<<"\n"<<std::endl;
00146 }
00147 double t0,dedx;
00148 for(int i=0;i<35;i++){
00149 t0=test->get_t0(i);
00150 dedx=test->get_dedx(i);
00151 std::cout<<"t0 ="<<t0<<" dedx="<<dedx;
00152 if(i%5==0) std::cout<<"\n"<<std::endl;
00153 }
00154
00155 resol = test->getresol();
00156
00157
00158 for(i=0;i<43;i++)
00159 { layerg = test -> getlayerg(i);
00160
00161 }
00162
00163 for(i=0;i<60;i++)
00164 { wireg = test -> getwireg(i);
00165
00166 }
00167
00168 for(i=0;i<20;i++){
00169 std::cout<<"hadron "<<i<<" = "<<test->get_hadron(i);
00170 }
00171
00172
00173
00174 return StatusCode::SUCCESS;
00175 }