File indexing completed on 2025-08-06 08:12:37
0001 #include "sPhenixStyle.C"
0002
0003 int full_comparison()
0004 {
0005 SetsPhenixStyle();
0006
0007 std::string call_hist_name = "h2D_Mbin_NClus_Self";
0008
0009 std::string Data_Run54280_file_directory = "/sphenix/tg/tg01/commissioning/INTT/work/cwshih/seflgendata/run_54280_HR_Feb102025/lltoee_check/Data_MBin92_NClus.root";
0010 std::string MC_ZF_file_directory = "/sphenix/tg/tg01/commissioning/INTT/work/cwshih/sPH_dNdeta/Run24AuAuMC/Sim_HIJING_MDC2_ana472_20250307/lltoee_check/MC_MBin92_NClus.root";
0011
0012 std::string Data_Run54912_file_directory = "/sphenix/tg/tg01/commissioning/INTT/work/cwshih/seflgendata/run_54912/lltoee_check/Data_MBin92_NClus.root";
0013 std::string MC_Field_file_directory = "/sphenix/tg/tg01/commissioning/INTT/work/cwshih/sPH_dNdeta/Run24AuAuMC/WrongStrangeIncrease_FieldOn_Sim_HIJING_strangeness_MDC2_ana467_20250226/lltoee_check/MC_MBin92_NClus.root";
0014
0015
0016
0017 TFile *Data_Run54280_file = TFile::Open(Data_Run54280_file_directory.c_str());
0018 TFile *MC_ZF_file = TFile::Open(MC_ZF_file_directory.c_str());
0019
0020 TFile *Data_Run54912_file = TFile::Open(Data_Run54912_file_directory.c_str());
0021 TFile *MC_Field_file = TFile::Open(MC_Field_file_directory.c_str());
0022
0023
0024
0025 TH2D * h2D_Data_Run54280 = (TH2D*)Data_Run54280_file->Get(call_hist_name.c_str());
0026 TH2D * h2D_MC_ZF = (TH2D*)MC_ZF_file->Get(call_hist_name.c_str());
0027
0028 TH2D * h2D_Data_Run54912 = (TH2D*)Data_Run54912_file->Get(call_hist_name.c_str());
0029 TH2D * h2D_MC_Field = (TH2D*)MC_Field_file->Get(call_hist_name.c_str());
0030
0031 std::cout<<h2D_Data_Run54280->GetMean(2)<<std::endl;
0032 std::cout<<h2D_MC_ZF->GetMean(2)<<std::endl;
0033 std::cout<<h2D_Data_Run54912->GetMean(2)<<std::endl;
0034 std::cout<<h2D_MC_Field->GetMean(2)<<std::endl;
0035
0036
0037
0038 TProfile * profile_Data_Run54280 = h2D_Data_Run54280->ProfileX("profile_Data_Run54280");
0039 TProfile * profile_MC_ZF = h2D_MC_ZF->ProfileX("profile_MC_ZF");
0040
0041 TProfile * profile_Data_Run54912 = h2D_Data_Run54912->ProfileX("profile_Data_Run54912");
0042 TProfile * profile_MC_Field = h2D_MC_Field->ProfileX("profile_MC_Field");
0043
0044
0045
0046 TH1D * profile_Data_Run54280_hist = (TH1D*)profile_Data_Run54280->Clone("profile_Data_Run54280_hist");
0047 TH1D * profile_MC_ZF_hist = (TH1D*)profile_MC_ZF->Clone("profile_MC_ZF_hist");
0048
0049 TH1D * profile_Data_Run54912_hist = (TH1D*)profile_Data_Run54912->Clone("profile_Data_Run54912_hist");
0050 TH1D * profile_MC_Field_hist = (TH1D*)profile_MC_Field->Clone("profile_MC_Field_hist");
0051
0052
0053
0054
0055 TH1D * h1D_MC_ZFtoField = new TH1D("h1D_MC_ZFtoField", "h1D_MC_ZFtoField;centrality;MC_ZF / MC_FieldON", h2D_Data_Run54280->GetNbinsX(), h2D_Data_Run54280->GetXaxis()->GetXmin(), h2D_Data_Run54280->GetXaxis()->GetXmax());
0056 h1D_MC_ZFtoField->Divide(profile_MC_ZF, profile_MC_Field);
0057
0058 TH1D * h1D_Data54280to54912 = new TH1D("h1D_Data54280to54912", "h1D_Data54280to54912;centrality;Data_Run54280 / Data_Run54912", h2D_Data_Run54280->GetNbinsX(), h2D_Data_Run54280->GetXaxis()->GetXmin(), h2D_Data_Run54280->GetXaxis()->GetXmax());
0059 h1D_Data54280to54912->Divide(profile_Data_Run54280, profile_Data_Run54912);
0060
0061 TH1D * h1D_DATAtoMC = new TH1D("h1D_DATAtoMC", "h1D_DATAtoMC;centrality;Data / MC", h2D_Data_Run54280->GetNbinsX(), h2D_Data_Run54280->GetXaxis()->GetXmin(), h2D_Data_Run54280->GetXaxis()->GetXmax());
0062 h1D_DATAtoMC->Divide(h1D_Data54280to54912, h1D_MC_ZFtoField);
0063
0064 TH1D * h1D_DatatoMC_FieldON = new TH1D("h1D_DatatoMC_FieldON", "h1D_DatatoMC_FieldON;centrality;(FieldON) Data_54912 / MC", h2D_Data_Run54280->GetNbinsX(), h2D_Data_Run54280->GetXaxis()->GetXmin(), h2D_Data_Run54280->GetXaxis()->GetXmax());
0065 h1D_DatatoMC_FieldON->Divide(profile_Data_Run54912, profile_MC_Field);
0066
0067 TH1D * h1D_DatatoMC_FieldOff = new TH1D("h1D_DatatoMC_FieldOff", "h1D_DatatoMC_FieldOff;centrality;(FieldOff) Data_54280 / MC", h2D_Data_Run54280->GetNbinsX(), h2D_Data_Run54280->GetXaxis()->GetXmin(), h2D_Data_Run54280->GetXaxis()->GetXmax());
0068 h1D_DatatoMC_FieldOff->Divide(profile_Data_Run54280, profile_MC_ZF);
0069
0070 TCanvas * c1 = new TCanvas("c_Data54280_Profile", "c_Data54280_Profile", 950, 800);
0071 c1->cd();
0072 TPad * pad1 = (TPad*) c1 -> cd();
0073 pad1 -> SetLeftMargin(0.18);
0074 pad1 -> SetRightMargin(0.08);
0075 profile_Data_Run54280_hist->SetTitle("Data_Run54280;centrality interval;(run 54280) NClus");
0076 profile_Data_Run54280_hist->SetMarkerStyle(20);
0077 profile_Data_Run54280_hist->SetMarkerColor(kBlack);
0078 profile_Data_Run54280_hist->SetLineColor(kBlack);
0079 profile_Data_Run54280_hist->SetMarkerSize(0.5);
0080
0081 profile_Data_Run54280_hist->Draw("E1");
0082 c1 -> Print("c_Data54280_Profile.pdf");
0083
0084
0085 TCanvas * c2 = new TCanvas("c_MC_ZF_Profile", "c_MC_ZF_Profile", 950, 800);
0086 c2->cd();
0087 TPad * pad2 = (TPad*) c2 -> cd();
0088 pad2 -> SetLeftMargin(0.18);
0089 pad2 -> SetRightMargin(0.08);
0090 profile_MC_ZF_hist->SetTitle("MC_ZF;centrality interval; (MC zero field) NClus");
0091 profile_MC_ZF_hist->SetMarkerStyle(20);
0092 profile_MC_ZF_hist->SetMarkerColor(kRed);
0093 profile_MC_ZF_hist->SetLineColor(kRed);
0094 profile_MC_ZF_hist->SetMarkerSize(0.5);
0095
0096 profile_MC_ZF_hist->Draw("E1");
0097 c2 -> Print("c_MC_ZF_Profile.pdf");
0098
0099 TCanvas * c3 = new TCanvas("c_Data54912_Profile", "c_Data54912_Profile", 950, 800);
0100 c3->cd();
0101 TPad * pad3 = (TPad*) c3 -> cd();
0102 pad3 -> SetLeftMargin(0.18);
0103 pad3 -> SetRightMargin(0.08);
0104 profile_Data_Run54912_hist->SetTitle("Data_Run54912;centrality interval;(run 54912) NClus");
0105 profile_Data_Run54912_hist->SetMarkerStyle(20);
0106 profile_Data_Run54912_hist->SetMarkerColor(kBlack);
0107 profile_Data_Run54912_hist->SetLineColor(kBlack);
0108 profile_Data_Run54912_hist->SetMarkerSize(0.5);
0109
0110 profile_Data_Run54912_hist->Draw("E1");
0111 c3 -> Print("c_Data54912_Profile.pdf");
0112
0113 TCanvas * c4 = new TCanvas("c_MC_Field_Profile", "c_MC_Field_Profile", 950, 800);
0114 c4->cd();
0115 TPad * pad4 = (TPad*) c4 -> cd();
0116 pad4 -> SetLeftMargin(0.18);
0117 pad4 -> SetRightMargin(0.08);
0118 profile_MC_Field_hist->SetTitle("MC_FieldON;centrality interval; (MC field on) NClus");
0119 profile_MC_Field_hist->SetMarkerStyle(20);
0120 profile_MC_Field_hist->SetMarkerColor(kRed);
0121 profile_MC_Field_hist->SetLineColor(kRed);
0122 profile_MC_Field_hist->SetMarkerSize(0.5);
0123
0124 profile_MC_Field_hist->Draw("E1");
0125 c4 -> Print("c_MC_Field_Profile.pdf");
0126
0127 TCanvas * c5 = new TCanvas("c_MC_ZFtoField", "c_MC_ZFtoField", 950, 800);
0128 c5->cd();
0129 TPad * pad5 = (TPad*) c5 -> cd();
0130 pad5 -> SetLeftMargin(0.18);
0131 pad5 -> SetRightMargin(0.08);
0132 h1D_MC_ZFtoField->SetMinimum(0.5);
0133 h1D_MC_ZFtoField->SetMaximum(1.5);
0134 h1D_MC_ZFtoField->SetTitle("MC_ZFtoField;centrality interval;NClus ratio (MC_ZF / MC_FieldON)");
0135 h1D_MC_ZFtoField->SetMarkerStyle(20);
0136 h1D_MC_ZFtoField->SetMarkerColor(kBlack);
0137 h1D_MC_ZFtoField->SetLineColor(kBlack);
0138 h1D_MC_ZFtoField->SetMarkerSize(0.5);
0139
0140 h1D_MC_ZFtoField->Draw("E1");
0141 c5 -> Print("c_MC_ZFtoField.pdf");
0142
0143 TCanvas * c6 = new TCanvas("c_Data54280to54912", "c_Data54280to54912", 950, 800);
0144 c6->cd();
0145 TPad * pad6 = (TPad*) c6 -> cd();
0146 pad6 -> SetLeftMargin(0.18);
0147 pad6 -> SetRightMargin(0.08);
0148 h1D_Data54280to54912->SetMinimum(1.5);
0149 h1D_Data54280to54912->SetMaximum(2.5);
0150 h1D_Data54280to54912->SetTitle("Data54280to54912;centrality interval;Data NClus ratio (Run54280 / Run54912)");
0151 h1D_Data54280to54912->SetMarkerStyle(20);
0152 h1D_Data54280to54912->SetMarkerColor(kBlack);
0153 h1D_Data54280to54912->SetLineColor(kBlack);
0154 h1D_Data54280to54912->SetMarkerSize(0.5);
0155
0156 h1D_Data54280to54912->Draw("E1");
0157 c6 -> Print("c_Data54280to54912.pdf");
0158
0159
0160 TCanvas * c7 = new TCanvas("c_DATAtoMC", "c_DATAtoMC", 950, 800);
0161 c7->cd();
0162 TPad * pad7 = (TPad*) c7 -> cd();
0163 pad7 -> SetLeftMargin(0.18);
0164 pad7 -> SetRightMargin(0.08);
0165 h1D_DATAtoMC->SetMinimum(1.5);
0166 h1D_DATAtoMC->SetMaximum(2.5);
0167 h1D_DATAtoMC->SetTitle("DATAtoMC;centrality interval;Ratio of ratio (Data / MC)");
0168 h1D_DATAtoMC->SetMarkerStyle(20);
0169 h1D_DATAtoMC->SetMarkerColor(kBlack);
0170 h1D_DATAtoMC->SetLineColor(kBlack);
0171 h1D_DATAtoMC->SetMarkerSize(0.5);
0172
0173 h1D_DATAtoMC->Draw("E1");
0174 c7 -> Print("c_DATAtoMC.pdf");
0175
0176
0177 TCanvas * c8 = new TCanvas("c_DatatoMC_FieldON", "c_DatatoMC_FieldON", 950, 800);
0178 c8->cd();
0179 TPad * pad8 = (TPad*) c8 -> cd();
0180 pad8 -> SetLeftMargin(0.18);
0181 pad8 -> SetRightMargin(0.08);
0182 h1D_DatatoMC_FieldON->SetMinimum(0.5);
0183 h1D_DatatoMC_FieldON->SetMaximum(1.5);
0184 h1D_DatatoMC_FieldON->SetTitle("DatatoMC_FieldON;centrality interval;NClus Ratio (Data_54912 / MC_FieldON)");
0185 h1D_DatatoMC_FieldON->SetMarkerStyle(20);
0186 h1D_DatatoMC_FieldON->SetMarkerColor(kBlack);
0187 h1D_DatatoMC_FieldON->SetLineColor(kBlack);
0188 h1D_DatatoMC_FieldON->SetMarkerSize(0.5);
0189
0190 h1D_DatatoMC_FieldON->Draw("E1");
0191 c8 -> Print("c_DatatoMC_FieldON.pdf");
0192
0193 TCanvas * c9 = new TCanvas("c_DatatoMC_FieldOff", "c_DatatoMC_FieldOff", 950, 800);
0194 c9->cd();
0195 TPad * pad9 = (TPad*) c9 -> cd();
0196 pad9 -> SetLeftMargin(0.18);
0197 pad9 -> SetRightMargin(0.08);
0198 h1D_DatatoMC_FieldOff->SetMinimum(0.5);
0199 h1D_DatatoMC_FieldOff->SetMaximum(1.5);
0200 h1D_DatatoMC_FieldOff->SetTitle("DatatoMC_FieldOff;centrality interval;NClus Ratio (Data_54280 / MC_ZF)");
0201 h1D_DatatoMC_FieldOff->SetMarkerStyle(20);
0202 h1D_DatatoMC_FieldOff->SetMarkerColor(kBlack);
0203 h1D_DatatoMC_FieldOff->SetLineColor(kBlack);
0204 h1D_DatatoMC_FieldOff->SetMarkerSize(0.5);
0205
0206 h1D_DatatoMC_FieldOff->Draw("E1");
0207 c9 -> Print("c_DatatoMC_FieldOff.pdf");
0208
0209
0210 return 888;
0211 }
