From bff894ceb12295ffa73cf0aeeca8fc64bb112410 Mon Sep 17 00:00:00 2001
From: itayhorin <itayhorin@gmail.com>
Date: Tue, 25 Oct 2022 13:42:26 +0300
Subject: [PATCH] changed p2 calc and verified it by maxwell boltzman
 distirbution

---
 src/physics.cpp               | 4 +++-
 src/tallies/tally_scoring.cpp | 7 +++----
 2 files changed, 6 insertions(+), 5 deletions(-)

diff --git a/src/physics.cpp b/src/physics.cpp
index c59909f0998..0ceebe5b150 100644
--- a/src/physics.cpp
+++ b/src/physics.cpp
@@ -744,7 +744,9 @@ void elastic_scatter(int i_nuclide, const Reaction& rx, double kT, Particle& p)
     v_t = sample_target_velocity(*nuc, p.E(), p.u(), v_n,
       p.neutron_xs(i_nuclide).elastic, kT, p.current_seed());
   }
-  p.v_t() = C_LIGHT*std::sqrt(2/(awr*p.getMass()))*v_t;
+  //std::cout << "v_n = "<< p.speed() <<std::endl;
+  p.v_t() = C_LIGHT*std::sqrt(2/p.getMass())*v_t;
+  std::cout << "v_t = "<< p.v_t().norm() <<std::endl;
   // Velocity of center-of-mass
   Direction v_cm = (v_n + awr * v_t) / (awr + 1.0);
 
diff --git a/src/tallies/tally_scoring.cpp b/src/tallies/tally_scoring.cpp
index 5b140928b81..59e3d501f7a 100644
--- a/src/tallies/tally_scoring.cpp
+++ b/src/tallies/tally_scoring.cpp
@@ -2367,7 +2367,6 @@ void score_tracklength_tally(Particle& p, double distance)
 
 void score_collision_tally(Particle& p)
 {
-  fmt::print("collision tally called\n");
   // Determine the collision estimate of the flux
   double flux = 0.0;
   if (p.type() == ParticleType::neutron || p.type() == ParticleType::photon) {
@@ -3088,9 +3087,9 @@ void getMu_COM(double x_det , double y_det , double z_det ,Particle p_col , doub
   std::cout << "p1:  "<< p1[0] <<" "<< p1[1]<<" " <<p1[2] <<std::endl;
   std::cout << "p2:  "<< (p_col.v_t() * m2 / C_LIGHT).x <<" " <<(p_col.v_t() * m2 / C_LIGHT).y << " "<<(p_col.v_t() * m2 / C_LIGHT).z <<std::endl;
   //std::cout << (p.v_t() * m2 / C_LIGHT).x << (p.v_t() * m2 / C_LIGHT).y << (p.v_t() * m2 / C_LIGHT).z <<std::endl;
-  double p2_x = (p_col.v_t() * m2 / C_LIGHT).x;
-  double p2_y = (p_col.v_t() * m2 / C_LIGHT).y;
-  double p2_z = (p_col.v_t() * m2 / C_LIGHT).z;
+  double p2_x = (p_col.v_t() * m2 /C_LIGHT).x;
+  double p2_y = (p_col.v_t() * m2 /C_LIGHT).y;
+  double p2_z = (p_col.v_t() * m2 /C_LIGHT).z;
   double p2[3]={p2_x, p2_y, p2_z}; //3 momentum of target in lab  
  
   // calculate