GetFinalResult/macro/Run_54280_test.C

0001 #include "../GetFinalResult.h"
0002
0003 R__LOAD_LIBRARY(../libGetFinalResult.so)
0004
0005 int Run_54280_test(int condor_index)
0006 {
0007     std::vector<int> run_centrality_array = {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 70}; // note : 16 jobs in total
0008
0009     int SelectedMbin = run_centrality_array[condor_index];
0010
0011     int runnumber = 54280;
0012     // int SelectedMbin = 3; // note : 0, 1, ---- 10, 70, 100
0013     std::pair<double,double> vtxZ_range = {-10,10}; // note : cm
0014     bool isTypeA = false;
0015     bool ApplyGeoAccCorr = false;
0016     std::pair<bool, std::pair<double,double>> cut_EtaRange = {true, {-1.5, 1.5}}; // note : not used
0017
0018     std::string output_file_name_suffix = "";
0019
0020     std::string output_directory = "/sphenix/tg/tg01/commissioning/INTT/work/cwshih/seflgendata/run_54280_HR_Feb102025/Run6_EvtZFitWidthChange/EvtVtxZ/FinalResult_Nominal0Adc";
0021
0022     GetFinalResult * GFR = new GetFinalResult(
0023         runnumber,
0024         SelectedMbin,
0025         vtxZ_range,
0026         isTypeA,
0027         ApplyGeoAccCorr,
0028         cut_EtaRange,
0029
0030         output_file_name_suffix,
0031
0032         output_directory
0033     );
0034
0035     std::string data_input_directory = "/sphenix/tg/tg01/commissioning/INTT/work/cwshih/seflgendata/run_54280_HR_Feb102025/Run6_EvtZFitWidthChange/EvtVtxZ";
0036     std::string MC_input_directory = "/sphenix/user/ChengWei/sPH_dNdeta/Run24AuAuMC/Sim_HIJING_MDC2_ana467_20250225/Run7/EvtVtxZ";
0037
0038     std::string SetResultRangeFolderName = GFR -> GetSetResultRangeFolderName();
0039
0040     system(Form("if [ -d %s/completed/%s ]; then rm -r %s/completed/%s; fi;", output_directory.c_str(), SetResultRangeFolderName.c_str(), output_directory.c_str(), SetResultRangeFolderName.c_str()));
0041
0042     std::string baseline_subfolder = "TrackHist/noAdcCut/completed";
0043     std::string baseline_subfolder_MC_001 = "TrackHist_001/noAdcCut/completed";
0044     std::string baseline_subfolder_MC_002 = "TrackHist_002/noAdcCut/completed";
0045     // note : here are for preparing the reco. dNdEta
0046     GFR -> PrepareBaseLine(
0047         data_input_directory + "/" + baseline_subfolder,
0048         "Data_TrackHist_BcoFullDiffCut_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_00054280_merged_merged_001.root",
0049
0050         MC_input_directory + "/" + baseline_subfolder_MC_001,
0051         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_merged_001.root",
0052         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_merged_002.root"
0053     );
0054
0055
0056     std::string runseg_subfolder = "";
0057     GFR -> PrepareRunSegment(
0058         // data_input_directory + "/" + baseline_subfolder + "/" + runseg_subfolder,
0059         data_input_directory + "/" + baseline_subfolder,
0060         {
0061             "Data_TrackHist_BcoFullDiffCut_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_00054280_merged_merged_002.root"
0062             // "Data_TrackHist_BcoFullDiffCut_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_00054280_merged_002.root"
0063         },
0064
0065         MC_input_directory + "/" + baseline_subfolder_MC_002,
0066         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_merged_001.root",
0067         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_merged_002.root"
0068     ); // note : two times
0069
0070
0071     GFR -> PrepareDeltaPhiCut(
0072         {0.021, 0.031, 0.040, 0.060, 0.070},
0073         data_input_directory + "/" + baseline_subfolder,
0074         "Data_TrackHist_BcoFullDiffCut_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_00054280_merged_merged_001.root",
0075
0076         MC_input_directory + "/" + baseline_subfolder_MC_001,
0077         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_merged_001.root",
0078         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_merged_002.root"
0079     ); // note : two times
0080
0081
0082     // std::string ClusAdcCut_subfolder = "TrackHist_noAdcCut/completed";
0083     // GFR -> PrepareClusAdcCut(
0084     //     data_input_directory + "/" + ClusAdcCut_subfolder,
0085     //     {
0086     //         "Data_TrackHist_BcoFullDiffCut_VtxZQA_ClusQAAdc0PhiSize39_ColMulMask_00054280_merged.root"
0087     //     },
0088
0089     //     MC_input_directory + "/" + ClusAdcCut_subfolder,
0090     //     "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc0PhiSize39_ColMulMask_merged_001.root",
0091     //     "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc0PhiSize39_ColMulMask_merged_002.root"
0092     // );
0093
0094
0095     std::string NoClusAdcCut_subfolder = "TrackHist/baseline/completed";
0096     std::string NoClusAdcCut_subfolder_MC_001 = "TrackHist_001/noAdcCut/completed";
0097     GFR -> PrepareClusAdcCut(
0098         0,
0099         data_input_directory + "/" + NoClusAdcCut_subfolder,
0100         "Data_TrackHist_BcoFullDiffCut_VtxZQA_ClusQAAdc35PhiSize40_ColMulMask_00054280_merged_merged_001.root",
0101
0102         MC_input_directory + "/" + NoClusAdcCut_subfolder_MC_001,
0103         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_merged_001.root",
0104         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc0PhiSize40_ColMulMask_merged_002.root"
0105     );
0106     std::string ClusAdcCut50_subfolder = "TrackHist/baseline/completed";
0107     std::string ClusAdcCut50_subfolder_MC_001 = "TrackHist_001/baseline/completed";
0108     GFR -> PrepareClusAdcCut(
0109         1,
0110         data_input_directory + "/" + ClusAdcCut50_subfolder,
0111         "Data_TrackHist_BcoFullDiffCut_VtxZQA_ClusQAAdc35PhiSize40_ColMulMask_00054280_merged_merged_001.root",
0112
0113         MC_input_directory + "/" + ClusAdcCut50_subfolder_MC_001,
0114         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc35PhiSize40_ColMulMask_merged_001.root",
0115         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc35PhiSize40_ColMulMask_merged_002.root"
0116     );
0117
0118
0119
0120     std::string ClusPhiCut_subfolder = "TrackHist/noPhiCut/completed";
0121     std::string ClusPhiCut_subfolder_MC_001 = "TrackHist_001/noPhiCut/completed";
0122     GFR -> PrepareClusPhiCut(
0123         data_input_directory + "/" + ClusPhiCut_subfolder,
0124         {
0125             "Data_TrackHist_BcoFullDiffCut_VtxZQA_ClusQAAdc35PhiSize350_ColMulMask_00054280_merged_merged_001.root"
0126         },
0127
0128         MC_input_directory + "/" + ClusPhiCut_subfolder_MC_001,
0129         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc35PhiSize350_ColMulMask_merged_001.root",
0130         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc35PhiSize350_ColMulMask_merged_002.root"
0131     );
0132
0133     system(Form("mv %s/%s %s/completed", output_directory.c_str(), SetResultRangeFolderName.c_str(), output_directory.c_str()));
0134
0135     return 777;
0136 }