diff --git a/examples/xml/pincell/settings.xml b/examples/xml/pincell/settings.xml
index 2908c2aab..e979a331f 100644
--- a/examples/xml/pincell/settings.xml
+++ b/examples/xml/pincell/settings.xml
@@ -5,7 +5,7 @@
   <run_mode>eigenvalue</run_mode>
   <batches>2</batches>
   <inactive>1</inactive>
-  <particles>1</particles>
+  <particles>2</particles>
   <seed>1337</seed>
 
   <!-- The starting source is a uniform distribution over the entire pin
diff --git a/src/simulation.cpp b/src/simulation.cpp
index f8ff9023c..6e2d63c97 100644
--- a/src/simulation.cpp
+++ b/src/simulation.cpp
@@ -220,6 +220,7 @@ void process_calculate_xs_events(int * queue, int n)
           // If particle energy is greater than the highest energy for the
           // S(a,b) table, then don't use the S(a,b) table
           //if (p->E_ > data::thermal_scatt[i_sab]->threshold()) i_sab = C_NONE;
+			std::cout << "S(alpha, beta) max energy = " << data::thermal_scatt[i_sab]->energy_max_ <<std::endl;
           if (p->E_ > data::thermal_scatt[i_sab]->energy_max_) i_sab = C_NONE;
         }
       }
@@ -584,7 +585,7 @@ void transport()
 			std::cout << "Surface Crossings = " << surface_crossing_queue_length << std::endl;
 			std::cout << "Collisions = " << collision_queue_length << std::endl;
 			*/
-			Particle * p = particles;
+			Particle * p = particles +  1;
 		   std::cout << "E = " << p->E_ << " and Position {" <<
 			   p->r().x << ", " <<
 			   p->r().y << ", " <<