Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 // G4Trajectory << 27 // 23 // 28 // Class description: << 24 // $Id: G4Trajectory.hh,v 1.17 2002/10/28 11:10:57 johna Exp $ >> 25 // GEANT4 tag $Name: geant4-05-01 $ >> 26 // >> 27 // >> 28 //--------------------------------------------------------------- >> 29 // >> 30 // G4Trajectory.hh >> 31 // >> 32 // class description: >> 33 // This class represents the trajectory of a particle tracked. >> 34 // It includes information of >> 35 // 1) List of trajectory points which compose the trajectory, >> 36 // 2) static information of particle which generated the >> 37 // trajectory, >> 38 // 3) trackID and parent particle ID of the trajectory, 29 // 39 // 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: 40 // Contact: 38 // Questions and comments to this code shoul 41 // Questions and comments to this code should be sent to 39 // Katsuya Amako (e-mail: Katsuya.Amako@k 42 // Katsuya Amako (e-mail: Katsuya.Amako@kek.jp) 40 // Makoto Asai (e-mail: asai@slac.stanf << 43 // Makoto Asai (e-mail: asai@kekvax.kek.jp) 41 // Takashi Sasaki (e-mail: Takashi.Sasaki@ 44 // Takashi Sasaki (e-mail: Takashi.Sasaki@kek.jp) 42 // ------------------------------------------- << 45 // 43 #ifndef G4Trajectory_hh << 46 // --------------------------------------------------------------- 44 #define G4Trajectory_hh 1 << 47 >> 48 class G4Trajectory; 45 49 >> 50 #ifndef G4Trajectory_h >> 51 #define G4Trajectory_h 1 >> 52 >> 53 #include "G4VTrajectory.hh" 46 #include "G4Allocator.hh" 54 #include "G4Allocator.hh" >> 55 #include <stdlib.h> // Include from 'system' >> 56 #include "G4ios.hh" // Include from 'system' >> 57 #include "g4std/vector" // G4RWTValOrderedVector >> 58 #include "globals.hh" // Include from 'global' 47 #include "G4ParticleDefinition.hh" // Include 59 #include "G4ParticleDefinition.hh" // Include from 'particle+matter' 48 #include "G4Step.hh" << 60 #include "G4TrajectoryPoint.hh" // Include from 'tracking' 49 #include "G4Track.hh" 61 #include "G4Track.hh" 50 #include "G4TrajectoryPoint.hh" // Include fr << 62 #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 63 61 class G4Polyline; << 64 class G4Polyline; // Forward declaration. 62 class G4ClonedTrajectory; << 63 65 >> 66 typedef G4std::vector<G4VTrajectoryPoint*> TrajectoryPointContainer; >> 67 /////////////////// 64 class G4Trajectory : public G4VTrajectory 68 class G4Trajectory : public G4VTrajectory >> 69 /////////////////// 65 { 70 { 66 using G4TrajectoryPointContainer = std::vect << 67 << 68 friend class G4ClonedTrajectory; << 69 71 70 public: << 72 //-------- 71 // Constructors/Destructor << 73 public: // with description >> 74 //-------- >> 75 >> 76 // Constructor/Destrcutor >> 77 >> 78 G4Trajectory(); >> 79 >> 80 G4Trajectory(const G4Track* aTrack); >> 81 G4Trajectory(G4Trajectory &); >> 82 virtual ~G4Trajectory(); >> 83 >> 84 // Operators >> 85 inline void* operator new(size_t); >> 86 inline void operator delete(void*); >> 87 inline int operator == (const G4Trajectory& right) const >> 88 {return (this==&right);} >> 89 >> 90 // Get/Set functions >> 91 inline G4int GetTrackID() const >> 92 { return fTrackID; } >> 93 inline G4int GetParentID() const >> 94 { return fParentID; } >> 95 inline G4String GetParticleName() const >> 96 { return ParticleName; } >> 97 inline G4double GetCharge() const >> 98 { return PDGCharge; } >> 99 inline G4int GetPDGEncoding() const >> 100 { return PDGEncoding; } >> 101 inline G4ThreeVector GetInitialMomentum() const >> 102 { return initialMomentum; } >> 103 >> 104 // Other member functions >> 105 virtual void ShowTrajectory(G4std::ostream& os=G4cout) const; >> 106 virtual void DrawTrajectory(G4int i_mode=0) const; >> 107 virtual void AppendStep(const G4Step* aStep); >> 108 virtual int GetPointEntries() const { return positionRecord->size(); } >> 109 virtual G4VTrajectoryPoint* GetPoint(G4int i) const >> 110 { return (*positionRecord)[i]; } >> 111 virtual void MergeTrajectory(G4VTrajectory* secondTrajectory); >> 112 >> 113 G4ParticleDefinition* GetParticleDefinition(); >> 114 >> 115 virtual const G4std::map<G4String,G4AttDef>* GetAttDefs() const; >> 116 virtual G4std::vector<G4AttValue>* CreateAttValues() const; >> 117 >> 118 //--------- >> 119 private: >> 120 //--------- >> 121 >> 122 TrajectoryPointContainer* positionRecord; >> 123 G4int fTrackID; >> 124 G4int fParentID; >> 125 G4int PDGEncoding; >> 126 G4double PDGCharge; >> 127 G4String ParticleName; >> 128 G4ThreeVector initialMomentum; 72 129 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 }; 130 }; 121 131 122 extern G4TRACKING_DLL G4Allocator<G4Trajectory << 132 extern G4Allocator<G4Trajectory> aTrajectoryAllocator; 123 133 124 inline void* G4Trajectory::operator new(size_t 134 inline void* G4Trajectory::operator new(size_t) 125 { 135 { 126 if (aTrajectoryAllocator() == nullptr) { << 136 void* aTrajectory; 127 aTrajectoryAllocator() = new G4Allocator<G << 137 aTrajectory = (void*)aTrajectoryAllocator.MallocSingle(); 128 } << 138 return aTrajectory; 129 return (void*)aTrajectoryAllocator()->Malloc << 130 } 139 } 131 140 132 inline void G4Trajectory::operator delete(void 141 inline void G4Trajectory::operator delete(void* aTrajectory) 133 { 142 { 134 aTrajectoryAllocator()->FreeSingle((G4Trajec << 143 aTrajectoryAllocator.FreeSingle((G4Trajectory*)aTrajectory); 135 } 144 } 136 145 137 inline G4bool G4Trajectory::operator==(const G << 138 << 139 #endif 146 #endif >> 147 >> 148 >> 149 >> 150 >> 151 >> 152 >> 153 >> 154 >> 155 >> 156 140 157