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 // G4SmoothTrajectoryPoint class implementatio << 27 // 26 // 28 // Contact: << 27 // $Id: G4SmoothTrajectoryPoint.cc,v 1.15 2006/06/29 21:16:01 gunter Exp $ 29 // Questions and comments to this code shoul << 28 // GEANT4 tag $Name: geant4-09-01-patch-03 $ 30 // Katsuya Amako (e-mail: Katsuya.Amako@k << 29 // 31 // Makoto Asai (e-mail: asai@slac.stanf << 30 // 32 // Takashi Sasaki (e-mail: Takashi.Sasaki@ << 31 // --------------------------------------------------------------- 33 // ------------------------------------------- << 32 // >> 33 // G4SmoothTrajectoryPoint.cc >> 34 // >> 35 // --------------------------------------------------------------- 34 36 35 #include "G4SmoothTrajectoryPoint.hh" 37 #include "G4SmoothTrajectoryPoint.hh" 36 38 37 #include "G4AttDef.hh" << 38 #include "G4AttDefStore.hh" 39 #include "G4AttDefStore.hh" >> 40 #include "G4AttDef.hh" 39 #include "G4AttValue.hh" 41 #include "G4AttValue.hh" 40 #include "G4UnitsTable.hh" 42 #include "G4UnitsTable.hh" 41 43 42 // #define G4ATTDEBUG << 44 //#define G4ATTDEBUG 43 #ifdef G4ATTDEBUG 45 #ifdef G4ATTDEBUG 44 # include "G4AttCheck.hh" << 46 #include "G4AttCheck.hh" 45 #endif 47 #endif 46 48 47 G4Allocator<G4SmoothTrajectoryPoint>*& aSmooth << 49 G4Allocator<G4SmoothTrajectoryPoint> aSmoothTrajectoryPointAllocator; >> 50 >> 51 G4SmoothTrajectoryPoint::G4SmoothTrajectoryPoint() >> 52 : fAuxiliaryPointVector(0) 48 { 53 { 49 G4ThreadLocalStatic G4Allocator<G4SmoothTraj << 54 fPosition = G4ThreeVector(0.,0.,0.); 50 return _instance; << 51 } 55 } 52 56 53 G4SmoothTrajectoryPoint::G4SmoothTrajectoryPoi << 57 G4SmoothTrajectoryPoint::G4SmoothTrajectoryPoint(G4ThreeVector pos) >> 58 : fAuxiliaryPointVector(0) >> 59 { >> 60 fPosition = pos; >> 61 } 54 62 55 G4SmoothTrajectoryPoint::G4SmoothTrajectoryPoi << 63 G4SmoothTrajectoryPoint::G4SmoothTrajectoryPoint(G4ThreeVector pos, 56 G4ThreeVector pos, std::vector<G4ThreeVector << 64 std::vector<G4ThreeVector>* auxiliaryPoints) 57 : fPosition(pos), fAuxiliaryPointVector(auxi << 65 : fPosition(pos), >> 66 fAuxiliaryPointVector(auxiliaryPoints) 58 {} 67 {} 59 68 60 G4SmoothTrajectoryPoint::G4SmoothTrajectoryPoi << 69 G4SmoothTrajectoryPoint::G4SmoothTrajectoryPoint(const G4SmoothTrajectoryPoint &right) 61 : fPosition(right.fPosition), fAuxiliaryPoin << 70 : G4VTrajectoryPoint(), 62 {} << 71 fPosition(right.fPosition),fAuxiliaryPointVector(right.fAuxiliaryPointVector) >> 72 { >> 73 } >> 74 >> 75 G4SmoothTrajectoryPoint::~G4SmoothTrajectoryPoint() >> 76 { >> 77 if(fAuxiliaryPointVector) { >> 78 delete fAuxiliaryPointVector; >> 79 } >> 80 } 63 81 64 G4SmoothTrajectoryPoint::~G4SmoothTrajectoryPo << 65 82 66 const std::map<G4String, G4AttDef>* G4SmoothTr << 83 const std::map<G4String,G4AttDef>* >> 84 G4SmoothTrajectoryPoint::GetAttDefs() const 67 { 85 { 68 G4bool isNew; 86 G4bool isNew; 69 std::map<G4String, G4AttDef>* store = << 87 std::map<G4String,G4AttDef>* store 70 G4AttDefStore::GetInstance("G4SmoothTrajec << 88 = G4AttDefStore::GetInstance("G4SmoothTrajectoryPoint",isNew); 71 if (isNew) { 89 if (isNew) { 72 G4String Pos("Pos"); 90 G4String Pos("Pos"); 73 (*store)[Pos] = G4AttDef(Pos, "Step Positi << 91 (*store)[Pos] = G4AttDef(Pos, "Step Position", >> 92 "Physics","G4BestUnit","G4ThreeVector"); 74 G4String Aux("Aux"); 93 G4String Aux("Aux"); 75 (*store)[Aux] = << 94 (*store)[Aux] = G4AttDef(Aux, "Auxiliary Point Position", 76 G4AttDef(Aux, "Auxiliary Point Position" << 95 "Physics","G4BestUnit","G4ThreeVector"); 77 } 96 } 78 return store; 97 return store; 79 } 98 } 80 99 81 std::vector<G4AttValue>* G4SmoothTrajectoryPoi 100 std::vector<G4AttValue>* G4SmoothTrajectoryPoint::CreateAttValues() const 82 { 101 { 83 auto values = new std::vector<G4AttValue>; << 102 std::vector<G4AttValue>* values = new std::vector<G4AttValue>; 84 103 85 if (fAuxiliaryPointVector != nullptr) { << 104 if (fAuxiliaryPointVector) { 86 for (const auto& iAux : *fAuxiliaryPointVe << 105 std::vector<G4ThreeVector>::iterator iAux; 87 values->push_back(G4AttValue("Aux", G4Be << 106 for (iAux = fAuxiliaryPointVector->begin(); >> 107 iAux != fAuxiliaryPointVector->end(); ++iAux) { >> 108 values->push_back(G4AttValue("Aux",G4BestUnit(*iAux,"Length"),"")); 88 } 109 } 89 } 110 } 90 111 91 values->push_back(G4AttValue("Pos", G4BestUn << 112 values->push_back(G4AttValue("Pos",G4BestUnit(fPosition,"Length"),"")); 92 113 93 #ifdef G4ATTDEBUG 114 #ifdef G4ATTDEBUG 94 G4cout << G4AttCheck(values, GetAttDefs()); << 115 G4cout << G4AttCheck(values,GetAttDefs()); 95 #endif 116 #endif 96 117 97 return values; 118 return values; 98 } 119 } 99 120