From 27388635d095fc266829e39acedc3c0246bb03a0 Mon Sep 17 00:00:00 2001
From: Whitney Armstrong <warmstrong@anl.gov>
Date: Mon, 19 Oct 2020 20:46:38 -0500
Subject: [PATCH] Removed dummies and moved scripts around.
- Each benchmark or set ofcommon benchmarks will have their own
directory.
modified: emcal_electrons.sh
new file: scripts/emcal_electrons.cxx
new file: scripts/emcal_electrons_reader.cxx
modified: emcal_electrons.sh
modified: scripts/emcal_electrons.cxx
---
ecal/emcal_electrons.sh | 2 +-
ecal/scripts/emcal_electrons.cxx | 4 +++-
2 files changed, 4 insertions(+), 2 deletions(-)
diff --git a/ecal/emcal_electrons.sh b/ecal/emcal_electrons.sh
index eefa87ad..3173c35d 100644
--- a/ecal/emcal_electrons.sh
+++ b/ecal/emcal_electrons.sh
@@ -5,7 +5,7 @@ if [[ ! -n "${JUGGLER_DETECTOR}" ]] ; then
fi
if [[ ! -n "${JUGGLER_N_EVENTS}" ]] ; then
- export JUGGLER_N_EVENTS=1000
+ export JUGGLER_N_EVENTS=100
fi
export JUGGLER_FILE_NAME_TAG="emcal_uniform_electrons"
diff --git a/ecal/scripts/emcal_electrons.cxx b/ecal/scripts/emcal_electrons.cxx
index ca45a209..02ac3fe0 100644
--- a/ecal/scripts/emcal_electrons.cxx
+++ b/ecal/scripts/emcal_electrons.cxx
@@ -19,7 +19,7 @@
using namespace HepMC3;
-void emcal_electrons(int n_events = 1e6, const char* out_fname = "./data/emcal_electron_0GeVto30GeV_100kEvt.hepmc")
+void emcal_electrons(int n_events = 1e3, const char* out_fname = "./data/emcal_electron_0GeVto30GeV_100kEvt.hepmc")
{
WriterAscii hepmc_output(out_fname);
int events_parsed = 0;
@@ -48,6 +48,8 @@ void emcal_electrons(int n_events = 1e6, const char* out_fname = "./data/emcal_e
// sphere of given radius, in this case momentum.
r1->Sphere(px, py, pz, p);
+ //std::cout << std::sqrt(px*px + py*py + pz*pz) - p << " is zero? \n";
+
// type 1 is final state
// pdgid 11 - electron 0.510 MeV/c^2
GenParticlePtr p3 = std::make_shared<GenParticle>(
--
GitLab