Skip to content
Snippets Groups Projects

Draft: Resolve "Add energy scan for Barrel Ecal"

Closed Maria Zurek requested to merge 45-add-energy-scan-for-barrel-ecal into master
1 file
+ 3
4
Compare changes
  • Side-by-side
  • Inline
@@ -87,14 +87,13 @@ void energy_scan_analysis(const char* input_fname = "sim_output/rec_emcal_barrel
@@ -87,14 +87,13 @@ void energy_scan_analysis(const char* input_fname = "sim_output/rec_emcal_barrel
auto hEthr = d1.Histo1D({"hEthr", "Thrown Energy; Thrown Energy [GeV]; Events", 100, 0.0, 25.0},"Ethr");
auto hEthr = d1.Histo1D({"hEthr", "Thrown Energy; Thrown Energy [GeV]; Events", 100, 0.0, 25.0},"Ethr");
//auto hNhits =d1.Histo1D({"hNhits", "Number of hits per events; Number of hits; Events", 100, 0.0, 2000.0},"nhits");
//auto hNhits =d1.Histo1D({"hNhits", "Number of hits per events; Number of hits; Events", 100, 0.0, 2000.0},"nhits");
auto hErecImg = d1.Histo1D({"hErecImg", "Reconstructed energy deposit; Energy Deposit [GeV]; Events", 500, 0.0, 0.5},"ErecImg");
auto hErecImg = d1.Histo1D({"hErecImg", "Reconstructed energy deposit; Energy Deposit [GeV]; Events", 500, 0.0, 0.5},"ErecImg");
auto hEClusterImg = d1.Histo1D({"hEClusterImg", "Cluster Energy; Cluster Energy [GeV]; Events", 500, 0.0, 25.0},"EClusterImg");
auto hEClusterImg = d1.Histo1D({"hEClusterImg", "Cluster Energy; Cluster Energy [GeV]; Events", 200, 0.0, 25.0},"EClusterImg");
auto hNClusterImg = d1.Histo1D({"hNClusterImg", "Number of Clusters; # of Clusters; Events", 100, 0.0, 100.0}, "NClusterImg");
auto hNClusterImg = d1.Histo1D({"hNClusterImg", "Number of Clusters; # of Clusters; Events", 100, 0.0, 100.0}, "NClusterImg");
auto hfsamImg = d1.Histo1D({"hfsamImg", "Sampling Fraction; Sampling Fraction; Events", 800, 0.5, 1.5},"fsamImg");
auto hfsamImg = d1.Histo1D({"hfsamImg", "Sampling Fraction; Sampling Fraction; Events", 100, 0.0, 1.5},"fsamImg");
auto hErecScFi = d1.Histo1D({"hErecScFi", "Reconstructed energy deposit; Energy Deposit [GeV]; Events", 500, 0.0, 10.0},"ErecScFi");
auto hErecScFi = d1.Histo1D({"hErecScFi", "Reconstructed energy deposit; Energy Deposit [GeV]; Events", 500, 0.0, 10.0},"ErecScFi");
auto hEClusterScFi = d1.Histo1D({"hEClusterScFi", "Cluster Energy; Cluster Energy [GeV]; Events", 500, 0.0, 25.0},"EClusterScFi");
auto hEClusterScFi = d1.Histo1D({"hEClusterScFi", "Cluster Energy; Cluster Energy [GeV]; Events", 500, 0.0, 25.0},"EClusterScFi");
auto hNClusterScFi = d1.Histo1D({"hNClusterScFi", "Number of Clusters; # of Clusters; Events", 100, 0.0, 100.0}, "NClusterScFi");
auto hNClusterScFi = d1.Histo1D({"hNClusterScFi", "Number of Clusters; # of Clusters; Events", 100, 0.0, 100.0}, "NClusterScFi");
auto hfsamScFi = d1.Histo1D({"hfsamScFi", "Sampling Fraction; Sampling Fraction; Events", 800, 0.8, 1.2},"fsamScFi");
auto hfsamScFi = d1.Histo1D({"hfsamScFi", "Sampling Fraction; Sampling Fraction; Events", 100, 0.8, 1.2},"fsamScFi");
// Event Counts
// Event Counts
auto nevents_thrown = d1.Count();
auto nevents_thrown = d1.Count();
Loading