diff --git a/benchmarks/far_forward/scripts/hits_far_forward_protons.cxx b/benchmarks/far_forward/scripts/hits_far_forward_protons.cxx index 29ee91f1964ad726b740f3ffcf2bb6dc8081badf..5333ab7ab6697c733054fa5d6a5e233b4f2ad265 100644 --- a/benchmarks/far_forward/scripts/hits_far_forward_protons.cxx +++ b/benchmarks/far_forward/scripts/hits_far_forward_protons.cxx @@ -117,21 +117,21 @@ int hits_far_forward_protons(const char* fname = "sim_far_forward_protons.root") auto hEndcap_Nhits = df0.Histo1D({"hEndcap_Nhits", "; #theta [deg.]", 20, 0, 20 }, "N_EndcapHits"); //auto hVtxBarrel_Nhits = df0.Histo1D({"hVtxBarrel_Nhits", "; #theta [deg.]", 20, 0, 20 }, "N_VtxBarrelHits"); - auto hBarrel_Ntheta = df0.Histo1D({"hBarrel_Ntheta", "; #theta [deg.]", 20, 0, 180 }, "theta0"); - auto hEndcap_Ntheta = df0.Histo1D({"hEndcap_Ntheta", "; #theta [deg.]", 20, 0, 180 }, "theta0"); + auto hB0_theta = df0.Histo1D({"theta_b0", "; #theta [deg.]", 100, 0, 1.5 }, "theta0"); + //auto hEndcap_Ntheta = df0.Histo1D({"hEndcap_Ntheta", "; #theta [deg.]", 20, 0, 180 }, "theta0"); //auto hVtxBarrel_Ntheta = df0.Histo1D({"hVtxBarrel_Ntheta", "; #theta [deg.]", 20, 0, 180 }, "theta0"); auto c = new TCanvas(); auto hs = new THStack("n_hits","; #theta "); auto h1 = (TH1D*) hBarrel_N_vs_theta->Clone(); - auto h2 = (TH1D*) hBarrel_Ntheta->Clone(); + auto h2 = (TH1D*) hB0_theta->Clone(); h1->Divide(h2); hs->Add(h1); - h1 = (TH1D*) hEndcap_N_vs_theta->Clone(); - h2 = (TH1D*) hEndcap_Ntheta->Clone(); - h1->Divide(h2); + //h1 = (TH1D*) hEndcap_N_vs_theta->Clone(); + //h2 = (TH1D*) hEndcap_Ntheta->Clone(); + //h1->Divide(h2); h1->SetLineColor(2); - hs->Add(h1); + //hs->Add(h1); //h1 = (TH1D*) hVtxBarrel_vs_theta->Clone(); //h1->SetLineColor(4); //h1->SetFillStyle(3001); @@ -144,16 +144,16 @@ int hits_far_forward_protons(const char* fname = "sim_far_forward_protons.root") c = new TCanvas(); hs = new THStack("theta","; #theta "); - h1 = (TH1D*) hBarrel_N_vs_theta->Clone(); - h2 = (TH1D*) hBarrel_Ntheta->Clone(); + //h1 = (TH1D*) hBarrel_N_vs_theta->Clone(); + h2 = (TH1D*) hB0_theta->Clone(); //h1->Divide(h2); hs->Add(h2); - h1 = (TH1D*) hEndcap_N_vs_theta->Clone(); - h2 = (TH1D*) hEndcap_Ntheta->Clone(); + //h1 = (TH1D*) hEndcap_N_vs_theta->Clone(); + //h2 = (TH1D*) hEndcap_Ntheta->Clone(); //h1->Divide(h2); h1->SetLineColor(2); - h2->SetLineColor(2); - hs->Add(h2); + //h2->SetLineColor(2); + //hs->Add(h2); //h1 = (TH1D*) hVtxBarrel_vs_theta->Clone(); //h1->SetLineColor(4); //h1->SetFillStyle(3001); @@ -166,12 +166,12 @@ int hits_far_forward_protons(const char* fname = "sim_far_forward_protons.root") c = new TCanvas(); hs = new THStack("hits","; hits "); - h1 = (TH1D*) hBarrel_Nhits->Clone(); + h1 = (TH1D*) hB0_theta->Clone(); hs->Add(h1); - h1 = (TH1D*) hEndcap_Nhits->Clone(); + //h1 = (TH1D*) hEndcap_Nhits->Clone(); h1->SetLineColor(2); - h2->SetLineColor(2); - hs->Add(h2); + //h2->SetLineColor(2); + //hs->Add(h2); //h1 = (TH1D*) hVtxBarrel_Nhits->Clone(); //h1->SetLineColor(4); //h1->SetFillStyle(3001);