GetFinalResult/macro/Run_54280_EPOS.C

0001 #include "../GetFinalResult.h"
0002
0003 R__LOAD_LIBRARY(../libGetFinalResult.so)
0004
0005 int Run_54280_EPOS(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/user/ChengWei/sPH_dNdeta/Run24AuAuMC/Sim_EPOS_MDC2_ana472_20250310/Run7/EvtVtxZ/FinalResult_10cm_Pol2BkgFit_DeltaPhi0p026";
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     GFR -> SetBaseLineDeltaPhiCut({-0.026,0.026}); // todo : change here
0036
0037     std::string data_input_directory = "/sphenix/tg/tg01/commissioning/INTT/work/cwshih/seflgendata/run_54280_HR_Feb102025/Run6_EvtZFitWidthChange/EvtVtxZ";
0038     std::string MC_input_directory = "/sphenix/user/ChengWei/sPH_dNdeta/Run24AuAuMC/Sim_EPOS_MDC2_ana472_20250310/Run7/EvtVtxZ";
0039
0040     std::string SetResultRangeFolderName = GFR -> GetSetResultRangeFolderName();
0041
0042     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()));
0043
0044     std::string baseline_subfolder = "TrackHist/baseline/completed";
0045     std::string baseline_subfolder_MC_001 = "TrackHist_DeltaPhi0p04_001/baseline/completed"; // note : for the first half of data
0046     // std::string baseline_subfolder_MC_002 = "TrackHist_002/baseline/completed"; // note : for the second half of data
0047     // note : here are for preparing the reco. dNdEta
0048     GFR -> PrepareBaseLine(
0049         data_input_directory + "/" + baseline_subfolder,
0050         "Data_TrackHist_BcoFullDiffCut_VtxZQA_ClusQAAdc35PhiSize40_ColMulMask_00054280_merged_merged_001.root",
0051
0052         MC_input_directory + "/" + baseline_subfolder_MC_001,
0053         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc35PhiSize40_ColMulMask_merged_001.root", // note : derive the alpha correction
0054         "MC_TrackHist_vtxZReweight_VtxZQA_ClusQAAdc35PhiSize40_ColMulMask_merged_002.root"  // note : closure test
0055     );
0056
0057
0058     system(Form("mv %s/%s %s/completed", output_directory.c_str(), SetResultRangeFolderName.c_str(), output_directory.c_str()));
0059
0060     return 777;
0061 }