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