Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // G4StackChecker 27 // 28 // Author: Makoto Asai, 2003 29 // -------------------------------------------------------------------- 30 #include "G4StackChecker.hh" 31 #include "G4Track.hh" 32 #include "G4EventManager.hh" 33 #include "G4Event.hh" 34 #include "G4ios.hh" 35 36 G4ClassificationOfNewTrack 37 G4StackChecker::ClassifyNewTrack(const G4Track* track) 38 { 39 G4ClassificationOfNewTrack result = fUrgent; 40 G4double e = track->GetKineticEnergy(); 41 if ( (!(e < 0.0) && !(e > 0.0) && !(e == 0.0)) 42 || track->GetMomentumDirection() == nullDirection) 43 { 44 result = fKill; 45 G4String nam = track->GetDefinition()->GetParticleName(); 46 G4cout << "### G4StackChecker: event# " 47 << (G4EventManager::GetEventManager())->GetConstCurrentEvent()->GetEventID() 48 << " unacceptable " << nam << " is killed in the stack" << G4endl; 49 G4cout << "### " << nam << " have been produced by the process " 50 << track->GetCreatorProcess()->GetProcessName() 51 << " trackID= " << track->GetTrackID() 52 << " parentID= " << track->GetParentID() 53 << G4endl; 54 G4cout << "### E= " << track->GetKineticEnergy() 55 << " position= " << track->GetPosition() 56 << " direction= " << track->GetMomentumDirection() 57 << " time= " << track->GetGlobalTime() 58 << G4endl; 59 } 60 return result; 61 } 62