| 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 // This example is provided by the Geant4-DNA 26 // This example is provided by the Geant4-DNA collaboration
27 // Any report or published results obtained us << 27 // Any report or published results obtained using the Geant4-DNA software
28 // shall cite the following Geant4-DNA collabo 28 // shall cite the following Geant4-DNA collaboration publications:
29 // Med. Phys. 45 (2018) e722-e739 <<
30 // Phys. Med. 31 (2015) 861-874 29 // Phys. Med. 31 (2015) 861-874
31 // Med. Phys. 37 (2010) 4692-4708 30 // Med. Phys. 37 (2010) 4692-4708
32 // Int. J. Model. Simul. Sci. Comput. 1 (2010) <<
33 // <<
34 // The Geant4-DNA web site is available at htt 31 // The Geant4-DNA web site is available at http://geant4-dna.org
35 // 32 //
>> 33 //
36 /// \file medical/dna/slowing/src/DetectorCons 34 /// \file medical/dna/slowing/src/DetectorConstruction.cc
37 /// \brief Implementation of the DetectorConst 35 /// \brief Implementation of the DetectorConstruction class
38 36
39 #include "DetectorConstruction.hh" 37 #include "DetectorConstruction.hh"
40 <<
41 #include "DetectorMessenger.hh" 38 #include "DetectorMessenger.hh"
42 39
43 #include "G4NistManager.hh" <<
44 #include "G4RunManager.hh" <<
45 #include "G4SystemOfUnits.hh" 40 #include "G4SystemOfUnits.hh"
46 #include "G4UnitsTable.hh" <<
47 #include "G4UserLimits.hh" 41 #include "G4UserLimits.hh"
>> 42 #include "G4NistManager.hh"
>> 43 #include "G4UnitsTable.hh"
>> 44 #include "G4LogicalVolumeStore.hh"
>> 45 #include "G4RunManager.hh"
48 46
49 //....oooOO0OOooo........oooOO0OOooo........oo 47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
50 48
51 DetectorConstruction::DetectorConstruction() : << 49 DetectorConstruction::DetectorConstruction() :
>> 50 G4VUserDetectorConstruction(), fWaterMaterial(0)
52 { 51 {
53 // Create commands for interactive definitio << 52 // create commands for interactive definition of the detector
54 fDetectorMessenger = new DetectorMessenger(t 53 fDetectorMessenger = new DetectorMessenger(this);
55 54
56 // Default tracking cut << 55 //default tracking cut
57 fTrackingCut = 7.4 * eV; << 56 fTrackingCut = 7.4*eV;
58 } 57 }
59 58
60 //....oooOO0OOooo........oooOO0OOooo........oo 59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
61 60
62 DetectorConstruction::~DetectorConstruction() 61 DetectorConstruction::~DetectorConstruction()
63 { 62 {
64 delete fDetectorMessenger; 63 delete fDetectorMessenger;
65 } 64 }
66 65
67 //....oooOO0OOooo........oooOO0OOooo........oo 66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
68 67
69 G4VPhysicalVolume* DetectorConstruction::Const 68 G4VPhysicalVolume* DetectorConstruction::Construct()
>> 69
70 { 70 {
71 DefineMaterials(); 71 DefineMaterials();
72 return ConstructDetector(); 72 return ConstructDetector();
73 } 73 }
74 74
75 //....oooOO0OOooo........oooOO0OOooo........oo 75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
76 76
77 void DetectorConstruction::DefineMaterials() 77 void DetectorConstruction::DefineMaterials()
78 { 78 {
>> 79
79 // Water is defined from NIST material datab 80 // Water is defined from NIST material database
80 G4NistManager* man = G4NistManager::Instance << 81 G4NistManager * man = G4NistManager::Instance();
81 82
82 G4Material* H2O = man->FindOrBuildMaterial(" << 83 G4Material * H2O = man->FindOrBuildMaterial("G4_WATER");
83 84
84 fWaterMaterial = H2O; 85 fWaterMaterial = H2O;
85 86
86 G4cout << *(G4Material::GetMaterialTable()) 87 G4cout << *(G4Material::GetMaterialTable()) << G4endl;
87 } 88 }
88 89
89 //....oooOO0OOooo........oooOO0OOooo........oo 90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
90 91
91 G4VPhysicalVolume* DetectorConstruction::Const 92 G4VPhysicalVolume* DetectorConstruction::ConstructDetector()
92 { 93 {
>> 94
93 // WORLD VOLUME 95 // WORLD VOLUME
94 96
95 G4double worldSizeX = 1 * m; 97 G4double worldSizeX = 1 * m;
96 G4double worldSizeY = worldSizeX; 98 G4double worldSizeY = worldSizeX;
97 G4double worldSizeZ = worldSizeX; 99 G4double worldSizeZ = worldSizeX;
98 100
99 G4Box* solidWorld = new G4Box("World", // i << 101 G4Box* solidWorld = new G4Box("World", //its name
100 worldSizeX / 2 << 102 worldSizeX / 2, worldSizeY / 2, worldSizeZ / 2); //its size
101 103
102 fLogicWorld = new G4LogicalVolume(solidWorld << 104 fLogicWorld = new G4LogicalVolume(solidWorld, //its solid
103 fWaterMate << 105 fWaterMaterial, //its material
104 "World"); << 106 "World"); //its name
105 << 107
106 G4PVPlacement* physiWorld = new G4PVPlacemen << 108 G4PVPlacement* physiWorld = new G4PVPlacement(
107 << 109 0, //no rotation
108 << 110 G4ThreeVector(), //at (0,0,0)
109 << 111 fLogicWorld,
110 << 112 "World", //its name
111 << 113 0, //its mother volume
112 << 114 false, //no boolean operation
>> 115 0); //copy number
113 116
114 // Visualization attributes 117 // Visualization attributes
115 G4VisAttributes* worldVisAtt = new G4VisAttr 118 G4VisAttributes* worldVisAtt = new G4VisAttributes(G4Colour(1.0, 1.0, 1.0));
116 worldVisAtt->SetVisibility(true); 119 worldVisAtt->SetVisibility(true);
117 fLogicWorld->SetVisAttributes(worldVisAtt); 120 fLogicWorld->SetVisAttributes(worldVisAtt);
118 121
119 G4VisAttributes* worldVisAtt1 = new G4VisAtt 122 G4VisAttributes* worldVisAtt1 = new G4VisAttributes(G4Colour(1.0, 0.0, 0.0));
120 worldVisAtt1->SetVisibility(true); 123 worldVisAtt1->SetVisibility(true);
121 124
122 // Tracking cut 125 // Tracking cut
123 fLogicWorld->SetUserLimits(new G4UserLimits( << 126 fLogicWorld->SetUserLimits(new G4UserLimits(DBL_MAX,DBL_MAX,DBL_MAX,
>> 127 fTrackingCut));
124 128
125 PrintParameters(); 129 PrintParameters();
126 130
127 return physiWorld; 131 return physiWorld;
128 } 132 }
129 133
130 //....oooOO0OOooo........oooOO0OOooo........oo 134 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
131 135
132 void DetectorConstruction::SetMaterial(G4Strin 136 void DetectorConstruction::SetMaterial(G4String materialChoice)
133 { 137 {
134 // Search the material by its name << 138 // search the material by its name
135 G4Material* pttoMaterial = G4Material::GetMa << 139 G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);
136 if (pttoMaterial) fWaterMaterial = pttoMater 140 if (pttoMaterial) fWaterMaterial = pttoMaterial;
137 G4RunManager::GetRunManager()->PhysicsHasBee << 141 G4RunManager::GetRunManager()->PhysicsHasBeenModified();
138 } 142 }
139 143
140 //....oooOO0OOooo........oooOO0OOooo........oo 144 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
141 145
142 void DetectorConstruction::SetTrackingCut(G4do 146 void DetectorConstruction::SetTrackingCut(G4double value)
143 { 147 {
144 fTrackingCut = value; 148 fTrackingCut = value;
145 G4RunManager::GetRunManager()->ReinitializeG << 149 G4RunManager::GetRunManager()->ReinitializeGeometry();
146 } 150 }
147 151
148 //....oooOO0OOooo........oooOO0OOooo........oo 152 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
149 153
150 void DetectorConstruction::PrintParameters() c 154 void DetectorConstruction::PrintParameters() const
151 { 155 {
152 G4cout << "\n------------------------------- 156 G4cout << "\n---------------------------------------------------------\n";
153 G4cout << "---> The tracking cut is set to " << 157 G4cout << "---> The tracking cut is set to "
>> 158 << G4BestUnit(fTrackingCut,"Energy") << G4endl;
154 G4cout << "\n------------------------------- 159 G4cout << "\n---------------------------------------------------------\n";
155 } 160 }
>> 161
>> 162 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
156 163