solving rad decay bug

icecube · Sep 5, 2024 · 5a59777 · 5a59777
1 parent 64ea9b3
commit 5a59777
Show file tree

Hide file tree

Showing 3 changed files with 4 additions and 3 deletions.
diff --git a/simulations/CMakeLists.txt b/simulations/CMakeLists.txt
@@ -1,4 +1,4 @@
-add_subdirectory(effective_area)
-#add_subdirectory(radioactive_decays)
+#add_subdirectory(effective_area)
+add_subdirectory(radioactive_decays)
 #add_subdirectory(supernova)
 #add_subdirectory(efficiency_calibration)
diff --git a/simulations/radioactive_decays/src/OMSimDecaysAnalysis.cc b/simulations/radioactive_decays/src/OMSimDecaysAnalysis.cc
@@ -74,6 +74,7 @@ void OMSimDecaysAnalysis::writeMultiplicity(G4double p_timeWindow)
  */
 void OMSimDecaysAnalysis::writeThreadDecayInformation()
 {
+ if (!m_threadDecayStats) return;
  log_trace("Writing decay information of {} decays", m_threadDecayStats->eventId.size());
 
  G4String outputSufix = OMSimCommandArgsTable::getInstance().get<std::string>("output_file");

diff --git a/simulations/radioactive_decays/src/OMSimG4RadioactiveDecay.cc b/simulations/radioactive_decays/src/OMSimG4RadioactiveDecay.cc
@@ -1076,7 +1076,7 @@ void G4RadioactiveDecay::DecayAnalog(const G4Track& theTrack,
  } 
  }
 
-  
+ log_info("DecayAnalog");
  OMSimDecaysAnalysis::getInstance().appendDecay(
  theParticleDef->GetParticleName(),
  finalGlobalTime / s, theTrack.GetPosition());