| 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 // $Id$
>> 34 //
36 /// \file medical/dna/slowing/src/DetectorCons 35 /// \file medical/dna/slowing/src/DetectorConstruction.cc
37 /// \brief Implementation of the DetectorConst 36 /// \brief Implementation of the DetectorConstruction class
38 37
39 #include "DetectorConstruction.hh" 38 #include "DetectorConstruction.hh"
40 <<
41 #include "DetectorMessenger.hh" 39 #include "DetectorMessenger.hh"
42 40
43 #include "G4NistManager.hh" <<
44 #include "G4RunManager.hh" <<
45 #include "G4SystemOfUnits.hh" 41 #include "G4SystemOfUnits.hh"
46 #include "G4UnitsTable.hh" <<
47 #include "G4UserLimits.hh" 42 #include "G4UserLimits.hh"
>> 43 #include "G4NistManager.hh"
>> 44 #include "G4UnitsTable.hh"
>> 45 #include "G4LogicalVolumeStore.hh"
>> 46 #include "G4RunManager.hh"
48 47
49 //....oooOO0OOooo........oooOO0OOooo........oo 48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
50 49
51 DetectorConstruction::DetectorConstruction() : << 50 DetectorConstruction::DetectorConstruction() :
>> 51 G4VUserDetectorConstruction(), fWaterMaterial(0)
52 { 52 {
53 // Create commands for interactive definitio << 53 // create commands for interactive definition of the detector
54 fDetectorMessenger = new DetectorMessenger(t 54 fDetectorMessenger = new DetectorMessenger(this);
55 55
56 // Default tracking cut << 56 //default tracking cut
57 fTrackingCut = 7.4 * eV; << 57 fTrackingCut = 7.4*eV;
58 } 58 }
59 59
60 //....oooOO0OOooo........oooOO0OOooo........oo 60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
61 61
62 DetectorConstruction::~DetectorConstruction() 62 DetectorConstruction::~DetectorConstruction()
63 { 63 {
64 delete fDetectorMessenger; 64 delete fDetectorMessenger;
65 } 65 }
66 66
67 //....oooOO0OOooo........oooOO0OOooo........oo 67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
68 68
69 G4VPhysicalVolume* DetectorConstruction::Const 69 G4VPhysicalVolume* DetectorConstruction::Construct()
>> 70
70 { 71 {
71 DefineMaterials(); 72 DefineMaterials();
72 return ConstructDetector(); 73 return ConstructDetector();
73 } 74 }
74 75
75 //....oooOO0OOooo........oooOO0OOooo........oo 76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
76 77
77 void DetectorConstruction::DefineMaterials() 78 void DetectorConstruction::DefineMaterials()
78 { 79 {
>> 80
79 // Water is defined from NIST material datab 81 // Water is defined from NIST material database
80 G4NistManager* man = G4NistManager::Instance << 82 G4NistManager * man = G4NistManager::Instance();
81 83
82 G4Material* H2O = man->FindOrBuildMaterial(" << 84 G4Material * H2O = man->FindOrBuildMaterial("G4_WATER");
83 85
>> 86 /*
>> 87 If one wishes to test other density value for water material,
>> 88 one should use instead:
>> 89 G4Material * H2O = man->BuildMaterialWithNewDensity("G4_WATER_MODIFIED",
>> 90 "G4_WATER",1.100*g/cm3);
>> 91
>> 92 Note: any string for "G4_WATER_MODIFIED" parameter is accepted
>> 93 and "G4_WATER" parameter should not be changed
>> 94 Both materials are created and can be selected from dna.mac
>> 95 */
84 fWaterMaterial = H2O; 96 fWaterMaterial = H2O;
85 97
>> 98 //G4cout << "-> Density of water material (g/cm3)="
>> 99 // << fWaterMaterial->GetDensity()/(g/cm/cm/cm) << G4endl;
>> 100
86 G4cout << *(G4Material::GetMaterialTable()) 101 G4cout << *(G4Material::GetMaterialTable()) << G4endl;
87 } 102 }
88 103
89 //....oooOO0OOooo........oooOO0OOooo........oo 104 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
90 105
91 G4VPhysicalVolume* DetectorConstruction::Const 106 G4VPhysicalVolume* DetectorConstruction::ConstructDetector()
92 { 107 {
>> 108
93 // WORLD VOLUME 109 // WORLD VOLUME
94 110
95 G4double worldSizeX = 1 * m; 111 G4double worldSizeX = 1 * m;
96 G4double worldSizeY = worldSizeX; 112 G4double worldSizeY = worldSizeX;
97 G4double worldSizeZ = worldSizeX; 113 G4double worldSizeZ = worldSizeX;
98 114
99 G4Box* solidWorld = new G4Box("World", // i << 115 G4Box* solidWorld = new G4Box("World", //its name
100 worldSizeX / 2 << 116 worldSizeX / 2, worldSizeY / 2, worldSizeZ / 2); //its size
101 117
102 fLogicWorld = new G4LogicalVolume(solidWorld << 118 fLogicWorld = new G4LogicalVolume(solidWorld, //its solid
103 fWaterMate << 119 fWaterMaterial, //its material
104 "World"); << 120 "World"); //its name
105 << 121
106 G4PVPlacement* physiWorld = new G4PVPlacemen << 122 G4PVPlacement* physiWorld = new G4PVPlacement(
107 << 123 0, //no rotation
108 << 124 G4ThreeVector(), //at (0,0,0)
109 << 125 fLogicWorld,
110 << 126 "World", //its name
111 << 127 0, //its mother volume
112 << 128 false, //no boolean operation
>> 129 0); //copy number
113 130
114 // Visualization attributes 131 // Visualization attributes
115 G4VisAttributes* worldVisAtt = new G4VisAttr << 132 G4VisAttributes* worldVisAtt = new G4VisAttributes(G4Colour(1.0, 1.0, 1.0)); //White
116 worldVisAtt->SetVisibility(true); 133 worldVisAtt->SetVisibility(true);
117 fLogicWorld->SetVisAttributes(worldVisAtt); 134 fLogicWorld->SetVisAttributes(worldVisAtt);
118 135
119 G4VisAttributes* worldVisAtt1 = new G4VisAtt 136 G4VisAttributes* worldVisAtt1 = new G4VisAttributes(G4Colour(1.0, 0.0, 0.0));
120 worldVisAtt1->SetVisibility(true); 137 worldVisAtt1->SetVisibility(true);
121 138
122 // Tracking cut 139 // Tracking cut
123 fLogicWorld->SetUserLimits(new G4UserLimits( << 140 fLogicWorld->SetUserLimits(new G4UserLimits(DBL_MAX,DBL_MAX,DBL_MAX,
>> 141 fTrackingCut));
124 142
125 PrintParameters(); 143 PrintParameters();
126 144
127 return physiWorld; 145 return physiWorld;
128 } 146 }
129 147
130 //....oooOO0OOooo........oooOO0OOooo........oo 148 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
131 149
132 void DetectorConstruction::SetMaterial(G4Strin 150 void DetectorConstruction::SetMaterial(G4String materialChoice)
133 { 151 {
134 // Search the material by its name << 152 // Search the material by its name
135 G4Material* pttoMaterial = G4Material::GetMa << 153 G4Material* pttoMaterial = G4NistManager::Instance()->FindOrBuildMaterial
136 if (pttoMaterial) fWaterMaterial = pttoMater << 154 (materialChoice);
137 G4RunManager::GetRunManager()->PhysicsHasBee << 155
>> 156 if (pttoMaterial)
>> 157 {
>> 158 fWaterMaterial = pttoMaterial;
>> 159 G4LogicalVolume* logicWorld =
>> 160 G4LogicalVolumeStore::GetInstance()->GetVolume("World");
>> 161 logicWorld->SetMaterial(fWaterMaterial);
>> 162 G4RunManager::GetRunManager()->GeometryHasBeenModified();
>> 163 }
>> 164 }
>> 165
>> 166 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 167
>> 168 void DetectorConstruction::UpdateGeometry()
>> 169 {
>> 170 G4RunManager::GetRunManager()->GeometryHasBeenModified();
>> 171 G4RunManager::GetRunManager()->DefineWorldVolume(ConstructDetector());
138 } 172 }
139 173
140 //....oooOO0OOooo........oooOO0OOooo........oo 174 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
141 175
142 void DetectorConstruction::SetTrackingCut(G4do 176 void DetectorConstruction::SetTrackingCut(G4double value)
143 { 177 {
144 fTrackingCut = value; 178 fTrackingCut = value;
145 G4RunManager::GetRunManager()->ReinitializeG << 179 G4RunManager::GetRunManager()->ReinitializeGeometry();
146 } 180 }
147 181
148 //....oooOO0OOooo........oooOO0OOooo........oo 182 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
149 183
150 void DetectorConstruction::PrintParameters() c 184 void DetectorConstruction::PrintParameters() const
151 { 185 {
152 G4cout << "\n------------------------------- 186 G4cout << "\n---------------------------------------------------------\n";
153 G4cout << "---> The tracking cut is set to " << 187 G4cout << "---> The tracking cut is set to "
>> 188 << G4BestUnit(fTrackingCut,"Energy") << G4endl;
154 G4cout << "\n------------------------------- 189 G4cout << "\n---------------------------------------------------------\n";
155 } 190 }
>> 191
>> 192 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
156 193