00001 // 00002 // ******************************************************************** 00003 // * License and Disclaimer * 00004 // * * 00005 // * The Geant4 software is copyright of the Copyright Holders of * 00006 // * the Geant4 Collaboration. It is provided under the terms and * 00007 // * conditions of the Geant4 Software License, included in the file * 00008 // * LICENSE and available at http://cern.ch/geant4/license . These * 00009 // * include a list of copyright holders. * 00010 // * * 00011 // * Neither the authors of this software system, nor their employing * 00012 // * institutes,nor the agencies providing financial support for this * 00013 // * work make any representation or warranty, express or implied, * 00014 // * regarding this software system or assume any liability for its * 00015 // * use. Please see the license in the file LICENSE and URL above * 00016 // * for the full disclaimer and the limitation of liability. * 00017 // * * 00018 // * This code implementation is the result of the scientific and * 00019 // * technical work of the GEANT4 collaboration. * 00020 // * By using, copying, modifying or distributing the software (or * 00021 // * any work based on the software) you agree to acknowledge its * 00022 // * use in resulting scientific publications, and indicate your * 00023 // * acceptance of all terms of the Geant4 Software license. * 00024 // ******************************************************************** 00025 // 00026 // $Id$ 00027 // 00028 // ------------------------------------------------------------------- 00029 // 00030 // GEANT4 Class file 00031 // 00032 // File name: G4DecayStrongResonances 00033 // 00034 // Modified: 00035 // 02.11.2010 V.Ivanchenko moved constructor and destructor to source 00036 // 07.27.2011 M.Kelsey -- Use new decay utility to process input list 00037 00038 #include "G4DecayStrongResonances.hh" 00039 00040 #include "G4DecayKineticTracks.hh" 00041 #include "G4HadTmpUtil.hh" 00042 #include <algorithm> 00043 00044 00045 G4DecayStrongResonances::G4DecayStrongResonances() {} 00046 00047 G4DecayStrongResonances::~G4DecayStrongResonances() {} 00048 00049 G4ReactionProductVector* 00050 G4DecayStrongResonances::Propagate(G4KineticTrackVector* theSecondaries, 00051 G4V3DNucleus* ) { 00052 G4DecayKineticTracks decay(theSecondaries); // Changes input list in situ 00053 00054 // translate to ReactionProducts 00055 G4ReactionProductVector * theResult; 00056 try { theResult = new G4ReactionProductVector; } 00057 catch(...) { 00058 throw G4HadronicException(__FILE__, __LINE__, "DecayStrongRes: out of memory "); 00059 } 00060 00061 G4ReactionProduct * it = NULL; 00062 00063 G4KineticTrackVector::iterator secIter = theSecondaries->begin(); 00064 for(; secIter != theSecondaries->end(); ++secIter) { 00065 G4KineticTrack* aSecondary = *secIter; 00066 if (!aSecondary) continue; // Skip null pointers 00067 00068 try { it = new G4ReactionProduct(); } 00069 catch(...) { 00070 throw G4HadronicException(__FILE__, __LINE__, "DecayStrongRes: out of memory "); 00071 } 00072 00073 it->SetDefinition(aSecondary->GetDefinition()); 00074 it->SetMass(aSecondary->GetDefinition()->GetPDGMass()); 00075 it->SetTotalEnergy(aSecondary->Get4Momentum().t()); 00076 it->SetMomentum(aSecondary->Get4Momentum().vect()); 00077 delete aSecondary; 00078 try { theResult->push_back(it); } 00079 catch(...){ 00080 throw G4HadronicException(__FILE__, __LINE__, "DecayStrongRes: push to result failed - out of mem."); 00081 } 00082 } 00083 delete theSecondaries; 00084 00085 return theResult; 00086 } 00087 00088