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 /// \file runAndEvent/RE01/src/RE01TrackInformation.cc 27 /// \brief Implementation of the RE01TrackInformation class 28 // 29 // 30 // 31 32 #include "RE01TrackInformation.hh" 33 34 #include "G4SystemOfUnits.hh" 35 #include "G4ios.hh" 36 37 G4ThreadLocal G4Allocator<RE01TrackInformation>* aTrackInformationAllocator = 0; 38 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 40 RE01TrackInformation::RE01TrackInformation() : G4VUserTrackInformation() 41 { 42 fOriginalTrackID = 0; 43 fParticleDefinition = 0; 44 fOriginalPosition = G4ThreeVector(0., 0., 0.); 45 fOriginalMomentum = G4ThreeVector(0., 0., 0.); 46 fOriginalEnergy = 0.; 47 fOriginalTime = 0.; 48 fTrackingStatus = 1; 49 fSourceTrackID = -1; 50 fSourceDefinition = 0; 51 fSourcePosition = G4ThreeVector(0., 0., 0.); 52 fSourceMomentum = G4ThreeVector(0., 0., 0.); 53 fSourceEnergy = 0.; 54 fSourceTime = 0.; 55 fSuspendedStepID = -1; 56 } 57 58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 59 RE01TrackInformation::RE01TrackInformation(const G4Track* aTrack) : G4VUserTrackInformation() 60 { 61 fOriginalTrackID = aTrack->GetTrackID(); 62 fParticleDefinition = aTrack->GetDefinition(); 63 fOriginalPosition = aTrack->GetPosition(); 64 fOriginalMomentum = aTrack->GetMomentum(); 65 fOriginalEnergy = aTrack->GetTotalEnergy(); 66 fOriginalTime = aTrack->GetGlobalTime(); 67 fTrackingStatus = 1; 68 fSourceTrackID = -1; 69 fSourceDefinition = 0; 70 fSourcePosition = G4ThreeVector(0., 0., 0.); 71 fSourceMomentum = G4ThreeVector(0., 0., 0.); 72 fSourceEnergy = 0.; 73 fSourceTime = 0.; 74 fSuspendedStepID = -1; 75 } 76 77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 78 RE01TrackInformation ::RE01TrackInformation(const RE01TrackInformation* aTrackInfo) 79 : G4VUserTrackInformation() 80 { 81 fOriginalTrackID = aTrackInfo->fOriginalTrackID; 82 fParticleDefinition = aTrackInfo->fParticleDefinition; 83 fOriginalPosition = aTrackInfo->fOriginalPosition; 84 fOriginalMomentum = aTrackInfo->fOriginalMomentum; 85 fOriginalEnergy = aTrackInfo->fOriginalEnergy; 86 fOriginalTime = aTrackInfo->fOriginalTime; 87 fTrackingStatus = aTrackInfo->fTrackingStatus; 88 fSourceTrackID = aTrackInfo->fSourceTrackID; 89 fSourceDefinition = aTrackInfo->fSourceDefinition; 90 fSourcePosition = aTrackInfo->fSourcePosition; 91 fSourceMomentum = aTrackInfo->fSourceMomentum; 92 fSourceEnergy = aTrackInfo->fSourceEnergy; 93 fSourceTime = aTrackInfo->fSourceTime; 94 fSuspendedStepID = -1; 95 } 96 97 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 98 RE01TrackInformation::~RE01TrackInformation() 99 { 100 ; 101 } 102 103 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 104 RE01TrackInformation& RE01TrackInformation ::operator=(const RE01TrackInformation& aTrackInfo) 105 { 106 fOriginalTrackID = aTrackInfo.fOriginalTrackID; 107 fParticleDefinition = aTrackInfo.fParticleDefinition; 108 fOriginalPosition = aTrackInfo.fOriginalPosition; 109 fOriginalMomentum = aTrackInfo.fOriginalMomentum; 110 fOriginalEnergy = aTrackInfo.fOriginalEnergy; 111 fOriginalTime = aTrackInfo.fOriginalTime; 112 fTrackingStatus = aTrackInfo.fTrackingStatus; 113 fSourceTrackID = aTrackInfo.fSourceTrackID; 114 fSourceDefinition = aTrackInfo.fSourceDefinition; 115 fSourcePosition = aTrackInfo.fSourcePosition; 116 fSourceMomentum = aTrackInfo.fSourceMomentum; 117 fSourceEnergy = aTrackInfo.fSourceEnergy; 118 fSourceTime = aTrackInfo.fSourceTime; 119 fSuspendedStepID = -1; 120 121 return *this; 122 } 123 124 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 125 void RE01TrackInformation::SetSourceTrackInformation(const G4Track* aTrack) 126 { 127 fSourceTrackID = aTrack->GetTrackID(); 128 fSourceDefinition = aTrack->GetDefinition(); 129 fSourcePosition = aTrack->GetPosition(); 130 fSourceMomentum = aTrack->GetMomentum(); 131 fSourceEnergy = aTrack->GetTotalEnergy(); 132 fSourceTime = aTrack->GetGlobalTime(); 133 } 134 135 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 136 void RE01TrackInformation::Print() const 137 { 138 G4cout << "Source track ID " << fSourceTrackID << " (" << fSourceDefinition->GetParticleName() 139 << "," << fSourceEnergy / GeV << "[GeV]) at " << fSourcePosition << G4endl; 140 G4cout << "Original primary track ID " << fOriginalTrackID << " (" 141 << fParticleDefinition->GetParticleName() << "," << fOriginalEnergy / GeV << "[GeV])" 142 << G4endl; 143 } 144