File indexing completed on 2025-08-10 08:12:57
0001 int DataMcComp()
0002 {
0003 TFile * file_in_data = TFile::Open("/sphenix/tg/tg01/commissioning/INTT/work/cwshih/seflgendata/run_54280_HR_Dec042024/completed/Run24NewCode_TrackHist/completed/Run24NewCode_dNdEta/completed/Data_PreparedNdEtaEach_ApplyAlphaCorr_AllSensor_VtxZ10_Mbin70_SecondRun_00054280_00000_dNdEta.root");
0004 TFile * file_in_mc = TFile::Open("/sphenix/user/ChengWei/sPH_dNdeta/Run24AuAuMC/Sim_Ntuple_HIJING_ana443_20241102/Run24NewCode_dNdEta/completed/MC_PreparedNdEtaEach_ApplyAlphaCorr_AllSensor_VtxZ10_Mbin70_SecondRun_00002_dNdEta.root");
0005
0006 TH1D * data_h1D_FitBkg_RecoTrackletEtaPerEvt = (TH1D*)file_in_data->Get("h1D_FitBkg_RecoTrackletEtaPerEvt");
0007 TH1D * data_h1D_RotatedBkg_RecoTrackletEtaPerEvt = (TH1D*)file_in_data->Get("h1D_RotatedBkg_RecoTrackletEtaPerEvt");
0008 TH1D * data_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC = (TH1D*)file_in_data->Get("h1D_FitBkg_RecoTrackletEtaPerEvtPostAC");
0009 TH1D * data_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC = (TH1D*)file_in_data->Get("h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC");
0010
0011 data_h1D_FitBkg_RecoTrackletEtaPerEvt -> SetMarkerColor(1);
0012 data_h1D_FitBkg_RecoTrackletEtaPerEvt -> SetLineColor(1);
0013
0014 data_h1D_RotatedBkg_RecoTrackletEtaPerEvt -> SetMarkerColor(1);
0015 data_h1D_RotatedBkg_RecoTrackletEtaPerEvt -> SetLineColor(1);
0016
0017 data_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC -> SetMarkerColor(1);
0018 data_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC -> SetLineColor(1);
0019
0020 data_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC -> SetMarkerColor(1);
0021 data_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC -> SetLineColor(1);
0022
0023
0024
0025 TH1D * MC_h1D_FitBkg_RecoTrackletEtaPerEvt = (TH1D*)file_in_mc->Get("h1D_FitBkg_RecoTrackletEtaPerEvt");
0026 TH1D * MC_h1D_RotatedBkg_RecoTrackletEtaPerEvt = (TH1D*)file_in_mc->Get("h1D_RotatedBkg_RecoTrackletEtaPerEvt");
0027 TH1D * MC_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC = (TH1D*)file_in_mc->Get("h1D_FitBkg_RecoTrackletEtaPerEvtPostAC");
0028 TH1D * MC_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC = (TH1D*)file_in_mc->Get("h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC");
0029
0030 MC_h1D_FitBkg_RecoTrackletEtaPerEvt -> SetMarkerColor(2);
0031 MC_h1D_FitBkg_RecoTrackletEtaPerEvt -> SetLineColor(2);
0032
0033 MC_h1D_RotatedBkg_RecoTrackletEtaPerEvt -> SetMarkerColor(2);
0034 MC_h1D_RotatedBkg_RecoTrackletEtaPerEvt -> SetLineColor(2);
0035
0036 MC_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC -> SetMarkerColor(2);
0037 MC_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC -> SetLineColor(2);
0038
0039 MC_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC -> SetMarkerColor(2);
0040 MC_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC -> SetLineColor(2);
0041
0042
0043
0044 TH1D * h1D_TruedNdEta = (TH1D*)file_in_mc->Get("h1D_TruedNdEta");
0045 h1D_TruedNdEta -> SetLineColor(3);
0046
0047
0048 TFile * file_out = new TFile("SecondRun_Mbin70.root", "RECREATE");
0049 TCanvas * c1 = new TCanvas("c1", "c1", 800, 600);
0050
0051 c1 -> cd();
0052 data_h1D_FitBkg_RecoTrackletEtaPerEvt -> Draw("ep");
0053 MC_h1D_FitBkg_RecoTrackletEtaPerEvt -> Draw("hist same");
0054 c1 -> Write("h1D_FitBkg_RecoTrackletEtaPerEvt");
0055 c1 -> Clear();
0056
0057 c1 -> cd();
0058 data_h1D_RotatedBkg_RecoTrackletEtaPerEvt -> Draw("ep");
0059 MC_h1D_RotatedBkg_RecoTrackletEtaPerEvt -> Draw("hist same");
0060 c1 -> Write("h1D_RotatedBkg_RecoTrackletEtaPerEvt");
0061 c1 -> Clear();
0062
0063 c1 -> cd();
0064 data_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC -> Draw("ep");
0065 MC_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC -> Draw("hist same");
0066 h1D_TruedNdEta -> Draw("hist same");
0067 c1 -> Write("h1D_FitBkg_RecoTrackletEtaPerEvtPostAC");
0068 c1 -> Clear();
0069
0070 c1 -> cd();
0071 data_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC -> Draw("ep");
0072 MC_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC -> Draw("hist same");
0073 h1D_TruedNdEta -> Draw("hist same");
0074 c1 -> Write("h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC");
0075 c1 -> Clear();
0076
0077 file_out -> Close();
0078
0079
0080
0081 return 0;
0082 }
