From e4f377ebb64bea0dff259e6cd242864694e7d8d5 Mon Sep 17 00:00:00 2001
From: John Tramm <john.tramm@gmail.com>
Date: Wed, 27 Nov 2019 18:18:54 +0000
Subject: [PATCH] added particle trace print

---
 examples/xml/pincell/settings.xml | 3 ++-
 src/simulation.cpp                | 5 +++++
 2 files changed, 7 insertions(+), 1 deletion(-)

diff --git a/examples/xml/pincell/settings.xml b/examples/xml/pincell/settings.xml
index 48da5cc79..2908c2aab 100644
--- a/examples/xml/pincell/settings.xml
+++ b/examples/xml/pincell/settings.xml
@@ -5,7 +5,8 @@
   <run_mode>eigenvalue</run_mode>
   <batches>2</batches>
   <inactive>1</inactive>
-  <particles>10</particles>
+  <particles>1</particles>
+  <seed>1337</seed>
 
   <!-- The starting source is a uniform distribution over the entire pin
        cell. Note that since this is effectively a 2D model, the z coordinates
diff --git a/src/simulation.cpp b/src/simulation.cpp
index fd7c05df5..f8ff9023c 100644
--- a/src/simulation.cpp
+++ b/src/simulation.cpp
@@ -584,6 +584,11 @@ void transport()
 			std::cout << "Surface Crossings = " << surface_crossing_queue_length << std::endl;
 			std::cout << "Collisions = " << collision_queue_length << std::endl;
 			*/
+			Particle * p = particles;
+		   std::cout << "E = " << p->E_ << " and Position {" <<
+			   p->r().x << ", " <<
+			   p->r().y << ", " <<
+			   p->r().z << "}" << std::endl;
 			int max = std::max({calculate_fuel_xs_queue_length, calculate_nonfuel_xs_queue_length, advance_particle_queue_length, surface_crossing_queue_length, collision_queue_length});
 			//check_energies();
 			if (max == 0) {