SystUnc/macro/Run_Final.C

0001 #include "../FinalResult.h"
0002 R__LOAD_LIBRARY(../libFinalResult.so)
0003
0004 int Run_Final()
0005 {
0006     const string CW_folder = "/sphenix/tg/tg01/commissioning/INTT/work/cwshih";
0007     const string data_folder_mother = "/seflgendata/run_54280_HR_Dec042024/completed/Run3";
0008
0009     int runnumber = 54280;
0010     int Mbin = 70;
0011     std::string StandardData_directory = CW_folder + data_folder_mother + "/EvtVtxZ/TrackHist/completed/dNdEta/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed";
0012     std::string StandardData_file_name = "Data_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00054280_00000_dNdEta.root";
0013     std::string StandardMC_directory   = CW_folder + data_folder_mother + "/EvtVtxZ/TrackHist/completed/dNdEta/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed";
0014     std::string StandardMC_file_name   = "MC_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00002_dNdEta.root"; // note : for generator level
0015     std::string sPH_label = "Internal";
0016
0017     // std::string Output_directory = "Not_given";
0018
0019     FinalResult * final = new FinalResult(
0020         runnumber,
0021         Mbin,
0022         StandardData_directory,
0023         StandardData_file_name,
0024         StandardMC_directory,
0025         StandardMC_file_name,
0026         sPH_label
0027     );
0028
0029     final -> SetEtaRange({-1.5, 1.5});
0030
0031     final -> SetCollisionStr({{0.2, 0.96}, "Au+Au #sqrt{s_{NN}} = 200 GeV"});
0032     final -> SetAnaDescription({{0.21, 0.9}, "Centrality [0-70]%, VtxZ [-10, 10] cm"});
0033
0034
0035     final -> SetFinal_Data_MC_text(
0036         {
0037             "Data (PHOBOS approach)",
0038             "HIJING (generator)"
0039         }
0040     );
0041
0042     std::string output_folder_name = final -> GetOutputFileName();
0043     std::cout << "output_folder_name = " << output_folder_name << std::endl;
0044
0045     final -> PrepareStatisticalError();
0046
0047     // Division : -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
0048     std::string RunSegMother_directory = CW_folder + data_folder_mother + "/EvtVtxZ/TrackHist/completed/merged_files_Data_TrackHist_BcoFullDiffCut_VtxZQA_ClusQAAdc35PhiSize500_ColMulMask_00054280_";
0049     final -> PrepareRunSegmentError(
0050         {
0051             RunSegMother_directory + "/dNdEta_001/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed/Data_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00054280_00000_dNdEta.root",
0052             RunSegMother_directory + "/dNdEta_002/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed/Data_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00054280_00000_dNdEta.root"
0053         }
0054     ); // note : run1, run2
0055
0056     // Division : -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
0057     std::string ClusAdcMother_directory = CW_folder + data_folder_mother + "/EvtVtxZ/TrackHist_NoClusQA/completed";
0058     final -> PrepareClusAdcError(
0059         {
0060             ClusAdcMother_directory + "/dNdEta/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed/Data_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00054280_00000_dNdEta.root"
0061         }
0062     );
0063
0064     // Division : -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
0065     final -> PrepareGeoOffsetError(
0066         "/sphenix/user/ChengWei/sPH_dNdeta/Run24AuAuMC/Sim_Ntuple_HIJING_ana443_20241102/GeoOffset_v1/completed/merged_result/FromdNdEta.root",
0067         CW_folder + data_folder_mother + "/EvtVtxZ/TrackHist/completed/dNdEta/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed/MC_PreparedNdEtaEach_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00001_dNdEta.root" // note : for the alpha correction
0068     );
0069
0070     // Division : -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
0071     std::string DeltaPhiMother_directory = CW_folder + data_folder_mother + "/EvtVtxZ/TrackHist/completed";
0072     final -> PrepareDeltaPhiError(
0073         {
0074             DeltaPhiMother_directory + "/dNdEta_DeltaPhi0p018/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed/Data_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00054280_00000_dNdEta.root",
0075             DeltaPhiMother_directory + "/dNdEta_DeltaPhi0p024/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed/Data_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00054280_00000_dNdEta.root"
0076         }
0077
0078     ); // note : 0.018 and 0.024
0079
0080     // Division : -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
0081     std::string ClusPhiSizeMother_directory = CW_folder + data_folder_mother + "/EvtVtxZ/TrackHist_PhiCut/completed";
0082     final -> PrepareClusPhiSizeError(
0083         {
0084             ClusPhiSizeMother_directory + "/dNdEta/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed/Data_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00054280_00000_dNdEta.root"
0085         }
0086     );
0087
0088     // Division : -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
0089     // std::string MCMergedMother_directory = CW_folder + data_folder_mother + "/EvtVtxZ/TrackHist/completed";
0090     // final -> PrepareMCMergedError(
0091     //     {
0092     //         MCMergedMother_directory + "/dNdEta_MCMergeError1/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed/Data_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00054280_00000_dNdEta.root",
0093     //         MCMergedMother_directory + "/dNdEta_MCMergeError2/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed/Data_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00054280_00000_dNdEta.root",
0094     //         MCMergedMother_directory + "/dNdEta_MCMergeError3/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed/Data_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00054280_00000_dNdEta.root",
0095     //         MCMergedMother_directory + "/dNdEta_MCMergeError4/dNdEta_AllSensor_GeoAccCorr_VtxZ10_Mbin70/completed/Data_PreparedNdEtaEach_AlphaCorr_GeoAccCorr_AllSensor_VtxZ10_Mbin70_00054280_00000_dNdEta.root"
0096     //     }
0097     // );
0098
0099     // Division : -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
0100     final -> PrepareFinalError();
0101     final -> PrepareFinalResult();
0102     final -> EndRun();
0103
0104     return 9;
0105 }