| 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 medical/electronScattering2/src/Elec
28 /// \brief Implementation of the ElectronBench
29
30 #include "ElectronBenchmarkDetector.hh"
31
32 #include "ElectronBenchmarkDetectorMessenger.h
33
34 #include "G4Colour.hh"
35 #include "G4GeometryManager.hh"
36 #include "G4LogicalVolume.hh"
37 #include "G4LogicalVolumeStore.hh"
38 #include "G4Material.hh"
39 #include "G4MultiFunctionalDetector.hh"
40 #include "G4NistManager.hh"
41 #include "G4PSCellFlux.hh"
42 #include "G4PSPopulation.hh"
43 #include "G4PVPlacement.hh"
44 #include "G4PVReplica.hh"
45 #include "G4PhysicalVolumeStore.hh"
46 #include "G4RunManager.hh"
47 #include "G4SDManager.hh"
48 #include "G4SDParticleFilter.hh"
49 #include "G4SolidStore.hh"
50 #include "G4SystemOfUnits.hh"
51 #include "G4Tubs.hh"
52 #include "G4UImanager.hh"
53 #include "G4VPrimitiveScorer.hh"
54 #include "G4VisAttributes.hh"
55
56 //....oooOO0OOooo........oooOO0OOooo........oo
57
58 ElectronBenchmarkDetector::ElectronBenchmarkDe
59 {
60 // Exit Window
61 fPosWindow0 = 0.000000 * cm;
62 fPosWindow1 = 0.004120 * cm;
63
64 // Scattering Foil
65 fPosPrimFoil = 2.650000 * cm;
66 fHalfThicknessPrimFoil = 0.0 * cm;
67
68 // Monitor Chamber
69 fPosMon0 = 5.000000 * cm;
70 fPosMon1 = 5.011270 * cm;
71
72 // Helium Bag
73 fPosBag0 = 6.497500 * cm;
74 fPosHelium0 = 6.500000 * cm;
75 fPosHelium1 = 116.500000 * cm;
76 fPosBag1 = 116.502500 * cm;
77 fThicknessRing = 1.4 * cm;
78
79 // Scoring Plane
80 fPosScorer = 118.200000 * cm;
81 fThicknessScorer = 0.001 * cm;
82 fWidthScorerRing = 0.1 * cm;
83
84 // Radii
85 fRadOverall = 23.3 * cm;
86 fRadRingInner = 20.0 * cm;
87
88 // Extra space remaining in world volume aro
89 fPosDelta = 1. * cm;
90 fRadDelta = 0.1 * cm;
91
92 fMessenger = new ElectronBenchmarkDetectorMe
93 DefineMaterials();
94 }
95
96 //....oooOO0OOooo........oooOO0OOooo........oo
97
98 ElectronBenchmarkDetector::~ElectronBenchmarkD
99 {
100 delete fMessenger;
101
102 delete fWorldVisAtt;
103 delete fWindowVisAtt;
104 delete fPrimFoilVisAtt;
105 delete fMonVisAtt;
106 delete fBagVisAtt;
107 delete fHeliumVisAtt;
108 delete fRingVisAtt;
109 delete fScorerVisAtt;
110 }
111
112 //....oooOO0OOooo........oooOO0OOooo........oo
113
114 G4VPhysicalVolume* ElectronBenchmarkDetector::
115 {
116 return CreateGeometry();
117 }
118
119 //....oooOO0OOooo........oooOO0OOooo........oo
120
121 void ElectronBenchmarkDetector::DefineMaterial
122 {
123 // Use NIST database for elements and materi
124 G4NistManager* man = G4NistManager::Instance
125 man->SetVerbose(1);
126
127 // Take all elements and materials from NIST
128 man->FindOrBuildMaterial("G4_He");
129 man->FindOrBuildMaterial("G4_Be");
130 man->FindOrBuildMaterial("G4_Al");
131 man->FindOrBuildMaterial("G4_Ti");
132 man->FindOrBuildMaterial("G4_Ta");
133 man->FindOrBuildMaterial("G4_AIR");
134 man->FindOrBuildMaterial("G4_MYLAR");
135
136 G4Element* C = man->FindOrBuildElement("C");
137 G4Element* Cu = man->FindOrBuildElement("Cu"
138 G4Element* Au = man->FindOrBuildElement("Au"
139 G4Element* Ti = man->FindOrBuildElement("Ti"
140 G4Element* Al = man->FindOrBuildElement("Al"
141 G4Element* V = man->FindOrBuildElement("V");
142
143 // Define materials not in NIST.
144 // While the NIST database does contain defa
145 // those defaults have different densities t
146 // benchmark specification.
147 G4double density;
148 G4int ncomponents;
149 G4double fractionmass;
150
151 G4Material* G4_C = new G4Material("G4_C", de
152 G4_C->AddElement(C, fractionmass = 1.00);
153
154 G4Material* G4_Cu = new G4Material("G4_Cu",
155 G4_Cu->AddElement(Cu, fractionmass = 1.00);
156
157 G4Material* G4_Au = new G4Material("G4_Au",
158 G4_Au->AddElement(Au, fractionmass = 1.00);
159
160 G4Material* TiAlloy = new G4Material("TiAllo
161 TiAlloy->AddElement(Ti, fractionmass = 0.90)
162 TiAlloy->AddElement(Al, fractionmass = 0.06)
163 TiAlloy->AddElement(V, fractionmass = 0.04);
164
165 // Print materials table
166 G4cout << *(G4Material::GetMaterialTable())
167 }
168
169 //....oooOO0OOooo........oooOO0OOooo........oo
170
171 G4VPhysicalVolume* ElectronBenchmarkDetector::
172 {
173 if (fPhysiWorld) return fPhysiWorld;
174
175 // Instantiate the world
176 fPhysiWorld = CreateWorld();
177 fLogWorld = fPhysiWorld->GetLogicalVolume();
178
179 // Instantiate the geometry
180 CreateExitWindow(fLogWorld);
181 CreatePrimaryFoil(fLogWorld);
182 CreateMonitor(fLogWorld);
183 CreateHeliumBag(fLogWorld);
184
185 // Create the scorers
186 CreateScorer(fLogWorld);
187
188 return fPhysiWorld;
189 }
190
191 //....oooOO0OOooo........oooOO0OOooo........oo
192
193 G4VPhysicalVolume* ElectronBenchmarkDetector::
194 {
195 G4double halfLengthWorld = fPosScorer / 2. +
196 G4double radWorld = fRadOverall + fRadDelta;
197 G4VSolid* worldSolid =
198 new G4Tubs("WorldSolid", 0. * cm, radWorld
199 G4LogicalVolume* worldLog =
200 new G4LogicalVolume(worldSolid, G4Material
201
202 fWorldVisAtt = new G4VisAttributes(G4Colour(
203 worldLog->SetVisAttributes(fWorldVisAtt);
204
205 G4VPhysicalVolume* worldPhys =
206 new G4PVPlacement(0, G4ThreeVector(0., 0.,
207
208 return worldPhys;
209 }
210
211 //....oooOO0OOooo........oooOO0OOooo........oo
212
213 void ElectronBenchmarkDetector::CreateExitWind
214 {
215 G4double halfLengthWorld = fPosScorer / 2.;
216 G4double halfThicknessWindow = fPosWindow1 /
217 G4VSolid* windowSolid =
218 new G4Tubs("windowSolid", 0. * cm, fRadOve
219 G4LogicalVolume* windowLog =
220 new G4LogicalVolume(windowSolid, G4Materia
221
222 fWindowVisAtt = new G4VisAttributes(G4Colour
223 windowLog->SetVisAttributes(fWindowVisAtt);
224
225 new G4PVPlacement(0, G4ThreeVector(0., 0., h
226 "ExitWindow", worldLog, fa
227 }
228
229 //....oooOO0OOooo........oooOO0OOooo........oo
230
231 void ElectronBenchmarkDetector::CreatePrimaryF
232 {
233 G4double halfLengthWorld = fPosScorer / 2.;
234
235 // For some energies, we have no Primary Foi
236 if (fHalfThicknessPrimFoil == 0.) return;
237
238 fSolidPrimFoil =
239 new G4Tubs("PrimFoilSolid", 0. * cm, fRadO
240 fLogPrimFoil = new G4LogicalVolume(fSolidPri
241
242 fPrimFoilVisAtt = new G4VisAttributes(G4Colo
243 fLogPrimFoil->SetVisAttributes(fPrimFoilVisA
244
245 new G4PVPlacement(0,
246 G4ThreeVector(0., 0., fPos
247 fLogPrimFoil, "ScatteringF
248 }
249
250 //....oooOO0OOooo........oooOO0OOooo........oo
251
252 void ElectronBenchmarkDetector::CreateMonitor(
253 {
254 G4double halfLengthWorld = fPosScorer / 2.;
255 G4double halfThicknessMon = (fPosMon1 - fPos
256 G4VSolid* monSolid =
257 new G4Tubs("monSolid", 0. * cm, fRadOveral
258 G4LogicalVolume* monLog =
259 new G4LogicalVolume(monSolid, G4Material::
260
261 fMonVisAtt = new G4VisAttributes(G4Colour(0.
262 monLog->SetVisAttributes(fMonVisAtt);
263
264 new G4PVPlacement(0, G4ThreeVector(0., 0., f
265 "MonitorChamber", worldLog
266 }
267
268 //....oooOO0OOooo........oooOO0OOooo........oo
269
270 void ElectronBenchmarkDetector::CreateHeliumBa
271 {
272 G4double halfLengthWorld = fPosScorer / 2.;
273
274 // Construct cylinder of Mylar
275 G4double halfThicknessBag = (fPosBag1 - fPos
276 G4VSolid* bagSolid =
277 new G4Tubs("bagSolid", 0. * cm, fRadOveral
278 G4LogicalVolume* bagLog =
279 new G4LogicalVolume(bagSolid, G4Material::
280
281 fBagVisAtt = new G4VisAttributes(G4Colour(0.
282 bagLog->SetVisAttributes(fBagVisAtt);
283
284 new G4PVPlacement(0, G4ThreeVector(0., 0., f
285 "HeliumBag", worldLog, fal
286
287 // Insert cylinder of Helium into the Cylind
288 G4double halfThicknessHelium = (fPosHelium1
289 G4VSolid* heliumSolid =
290 new G4Tubs("heliumSolid", 0. * cm, fRadOve
291 G4LogicalVolume* heliumLog =
292 new G4LogicalVolume(heliumSolid, G4Materia
293
294 fHeliumVisAtt = new G4VisAttributes(G4Colour
295 heliumLog->SetVisAttributes(fHeliumVisAtt);
296
297 new G4PVPlacement(0, G4ThreeVector(0., 0., 0
298
299 // Insert two rings of Aluminum into the Cyl
300 G4double halfThicknessRing = fThicknessRing
301 G4VSolid* ringSolid =
302 new G4Tubs("ringSolid", fRadRingInner, fRa
303 G4LogicalVolume* ring0Log =
304 new G4LogicalVolume(ringSolid, G4Material:
305 G4LogicalVolume* ring1Log =
306 new G4LogicalVolume(ringSolid, G4Material:
307
308 fRingVisAtt = new G4VisAttributes(G4Colour(0
309 ring0Log->SetVisAttributes(fRingVisAtt);
310 ring1Log->SetVisAttributes(fRingVisAtt);
311
312 new G4PVPlacement(0, G4ThreeVector(0., 0., -
313 "Ring0", heliumLog, false,
314
315 new G4PVPlacement(0, G4ThreeVector(0., 0., h
316 "Ring1", heliumLog, false,
317 }
318
319 //....oooOO0OOooo........oooOO0OOooo........oo
320
321 void ElectronBenchmarkDetector::CreateScorer(G
322 {
323 G4double halfLengthWorld = fPosScorer / 2.;
324 G4double halfThicknessScorer = fThicknessSco
325
326 G4VSolid* scorerSolid =
327 new G4Tubs("scorerSolid", 0. * cm, fRadOve
328 G4LogicalVolume* scorerLog =
329 new G4LogicalVolume(scorerSolid, G4Materia
330
331 fScorerVisAtt = new G4VisAttributes(G4Colour
332 scorerLog->SetVisAttributes(fScorerVisAtt);
333 new G4PVPlacement(0, G4ThreeVector(0., 0., h
334 "Scorer", worldLog, false,
335
336 G4VSolid* scorerRingSolid =
337 new G4Tubs("scorerRingSolid", 0. * cm, fRa
338 fScorerRingLog =
339 new G4LogicalVolume(scorerRingSolid, G4Mat
340 new G4PVReplica("ScorerRing", fScorerRingLog
341 G4int(fRadOverall / fWidthSc
342 }
343
344 //....oooOO0OOooo........oooOO0OOooo........oo
345
346 // Note that this method is called both at sta
347 // any command causes a change to detector geo
348 void ElectronBenchmarkDetector::ConstructSDand
349 {
350 G4SDManager::GetSDMpointer()->SetVerboseLeve
351
352 // G4Cache mechanism is necessary for multi-
353 // as it allows us to store separate detecto
354 G4MultiFunctionalDetector*& sensitiveDetecto
355
356 if (!sensitiveDetector) {
357 sensitiveDetector = new G4MultiFunctionalD
358
359 G4VPrimitiveScorer* primitive;
360
361 G4SDParticleFilter* electronFilter = new G
362
363 primitive = new G4PSCellFlux("cell flux");
364 sensitiveDetector->RegisterPrimitive(primi
365
366 primitive = new G4PSCellFlux("e cell flux"
367 primitive->SetFilter(electronFilter);
368 sensitiveDetector->RegisterPrimitive(primi
369
370 primitive = new G4PSPopulation("population
371 sensitiveDetector->RegisterPrimitive(primi
372
373 primitive = new G4PSPopulation("e populati
374 primitive->SetFilter(electronFilter);
375 sensitiveDetector->RegisterPrimitive(primi
376 }
377 G4SDManager::GetSDMpointer()->AddNewDetector
378 fScorerRingLog->SetSensitiveDetector(sensiti
379 }
380
381 //....oooOO0OOooo........oooOO0OOooo........oo
382
383 void ElectronBenchmarkDetector::SetPrimFoilMat
384 {
385 G4Material* material = G4NistManager::Instan
386
387 if (material && material != fMaterialPrimFoi
388 fMaterialPrimFoil = material;
389 if (fLogPrimFoil) {
390 fLogPrimFoil->SetMaterial(fMaterialPrimF
391 }
392 G4RunManager::GetRunManager()->PhysicsHasB
393 }
394 }
395
396 //....oooOO0OOooo........oooOO0OOooo........oo
397
398 void ElectronBenchmarkDetector::SetPrimFoilThi
399 {
400 fHalfThicknessPrimFoil = thicknessPrimFoil /
401 }
402
403 //....oooOO0OOooo........oooOO0OOooo........oo
404