Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 /// \file field/field04/src/F04ElementField.cc 28 /// \brief Implementation of the F04ElementFie 29 // 30 31 #include "F04ElementField.hh" 32 33 #include "F04GlobalField.hh" 34 35 #include "G4SystemOfUnits.hh" 36 #include "G4TouchableHandle.hh" 37 38 //....oooOO0OOooo........oooOO0OOooo........oo 39 40 G4ThreadLocal G4Navigator* F04ElementField::fN 41 42 //....oooOO0OOooo........oooOO0OOooo........oo 43 44 F04ElementField::F04ElementField(G4ThreeVector 45 { 46 F04GlobalField::GetObject()->AddElementField 47 48 fUserLimits = new G4UserLimits(); 49 fUserLimits->SetMaxAllowedStep(fMaxStep); 50 fUserLimits->SetUserMaxTrackLength(500. * m) 51 fUserLimits->SetUserMaxTime(10 * ms); 52 fUserLimits->SetUserMinEkine(0.1 * MeV); 53 // fUserLimits->SetUserMinRange(1*mm); 54 55 fVolume->SetVisAttributes(GetVisAttribute(fC 56 fVolume->SetUserLimits(fUserLimits); 57 } 58 59 //....oooOO0OOooo........oooOO0OOooo........oo 60 61 void F04ElementField::Construct() 62 { 63 G4Navigator* theNavigator = 64 G4TransportationManager::GetTransportation 65 if (!fNavigator) { 66 fNavigator = new G4Navigator(); 67 if (theNavigator->GetWorldVolume()) fNavig 68 } 69 70 fNavigator->LocateGlobalPointAndSetup(fCente 71 72 G4TouchableHandle touchable = fNavigator->Cr 73 74 G4int depth = touchable->GetHistoryDepth(); 75 for (G4int i = 0; i < depth; ++i) { 76 if (touchable->GetVolume()->GetLogicalVolu 77 touchable->MoveUpHistory(); 78 } 79 80 // set fGlobal2local transform 81 fGlobal2local = touchable->GetHistory()->Get 82 83 // set global bounding box 84 G4double local[4], global[4]; 85 86 G4ThreeVector globalPosition; 87 local[3] = 0.0; 88 for (G4int i = 0; i < 2; ++i) { 89 local[0] = (i == 0 ? -1.0 : 1.0) * GetWidt 90 for (G4int j = 0; j < 2; ++j) { 91 local[1] = (j == 0 ? -1.0 : 1.0) * GetHe 92 for (G4int k = 0; k < 2; ++k) { 93 local[2] = (k == 0 ? -1.0 : 1.0) * Get 94 G4ThreeVector localPosition(local[0], 95 globalPosition = fGlobal2local.Inverse 96 global[0] = globalPosition.x(); 97 global[1] = globalPosition.y(); 98 global[2] = globalPosition.z(); 99 SetGlobalPoint(global); 100 } 101 } 102 } 103 } 104 105 //....oooOO0OOooo........oooOO0OOooo........oo 106 107 G4VisAttributes* F04ElementField::GetVisAttrib 108 { 109 G4VisAttributes* p = nullptr; 110 if (color.size() > 0 && (isdigit(color.c_str 111 G4double red = 0.0, green = 0.0, blue = 0. 112 if (sscanf(color.c_str(), "%lf,%lf,%lf", & 113 p = new G4VisAttributes(true, G4Color(re 114 } 115 else { 116 G4cout << " Invalid color " << color << 117 } 118 } 119 120 if (!p) p = new G4VisAttributes(G4VisAttribu 121 p->SetDaughtersInvisible(false); 122 123 return p; 124 } 125