93 extTrkNum = extTrkNum + ntrk;
112 TofTrackVec::iterator iter_tof = tofTrackVec->begin();
113 for( ; iter_tof!=tofTrackVec->end(); iter_tof++ ) {
183 std::cout<<
"############################# TofRec Summary #####################################" << std::endl;
184 std::cout<<
" ===> Number of Extrapolated tracks: " << extTrkNum <<std::endl;
185 std::cout<<
" ===> Number of valid Extraploated tracks: " << (inner1+outer1+double1+east1+west1) <<
" Barrel: " << (inner1+outer1+double1) <<
" Endcap: " << (east1+west1) << std::endl;
186 std::cout<<
" ===> TofRec: number of barrel tracks: " << (inner3+outer3+double3) <<
" Efficiency: " << setprecision(4) << (inner3+outer3+double3)*100.0/((inner1+outer1+double1)*1.0) <<
"%" << std::endl;
187 std::cout<<
" ===> Inner Layer hit only in total: " << inner1 <<
" \t" << inner2 <<
" \t" << inner3 <<
" \t" << inner4 << std::endl;
188 std::cout<<
" ===> Outer Layer hit only in total: " << outer1 <<
"\t" << outer2 <<
" \t" << outer3 <<
" \t" << outer4 << std::endl;
189 std::cout<<
" ===> Double Layer hit in total: " << double1 <<
" \t" << double2 <<
" \t" << double3 <<
" \t" << double4 << std::endl;
190 std::cout<<
" ===> TofRec: number of endcap tracks: " << (east3+west3) <<
" Efficiency: " << setprecision(4) << (east3+west3)*100.0/((east1+west1)*1.0) <<
"%" << std::endl;
191 std::cout<<
" ===> East endcap tracks in total: " << east1 <<
" \t" << east2 <<
" \t" << east3 <<
" \t" << east4 << std::endl;
192 std::cout<<
" ===> West endcap tracks in total: " << west1 <<
" \t" << west2 <<
" \t" << west3 <<
" \t" << west4 << std::endl;
193 std::cout<<
" ===> The invalid Extend tracks in total: "<< nohit1 <<
" \t" << nohit2 <<
" \t" << nohit3 <<
" \t" << nohit4 <<std::endl;
194 std::cout<<
"############################# End Summary #######################################" <<std::endl;
216 std::cout<<
"$$$$$$$ Summary of Bhabha Selection Criteria $$$$$$$$$$$$$$$$$$$$$"<<std::endl;
217 std::cout <<
"0 Total umber of events: "<<nbb[0]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[0])*100.0/double(nbb[0])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[0])*100.0/double(nbb[0])<<
"%"<<std::endl;
218 std::cout <<
"1 Event Start Time: "<<nbb[1]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[1])*100.0/double(nbb[0])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[1])*100.0/double(nbb[0])<<
"%"<<std::endl;
219 std::cout <<
"2 N_ext_track = 2: "<<nbb[2]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[2])*100.0/double(nbb[1])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[2])*100.0/double(nbb[0])<<
"%"<<std::endl;
220 std::cout <<
"3 N_mdc_track = 2: "<<nbb[3]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[3])*100.0/double(nbb[2])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[3])*100.0/double(nbb[0])<<
"%"<<std::endl;
221 std::cout <<
"4 EMC reconstructed: "<<nbb[4]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[4])*100.0/double(nbb[3])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[4])*100.0/double(nbb[0])<<
"%"<<std::endl;
222 std::cout <<
"5 N_shower >= 2: "<<nbb[5]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[5])*100.0/double(nbb[4])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[5])*100.0/double(nbb[0])<<
"%"<<std::endl;
223 std::cout <<
"6 Total charge = 0: "<<nbb[6]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[6])*100.0/double(nbb[5])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[6])*100.0/double(nbb[0])<<
"%"<<std::endl;
224 std::cout <<
"7 Back to back: "<<nbb[7]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[7])*100.0/double(nbb[6])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[7])*100.0/double(nbb[0])<<
"%"<<std::endl;
225 std::cout <<
"8 Kalman Filter: "<<nbb[8]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[8])*100.0/double(nbb[7])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[8])*100.0/double(nbb[0])<<
"%"<<std::endl;
226 std::cout <<
"9 Distance of ext-shower: "<<nbb[9]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[9])*100.0/double(nbb[8])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[9])*100.0/double(nbb[0])<<
"%"<<std::endl;
227 std::cout <<
"10 Vertex cut of 1st trk: "<<nbb[10]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[10])*100.0/double(nbb[9])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[10])*100.0/double(nbb[0])<<
"%"<<std::endl;
228 std::cout <<
"11 Vertex cut of 2nd trk: "<<nbb[11]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[11])*100.0/double(nbb[10])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[11])*100.0/double(nbb[0])<<
"%"<<std::endl;
229 std::cout <<
"12 Delta phi: "<<nbb[12]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[12])*100.0/double(nbb[11])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[12])*100.0/double(nbb[0])<<
"%"<<std::endl;
230 std::cout <<
"13 Energy of Shower: "<<nbb[13]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[13])*100.0/double(nbb[12])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[13])*100.0/double(nbb[0])<<
"%"<<std::endl;
231 if( calibData==
"Bhabha" ) {
232 std::cout <<
"14 Energy of others: "<<nbb[14]<<
"\t"<<
"Ratio1: "<<setprecision(4)<<double(nbb[14])*100.0/double(nbb[13])<<
"%"<<
"\t"<<
"Ratio2: "<<setprecision(4)<<double(nbb[14])*100.0/double(nbb[0])<<
"%"<<std::endl;
234 else if( calibData==
"Dimu") {
236 std::cout<<
"$$$$$$$ End of Summary $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$"<<std::endl;
const unsigned int nParBhabha
void finalBhabha(std::string calibData)
void setExtTrackNum(unsigned int ntrk)
void setTrack1Col(std::vector< TofTrack * > *&tofTrackVec)
void setTrack1(TofTrack *&tof)
void addNumber(unsigned int i)
void setTrack3(TofTrack *&tof)
void setTrack2(TofTrack *&tof)
void setTrack4(TofTrack *&tof)
ExtTrackCase hitCase() const