analysis_INTT_CW/old_code/INTTDSTchain_test.C

0001 #include "INTTDSTchain.C"
0002
0003 void INTTDSTchain_test()
0004 {
0005     string folder_direction = "/sphenix/user/ChengWei/sPH_dNdeta/dNdEta_INTT_MC/";
0006     vector <string> file_list = {"INTTRecoClusters_test_400_0.root","INTTRecoClusters_test_400_2.root"};
0007
0008     TChain * chain_in = new TChain("EventTree");
0009     INTTDSTchain inttDSTread(chain_in,folder_direction,file_list);
0010     std::printf("inttDSTread N event : %lli \n", chain_in->GetEntries());
0011
0012     for (int i = 399; i < 401; i++) {
0013         inttDSTread.LoadTree(i);
0014         inttDSTread.GetEntry(i);
0015
0016         cout<<"----------------------------------------------- -----------------------------------------------"<<endl;
0017         printf(" event : %i, NTruthVtx : %i, NClus : %i, NTrkrhits : %i \n", inttDSTread.event, inttDSTread.NTruthVtx, inttDSTread.NClus, inttDSTread.NTrkrhits);
0018         printf(" size of ClusX : %i \n", inttDSTread.ClusX -> size());
0019         printf("triggered VTX ? : %.2f, %.2f, %.2f \n", inttDSTread.TruthPV_trig_x, inttDSTread.TruthPV_trig_y, inttDSTread.TruthPV_trig_z);
0020         for (int vtx_i = 0; vtx_i < inttDSTread.TruthPV_x -> size(); vtx_i++){
0021             printf("VTX vec ID-%i : %.2f, %.2f, %.2f \n", vtx_i, inttDSTread.TruthPV_x -> at(vtx_i), inttDSTread.TruthPV_y -> at(vtx_i), inttDSTread.TruthPV_z -> at(vtx_i));
0022         }
0023         // cout<<"~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ "<<endl;
0024         // for (int clu_i = 0; clu_i < ClusX -> size(); clu_i++){
0025         //     printf("cluster %zu : X %.2f, Y %.2f, Z %.2f, sizePhi : %f, sizeZ : %f,  \n", clu_i, ClusX -> at(clu_i), ClusY -> at(clu_i), ClusZ -> at(clu_i), ClusPhiSize -> at(clu_i), ClusZSize -> at(clu_i));
0026         // }
0027     }
0028
0029 }