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