254 {
255 std::cout<<"############################# TofRec Summary #####################################" << std::endl;
256 std::cout<<" ===> Number of Extrapolated tracks: " << extTrkNum <<std::endl;
257 std::cout<<" ===> Number of valid Extraploated tracks: " << (inner1+outer1+double1+east1+west1+mrpc_east1_1+mrpc_east2_1+mrpc_west1_1+mrpc_west2_1) << " Barrel: " << (inner1+outer1+double1) << " Endcap: " << (east1+west1+mrpc_east1_1+mrpc_east2_1+mrpc_west1_1+mrpc_west2_1) << std::endl;
258 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;
259 std::cout<<" ===> Inner Layer hit only in total: " << inner1 << " \t" << inner2 << " \t" << inner3 << " \t" << inner4 << std::endl;
260 std::cout<<" ===> Outer Layer hit only in total: " << outer1 << "\t" << outer2 << " \t" << outer3 << " \t" << outer4 << std::endl;
261 std::cout<<" ===> Double Layer hit in total: " << double1 << " \t" << double2 << " \t" << double3 << " \t" << double4 << std::endl;
262 std::cout<<" ===> TofRec: number of endcap tracks: " << (east3+west3+mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3) << " Efficiency: " << setprecision(4) << (east3+west3+mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3)*100.0/((east1+west1+mrpc_east1_1+mrpc_east2_1+mrpc_west1_1+mrpc_west2_1)*1.0) << "%" << std::endl;
263 std::cout<<" ===> East endcap tracks in total: " << east1 << " \t" << east2 << " \t" << east3 << " \t" << east4 << std::endl;
264 std::cout<<" ===> West endcap tracks in total: " << west1 << " \t" << west2 << " \t" << west3 << " \t" << west4 << std::endl;
265 std::cout<<" ===> East 1 mrpc tracks in total: " << mrpc_east1_1 << " \t" << mrpc_east1_2 << " \t" << mrpc_east1_3 << " \t" << mrpc_east1_4 << std::endl;
266 std::cout<<" ===> West 1 mrpc tracks in total: " << mrpc_west1_1 << " \t" << mrpc_west1_2 << " \t" << mrpc_west1_3 << " \t" << mrpc_west1_4 << std::endl;
267 std::cout<<" ===> East 2 mrpc tracks in total: " << mrpc_east2_1 << " \t" << mrpc_east2_2 << " \t" << mrpc_east2_3 << " \t" << mrpc_east2_4 << std::endl;
268 std::cout<<" ===> West 2 mrpc tracks in total: " << mrpc_west2_1 << " \t" << mrpc_west2_2 << " \t" << mrpc_west2_3 << " \t" << mrpc_west2_4 << std::endl;
269
270 std::cout<<" ===> The invalid Extend tracks in total: "<< nohit1 << " \t" << nohit2 << " \t" << nohit3 << " \t" << nohit4 <<std::endl;
271
272 if((mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3) >0)
273 {
274 std::cout<<" ==> MRPC extrapolation (before calibration): " << endl;
275 std::cout<<" Determine transmision time without extrapolation : "<< setprecision(3) << (double)mrpc_trans_with_dbs_true_3*100.0/((double)(mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3)*1.0)<< " % "<<endl ;
276 std::cout<<" Determine transmision time with extrapolation : "<< setprecision(3) << (double)mrpc_trans_with_ext_true_3*100.0/((double)(mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3)*1.0)<< " % "<<endl ;
277 std::cout<<" Determine transmision time with average value : "<< setprecision(3) <<
278 ((double)(mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3)-mrpc_trans_with_dbs_true_3-mrpc_trans_with_ext_true_3)*100.0/((double)(mrpc_east1_3+mrpc_east2_3+mrpc_west1_3+mrpc_west2_3)*1.0)<< " % "<<endl ;
279 }
280 if((mrpc_east1_4+mrpc_east2_4+mrpc_west1_4+mrpc_west2_4) >0)
281 {
282 std::cout<<" ==> MRPC extrapolation (after calibration): " << endl;
283 std::cout<<" Determine transmision time without extrapolation : "<< setprecision(3) << (double)mrpc_trans_with_dbs_true_4*100.0/((double)(mrpc_east1_4+mrpc_east2_4+mrpc_west1_4+mrpc_west2_4)*1.0)<< " % "<<endl ;
284 std::cout<<" Determine transmision time with extrapolation : "<< setprecision(3) << (double)mrpc_trans_with_ext_true_4*100.0/((double)(mrpc_east1_4+mrpc_east2_4+mrpc_west1_4+mrpc_west2_4)*1.0)<< " % "<<endl ;
285 std::cout<<" Determine transmision time with average value : "<< setprecision(3) <<
286 ((double)(mrpc_east1_4+mrpc_east2_4+mrpc_west1_4+mrpc_west2_4)-mrpc_trans_with_dbs_true_4-mrpc_trans_with_ext_true_4)*100.0/((double)(mrpc_east1_4+mrpc_east2_4+mrpc_west1_4+mrpc_west2_4)*1.0)<< " % "<<endl ;
287 }
288
289
290
291 std::cout<<"############################# End Summary #######################################" <<std::endl;
292
293}