macro/Quick_check/McMcComp.C

0001 int McMcComp()
0002 {
0003     TFile * file_in_data = TFile::Open("/sphenix/user/ChengWei/sPH_dNdeta/Run24AuAuMC/Sim_Ntuple_HIJING_ana443_20241102/Run24NewCode_dNdEta/completed/MC_PreparedNdEtaEach_AllSensor_VtxZ10_Mbin100_SecondRun_00001_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_Mbin100_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     TH1D * data_h1D_TruedNdEta = (TH1D*)file_in_data->Get("h1D_TruedNdEta");
0024     data_h1D_TruedNdEta -> SetFillColorAlpha(1,0);
0025     data_h1D_TruedNdEta -> SetLineColor(4);
0026
0027
0028
0029     TH1D * MC_h1D_FitBkg_RecoTrackletEtaPerEvt = (TH1D*)file_in_mc->Get("h1D_FitBkg_RecoTrackletEtaPerEvt");
0030     TH1D * MC_h1D_RotatedBkg_RecoTrackletEtaPerEvt = (TH1D*)file_in_mc->Get("h1D_RotatedBkg_RecoTrackletEtaPerEvt");
0031     TH1D * MC_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC = (TH1D*)file_in_mc->Get("h1D_FitBkg_RecoTrackletEtaPerEvtPostAC");
0032     TH1D * MC_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC = (TH1D*)file_in_mc->Get("h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC");
0033
0034     MC_h1D_FitBkg_RecoTrackletEtaPerEvt -> SetMarkerColor(2);
0035     MC_h1D_FitBkg_RecoTrackletEtaPerEvt -> SetLineColor(2);
0036
0037     MC_h1D_RotatedBkg_RecoTrackletEtaPerEvt -> SetMarkerColor(2);
0038     MC_h1D_RotatedBkg_RecoTrackletEtaPerEvt -> SetLineColor(2);
0039
0040     MC_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC -> SetMarkerColor(2);
0041     MC_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC -> SetLineColor(2);
0042
0043     MC_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC -> SetMarkerColor(2);
0044     MC_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC -> SetLineColor(2);
0045
0046
0047
0048     TH1D * MC_h1D_TruedNdEta = (TH1D*)file_in_mc->Get("h1D_TruedNdEta");
0049     MC_h1D_TruedNdEta -> SetFillColorAlpha(2,0);
0050     MC_h1D_TruedNdEta -> SetLineColor(3);
0051
0052
0053     TFile * file_out = new TFile("SecondRun_Mbin70.root", "RECREATE");
0054     TCanvas * c1 = new TCanvas("c1", "c1", 800, 600);
0055
0056     c1 -> cd();
0057     data_h1D_FitBkg_RecoTrackletEtaPerEvt -> Draw("ep");
0058     MC_h1D_FitBkg_RecoTrackletEtaPerEvt -> Draw("hist same");
0059     c1 -> Write("h1D_FitBkg_RecoTrackletEtaPerEvt");
0060     c1 -> Clear();
0061
0062     c1 -> cd();
0063     data_h1D_RotatedBkg_RecoTrackletEtaPerEvt -> Draw("ep");
0064     MC_h1D_RotatedBkg_RecoTrackletEtaPerEvt -> Draw("hist same");
0065     c1 -> Write("h1D_RotatedBkg_RecoTrackletEtaPerEvt");
0066     c1 -> Clear();
0067
0068     c1 -> cd();
0069     data_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC -> Draw("ep");
0070     MC_h1D_FitBkg_RecoTrackletEtaPerEvtPostAC -> Draw("hist same");
0071     MC_h1D_TruedNdEta -> Draw("hist same");
0072     data_h1D_TruedNdEta -> Draw("hist same");
0073     c1 -> Write("h1D_FitBkg_RecoTrackletEtaPerEvtPostAC");
0074     c1 -> Clear();
0075
0076     c1 -> cd();
0077     data_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC -> Draw("ep");
0078     MC_h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC -> Draw("hist same");
0079     MC_h1D_TruedNdEta -> Draw("hist same");
0080     data_h1D_TruedNdEta -> Draw("hist same");
0081     c1 -> Write("h1D_RotatedBkg_RecoTrackletEtaPerEvtPostAC");
0082     c1 -> Clear();
0083
0084     c1 -> cd();
0085     MC_h1D_TruedNdEta -> Draw("hist");
0086     data_h1D_TruedNdEta -> Draw("ep same");
0087     c1 -> Write("h1D_TruedNdEta");
0088     c1 -> Clear();
0089
0090     file_out -> Close();
0091
0092
0093
0094     return 0;
0095 }