| 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 /// \file medical/fanoCavity2/src/DetectorCons << 26 // $Id: DetectorConstruction.cc,v 1.2 2007/11/05 13:19:16 maire Exp $
27 /// \brief Implementation of the DetectorConst << 27 // GEANT4 tag $Name: geant4-09-02-patch-04 $
28 // << 28
29 // 29 //
30 //....oooOO0OOooo........oooOO0OOooo........oo 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 32
33 #include "DetectorConstruction.hh" 33 #include "DetectorConstruction.hh"
34 <<
35 #include "DetectorMessenger.hh" 34 #include "DetectorMessenger.hh"
36 35
37 #include "G4GeometryManager.hh" <<
38 #include "G4LogicalVolume.hh" <<
39 #include "G4LogicalVolumeStore.hh" <<
40 #include "G4Material.hh" 36 #include "G4Material.hh"
41 #include "G4NistManager.hh" << 37 #include "G4Tubs.hh"
>> 38 #include "G4LogicalVolume.hh"
>> 39 #include "G4VPhysicalVolume.hh"
42 #include "G4PVPlacement.hh" 40 #include "G4PVPlacement.hh"
43 #include "G4PhysicalConstants.hh" << 41
>> 42 #include "G4GeometryManager.hh"
44 #include "G4PhysicalVolumeStore.hh" 43 #include "G4PhysicalVolumeStore.hh"
>> 44 #include "G4LogicalVolumeStore.hh"
45 #include "G4SolidStore.hh" 45 #include "G4SolidStore.hh"
46 #include "G4SystemOfUnits.hh" << 46
47 #include "G4Tubs.hh" <<
48 #include "G4UnitsTable.hh" 47 #include "G4UnitsTable.hh"
49 #include "G4VPhysicalVolume.hh" <<
50 48
51 //....oooOO0OOooo........oooOO0OOooo........oo 49 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
52 50
53 DetectorConstruction::DetectorConstruction() : << 51 DetectorConstruction::DetectorConstruction()
>> 52 :pWall(0), pCavity(0)
54 { 53 {
55 // default parameter values 54 // default parameter values
56 fWallThickness = 5 * cm; << 55 wallThickness = 5*cm;
57 fCavityThickness = 2 * mm; << 56 cavityThickness = 2*mm;
58 fCavityRadius = 10 * m; << 57 worldRadius = 10*m;
59 << 58
60 DefineMaterials(); 59 DefineMaterials();
61 SetWallMaterial("G4_WATER"); << 60 SetWallMaterial("Water");
62 SetCavityMaterial("g4Water_gas"); << 61
63 << 62 // create commands for interactive definition of the detector
64 // create commands for interactive definitio << 63 detectorMessenger = new DetectorMessenger(this);
65 fDetectorMessenger = new DetectorMessenger(t <<
66 } 64 }
67 65
68 //....oooOO0OOooo........oooOO0OOooo........oo 66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
69 67
70 DetectorConstruction::~DetectorConstruction() 68 DetectorConstruction::~DetectorConstruction()
>> 69 { delete detectorMessenger;}
>> 70
>> 71 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 72
>> 73 G4VPhysicalVolume* DetectorConstruction::Construct()
71 { 74 {
72 delete fDetectorMessenger; << 75 return ConstructVolumes();
73 } 76 }
74 77
75 //....oooOO0OOooo........oooOO0OOooo........oo 78 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
76 79
77 void DetectorConstruction::DefineMaterials() 80 void DetectorConstruction::DefineMaterials()
78 { << 81 {
79 G4double z, a; << 82 //
80 << 83 // define Elements
81 G4Element* H = new G4Element("Hydrogen", "H" << 84 //
82 G4Element* O = new G4Element("Oxygen", "O", << 85 G4double z,a;
>> 86
>> 87 G4Element* H = new G4Element("Hydrogen" ,"H" , z= 1., a= 1.01*g/mole);
>> 88 G4Element* O = new G4Element("Oxygen" ,"O" , z= 8., a= 16.00*g/mole);
83 89
84 G4Material* H2O = new G4Material("Water", 1. << 90 //
>> 91 // define materials
>> 92 //
>> 93 G4Material* H2O =
>> 94 new G4Material("Water", 1.0*g/cm3, 2);
85 H2O->AddElement(H, 2); 95 H2O->AddElement(H, 2);
86 H2O->AddElement(O, 1); 96 H2O->AddElement(O, 1);
87 H2O->GetIonisation()->SetMeanExcitationEnerg << 97 H2O->GetIonisation()->SetMeanExcitationEnergy(75.0*eV);
88 << 98
89 G4Material* gas = new G4Material("Water_gas" << 99 G4Material* gas =
>> 100 new G4Material("Water_gas", 1.0*mg/cm3, 2);
90 gas->AddElement(H, 2); 101 gas->AddElement(H, 2);
91 gas->AddElement(O, 1); 102 gas->AddElement(O, 1);
92 gas->GetIonisation()->SetMeanExcitationEnerg << 103 gas->GetIonisation()->SetMeanExcitationEnergy(75.0*eV);
93 << 104
94 new G4Material("Graphite", 6, 12.01 * g / mo << 105 new G4Material("Graphite", 6, 12.01*g/mole, 2.265*g/cm3);
95 new G4Material("Graphite_gas", 6, 12.01 * g << 106 new G4Material("Graphite_gas", 6, 12.01*g/mole, 2.265*mg/cm3);
96 << 107
97 new G4Material("Aluminium", 13, 26.98 * g / << 108 new G4Material("Aluminium", 13, 26.98*g/mole, 2.700*g/cm3);
98 new G4Material("Aluminium_gas", 13, 26.98 * << 109 new G4Material("Aluminium_gas", 13, 26.98*g/mole, 2.700*mg/cm3);
99 << 110
100 // alternatively, use G4 data base << 111 G4cout << *(G4Material::GetMaterialTable()) << G4endl;
101 // <<
102 G4NistManager* nist = G4NistManager::Instanc <<
103 <<
104 nist->FindOrBuildMaterial("G4_WATER"); <<
105 nist->BuildMaterialWithNewDensity("g4Water_g <<
106 <<
107 // printout <<
108 G4cout << *(G4Material::GetMaterialTable()) <<
109 } 112 }
110 113
111 //....oooOO0OOooo........oooOO0OOooo........oo 114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
112 << 115
113 G4VPhysicalVolume* DetectorConstruction::Const << 116 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
114 { << 117 {
115 if (fWall) { << 118 G4GeometryManager::GetInstance()->OpenGeometry();
116 return fWall; << 119 G4PhysicalVolumeStore::GetInstance()->Clean();
117 } << 120 G4LogicalVolumeStore::GetInstance()->Clean();
118 << 121 G4SolidStore::GetInstance()->Clean();
119 // Chamber << 122
>> 123 // Wall
120 // 124 //
121 fTotalThickness = fCavityThickness + 2 * fWa << 125 worldThickness = cavityThickness + 2*wallThickness;
122 fWallRadius = fCavityRadius; << 126
123 << 127 G4Tubs*
124 G4Tubs* sChamber = new G4Tubs("Chamber", // << 128 sWall = new G4Tubs("Wall", //name
125 0., fWallRadiu << 129 0.,worldRadius,0.5*worldThickness,0.,twopi); //size
126 << 130
127 G4LogicalVolume* lChamber = new G4LogicalVol << 131 G4LogicalVolume*
128 << 132 lWall = new G4LogicalVolume(sWall, //solid
129 << 133 wallMaterial, //material
130 << 134 "Wall"); //name
131 fWall = new G4PVPlacement(0, // no rotation << 135
132 G4ThreeVector(), << 136 pWall = new G4PVPlacement(0, //no rotation
133 lChamber, // logi << 137 G4ThreeVector(), //at (0,0,0)
134 "Wall", // name << 138 lWall, //logical volume
135 0, // mother vol << 139 "Wall", //name
136 false, // no bool << 140 0, //mother volume
137 0); // copy numbe << 141 false, //no boolean operation
>> 142 0); //copy number
138 143
139 // Cavity 144 // Cavity
140 // << 145 //
141 G4Tubs* sCavity = new G4Tubs("Cavity", 0., f << 146 G4Tubs*
142 << 147 sCavity = new G4Tubs("Cavity",
143 G4LogicalVolume* lCavity = new G4LogicalVolu << 148 0.,worldRadius,0.5*cavityThickness,0.,twopi);
144 << 149
145 << 150 G4LogicalVolume*
146 << 151 lCavity = new G4LogicalVolume(sCavity, //shape
147 fCavity = new G4PVPlacement(0, // no rotati << 152 cavityMaterial, //material
148 G4ThreeVector(), << 153 "Cavity"); //name
149 lCavity, // log << 154
150 "Cavity", // na << 155 pCavity = new G4PVPlacement(0, //no rotation
151 lChamber, // mo << 156 G4ThreeVector(), //at (0,0,0)
152 false, // no bo << 157 lCavity, //logical volume
153 1); // copy num << 158 "Cavity", //name
154 << 159 lWall, //mother volume
>> 160 false, //no boolean operation
>> 161 1); //copy number
>> 162
155 PrintParameters(); 163 PrintParameters();
156 << 164
157 // 165 //
158 // always return the root volume << 166 //always return the root volume
159 // << 167 //
160 return fWall; << 168 return pWall;
161 } 169 }
162 170
163 //....oooOO0OOooo........oooOO0OOooo........oo 171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
164 172
165 void DetectorConstruction::PrintParameters() 173 void DetectorConstruction::PrintParameters()
166 { 174 {
167 G4cout << "\n------------------------------- 175 G4cout << "\n---------------------------------------------------------\n";
168 G4cout << "---> The Wall is " << G4BestUnit( << 176 G4cout << "---> The Wall is " << G4BestUnit(wallThickness,"Length")
169 << fWallMaterial->GetName() << " ( " << 177 << " of " << wallMaterial->GetName() << " ( "
170 << G4BestUnit(fWallMaterial->GetDensi << 178 << G4BestUnit(wallMaterial->GetDensity(),"Volumic Mass") << " )\n";
171 G4cout << " The Cavity is " << G4BestUni << 179 G4cout << " The Cavity is " << G4BestUnit(cavityThickness,"Length")
172 << fCavityMaterial->GetName() << " ( << 180 << " of " << cavityMaterial->GetName() << " ( "
173 << G4BestUnit(fCavityMaterial->GetDen << 181 << G4BestUnit(cavityMaterial->GetDensity(),"Volumic Mass") << " )";
174 G4cout << "\n------------------------------- 182 G4cout << "\n---------------------------------------------------------\n";
175 G4cout << G4endl; 183 G4cout << G4endl;
176 } 184 }
177 << 185
178 //....oooOO0OOooo........oooOO0OOooo........oo 186 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
179 187
180 void DetectorConstruction::SetWallThickness(G4 << 188 void DetectorConstruction::SetWallMaterial(G4String materialChoice)
181 { 189 {
182 fWallThickness = value; << 190 // search the material by its name
>> 191 G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);
>> 192 if (pttoMaterial) wallMaterial = pttoMaterial;
>> 193
>> 194 pttoMaterial = G4Material::GetMaterial(materialChoice + "_gas");
>> 195 if (pttoMaterial) cavityMaterial = pttoMaterial;
183 } 196 }
184 197
185 //....oooOO0OOooo........oooOO0OOooo........oo 198 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
186 199
187 void DetectorConstruction::SetWallMaterial(con << 200 void DetectorConstruction::SetWallThickness(G4double value)
188 { 201 {
189 // search the material by its name << 202 wallThickness = value;
190 G4Material* pttoMaterial = G4Material::GetMa <<
191 if (pttoMaterial) fWallMaterial = pttoMateri <<
192 } 203 }
193 204
194 //....oooOO0OOooo........oooOO0OOooo........oo 205 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
195 206
196 void DetectorConstruction::SetCavityThickness( 207 void DetectorConstruction::SetCavityThickness(G4double value)
197 { 208 {
198 fCavityThickness = value; << 209 cavityThickness = value;
>> 210
199 } 211 }
200 <<
201 //....oooOO0OOooo........oooOO0OOooo........oo 212 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
202 213
203 void DetectorConstruction::SetCavityRadius(G4d << 214 void DetectorConstruction::SetWorldRadius(G4double value)
204 { 215 {
205 fCavityRadius = value; << 216 worldRadius = value;
206 } 217 }
207 << 218
208 //....oooOO0OOooo........oooOO0OOooo........oo 219 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
209 220
210 void DetectorConstruction::SetCavityMaterial(c << 221 #include "G4RunManager.hh"
>> 222
>> 223 void DetectorConstruction::UpdateGeometry()
211 { 224 {
212 // search the material by its name << 225 G4RunManager::GetRunManager()->DefineWorldVolume(ConstructVolumes());
213 G4Material* pttoMaterial = G4Material::GetMa <<
214 if (pttoMaterial) fCavityMaterial = pttoMate <<
215 } 226 }
216 227
217 //....oooOO0OOooo........oooOO0OOooo........oo 228 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
218 229