Resolve "Check ZDC particle gun acceptance"

Closes #64 (closed)

benchmarks/zdc/scripts/analysis_zdc_particles.cxx

@@ -69,6 +69,8 @@ void analysis_zdc_particles(
     .Define("nhits_Hcal", nhits, {"ZDCHcalHits"})
     .Define("Esim_Hcal", Esim, {"ZDCHcalHits"})
     .Define("fsam_Hcal", fsam, {"Esim_Hcal", "Ethr"})
+    .Define("hit_x_position", hit_x_position, {"ZDCHcalHits"})
+	  .Define("hit_y_position", hit_y_position, {"ZDCHcalHits"})
   ;
 
   // Define Histograms
@@ -96,6 +98,9 @@ void analysis_zdc_particles(
       {"hfsam_Hcal", "Hcal Sampling Fraction; Hcal Sampling Fraction; Events", 100, 0.0, 0.1},
       "fsam_Hcal");
 
+  auto hXY_HitMap_Hcal = d1.Histo2D({"hXY_HitMap_Hcal", "hit position Y vs. X histogram; hit position X [mm]; hit position Y [mm]", 100,-30,30, 100, -30, 30}, "hit_x_position", "hit_y_position");
+
+
   // Event Counts
   auto nevents_thrown = d1.Count();
   std::cout << "Number of Thrown Events: " << (*nevents_thrown) << "\n";
@@ -187,4 +192,10 @@ void analysis_zdc_particles(
     c4->SaveAs(TString(results_path + "/zdc_fsam_Hcal.png"));
     c4->SaveAs(TString(results_path + "/zdc_fsam_Hcal.pdf"));
   }
+  {
+    TCanvas* c5 = new TCanvas("c5", "c5", 600, 500);
+    hXY_HitMap_Hcal->Draw("COLZ");
+    c5->SaveAs(TString(results_path + "/zdc_xy_map_Hcal.png"));
+    c5->SaveAs(TString(results_path + "/zdc_xy_map_Hcal.pdf"));
+  }
 }