From ed6a8a42341f2107775d44f200f5ef9bacf7ec73 Mon Sep 17 00:00:00 2001
From: Marshall Scott <mbscott@anl.gov>
Date: Mon, 3 May 2021 13:05:19 -0400
Subject: [PATCH] see sampling fraction mean 3
---
benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx b/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx
index de2a14ca..eef7ac6b 100644
--- a/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx
+++ b/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx
@@ -56,7 +56,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
ROOT::RDataFrame d0("events", input_fname);
// Sampling Fraction
- double samp_frac = 1.0;
+ double samp_frac = 0.0136;
// Thrown Energy [GeV]
auto Ethr = [](std::vector<dd4pod::Geant4ParticleData> const& input) {
@@ -132,6 +132,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
auto hpid = d1.Histo1D({"hpid", "PID; PID; Count", 100, -220, 220}, "pid");
// Set sampling Fraction, ideally this will be taken from a json file
+ cout << "The sampling fraction mean is " << samp_frac << endl;
samp_frac = hfsam -> GetMean();
cout << "The sampling fraction mean is " << samp_frac << endl;
--
GitLab