Add helper function to write out details

Author: EinarElen

include/SimCore/KaonPhysics.h

@@ -9,6 +9,8 @@
 #include <G4ParticleDefinition.hh>
 #include <G4VDecayChannel.hh>
 #include <G4VPhysicsConstructor.hh>
+#include <iomanip>
+#include <iostream>
 #include <numeric>
 #include <vector>
 
@@ -127,6 +129,8 @@ class KaonPhysics : public G4VPhysicsConstructor {
 
   void ConstructParticle() override;
 
+  void DumpDecayDetails(const G4ParticleDefinition* kaon) const;
+
   /**
    * Construct processes
    *

src/SimCore/KaonPhysics.cxx

@@ -79,4 +79,26 @@ void KaonPhysics::ConstructParticle() {
   setDecayProperties(kaonShort, k0s_branching_ratios, k0s_lifetime_factor);
 }
 
+void KaonPhysics::DumpDecayDetails(const G4ParticleDefinition* kaon) const {
+  std::cout << "Decay table details for " << kaon->GetParticleName()
+            << std::endl
+            << std::scientific << std::setprecision(15);
+  std::cout << "PDG Lifetime " << kaon->GetPDGLifeTime() << std::endl;
+  const auto table{kaon->GetDecayTable()};
+  const int entries{table->entries()};
+  for (auto i{0}; i < entries; ++i) {
+    const auto channel{(*table)[i]};
+    std::cout << "Channel " << i << " Kinematics type "
+              << channel->GetKinematicsName() << " with BR " << channel->GetBR()
+              << std::endl;
+    std::cout << kaon->GetParticleName() << " -> ";
+    const auto daughters{channel->GetNumberOfDaughters()};
+    for (auto j{0}; j < daughters - 1; ++j) {
+      std::cout << channel->GetDaughter(j)->GetParticleName() << " + ";
+    }
+    // Special formatting for last one :)
+    std::cout << channel->GetDaughter(daughters - 1)->GetParticleName()
+              << std::endl;
+  }
+}
 }  // namespace simcore