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

Last change on this file was 1337, checked in by garnier, 14 years ago

tag geant4.9.4 beta 1 + modifs locales
File size: 13.5 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.14 2009/11/21 17:28:16 vnivanch Exp $
28// GEANT4 tag $Name: geant4-09-04-beta-01 $
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 "G4ProductionCuts.hh"
61#include "G4PhysicalVolumeStore.hh"
62#include "G4LogicalVolumeStore.hh"
63#include "G4SolidStore.hh"
64#include "G4NistManager.hh"
65
66#include "G4VisAttributes.hh"
67#include "G4Colour.hh"
68
69#include "G4UnitsTable.hh"
70#include "G4ios.hh"
71
72//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
73
74DetectorConstruction::DetectorConstruction()
75  : G4VUserDetectorConstruction()
76{
77  detectorMessenger = new DetectorMessenger(this);
78
79  ecalLength   = 36.*cm;
80  ecalWidth    = 6.*cm;
81  vertexLength = 3.*cm;
82  padLength    = 0.1*mm;
83  padWidth     = 0.02*mm;
84  absLength    = 2.*mm;
85  vertexRegion = 0;
86  muonRegion   = 0;
87  logicWorld   = 0;
88  logicCal     = 0;
89  logicA1      = 0;
90  DefineMaterials();
91  vertexDetectorCuts = new G4ProductionCuts();
92  muonDetectorCuts   = new G4ProductionCuts();
93}
94
95//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
96
97DetectorConstruction::~DetectorConstruction()
98{ 
99  delete detectorMessenger;
100  delete vertexDetectorCuts;
101  delete muonDetectorCuts;
102}
103
104//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
105
106G4VPhysicalVolume* DetectorConstruction::Construct()
107{
108  return ConstructVolumes();
109}
110
111//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
112
113void DetectorConstruction::DefineMaterials()
114{
115  // Default materials
116
117  G4NistManager* man = G4NistManager::Instance();
118  //  man->SetVerbose(1);
119  worldMaterial = man->FindOrBuildMaterial("G4_AIR");
120  absMaterial   = man->FindOrBuildMaterial("G4_Al");
121  vertMaterial  = man->FindOrBuildMaterial("G4_Si");
122  yorkMaterial  = man->FindOrBuildMaterial("G4_Fe");
123  calMaterial   = man->FindOrBuildMaterial("G4_CESIUM_IODIDE");
124}
125
126//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
127
128G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
129{
130  // Cleanup old geometry
131
132  G4GeometryManager::GetInstance()->OpenGeometry();
133
134  if(G4NistManager::Instance()->GetVerbose() > 0)
135    G4cout << *(G4Material::GetMaterialTable()) << G4endl;
136
137  if(vertexRegion) {
138    delete vertexRegion;
139    delete muonRegion;
140  }
141  vertexRegion = new G4Region("VertexDetector");
142  vertexRegion->SetProductionCuts(vertexDetectorCuts);
143
144  muonRegion   = new G4Region("MuonDetector");
145  muonRegion->SetProductionCuts(muonDetectorCuts);
146
147  G4SolidStore::GetInstance()->Clean();
148  G4LogicalVolumeStore::GetInstance()->Clean();
149  G4PhysicalVolumeStore::GetInstance()->Clean();
150
151  if(vertexLength < padLength*5.0) vertexLength = padLength*5.0;
152  G4double gap    = 0.01*mm;
153  G4double biggap = 2.*cm;
154  G4double york   = 10.*cm;
155
156  worldZ = 2.*vertexLength + 3.*absLength + 0.5*(ecalLength + york) + biggap*2.;
157
158  G4double worldX = ecalWidth*3.0;
159  G4double vertexZ= -worldZ + vertexLength*2.0 + absLength     + biggap;
160  G4double absZ2  = -worldZ + vertexLength*4.0 + absLength*3.5 + biggap;
161  G4double ecalZ  = -worldZ + vertexLength*4.0 + absLength*4.0 + ecalLength*0.5 
162    + 2.*biggap;
163  G4double yorkZ  = -worldZ + vertexLength*4.0 + absLength*5.0 + ecalLength
164    + york*0.5 + 3.*biggap;
165
166  //
167  // World
168  //
169  G4Box* solidW = new G4Box("World",worldX,worldX,worldZ);
170  logicWorld = new G4LogicalVolume( solidW,worldMaterial,"World");
171  G4VPhysicalVolume* world = new G4PVPlacement(0,G4ThreeVector(),
172                                               "World",logicWorld,0,false,0);
173
174  //
175  // Ecal
176  //
177  G4Box* solidE = new G4Box("VolE",worldX,worldX,ecalLength*0.5 + gap);
178  logicECal = new G4LogicalVolume( solidE,worldMaterial,"VolE");
179  G4VPhysicalVolume* physE = new G4PVPlacement(0,G4ThreeVector(0.,0.,ecalZ),
180                                               "VolE",logicECal,world,false,0);
181
182  G4Box* solidC = new G4Box("Ecal",ecalWidth*0.5,ecalWidth*0.5,ecalLength*0.5);
183  logicCal = new G4LogicalVolume( solidC,calMaterial,"Ecal");
184
185  G4cout << "Ecal is " << G4BestUnit(ecalLength,"Length")
186         << " of " << calMaterial->GetName() << G4endl;
187
188  // Crystals
189
190  G4double x0 = -(ecalWidth + gap)*2.0;
191  G4double y  = x0;
192  G4double x;
193  G4int k = 0;
194  G4VPhysicalVolume* pv;
195  G4int i,j;
196
197  for (i=0; i<5; i++) {
198    x  = x0;
199    for (j=0; j<5; j++) {
200
201      pv = new G4PVPlacement(0,G4ThreeVector(x,y,0.),"Ecal",logicCal,
202                                    physE,false,k);
203      k++;
204      x += ecalWidth + gap;
205    }
206    y += ecalWidth + gap;
207  }
208
209  //Absorber
210
211  G4Box* solidA = new G4Box("Abso",worldX,worldX,absLength*0.5);
212  logicA2 = new G4LogicalVolume( solidA,absMaterial,"Abs2");
213  pv = new G4PVPlacement(0,G4ThreeVector(0.,0.,absZ2),
214                         "Abs2",logicA2,world,false,0);
215
216  G4cout << "Absorber is " << G4BestUnit(absLength,"Length")
217         << " of " << absMaterial->GetName() << G4endl;
218
219  //York
220
221  G4Box* solidYV = new G4Box("VolY",worldX,worldX,york*0.5+absLength);
222  logicYV = new G4LogicalVolume( solidYV,yorkMaterial,"VolY");
223  G4VPhysicalVolume* physYV = new G4PVPlacement(0,G4ThreeVector(0.,0.,yorkZ),
224                                                "VolY",logicYV,world,false,0);
225
226  G4Box* solidY = new G4Box("York",worldX,worldX,york*0.5);
227  logicY = new G4LogicalVolume( solidY,yorkMaterial,"York");
228  pv = new G4PVPlacement(0,G4ThreeVector(),
229                         "York",logicY,physYV,false,0);
230
231  logicA3 = new G4LogicalVolume( solidA,absMaterial,"Abs3");
232  logicA4 = new G4LogicalVolume( solidA,absMaterial,"Abs4");
233
234  pv = new G4PVPlacement(0,G4ThreeVector(0.,0.,-(york+absLength)*0.5),
235                         "Abs3",logicA3,physYV,false,0);
236  pv = new G4PVPlacement(0,G4ThreeVector(0.,0.,(york+absLength)*0.5),
237                         "Abs4",logicA4,physYV,false,0);
238
239  //Vertex volume
240  G4Box* solidVV = new G4Box("VolV",worldX,worldX,vertexLength*2.+absLength+gap);
241  logicVV = new G4LogicalVolume( solidVV,worldMaterial,"VolV");
242  G4VPhysicalVolume* physVV = new G4PVPlacement(0,G4ThreeVector(0.,0.,vertexZ),
243                                                "VolV",logicVV,world,false,0);
244
245  //Absorber
246  logicA1 = new G4LogicalVolume( solidA,absMaterial,"Abs1");
247  pv = new G4PVPlacement(0,G4ThreeVector(0.,0.,vertexLength*2.-absLength*0.5),
248                         "Abs1",logicA1,physVV,false,0);
249
250  //Vertex
251  G4double vertWidth = ecalWidth/5.;
252  G4int npads = (G4int)(vertWidth/padWidth);
253  //G4cout << " vertWidth= " << vertWidth << " padWidth= " << padWidth
254  //     << " npads= " << npads << G4endl;
255  // insure beam to hit a middle of central pad
256  npads = (npads/2)*2 + 1;
257  x0 = -0.5*(padWidth + vertWidth);
258  G4double x1 = 0.5*vertWidth + gap; 
259  G4double z  = -(vertexLength+absLength);
260
261  G4Box* solidVD = new G4Box("VertDet",x1,ecalWidth*0.5+gap,padLength*0.5);
262  logicVD = new G4LogicalVolume( solidVD,vertMaterial,"VertDet");
263  logicVD->SetSolid(solidVD);
264
265  G4Box* solidV = new G4Box("Vert",padWidth*0.5,ecalWidth*0.5,padLength*0.5);
266  logicV = new G4LogicalVolume( solidV,vertMaterial,"Vert");
267
268  for (i=0; i<3; i++) {
269    pv = new G4PVPlacement(0,G4ThreeVector(0.,0.,z),"VertDet",logicVD,
270                                    physVV,false,i);
271    z += vertexLength;
272  }
273  x = x0;
274
275  for (j=0; j<npads; j++) {
276
277    pv = new G4PVPlacement(0,G4ThreeVector(x,0.,0.),logicV,"Vert",logicVD,
278                                    false,k);
279    x += padWidth;
280  }
281
282  G4cout << "Vertex is " << G4BestUnit(vertexLength,"Length")
283         << " of 3 layers of Si of " << G4BestUnit(padLength,"Length")
284         << " npads= " << npads
285         << G4endl;
286
287  // Define region for the vertex detector
288  vertexRegion->AddRootLogicalVolume(logicVV);
289  vertexRegion->AddRootLogicalVolume(logicA3);
290
291  // Define region for the muon detector
292  muonRegion->AddRootLogicalVolume(logicYV);
293
294  // color regions
295  logicVV-> SetVisAttributes(G4VisAttributes::Invisible);
296  logicV-> SetVisAttributes(G4VisAttributes::Invisible);
297  logicECal-> SetVisAttributes(G4VisAttributes::Invisible);
298  logicYV-> SetVisAttributes(G4VisAttributes::Invisible);
299
300  G4VisAttributes* regWcolor = new G4VisAttributes(G4Colour(0.3, 0.3, 0.3));
301  logicWorld->SetVisAttributes(regWcolor);
302
303  G4VisAttributes* regVcolor = new G4VisAttributes(G4Colour(0., 0.3, 0.7));
304  logicVD->SetVisAttributes(regVcolor);
305
306  G4VisAttributes* regCcolor = new G4VisAttributes(G4Colour(0., 0.7, 0.3));
307  logicCal->SetVisAttributes(regCcolor);
308
309  G4VisAttributes* regAcolor = new G4VisAttributes(G4Colour(1., 0.5, 0.5));
310  logicA1->SetVisAttributes(regAcolor);
311  logicA2->SetVisAttributes(regAcolor);
312  logicA3->SetVisAttributes(regAcolor);
313  logicA4->SetVisAttributes(regAcolor);
314
315  G4VisAttributes* regMcolor = new G4VisAttributes(G4Colour(1., 1., 0.));
316  logicY->SetVisAttributes(regMcolor);
317
318  // always return world
319  G4cout << "### New geometry is constructed" << G4endl;
320 
321  return world;
322}
323
324//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
325
326void DetectorConstruction::SetEcalMaterial(const G4String& mat)
327{
328  // search the material by its name
329  G4Material* pttoMaterial = 
330    G4NistManager::Instance()->FindOrBuildMaterial(mat, false);
331  if (pttoMaterial) {
332    calMaterial = pttoMaterial;
333    if(logicCal) {
334      logicCal->SetMaterial(calMaterial);
335      G4RunManager::GetRunManager()->PhysicsHasBeenModified();
336    }
337  }
338}
339
340//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
341
342void DetectorConstruction::SetAbsMaterial(const G4String& mat)
343{
344  // search the material by its name
345  G4Material* pttoMaterial = 
346    G4NistManager::Instance()->FindOrBuildMaterial(mat, false);
347  if (pttoMaterial) {
348    absMaterial = pttoMaterial;
349    if(logicA1) {
350      logicA1->SetMaterial(absMaterial);
351      logicA2->SetMaterial(absMaterial);
352      logicA3->SetMaterial(absMaterial);
353      logicA4->SetMaterial(absMaterial);
354      G4RunManager::GetRunManager()->PhysicsHasBeenModified();
355    }
356  }
357}
358
359//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
360
361void DetectorConstruction::UpdateGeometry()
362{
363  G4RunManager::GetRunManager()->PhysicsHasBeenModified();
364  G4RunManager::GetRunManager()->DefineWorldVolume(ConstructVolumes());
365}
366
367//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
368
369void DetectorConstruction::SetEcalLength (G4double val)   
370{
371  if(val > 0.0) {
372    ecalLength = val;
373    G4RunManager::GetRunManager()->GeometryHasBeenModified();
374  }
375}
376
377//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
378
379void DetectorConstruction::SetEcalWidth  (G4double val)   
380{
381  if(val > 0.0) {
382    ecalWidth = val;
383    G4RunManager::GetRunManager()->GeometryHasBeenModified();
384  }
385}
386
387//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
388
389void DetectorConstruction::SetVertexLength (G4double val) 
390{
391  if(val > 0.0) {
392    vertexLength = val;
393    G4RunManager::GetRunManager()->GeometryHasBeenModified();
394  }
395}
396
397//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
398
399void DetectorConstruction::SetPadLength  (G4double val)   
400{
401  if(val > 0.0) {
402    padLength = val;
403    G4RunManager::GetRunManager()->GeometryHasBeenModified();
404  }
405}
406
407//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
408
409void DetectorConstruction::SetPadWidth  (G4double val)   
410{
411  if(val > 0.0) {
412    padWidth = val;
413    G4RunManager::GetRunManager()->GeometryHasBeenModified();
414  }
415}
416
417//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
418
419void DetectorConstruction::SetAbsLength(G4double val)     
420{
421  if(val > 0.0) {
422    absLength = val;
423    G4RunManager::GetRunManager()->GeometryHasBeenModified();
424  }
425}
426
427//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Note: See TracBrowser for help on using the repository browser.