| 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 // <<
27 /// \file field/field04/src/F04ElementField.cc 26 /// \file field/field04/src/F04ElementField.cc
28 /// \brief Implementation of the F04ElementFie 27 /// \brief Implementation of the F04ElementField class
29 // 28 //
>> 29 //
>> 30 #include "G4GeometryManager.hh"
30 31
31 #include "F04ElementField.hh" 32 #include "F04ElementField.hh"
32 33
33 #include "F04GlobalField.hh" 34 #include "F04GlobalField.hh"
34 35
35 #include "G4SystemOfUnits.hh" 36 #include "G4SystemOfUnits.hh"
36 #include "G4TouchableHandle.hh" <<
37 <<
38 //....oooOO0OOooo........oooOO0OOooo........oo <<
39 37
40 G4ThreadLocal G4Navigator* F04ElementField::fN << 38 G4Navigator* F04ElementField::fNavigator;
41 39
42 //....oooOO0OOooo........oooOO0OOooo........oo 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
43 41
44 F04ElementField::F04ElementField(G4ThreeVector << 42 F04ElementField::F04ElementField(G4ThreeVector c, G4LogicalVolume* lv)
45 { 43 {
>> 44 fCenter = c;
>> 45
>> 46 fMinX = fMinY = fMinZ = -DBL_MAX;
>> 47 fMaxX = fMaxY = fMaxZ = DBL_MAX;
>> 48
46 F04GlobalField::GetObject()->AddElementField 49 F04GlobalField::GetObject()->AddElementField(this);
47 50
>> 51 fColor = "1,1,1";
>> 52
48 fUserLimits = new G4UserLimits(); 53 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 54
>> 55 fVolume = lv;
55 fVolume->SetVisAttributes(GetVisAttribute(fC 56 fVolume->SetVisAttributes(GetVisAttribute(fColor));
>> 57
>> 58 fMaxStep = 1*mm;
>> 59
>> 60 fUserLimits->SetMaxAllowedStep(fMaxStep);
>> 61
>> 62 fUserLimits->SetUserMaxTrackLength(500.*m);
>> 63 fUserLimits->SetUserMaxTime(10*ms);
>> 64 fUserLimits->SetUserMinEkine(0.1*MeV);
>> 65 // fUserLimits->SetUserMinRange(1*mm);
>> 66
56 fVolume->SetUserLimits(fUserLimits); 67 fVolume->SetUserLimits(fUserLimits);
57 } 68 }
58 69
59 //....oooOO0OOooo........oooOO0OOooo........oo 70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
60 71
61 void F04ElementField::Construct() 72 void F04ElementField::Construct()
62 { 73 {
63 G4Navigator* theNavigator = 74 G4Navigator* theNavigator =
64 G4TransportationManager::GetTransportation << 75 G4TransportationManager::GetTransportationManager()->
>> 76 GetNavigatorForTracking();
>> 77
65 if (!fNavigator) { 78 if (!fNavigator) {
66 fNavigator = new G4Navigator(); << 79 fNavigator = new G4Navigator();
67 if (theNavigator->GetWorldVolume()) fNavig << 80 if ( theNavigator->GetWorldVolume() )
>> 81 fNavigator->SetWorldVolume(theNavigator->GetWorldVolume());
>> 82 }
>> 83
>> 84 G4GeometryManager* geomManager = G4GeometryManager::GetInstance();
>> 85
>> 86 if (!geomManager->IsGeometryClosed()) {
>> 87 geomManager->OpenGeometry();
>> 88 geomManager->CloseGeometry(true);
68 } 89 }
69 90
70 fNavigator->LocateGlobalPointAndSetup(fCente << 91 fNavigator->LocateGlobalPointAndSetup(fCenter,0,false);
71 92
72 G4TouchableHandle touchable = fNavigator->Cr << 93 G4TouchableHistoryHandle touchable = fNavigator->
>> 94 CreateTouchableHistoryHandle();
73 95
74 G4int depth = touchable->GetHistoryDepth(); 96 G4int depth = touchable->GetHistoryDepth();
75 for (G4int i = 0; i < depth; ++i) { << 97 for (G4int i = 0; i<depth; ++i) {
76 if (touchable->GetVolume()->GetLogicalVolu << 98 if(touchable->GetVolume()->GetLogicalVolume() == fVolume)break;
77 touchable->MoveUpHistory(); << 99 touchable->MoveUpHistory();
78 } 100 }
79 101
80 // set fGlobal2local transform 102 // set fGlobal2local transform
81 fGlobal2local = touchable->GetHistory()->Get 103 fGlobal2local = touchable->GetHistory()->GetTopTransform();
82 104
83 // set global bounding box 105 // set global bounding box
84 G4double local[4], global[4]; 106 G4double local[4], global[4];
85 107
86 G4ThreeVector globalPosition; 108 G4ThreeVector globalPosition;
87 local[3] = 0.0; 109 local[3] = 0.0;
88 for (G4int i = 0; i < 2; ++i) { << 110 for (int i=0; i<2; ++i) {
89 local[0] = (i == 0 ? -1.0 : 1.0) * GetWidt << 111 local[0] = (i==0 ? -1.0 : 1.0) * GetWidth()/2.;
90 for (G4int j = 0; j < 2; ++j) { << 112 for (int j=0; j<2; ++j) {
91 local[1] = (j == 0 ? -1.0 : 1.0) * GetHe << 113 local[1] = (j==0 ? -1.0 : 1.0) * GetHeight()/2.;
92 for (G4int k = 0; k < 2; ++k) { << 114 for (int k=0; k<2; ++k) {
93 local[2] = (k == 0 ? -1.0 : 1.0) * Get << 115 local[2] = (k==0 ? -1.0 : 1.0) * GetLength()/2.;
94 G4ThreeVector localPosition(local[0], << 116 G4ThreeVector localPosition(local[0],local[1],local[2]);
95 globalPosition = fGlobal2local.Inverse << 117 globalPosition =
96 global[0] = globalPosition.x(); << 118 fGlobal2local.Inverse().TransformPoint(localPosition);
97 global[1] = globalPosition.y(); << 119 global[0] = globalPosition.x();
98 global[2] = globalPosition.z(); << 120 global[1] = globalPosition.y();
99 SetGlobalPoint(global); << 121 global[2] = globalPosition.z();
>> 122 SetGlobalPoint(global);
>> 123 }
100 } 124 }
101 } <<
102 } 125 }
103 } 126 }
104 127
105 //....oooOO0OOooo........oooOO0OOooo........oo 128 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
106 129
107 G4VisAttributes* F04ElementField::GetVisAttrib 130 G4VisAttributes* F04ElementField::GetVisAttribute(G4String color)
108 { 131 {
109 G4VisAttributes* p = nullptr; << 132 G4VisAttributes* p = NULL;
110 if (color.size() > 0 && (isdigit(color.c_str << 133 if(color.size() > 0 &&
111 G4double red = 0.0, green = 0.0, blue = 0. << 134 (isdigit(color.c_str()[0]) || color.c_str()[0] == '.')) {
112 if (sscanf(color.c_str(), "%lf,%lf,%lf", & << 135 G4double red=0.0, green=0.0, blue=0.0;
113 p = new G4VisAttributes(true, G4Color(re << 136 if (sscanf(color.c_str(),"%lf,%lf,%lf",&red,&green,&blue) == 3) {
114 } << 137 p = new G4VisAttributes(true,G4Color(red,green,blue));
115 else { << 138 } else {
116 G4cout << " Invalid color " << color << << 139 G4cout << " Invalid color " << color << G4endl;
117 } << 140 }
118 } << 141 }
119 142
120 if (!p) p = new G4VisAttributes(G4VisAttribu << 143 if (!p) p = new G4VisAttributes(G4VisAttributes::Invisible);
121 p->SetDaughtersInvisible(false); << 144 p->SetDaughtersInvisible(false);
122 145
123 return p; << 146 return p;
124 } 147 }
125 148