ColumnCheck/macro/Run_ColumnCheck_MC.C

0001 #include "../ColumnCheck.h"
0002
0003 R__LOAD_LIBRARY(../libColumnCheck.so)
0004
0005 TH2D * GetGoodColMap (std::string ColMulMask_map_dir_in, std::string ColMulMask_map_file_in, std::string map_name_in)
0006 {
0007   TFile * f = TFile::Open(Form("%s/%s", ColMulMask_map_dir_in.c_str(), ColMulMask_map_file_in.c_str()));
0008   TH2D * h = (TH2D*)f->Get(map_name_in.c_str());
0009   return h;
0010 }
0011
0012 void Run_ColumnCheck_MC(
0013   int process_id = 0,
0014   int run_num = 54280,
0015   int nevents = -1,
0016   string input_directory = "/sphenix/user/ChengWei/INTT/INTT/general_codes/CWShih/INTTBcoResolution/macro",
0017   string input_filename = "file_list_54280_intt.txt",
0018   string output_directory = "/sphenix/tg/tg01/commissioning/INTT/work/cwshih/seflgendata/run_54280/completed/BCO_check",
0019
0020   // todo : modify here
0021   std::string output_file_name_suffix = "",
0022   std::pair<double, double> vertexXYIncm = {-0.0217356, 0.223402},
0023   double SetMbinFloat = 70, // note : 0 - 100
0024
0025   std::pair<double, double> VtxZRange = {-30, 30},
0026   bool IsZClustering = false,
0027   bool BcoFullDiffCut = false,
0028   std::pair<bool, std::pair<double, double>> isClusQA = {true, {35,39}}, // note : {adc, phi size}
0029
0030   bool ColMulMask = false,
0031   std::string ColMulMask_map_dir = "/sphenix/tg/tg01/commissioning/INTT/work/cwshih/seflgendata/run_54280_HR_Dec042024/completed/Run3/EvtVtxZ/ColumnCheck/completed/MulMap/completed",
0032   std::string ColMulMask_map_file = "MulMap_BcoFullDiffCut_Mbin50_VtxZ-30to30cm_ClusQAAdc35PhiSize500_00054280.root"
0033 )
0034 {
0035
0036   ColumnCheck * CC = new ColumnCheck(
0037     process_id,
0038     run_num,
0039     nevents,
0040     input_directory,
0041     input_filename,
0042     output_directory,
0043
0044     output_file_name_suffix,
0045     vertexXYIncm,
0046     SetMbinFloat,
0047
0048     VtxZRange,
0049     IsZClustering,
0050     BcoFullDiffCut,
0051     isClusQA, // note : {adc, phi size}
0052
0053     ColMulMask
0054   );
0055
0056   string final_output_file_name = CC->GetOutputFileName();
0057   cout<<"final_output_file_name: "<<final_output_file_name<<endl;
0058   system(Form("if [ -f %s/completed/%s ]; then rm %s/completed/%s; fi;", output_directory.c_str(), final_output_file_name.c_str(), output_directory.c_str(), final_output_file_name.c_str()));
0059
0060   if (ColMulMask){
0061     CC -> SetGoodColMap(
0062       GetGoodColMap(ColMulMask_map_dir, ColMulMask_map_file, CC -> GetGoodColMapName())
0063     );
0064   }
0065
0066   CC -> MainProcess();
0067   CC -> EndRun();
0068
0069   system(Form("mv %s/%s %s/completed", output_directory.c_str(), final_output_file_name.c_str(), output_directory.c_str()));
0070
0071   return;
0072 }