Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // class G4TouchableHistory Implementation 27 // 28 // Created: Paul Kent, August 1996 29 // ---------------------------------------------------------------------- 30 31 #include "G4TouchableHistory.hh" 32 33 G4Allocator<G4TouchableHistory>*& aTouchableHistoryAllocator() 34 { 35 G4ThreadLocalStatic G4Allocator<G4TouchableHistory>* _instance = nullptr; 36 return _instance; 37 } 38 39 G4TouchableHistory::G4TouchableHistory() 40 : ftlate(G4ThreeVector(0.,0.,0.)) 41 42 { 43 fhistory.SetFirstEntry(nullptr); 44 } 45 46 G4TouchableHistory::G4TouchableHistory( const G4NavigationHistory& history ) 47 : fhistory(history) 48 { 49 const G4AffineTransform& tf = fhistory.GetTopTransform(); 50 ftlate = tf.InverseNetTranslation(); 51 frot = tf.InverseNetRotation(); 52 } 53 54 G4TouchableHistory::~G4TouchableHistory() = default; 55 56 const G4ThreeVector& 57 G4TouchableHistory::GetTranslation(G4int depth) const 58 { 59 // The value returned will change at the next call 60 // Copy it if you want to use it! 61 // 62 static G4ThreadLocal G4ThreeVector* ctrans = nullptr; 63 if ( ctrans == nullptr ) { ctrans = new G4ThreeVector; } 64 if(depth==0.0) 65 { 66 return ftlate; 67 } 68 else 69 { 70 *ctrans = 71 fhistory.GetTransform(CalculateHistoryIndex(depth)).NetTranslation(); 72 return *ctrans; 73 } 74 } 75 76 const G4RotationMatrix* 77 G4TouchableHistory::GetRotation(G4int depth) const 78 { 79 // The value returned will change at the next call 80 // Copy it if you want to use it! 81 // 82 static G4ThreadLocal G4RotationMatrix* rotM = nullptr; 83 if ( rotM == nullptr ) { rotM = new G4RotationMatrix; } 84 85 if(depth==0.0) 86 { 87 return &frot; 88 } 89 else 90 { 91 *rotM = fhistory.GetTransform(CalculateHistoryIndex(depth)).NetRotation(); 92 return rotM; 93 } 94 } 95