Skip to content
Snippets Groups Projects
Commit 3fd2e784 authored by Alex Jentsch's avatar Alex Jentsch
Browse files

Renaming a few histos and getting some familiarity with everything.

parent f29ca606
No related branches found
No related tags found
No related merge requests found
This commit is part of merge request !195. Comments created here will be created in the context of that merge request.
Hide whitespace changes
Inline Side-by-side
benchmarks/far_forward/scripts/hits_far_forward_protons.cxx
+ 17
− 17
Edit View file @ 3fd2e784
...@@ -117,21 +117,21 @@ int hits_far_forward_protons(const char* fname = "sim_far_forward_protons.root") ...@@ -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 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 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 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 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 hVtxBarrel_Ntheta = df0.Histo1D({"hVtxBarrel_Ntheta", "; #theta [deg.]", 20, 0, 180 }, "theta0");
auto c = new TCanvas(); auto c = new TCanvas();
auto hs = new THStack("n_hits","; #theta "); auto hs = new THStack("n_hits","; #theta ");
auto h1 = (TH1D*) hBarrel_N_vs_theta->Clone(); auto h1 = (TH1D*) hBarrel_N_vs_theta->Clone();
auto h2 = (TH1D*) hBarrel_Ntheta->Clone(); auto h2 = (TH1D*) hB0_theta->Clone();
h1->Divide(h2); h1->Divide(h2);
hs->Add(h1); hs->Add(h1);
h1 = (TH1D*) hEndcap_N_vs_theta->Clone(); //h1 = (TH1D*) hEndcap_N_vs_theta->Clone();
h2 = (TH1D*) hEndcap_Ntheta->Clone(); //h2 = (TH1D*) hEndcap_Ntheta->Clone();
h1->Divide(h2); //h1->Divide(h2);
h1->SetLineColor(2); h1->SetLineColor(2);
hs->Add(h1); //hs->Add(h1);
//h1 = (TH1D*) hVtxBarrel_vs_theta->Clone(); //h1 = (TH1D*) hVtxBarrel_vs_theta->Clone();
//h1->SetLineColor(4); //h1->SetLineColor(4);
//h1->SetFillStyle(3001); //h1->SetFillStyle(3001);
...@@ -144,16 +144,16 @@ int hits_far_forward_protons(const char* fname = "sim_far_forward_protons.root") ...@@ -144,16 +144,16 @@ int hits_far_forward_protons(const char* fname = "sim_far_forward_protons.root")
c = new TCanvas(); c = new TCanvas();
hs = new THStack("theta","; #theta "); hs = new THStack("theta","; #theta ");
h1 = (TH1D*) hBarrel_N_vs_theta->Clone(); //h1 = (TH1D*) hBarrel_N_vs_theta->Clone();
h2 = (TH1D*) hBarrel_Ntheta->Clone(); h2 = (TH1D*) hB0_theta->Clone();
//h1->Divide(h2); //h1->Divide(h2);
hs->Add(h2); hs->Add(h2);
h1 = (TH1D*) hEndcap_N_vs_theta->Clone(); //h1 = (TH1D*) hEndcap_N_vs_theta->Clone();
h2 = (TH1D*) hEndcap_Ntheta->Clone(); //h2 = (TH1D*) hEndcap_Ntheta->Clone();
//h1->Divide(h2); //h1->Divide(h2);
h1->SetLineColor(2); h1->SetLineColor(2);
h2->SetLineColor(2); //h2->SetLineColor(2);
hs->Add(h2); //hs->Add(h2);
//h1 = (TH1D*) hVtxBarrel_vs_theta->Clone(); //h1 = (TH1D*) hVtxBarrel_vs_theta->Clone();
//h1->SetLineColor(4); //h1->SetLineColor(4);
//h1->SetFillStyle(3001); //h1->SetFillStyle(3001);
...@@ -166,12 +166,12 @@ int hits_far_forward_protons(const char* fname = "sim_far_forward_protons.root") ...@@ -166,12 +166,12 @@ int hits_far_forward_protons(const char* fname = "sim_far_forward_protons.root")
c = new TCanvas(); c = new TCanvas();
hs = new THStack("hits","; hits "); hs = new THStack("hits","; hits ");
h1 = (TH1D*) hBarrel_Nhits->Clone(); h1 = (TH1D*) hB0_theta->Clone();
hs->Add(h1); hs->Add(h1);
h1 = (TH1D*) hEndcap_Nhits->Clone(); //h1 = (TH1D*) hEndcap_Nhits->Clone();
h1->SetLineColor(2); h1->SetLineColor(2);
h2->SetLineColor(2); //h2->SetLineColor(2);
hs->Add(h2); //hs->Add(h2);
//h1 = (TH1D*) hVtxBarrel_Nhits->Clone(); //h1 = (TH1D*) hVtxBarrel_Nhits->Clone();
//h1->SetLineColor(4); //h1->SetLineColor(4);
//h1->SetFillStyle(3001); //h1->SetFillStyle(3001);
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment