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 // G4Trajectory << 27 // 26 // 28 // Class description: << 27 // $Id: G4Trajectory.hh,v 1.24 2006/06/29 21:15:41 gunter Exp $ >> 28 // GEANT4 tag $Name: geant4-09-00 $ >> 29 // >> 30 //--------------------------------------------------------------- >> 31 // >> 32 // G4Trajectory.hh >> 33 // >> 34 // class description: >> 35 // This class represents the trajectory of a particle tracked. >> 36 // It includes information of >> 37 // 1) List of trajectory points which compose the trajectory, >> 38 // 2) static information of particle which generated the >> 39 // trajectory, >> 40 // 3) trackID and parent particle ID of the trajectory, 29 // 41 // 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: 42 // Contact: 38 // Questions and comments to this code shoul 43 // Questions and comments to this code should be sent to 39 // Katsuya Amako (e-mail: Katsuya.Amako@k 44 // Katsuya Amako (e-mail: Katsuya.Amako@kek.jp) 40 // Makoto Asai (e-mail: asai@slac.stanf << 45 // Makoto Asai (e-mail: asai@kekvax.kek.jp) 41 // Takashi Sasaki (e-mail: Takashi.Sasaki@ 46 // Takashi Sasaki (e-mail: Takashi.Sasaki@kek.jp) 42 // ------------------------------------------- << 47 // 43 #ifndef G4Trajectory_hh << 48 // --------------------------------------------------------------- 44 #define G4Trajectory_hh 1 << 49 >> 50 class G4Trajectory; >> 51 >> 52 #ifndef G4Trajectory_h >> 53 #define G4Trajectory_h 1 45 54 >> 55 #include "G4VTrajectory.hh" 46 #include "G4Allocator.hh" 56 #include "G4Allocator.hh" >> 57 #include <stdlib.h> // Include from 'system' >> 58 #include "G4ios.hh" // Include from 'system' >> 59 #include <vector> // G4RWTValOrderedVector >> 60 #include "globals.hh" // Include from 'global' 47 #include "G4ParticleDefinition.hh" // Include 61 #include "G4ParticleDefinition.hh" // Include from 'particle+matter' 48 #include "G4Step.hh" << 62 #include "G4TrajectoryPoint.hh" // Include from 'tracking' 49 #include "G4Track.hh" 63 #include "G4Track.hh" 50 #include "G4TrajectoryPoint.hh" // Include fr << 64 #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 65 61 class G4Polyline; << 66 class G4Polyline; // Forward declaration. 62 class G4ClonedTrajectory; << 63 67 >> 68 typedef std::vector<G4VTrajectoryPoint*> TrajectoryPointContainer; >> 69 /////////////////// 64 class G4Trajectory : public G4VTrajectory 70 class G4Trajectory : public G4VTrajectory >> 71 /////////////////// 65 { 72 { 66 using G4TrajectoryPointContainer = std::vect << 67 73 68 friend class G4ClonedTrajectory; << 74 //-------- >> 75 public: // with description >> 76 //-------- >> 77 >> 78 // Constructor/Destrcutor >> 79 >> 80 G4Trajectory(); >> 81 >> 82 G4Trajectory(const G4Track* aTrack); >> 83 G4Trajectory(G4Trajectory &); >> 84 virtual ~G4Trajectory(); >> 85 >> 86 // Operators >> 87 inline void* operator new(size_t); >> 88 inline void operator delete(void*); >> 89 inline int operator == (const G4Trajectory& right) const >> 90 {return (this==&right);} >> 91 >> 92 // Get/Set functions >> 93 inline G4int GetTrackID() const >> 94 { return fTrackID; } >> 95 inline G4int GetParentID() const >> 96 { return fParentID; } >> 97 inline G4String GetParticleName() const >> 98 { return ParticleName; } >> 99 inline G4double GetCharge() const >> 100 { return PDGCharge; } >> 101 inline G4int GetPDGEncoding() const >> 102 { return PDGEncoding; } >> 103 inline G4ThreeVector GetInitialMomentum() const >> 104 { return initialMomentum; } >> 105 >> 106 // Other member functions >> 107 virtual void ShowTrajectory(std::ostream& os=G4cout) const; >> 108 virtual void DrawTrajectory(G4int i_mode=0) const; >> 109 virtual void AppendStep(const G4Step* aStep); >> 110 virtual int GetPointEntries() const { return positionRecord->size(); } >> 111 virtual G4VTrajectoryPoint* GetPoint(G4int i) const >> 112 { return (*positionRecord)[i]; } >> 113 virtual void MergeTrajectory(G4VTrajectory* secondTrajectory); >> 114 >> 115 G4ParticleDefinition* GetParticleDefinition(); >> 116 >> 117 virtual const std::map<G4String,G4AttDef>* GetAttDefs() const; >> 118 virtual std::vector<G4AttValue>* CreateAttValues() const; >> 119 >> 120 //--------- >> 121 private: >> 122 //--------- >> 123 >> 124 TrajectoryPointContainer* positionRecord; >> 125 G4int fTrackID; >> 126 G4int fParentID; >> 127 G4int PDGEncoding; >> 128 G4double PDGCharge; >> 129 G4String ParticleName; >> 130 G4ThreeVector initialMomentum; 69 131 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 }; 132 }; 121 133 122 extern G4TRACKING_DLL G4Allocator<G4Trajectory << 134 #if defined G4TRACKING_ALLOC_EXPORT >> 135 extern G4DLLEXPORT G4Allocator<G4Trajectory> aTrajectoryAllocator; >> 136 #else >> 137 extern G4DLLIMPORT G4Allocator<G4Trajectory> aTrajectoryAllocator; >> 138 #endif 123 139 124 inline void* G4Trajectory::operator new(size_t 140 inline void* G4Trajectory::operator new(size_t) 125 { 141 { 126 if (aTrajectoryAllocator() == nullptr) { << 142 void* aTrajectory; 127 aTrajectoryAllocator() = new G4Allocator<G << 143 aTrajectory = (void*)aTrajectoryAllocator.MallocSingle(); 128 } << 144 return aTrajectory; 129 return (void*)aTrajectoryAllocator()->Malloc << 130 } 145 } 131 146 132 inline void G4Trajectory::operator delete(void 147 inline void G4Trajectory::operator delete(void* aTrajectory) 133 { 148 { 134 aTrajectoryAllocator()->FreeSingle((G4Trajec << 149 aTrajectoryAllocator.FreeSingle((G4Trajectory*)aTrajectory); 135 } 150 } 136 151 137 inline G4bool G4Trajectory::operator==(const G << 138 << 139 #endif 152 #endif >> 153 >> 154 >> 155 >> 156 >> 157 >> 158 >> 159 >> 160 >> 161 >> 162 140 163