Geant4 Cross Reference |
>> 1 // This code implementation is the intellectual property of >> 2 // the GEANT4 collaboration. >> 3 // >> 4 // By copying, distributing or modifying the Program (or any work >> 5 // based on the Program) you indicate your acceptance of this statement, >> 6 // and all its terms. >> 7 // >> 8 // $Id: G4Trajectory.hh,v 1.9 2000/01/26 04:20:31 asaim Exp $ >> 9 // GEANT4 tag $Name: geant4-02-00 $ >> 10 // >> 11 // >> 12 //--------------------------------------------------------------- >> 13 // >> 14 // G4Trajectory.hh >> 15 // >> 16 // Description: >> 17 // This class represents the trajectory of a particle tracked. >> 18 // It includes information of >> 19 // 1) List of trajectory points which compose the trajectory, >> 20 // 2) static information of particle which generated the >> 21 // trajectory, >> 22 // 3) trackID and parent particle ID of the trajectory, 1 // 23 // 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: 24 // Contact: 38 // Questions and comments to this code shoul 25 // Questions and comments to this code should be sent to 39 // Katsuya Amako (e-mail: Katsuya.Amako@k 26 // Katsuya Amako (e-mail: Katsuya.Amako@kek.jp) 40 // Makoto Asai (e-mail: asai@slac.stanf << 27 // Makoto Asai (e-mail: asai@kekvax.kek.jp) 41 // Takashi Sasaki (e-mail: Takashi.Sasaki@ 28 // Takashi Sasaki (e-mail: Takashi.Sasaki@kek.jp) 42 // ------------------------------------------- << 29 // 43 #ifndef G4Trajectory_hh << 30 // --------------------------------------------------------------- 44 #define G4Trajectory_hh 1 << 31 >> 32 class G4Trajectory; 45 33 >> 34 #ifndef G4Trajectory_h >> 35 #define G4Trajectory_h 1 >> 36 >> 37 #include "G4VTrajectory.hh" 46 #include "G4Allocator.hh" 38 #include "G4Allocator.hh" >> 39 #include <stdlib.h> // Include from 'system' >> 40 #include "G4ios.hh" // Include from 'system' >> 41 #include "g4rw/tpordvec.h" // G4RWTValOrderedVector >> 42 #include "globals.hh" // Include from 'global' 47 #include "G4ParticleDefinition.hh" // Include 43 #include "G4ParticleDefinition.hh" // Include from 'particle+matter' 48 #include "G4Step.hh" << 44 #include "G4TrajectoryPoint.hh" // Include from 'tracking' 49 #include "G4Track.hh" 45 #include "G4Track.hh" 50 #include "G4TrajectoryPoint.hh" // Include fr << 46 #include "G4Step.hh" 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 47 61 class G4Polyline; << 48 class G4Polyline; // Forward declaration. 62 class G4ClonedTrajectory; << 63 49 >> 50 /////////////////// 64 class G4Trajectory : public G4VTrajectory 51 class G4Trajectory : public G4VTrajectory >> 52 /////////////////// 65 { 53 { 66 using G4TrajectoryPointContainer = std::vect << 67 54 68 friend class G4ClonedTrajectory; << 55 //-------- >> 56 public: >> 57 //-------- >> 58 >> 59 // Constructor/Destrcutor >> 60 >> 61 G4Trajectory(); >> 62 >> 63 G4Trajectory(const G4Track* aTrack); >> 64 G4Trajectory(G4Trajectory &); >> 65 virtual ~G4Trajectory(); >> 66 >> 67 // Operators >> 68 inline void* operator new(size_t); >> 69 inline void operator delete(void*); >> 70 inline int operator == (const G4Trajectory& right) const >> 71 {return (this==&right);} >> 72 >> 73 // Get/Set functions >> 74 inline G4int GetTrackID() const >> 75 { return fTrackID; } >> 76 inline G4int GetParentID() const >> 77 { return fParentID; } >> 78 inline G4String GetParticleName() const >> 79 { return ParticleName; } >> 80 inline G4double GetCharge() const >> 81 { return PDGCharge; } >> 82 inline G4int GetPDGEncoding() const >> 83 { return PDGEncoding; } >> 84 >> 85 // Other member functions >> 86 virtual void ShowTrajectory() const; >> 87 virtual void DrawTrajectory(G4int i_mode=0) const; >> 88 virtual void AppendStep(const G4Step* aStep); >> 89 virtual int GetPointEntries() const { return positionRecord->entries(); } >> 90 virtual G4VTrajectoryPoint* GetPoint(G4int i) const >> 91 { return (*positionRecord)[i]; } >> 92 virtual void MergeTrajectory(G4VTrajectory* secondTrajectory); >> 93 >> 94 G4ParticleDefinition* GetParticleDefinition(); >> 95 >> 96 //--------- >> 97 private: >> 98 //--------- >> 99 >> 100 G4RWTPtrOrderedVector<G4VTrajectoryPoint>* positionRecord; >> 101 G4int fTrackID; >> 102 G4int fParentID; >> 103 G4String ParticleName; >> 104 G4double PDGCharge; >> 105 G4int PDGEncoding; 69 106 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 }; 107 }; 121 108 122 extern G4TRACKING_DLL G4Allocator<G4Trajectory << 109 extern G4Allocator<G4Trajectory> aTrajectoryAllocator; 123 110 124 inline void* G4Trajectory::operator new(size_t 111 inline void* G4Trajectory::operator new(size_t) 125 { 112 { 126 if (aTrajectoryAllocator() == nullptr) { << 113 void* aTrajectory; 127 aTrajectoryAllocator() = new G4Allocator<G << 114 aTrajectory = (void*)aTrajectoryAllocator.MallocSingle(); 128 } << 115 return aTrajectory; 129 return (void*)aTrajectoryAllocator()->Malloc << 130 } 116 } 131 117 132 inline void G4Trajectory::operator delete(void 118 inline void G4Trajectory::operator delete(void* aTrajectory) 133 { 119 { 134 aTrajectoryAllocator()->FreeSingle((G4Trajec << 120 aTrajectoryAllocator.FreeSingle((G4Trajectory*)aTrajectory); 135 } 121 } 136 122 137 inline G4bool G4Trajectory::operator==(const G << 138 << 139 #endif 123 #endif >> 124 >> 125 >> 126 >> 127 >> 128 >> 129 >> 130 >> 131 >> 132 >> 133 140 134