Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // G4Trajectory 27 // 28 // Class description: 29 // 30 // This class represents the trajectory of a p 31 // It includes information of: 32 // 1) List of trajectory points which comp 33 // 2) Static information of the particle w 34 // trajectory; 35 // 3) Track ID and parent particle ID of t 36 37 // Contact: 38 // Questions and comments to this code shoul 39 // Katsuya Amako (e-mail: Katsuya.Amako@k 40 // Makoto Asai (e-mail: asai@slac.stanf 41 // Takashi Sasaki (e-mail: Takashi.Sasaki@ 42 // ------------------------------------------- 43 #ifndef G4Trajectory_hh 44 #define G4Trajectory_hh 1 45 46 #include "G4Allocator.hh" 47 #include "G4ParticleDefinition.hh" // Include 48 #include "G4Step.hh" 49 #include "G4Track.hh" 50 #include "G4TrajectoryPoint.hh" // Include fr 51 #include "G4VTrajectory.hh" 52 #include "G4ios.hh" // Include from 'system' 53 #include "globals.hh" // Include from 'global 54 #include "G4Threading.hh" 55 56 #include "trkgdefs.hh" 57 #include <stdlib.h> // Include from 'system' 58 59 #include <vector> 60 61 class G4Polyline; 62 class G4ClonedTrajectory; 63 64 class G4Trajectory : public G4VTrajectory 65 { 66 using G4TrajectoryPointContainer = std::vect 67 68 friend class G4ClonedTrajectory; 69 70 public: 71 // Constructors/Destructor 72 73 G4Trajectory() = default; 74 G4Trajectory(const G4Track* aTrack); 75 G4Trajectory(G4Trajectory&); 76 ~G4Trajectory() override; 77 78 // Operators 79 80 inline void* operator new(size_t); 81 inline void operator delete(void*); 82 inline G4bool operator==(const G4Trajectory& 83 84 // cloning with the master thread allocator 85 G4VTrajectory* CloneForMaster() const overri 86 87 // Get/Set functions 88 89 inline G4int GetTrackID() const override { r 90 inline G4int GetParentID() const override { 91 inline G4String GetParticleName() const over 92 inline G4double GetCharge() const override { 93 inline G4int GetPDGEncoding() const override 94 inline G4double GetInitialKineticEnergy() co 95 inline G4ThreeVector GetInitialMomentum() co 96 97 // Other member functions 98 99 void ShowTrajectory(std::ostream& os = G4cou 100 void DrawTrajectory() const override; 101 void AppendStep(const G4Step* aStep) overrid 102 G4int GetPointEntries() const override { ret 103 G4VTrajectoryPoint* GetPoint(G4int i) const 104 void MergeTrajectory(G4VTrajectory* secondTr 105 106 G4ParticleDefinition* GetParticleDefinition( 107 108 const std::map<G4String, G4AttDef>* GetAttDe 109 std::vector<G4AttValue>* CreateAttValues() c 110 111 private: 112 G4TrajectoryPointContainer* positionRecord = 113 G4int fTrackID = 0; 114 G4int fParentID = 0; 115 G4int PDGEncoding = 0; 116 G4double PDGCharge = 0.0; 117 G4String ParticleName = "dummy"; 118 G4double initialKineticEnergy = 0.0; 119 G4ThreeVector initialMomentum; 120 }; 121 122 extern G4TRACKING_DLL G4Allocator<G4Trajectory 123 124 inline void* G4Trajectory::operator new(size_t 125 { 126 if (aTrajectoryAllocator() == nullptr) { 127 aTrajectoryAllocator() = new G4Allocator<G 128 } 129 return (void*)aTrajectoryAllocator()->Malloc 130 } 131 132 inline void G4Trajectory::operator delete(void 133 { 134 aTrajectoryAllocator()->FreeSingle((G4Trajec 135 } 136 137 inline G4bool G4Trajectory::operator==(const G 138 139 #endif 140