From 4ba47c2d9c94db812089ac1ae7df1c828831747c Mon Sep 17 00:00:00 2001 From: Patrick Green Date: Thu, 4 Dec 2025 19:31:51 -0600 Subject: [PATCH] fix g4rw daughters --- .../Calculators/Geant4/Geant4WeightCalc.cxx | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) diff --git a/sbncode/SBNEventWeight/Calculators/Geant4/Geant4WeightCalc.cxx b/sbncode/SBNEventWeight/Calculators/Geant4/Geant4WeightCalc.cxx index ec9d8cc6e..8677f7c55 100644 --- a/sbncode/SBNEventWeight/Calculators/Geant4/Geant4WeightCalc.cxx +++ b/sbncode/SBNEventWeight/Calculators/Geant4/Geant4WeightCalc.cxx @@ -389,10 +389,13 @@ std::vector Geant4WeightCalc::GetWeight(art::Event& e, size_t itruth ) { theTraj.AddStep( theStep ); - for( size_t k = 0; k < daughter_PDGs.size(); ++k ){ - theTraj.AddChild( new G4ReweightTraj(daughter_IDs[k], daughter_PDGs[k], mcpID, event_num, std::make_pair(0,0) ) ); - } } // end loop over nSteps (istep) + + // add daughters + for( size_t k = 0; k < daughter_PDGs.size(); ++k ){ + theTraj.AddChild( new G4ReweightTraj(daughter_IDs[k], daughter_PDGs[k], mcpID, event_num, std::make_pair(0,0) ) ); + } + p_track_length = theTraj.GetTotalLength(); p_init_momentum = sqrt( theTraj.GetEnergy()*theTraj.GetEnergy() - mass*mass );