source: trunk/examples/extended/electromagnetic/TestEm9/src/DetectorConstruction.cc@ 1036

Last change on this file since 1036 was 807, checked in by garnier, 17 years ago

update
File size: 11.6 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27// $Id: DetectorConstruction.cc,v 1.9 2006/06/29 17:02:59 gunter Exp $
28// GEANT4 tag $Name: $
29//
30//
31/////////////////////////////////////////////////////////////////////////
32//
33// TestEm9: Crystal calorimeter
34//
35// Created: 31.01.03 V.Ivanchenko
36//
37// Modified:
38//
39////////////////////////////////////////////////////////////////////////
40//
41
42
43//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
45
46#include "DetectorConstruction.hh"
47#include "DetectorMessenger.hh"
48
49#include "G4Box.hh"
50#include "G4LogicalVolume.hh"
51#include "G4PVPlacement.hh"
52#include "G4PVReplica.hh"
53
54#include "G4TransportationManager.hh"
55
56#include "G4GeometryManager.hh"
57#include "G4RunManager.hh"
58#include "G4Region.hh"
59#include "G4RegionStore.hh"
60#include "G4PhysicalVolumeStore.hh"
61#include "G4LogicalVolumeStore.hh"
62#include "G4SolidStore.hh"
63#include "G4NistManager.hh"
64
65#include "G4VisAttributes.hh"
66#include "G4Colour.hh"
67
68#include "G4UnitsTable.hh"
69#include "G4ios.hh"
70
71//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
72
73DetectorConstruction::DetectorConstruction()
74 :logicC(0),logicA1(0),logicA2(0)
75{
76 detectorMessenger = new DetectorMessenger(this);
77
78 ecalLength = 36.*cm;
79 ecalWidth = 6.*cm;
80 vertexLength = 3.*cm;
81 padLength = 0.1*mm;
82 padWidth = 0.02*mm;
83 absLength = 2.*mm;
84 vertexRegion = 0;
85 muonRegion = 0;
86 DefineMaterials();
87}
88
89//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
90
91DetectorConstruction::~DetectorConstruction()
92{ delete detectorMessenger;}
93
94//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
95
96G4VPhysicalVolume* DetectorConstruction::Construct()
97{
98 return ConstructVolumes();
99}
100
101//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
102
103void DetectorConstruction::DefineMaterials()
104{
105 // Default materials
106
107 G4NistManager* man = G4NistManager::Instance();
108 man->SetVerbose(1);
109 worldMaterial = man->FindOrBuildMaterial("G4_AIR");
110 absMaterial = man->FindOrBuildMaterial("G4_Al");
111 vertMaterial = man->FindOrBuildMaterial("G4_Si");
112 yorkMaterial = man->FindOrBuildMaterial("G4_Fe");
113 calMaterial = man->FindOrBuildMaterial("G4_CESIUM_IODIDE");
114}
115
116//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
117
118G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
119{
120 G4cout << *(G4Material::GetMaterialTable()) << G4endl;
121
122 // Cleanup old geometry
123
124 G4GeometryManager::GetInstance()->OpenGeometry();
125 G4PhysicalVolumeStore::GetInstance()->Clean();
126 G4LogicalVolumeStore::GetInstance()->Clean();
127 G4SolidStore::GetInstance()->Clean();
128
129 if(vertexRegion) delete vertexRegion;
130 if(muonRegion) delete muonRegion;
131 vertexRegion = new G4Region("VertexDetector");
132 muonRegion = new G4Region("MuonDetector");
133
134 if(vertexLength < padLength*5.0) vertexLength = padLength*5.0;
135 G4double gap = 0.01*mm;
136 G4double biggap = 2.*cm;
137 G4double york = 10.*cm;
138
139 worldZ = 2.*vertexLength + 3.*absLength + 0.5*(ecalLength + york) + biggap*2.;
140 G4double worldX = ecalWidth*3.0;
141 G4double vertexZ= -worldZ + vertexLength*2.0 + absLength + biggap;
142 G4double absZ2 = -worldZ + vertexLength*4.0 + absLength*3.5 + biggap;
143 G4double ecalZ = -worldZ + vertexLength*4.0 + absLength*4.0 + ecalLength*0.5 + 2.*biggap;
144 G4double yorkZ = -worldZ + vertexLength*4.0 + absLength*5.0 + ecalLength
145 + york*0.5 + 3.*biggap;
146
147 //
148 // World
149 //
150 G4Box* solidW = new G4Box("World",worldX,worldX,worldZ);
151 G4LogicalVolume* logicW = new G4LogicalVolume( solidW,worldMaterial,
152 "World");
153 G4VPhysicalVolume* world = new G4PVPlacement(0,G4ThreeVector(),
154 "World",logicW,0,false,0);
155
156 //
157 // Ecal
158 //
159 G4Box* solidE = new G4Box("VolE",worldX,worldX,ecalLength*0.5 + gap);
160 G4LogicalVolume* logicE = new G4LogicalVolume( solidE,worldMaterial,
161 "VolE");
162 G4VPhysicalVolume* physE = new G4PVPlacement(0,G4ThreeVector(0.,0.,ecalZ),
163 "VolE",logicE,world,false,0);
164
165 G4Box* solidC = new G4Box("Ecal",ecalWidth*0.5,ecalWidth*0.5,ecalLength*0.5);
166 logicC = new G4LogicalVolume( solidC,calMaterial,"Ecal");
167
168 G4cout << "Ecal is " << G4BestUnit(ecalLength,"Length")
169 << " of " << calMaterial->GetName() << G4endl;
170
171 // Crystals
172
173 G4double x0 = -(ecalWidth + gap)*2.0;
174 G4double y = x0;
175 G4double x;
176 G4int k = 0;
177 G4VPhysicalVolume* pv;
178 G4int i,j;
179
180 for (i=0; i<5; i++) {
181 x = x0;
182 for (j=0; j<5; j++) {
183
184 pv = new G4PVPlacement(0,G4ThreeVector(x,y,0.),"Ecal",logicC,
185 physE,false,k);
186 k++;
187 x += ecalWidth + gap;
188 }
189 y += ecalWidth + gap;
190 }
191
192 //Absorber
193
194 G4Box* solidA = new G4Box("Abso",worldX,worldX,absLength*0.5);
195 logicA2 = new G4LogicalVolume( solidA,absMaterial,"Abs2");
196 pv = new G4PVPlacement(0,G4ThreeVector(0.,0.,absZ2),
197 "Abs2",logicA2,world,false,0);
198
199 G4cout << "Absorber is " << G4BestUnit(absLength,"Length")
200 << " of " << absMaterial->GetName() << G4endl;
201
202 //York
203
204 G4Box* solidYV = new G4Box("VolY",worldX,worldX,york*0.5+absLength);
205 G4LogicalVolume* logicYV = new G4LogicalVolume( solidYV,yorkMaterial,"VolY");
206 G4VPhysicalVolume* physYV = new G4PVPlacement(0,G4ThreeVector(0.,0.,yorkZ),
207 "VolY",logicYV,world,false,0);
208
209 G4Box* solidY = new G4Box("York",worldX,worldX,york*0.5);
210 G4LogicalVolume* logicY = new G4LogicalVolume( solidY,yorkMaterial,"York");
211 pv = new G4PVPlacement(0,G4ThreeVector(),
212 "York",logicY,physYV,false,0);
213
214 logicA3 = new G4LogicalVolume( solidA,absMaterial,"Abs3");
215 logicA4 = new G4LogicalVolume( solidA,absMaterial,"Abs4");
216 pv = new G4PVPlacement(0,G4ThreeVector(0.,0.,-(york+absLength)*0.5),
217 "Abs3",logicA3,physYV,false,0);
218 pv = new G4PVPlacement(0,G4ThreeVector(0.,0.,(york+absLength)*0.5),
219 "Abs4",logicA4,physYV,false,0);
220
221 //Vertex volume
222
223 G4Box* solidVV = new G4Box("VolV",worldX,worldX,vertexLength*2.+absLength+gap);
224 G4LogicalVolume* logicVV = new G4LogicalVolume( solidVV,worldMaterial,"VolV");
225 G4VPhysicalVolume* physVV = new G4PVPlacement(0,G4ThreeVector(0.,0.,vertexZ),
226 "VolV",logicVV,world,false,0);
227
228 //Absorber
229
230 logicA1 = new G4LogicalVolume( solidA,absMaterial,"Abs1");
231 pv = new G4PVPlacement(0,G4ThreeVector(0.,0.,vertexLength*2.-absLength*0.5),
232 "Abs1",logicA1,physVV,false,0);
233
234 //Vertex
235
236 G4double vertWidth = ecalWidth/5.;
237 G4int npads = (G4int)(vertWidth/padWidth);
238 npads = (npads/2)*2 + 1;
239 x0 = -0.5*padWidth*((G4double)(npads-1));
240 G4double x1 = std::fabs(x0) + 0.5*padWidth + gap;
241 G4double z = -(vertexLength+absLength);
242
243 G4Box* solidVD = new G4Box("VertDet",x1,ecalWidth*0.5+gap,padLength*0.5);
244 G4LogicalVolume* logicVD = new G4LogicalVolume( solidVD,vertMaterial,"VertDet");
245
246 G4Box* solidV = new G4Box("Vert",padWidth*0.5,ecalWidth*0.5,padLength*0.5);
247 G4LogicalVolume* logicV = new G4LogicalVolume( solidV,vertMaterial,"Vert");
248
249 for (i=0; i<3; i++) {
250 pv = new G4PVPlacement(0,G4ThreeVector(0.,0.,z),"VertDet",logicVD,
251 physVV,false,i);
252 z += vertexLength;
253 }
254 x = x0;
255
256 for (j=0; j<npads; j++) {
257
258 pv = new G4PVPlacement(0,G4ThreeVector(x,0.,0.),logicV,"Vert",logicVD,
259 false,k);
260 x += padWidth;
261 }
262
263 G4cout << "Vertex is " << G4BestUnit(vertexLength,"Length")
264 << " of 3 layers of Si of " << G4BestUnit(padLength,"Length")
265 << " npads= " << npads
266 << G4endl;
267
268 // Define region for the vertex detector
269
270 vertexRegion->AddRootLogicalVolume(logicVV);
271 vertexRegion->AddRootLogicalVolume(logicA3);
272
273 // Define region for the muon detector
274
275 muonRegion->AddRootLogicalVolume(logicYV);
276
277 // color regions
278
279 logicVV-> SetVisAttributes(G4VisAttributes::Invisible);
280 logicV-> SetVisAttributes(G4VisAttributes::Invisible);
281 logicE-> SetVisAttributes(G4VisAttributes::Invisible);
282 logicYV-> SetVisAttributes(G4VisAttributes::Invisible);
283
284 G4VisAttributes* regWcolor = new G4VisAttributes(G4Colour(0.3, 0.3, 0.3));
285 logicW->SetVisAttributes(regWcolor);
286
287 G4VisAttributes* regVcolor = new G4VisAttributes(G4Colour(0., 0.3, 0.7));
288 logicVD->SetVisAttributes(regVcolor);
289
290 G4VisAttributes* regCcolor = new G4VisAttributes(G4Colour(0., 0.7, 0.3));
291 logicC->SetVisAttributes(regCcolor);
292
293 G4VisAttributes* regAcolor = new G4VisAttributes(G4Colour(1., 0.5, 0.5));
294 logicA1->SetVisAttributes(regAcolor);
295 logicA2->SetVisAttributes(regAcolor);
296 logicA3->SetVisAttributes(regAcolor);
297 logicA4->SetVisAttributes(regAcolor);
298
299 G4VisAttributes* regMcolor = new G4VisAttributes(G4Colour(1., 1., 0.));
300 logicY->SetVisAttributes(regMcolor);
301
302 // always return world
303 G4cout << "### New geometry is constructed" << G4endl;
304
305 return world;
306}
307
308//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
309
310void DetectorConstruction::SetEcalMaterial(const G4String& mat)
311{
312 // search the material by its name
313 G4Material* pttoMaterial =
314 G4NistManager::Instance()->FindOrBuildMaterial(mat, false);
315 if (pttoMaterial) calMaterial = pttoMaterial;
316}
317
318//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
319
320void DetectorConstruction::SetAbsMaterial(const G4String& mat)
321{
322 // search the material by its name
323 G4Material* pttoMaterial =
324 G4NistManager::Instance()->FindOrBuildMaterial(mat, false);
325 if (pttoMaterial) absMaterial = pttoMaterial;
326}
327
328//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
329
330void DetectorConstruction::UpdateGeometry()
331{
332 G4RunManager::GetRunManager()->DefineWorldVolume(ConstructVolumes());
333}
334
335//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Note: See TracBrowser for help on using the repository browser.