| 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 // 26 //
27 /// \file field/field06/src/F06DetectorConstru 27 /// \file field/field06/src/F06DetectorConstruction.cc
28 /// \brief Implementation of the F06DetectorCo 28 /// \brief Implementation of the F06DetectorConstruction class
29 // 29 //
30 30
31 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 33
34 #include "F06DetectorConstruction.hh" 34 #include "F06DetectorConstruction.hh"
35 35
36 #include "G4Box.hh" <<
37 #include "G4Colour.hh" <<
38 #include "G4FieldManager.hh" <<
39 #include "G4GeometryManager.hh" <<
40 #include "G4LogicalVolume.hh" <<
41 #include "G4LogicalVolumeStore.hh" <<
42 #include "G4Material.hh" 36 #include "G4Material.hh"
43 #include "G4NistManager.hh" 37 #include "G4NistManager.hh"
>> 38
>> 39 #include "G4Box.hh"
>> 40 #include "G4LogicalVolume.hh"
44 #include "G4PVPlacement.hh" 41 #include "G4PVPlacement.hh"
>> 42
>> 43 #include "G4GeometryManager.hh"
45 #include "G4PhysicalVolumeStore.hh" 44 #include "G4PhysicalVolumeStore.hh"
46 #include "G4RepleteEofM.hh" << 45 #include "G4LogicalVolumeStore.hh"
47 #include "G4SolidStore.hh" 46 #include "G4SolidStore.hh"
>> 47
>> 48 #include "G4VisAttributes.hh"
>> 49 #include "G4Colour.hh"
>> 50
>> 51 #include "G4UserLimits.hh"
48 #include "G4SystemOfUnits.hh" 52 #include "G4SystemOfUnits.hh"
49 #include "G4TransportationManager.hh" << 53
50 #include "G4UniformGravityField.hh" 54 #include "G4UniformGravityField.hh"
51 #include "G4UserLimits.hh" <<
52 #include "G4VisAttributes.hh" <<
53 // #include "G4EqGravityField.hh" <<
54 55
55 #include "G4ChordFinder.hh" << 56 #include "G4FieldManager.hh"
>> 57 #include "G4TransportationManager.hh"
>> 58
>> 59 #include "G4RepleteEofM.hh"
>> 60 //#include "G4EqGravityField.hh"
>> 61
56 #include "G4ClassicalRK4.hh" 62 #include "G4ClassicalRK4.hh"
57 #include "G4IntegrationDriver.hh" <<
58 #include "G4MagIntegratorStepper.hh" 63 #include "G4MagIntegratorStepper.hh"
>> 64 #include "G4IntegrationDriver.hh"
>> 65 #include "G4ChordFinder.hh"
59 #include "G4PropagatorInField.hh" 66 #include "G4PropagatorInField.hh"
60 67
61 //....oooOO0OOooo........oooOO0OOooo........oo 68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
62 69
63 F06DetectorConstruction::F06DetectorConstructi << 70 F06DetectorConstruction::F06DetectorConstruction()
>> 71 : fVacuum(0), fSolidWorld(0), fLogicWorld(0), fPhysiWorld(0)
64 { 72 {
65 // materials 73 // materials
66 DefineMaterials(); 74 DefineMaterials();
67 } 75 }
68 76
69 //....oooOO0OOooo........oooOO0OOooo........oo 77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
70 78
71 F06DetectorConstruction::~F06DetectorConstruct 79 F06DetectorConstruction::~F06DetectorConstruction()
72 { 80 {
73 delete fField; << 81 if (fField) delete fField;
74 } 82 }
75 83
76 //....oooOO0OOooo........oooOO0OOooo........oo 84 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
77 85
78 void F06DetectorConstruction::DefineMaterials( 86 void F06DetectorConstruction::DefineMaterials()
79 { 87 {
80 G4NistManager* nistMan = G4NistManager::Inst 88 G4NistManager* nistMan = G4NistManager::Instance();
81 89
82 fVacuum = nistMan->FindOrBuildMaterial("G4_G 90 fVacuum = nistMan->FindOrBuildMaterial("G4_Galactic");
83 91
84 G4cout << *(G4Material::GetMaterialTable()) 92 G4cout << *(G4Material::GetMaterialTable()) << G4endl;
85 } 93 }
86 94
87 //....oooOO0OOooo........oooOO0OOooo........oo 95 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
88 96
89 G4VPhysicalVolume* F06DetectorConstruction::Co 97 G4VPhysicalVolume* F06DetectorConstruction::Construct()
90 { 98 {
91 // 99 //
92 // World 100 // World
93 // 101 //
94 102
95 G4double expHall_x = 1.0 * m; << 103 G4double expHall_x = 1.0*m;
96 G4double expHall_y = 1.0 * m; << 104 G4double expHall_y = 1.0*m;
97 G4double expHall_z = 1.0 * m; << 105 G4double expHall_z = 1.0*m;
98 << 106
99 auto solidWorld = new G4Box("World", // its << 107 fSolidWorld = new G4Box("World", //its name
100 expHall_x, expHa << 108 expHall_x,expHall_y,expHall_z); //its size
101 << 109
102 auto logicWorld = new G4LogicalVolume(solidW << 110 fLogicWorld = new G4LogicalVolume(fSolidWorld, //its solid
103 fVacuu << 111 fVacuum, //its material
104 "World << 112 "World"); //its name
105 << 113
106 auto physiWorld = new G4PVPlacement(nullptr, << 114 fPhysiWorld = new G4PVPlacement(0, //no rotation
107 G4ThreeV << 115 G4ThreeVector(), //at (0,0,0)
108 logicWor << 116 fLogicWorld, //its logical volume
109 "World", << 117 "World", //its name
110 nullptr, << 118 0, //its mother volume
111 false, << 119 false, //no boolean operation
112 0); // << 120 0); //copy number
113 121
114 G4double maxStep = 1.0 * mm; << 122 G4double maxStep = 1.0*mm;
115 G4double maxTime = 41. * s; << 123 G4double maxTime = 41.*s;
116 124
117 auto stepLimit = new G4UserLimits(maxStep, D << 125 G4UserLimits* stepLimit = new G4UserLimits(maxStep,DBL_MAX,maxTime);
118 <<
119 logicWorld->SetUserLimits(stepLimit); <<
120 126
>> 127 fLogicWorld->SetUserLimits(stepLimit);
>> 128
121 // 129 //
122 // Visualization attributes 130 // Visualization attributes
123 // 131 //
124 // logicWorld->SetVisAttributes (G4VisAttrib << 132 // fLogicWorld->SetVisAttributes (G4VisAttributes::GetInvisible());
125 133
126 // 134 //
127 // always return the physical World << 135 //always return the physical World
128 // 136 //
129 return physiWorld; << 137 return fPhysiWorld;
130 } 138 }
131 139
132 //....oooOO0OOooo........oooOO0OOooo........oo 140 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
133 141
134 G4ThreadLocal G4UniformGravityField* F06Detect << 142 G4ThreadLocal G4UniformGravityField* F06DetectorConstruction::fField = 0;
135 143
136 void F06DetectorConstruction::ConstructSDandFi 144 void F06DetectorConstruction::ConstructSDandField()
137 { 145 {
138 using StepperType = G4ClassicalRK4; 146 using StepperType = G4ClassicalRK4;
139 << 147
140 if (!fField) { 148 if (!fField) {
141 fField = new G4UniformGravityField(); <<
142 149
143 auto equation = new G4RepleteEofM(fField); << 150 fField = new G4UniformGravityField();
144 // G4EqGravityField* equation = new G4 <<
145 151
146 G4TransportationManager* transportMgr = G4 << 152 G4RepleteEofM* equation = new G4RepleteEofM(fField);
>> 153 // G4EqGravityField* equation = new G4EqGravityField(fField);
147 154
148 G4FieldManager* fieldManager = transportMg << 155 G4TransportationManager* transportMgr =
149 fieldManager->SetDetectorField(fField); << 156 G4TransportationManager::GetTransportationManager();
150 <<
151 const int nVar = 8; // 12 for RepleteEofM <<
152 auto stepper = new StepperType(equation, n <<
153 <<
154 G4double minStep = 0.01 * mm; <<
155 G4ChordFinder* chordFinder = nullptr; <<
156 if (stepper) { <<
157 auto intgrDriver = <<
158 new G4IntegrationDriver<StepperType>(m <<
159 if (intgrDriver) { <<
160 chordFinder = new G4ChordFinder(intgrD <<
161 } <<
162 } <<
163 <<
164 // OLD -- and wrong <<
165 // new G4ChordFinder((G4MagneticField*)fFi <<
166 <<
167 // Set accuracy parameters <<
168 G4double deltaChord = 3.0 * mm; <<
169 chordFinder->SetDeltaChord(deltaChord); <<
170 <<
171 G4double deltaIntersection = 0.1 * mm; <<
172 fieldManager->SetDeltaIntersection(deltaIn <<
173 <<
174 // Control accuracy of integration <<
175 // <<
176 G4double deltaOneStep = 0.01 * mm; <<
177 fieldManager->SetAccuraciesWithDeltaOneSte <<
178 // <<
179 G4double epsMax = 1.0e-4; // Pure number <<
180 G4double epsMin = 2.5e-7; // <<
181 fieldManager->SetMinimumEpsilonStep(epsMin <<
182 fieldManager->SetMaximumEpsilonStep(epsMax <<
183 // The acceptable relative accuracy is cal <<
184 // but bounded to the interval between <<
185 157
186 fieldManager->SetChordFinder(chordFinder); << 158 G4FieldManager* fieldManager= transportMgr->GetFieldManager();
>> 159 fieldManager->SetDetectorField(fField);
>> 160
>> 161 const int nVar= 8; // 12 for RepleteEofM
>> 162 StepperType* stepper = new StepperType(equation,nVar);
>> 163
>> 164 G4double minStep = 0.01*mm;
>> 165 G4ChordFinder* chordFinder = nullptr;
>> 166 if( stepper )
>> 167 {
>> 168 auto intgrDriver = new G4IntegrationDriver<StepperType>(minStep,
>> 169 stepper,
>> 170 stepper->GetNumberOfVariables());
>> 171 if( intgrDriver ){
>> 172 chordFinder = new G4ChordFinder(intgrDriver);
>> 173 }
>> 174 }
>> 175
>> 176 // OLD -- and wrong
>> 177 // new G4ChordFinder((G4MagneticField*)fField,minStep,stepper);
>> 178
>> 179 // Set accuracy parameters
>> 180 G4double deltaChord = 3.0*mm;
>> 181 chordFinder->SetDeltaChord( deltaChord );
>> 182
>> 183 G4double deltaIntersection = 0.1*mm;
>> 184 fieldManager->SetDeltaIntersection(deltaIntersection);
>> 185
>> 186 // Control accuracy of integration
>> 187 //
>> 188 G4double deltaOneStep = 0.01*mm;
>> 189 fieldManager->SetAccuraciesWithDeltaOneStep(deltaOneStep);
>> 190 //
>> 191 G4double epsMax = 1.0e-4; // Pure number -- maximum relative integration error
>> 192 G4double epsMin = 2.5e-7; //
>> 193 fieldManager->SetMinimumEpsilonStep(epsMin);
>> 194 fieldManager->SetMaximumEpsilonStep(epsMax);
>> 195 // The acceptable relative accuracy is calculated as deltaOneStep / stepsize
>> 196 // but bounded to the interval between these values!
>> 197
>> 198 fieldManager->SetChordFinder(chordFinder);
187 } 199 }
188 } 200 }
189 201
190 //....oooOO0OOooo........oooOO0OOooo........oo 202 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
191 203