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 // 26 // 27 // 27 // 28 // 28 // 29 // Class Description: 29 // Class Description: 30 // 30 // 31 // Base class for the trajectory state 31 // Base class for the trajectory state 32 32 33 // History: 33 // History: 34 // - Created: P. Arce 34 // - Created: P. Arce 35 // ------------------------------------------- 35 // -------------------------------------------------------------------- 36 36 37 #ifndef G4ErrorTrajState_hh 37 #ifndef G4ErrorTrajState_hh 38 #define G4ErrorTrajState_hh 38 #define G4ErrorTrajState_hh 39 39 40 #include "globals.hh" 40 #include "globals.hh" 41 #include "G4Track.hh" 41 #include "G4Track.hh" 42 #include "G4Point3D.hh" 42 #include "G4Point3D.hh" 43 #include "G4Vector3D.hh" 43 #include "G4Vector3D.hh" 44 #include "G4ErrorTrajErr.hh" 44 #include "G4ErrorTrajErr.hh" 45 45 46 enum G4eTSType << 46 enum G4eTSType{ G4eTS_FREE, G4eTS_OS }; 47 { << 48 G4eTS_FREE, << 49 G4eTS_OS << 50 }; << 51 47 52 class G4ErrorTrajState 48 class G4ErrorTrajState 53 { 49 { 54 public: // with description 50 public: // with description 55 G4ErrorTrajState() << 51 56 : fCharge(0.) << 52 G4ErrorTrajState() : fCharge(0.), theG4Track(0), iverbose(0) {} 57 , theG4Track(0) << 53 58 , iverbose(0) << 54 G4ErrorTrajState( const G4String& partType, const G4Point3D& pos, 59 {} << 55 const G4Vector3D& mom, 60 << 56 const G4ErrorTrajErr& errmat = G4ErrorTrajErr(5,0) ); 61 G4ErrorTrajState(const G4String& partType, c << 57 // Constructor by providing particle, position and momentum 62 const G4Vector3D& mom, << 63 const G4ErrorTrajErr& errma << 64 // Constructor by providing particle, positi << 65 58 66 G4ErrorTrajState(const G4ErrorTrajState&); 59 G4ErrorTrajState(const G4ErrorTrajState&); 67 G4ErrorTrajState(G4ErrorTrajState&&); 60 G4ErrorTrajState(G4ErrorTrajState&&); 68 // The copy and move constructors << 61 // The copy and move constructors 69 62 70 virtual ~G4ErrorTrajState() {} << 63 virtual ~G4ErrorTrajState(){} 71 64 72 G4ErrorTrajState& operator=(const G4ErrorTra << 65 G4ErrorTrajState& operator = (const G4ErrorTrajState&); 73 G4ErrorTrajState& operator=(G4ErrorTrajState << 66 G4ErrorTrajState& operator = (G4ErrorTrajState&&); 74 // The copy and move assignment operators << 67 // The copy and move assignment operators 75 68 76 void SetData(const G4String& partType, const << 69 void SetData( const G4String& partType, const G4Point3D& pos, 77 const G4Vector3D& mom); << 70 const G4Vector3D& mom ); 78 // Set particle, position and momentum << 71 // Set particle, position and momentum 79 72 80 void BuildCharge(); 73 void BuildCharge(); 81 // Build charge based on particle type << 74 // Build charge based on particle type 82 << 83 friend std::ostream& operator<<(std::ostream << 84 << 85 virtual G4int PropagateError(const G4Track*) << 86 // Propagate the error along the step << 87 << 88 virtual G4int Update(const G4Track*) { retur << 89 // Update parameters from G4Track << 90 << 91 void UpdatePosMom(const G4Point3D& pos, cons << 92 // Update position and momentum << 93 75 94 void DumpPosMomError(std::ostream& out = G4c << 76 friend 95 // Dump position, momentum and error << 77 std::ostream& operator<<(std::ostream&, const G4ErrorTrajState& ts); 96 78 97 virtual void Dump(std::ostream& out = G4cout << 79 virtual G4int PropagateError( const G4Track* ); 98 // Abstract method to dump all TrajState par << 80 // Propagate the error along the step 99 81 100 // Set and Get methods << 82 virtual G4int Update( const G4Track* ){ return -1; } >> 83 // Update parameters from G4Track >> 84 >> 85 void UpdatePosMom( const G4Point3D& pos, const G4Vector3D& mom ); >> 86 // Update position and momentum >> 87 >> 88 void DumpPosMomError( std::ostream& out = G4cout ) const; >> 89 // Dump position, momentum and error >> 90 >> 91 virtual void Dump( std::ostream& out = G4cout ) const = 0; >> 92 // Abstract method to dump all TrajState parameters >> 93 >> 94 // Set and Get methods >> 95 >> 96 const G4String& GetParticleType() const >> 97 { return fParticleType;} >> 98 void SetParticleType( const G4String& partType ) >> 99 { fParticleType = partType;} >> 100 >> 101 G4Point3D GetPosition() const >> 102 { return fPosition; } >> 103 virtual void SetPosition( const G4Point3D pos ) >> 104 { fPosition = pos; } >> 105 >> 106 G4Vector3D GetMomentum() const >> 107 { return fMomentum; } >> 108 virtual void SetMomentum( const G4Vector3D& mom ) >> 109 { fMomentum = mom; } >> 110 >> 111 G4ErrorTrajErr GetError() const >> 112 { return fError; } >> 113 virtual void SetError( G4ErrorTrajErr em ) >> 114 { fError = em; } >> 115 >> 116 G4Track* GetG4Track() const >> 117 { return theG4Track; } >> 118 void SetG4Track( G4Track* trk ) >> 119 { theG4Track = trk; } >> 120 >> 121 G4double GetCharge() const >> 122 { return fCharge; } >> 123 void SetCharge( G4double ch ) >> 124 { fCharge = ch; } 101 125 102 const G4String& GetParticleType() const { re << 126 virtual G4eTSType GetTSType() const 103 void SetParticleType(const G4String& partTyp << 127 { return theTSType; } 104 << 105 G4Point3D GetPosition() const { return fPosi << 106 virtual void SetPosition(const G4Point3D pos << 107 << 108 G4Vector3D GetMomentum() const { return fMom << 109 virtual void SetMomentum(const G4Vector3D& m << 110 << 111 G4ErrorTrajErr GetError() const { return fEr << 112 virtual void SetError(G4ErrorTrajErr em) { f << 113 << 114 G4Track* GetG4Track() const { return theG4Tr << 115 void SetG4Track(G4Track* trk) { theG4Track = << 116 << 117 G4double GetCharge() const { return fCharge; << 118 void SetCharge(G4double ch) { fCharge = ch; << 119 << 120 virtual G4eTSType GetTSType() const { return << 121 128 122 protected: 129 protected: >> 130 123 G4String fParticleType; 131 G4String fParticleType; 124 G4Point3D fPosition; 132 G4Point3D fPosition; 125 G4Vector3D fMomentum; 133 G4Vector3D fMomentum; 126 G4double fCharge = 0.; 134 G4double fCharge = 0.; 127 G4ErrorTrajErr fError; 135 G4ErrorTrajErr fError; 128 136 129 G4eTSType theTSType; 137 G4eTSType theTSType; 130 138 131 G4Track* theG4Track = nullptr; 139 G4Track* theG4Track = nullptr; 132 140 133 G4int iverbose = 0; 141 G4int iverbose = 0; 134 }; 142 }; 135 143 136 #endif 144 #endif 137 145