Draft: Resolve "Add Roman Pots Reconstruction Matrix"

Closes #70

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);