source: trunk/examples/novice/N06/src/ExN06DetectorConstruction.cc @ 868

Last change on this file since 868 was 807, checked in by garnier, 16 years ago

update
File size: 10.0 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: ExN06DetectorConstruction.cc,v 1.15 2006/06/29 17:54:17 gunter Exp $
28// GEANT4 tag $Name:  $
29//
30//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32
33#include "ExN06DetectorConstruction.hh"
34
35#include "G4Material.hh"
36#include "G4MaterialTable.hh"
37#include "G4Element.hh"
38#include "G4ElementTable.hh"
39#include "G4LogicalBorderSurface.hh"
40#include "G4LogicalSkinSurface.hh"
41#include "G4Box.hh"
42#include "G4LogicalVolume.hh"
43#include "G4RotationMatrix.hh"
44#include "G4ThreeVector.hh"
45#include "G4Transform3D.hh"
46#include "G4PVPlacement.hh"
47#include "G4OpBoundaryProcess.hh"
48
49//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
50
51ExN06DetectorConstruction::ExN06DetectorConstruction()
52{
53  expHall_x = expHall_y = expHall_z = 10*m;
54  tank_x    = tank_y    = tank_z    = 5*m;
55  bubble_x  = bubble_y  = bubble_z  = 0.5*m;
56}
57
58//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
59
60ExN06DetectorConstruction::~ExN06DetectorConstruction(){;}
61
62//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
63
64G4VPhysicalVolume* ExN06DetectorConstruction::Construct()
65{
66
67//      ------------- Materials -------------
68
69  G4double a, z, density;
70  G4int nelements;
71
72// Air
73//
74  G4Element* N = new G4Element("Nitrogen", "N", z=7 , a=14.01*g/mole);
75  G4Element* O = new G4Element("Oxygen"  , "O", z=8 , a=16.00*g/mole);
76
77  G4Material* Air = new G4Material("Air", density=1.29*mg/cm3, nelements=2);
78  Air->AddElement(N, 70.*perCent);
79  Air->AddElement(O, 30.*perCent);
80
81// Water
82//
83  G4Element* H = new G4Element("Hydrogen", "H", z=1 , a=1.01*g/mole);
84
85  G4Material* Water = new G4Material("Water", density= 1.0*g/cm3, nelements=2);
86  Water->AddElement(H, 2);
87  Water->AddElement(O, 1);
88
89//
90// ------------ Generate & Add Material Properties Table ------------
91//
92  const G4int nEntries = 32;
93
94  G4double PhotonEnergy[nEntries] =
95            { 2.034*eV, 2.068*eV, 2.103*eV, 2.139*eV,
96              2.177*eV, 2.216*eV, 2.256*eV, 2.298*eV,
97              2.341*eV, 2.386*eV, 2.433*eV, 2.481*eV,
98              2.532*eV, 2.585*eV, 2.640*eV, 2.697*eV,
99              2.757*eV, 2.820*eV, 2.885*eV, 2.954*eV,
100              3.026*eV, 3.102*eV, 3.181*eV, 3.265*eV,
101              3.353*eV, 3.446*eV, 3.545*eV, 3.649*eV,
102              3.760*eV, 3.877*eV, 4.002*eV, 4.136*eV };
103//
104// Water
105//           
106  G4double RefractiveIndex1[nEntries] =
107            { 1.3435, 1.344,  1.3445, 1.345,  1.3455,
108              1.346,  1.3465, 1.347,  1.3475, 1.348,
109              1.3485, 1.3492, 1.35,   1.3505, 1.351,
110              1.3518, 1.3522, 1.3530, 1.3535, 1.354,
111              1.3545, 1.355,  1.3555, 1.356,  1.3568,
112              1.3572, 1.358,  1.3585, 1.359,  1.3595,
113              1.36,   1.3608};
114
115  G4double Absorption1[nEntries] =
116           {3.448*m,  4.082*m,  6.329*m,  9.174*m, 12.346*m, 13.889*m,
117           15.152*m, 17.241*m, 18.868*m, 20.000*m, 26.316*m, 35.714*m,
118           45.455*m, 47.619*m, 52.632*m, 52.632*m, 55.556*m, 52.632*m,
119           52.632*m, 47.619*m, 45.455*m, 41.667*m, 37.037*m, 33.333*m,
120           30.000*m, 28.500*m, 27.000*m, 24.500*m, 22.000*m, 19.500*m,
121           17.500*m, 14.500*m };
122
123  G4double ScintilFast[nEntries] =
124            { 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
125              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
126              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
127              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
128              1.00, 1.00, 1.00, 1.00 };
129  G4double ScintilSlow[nEntries] =
130            { 0.01, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00,
131              7.00, 8.00, 9.00, 8.00, 7.00, 6.00, 4.00,
132              3.00, 2.00, 1.00, 0.01, 1.00, 2.00, 3.00,
133              4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 8.00,
134              7.00, 6.00, 5.00, 4.00 };
135
136  G4MaterialPropertiesTable* myMPT1 = new G4MaterialPropertiesTable();
137  myMPT1->AddProperty("RINDEX",       PhotonEnergy, RefractiveIndex1,nEntries);
138  myMPT1->AddProperty("ABSLENGTH",    PhotonEnergy, Absorption1,     nEntries);
139  myMPT1->AddProperty("FASTCOMPONENT",PhotonEnergy, ScintilFast,     nEntries);
140  myMPT1->AddProperty("SLOWCOMPONENT",PhotonEnergy, ScintilSlow,     nEntries);
141 
142  myMPT1->AddConstProperty("SCINTILLATIONYIELD",50./MeV);
143  myMPT1->AddConstProperty("RESOLUTIONSCALE",1.0);
144  myMPT1->AddConstProperty("FASTTIMECONSTANT", 1.*ns);
145  myMPT1->AddConstProperty("SLOWTIMECONSTANT",10.*ns);
146  myMPT1->AddConstProperty("YIELDRATIO",0.8);
147 
148  Water->SetMaterialPropertiesTable(myMPT1);
149
150//
151// Air
152//
153  G4double RefractiveIndex2[nEntries] =
154            { 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
155              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
156              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
157              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
158              1.00, 1.00, 1.00, 1.00 };
159
160  G4MaterialPropertiesTable* myMPT2 = new G4MaterialPropertiesTable();
161  myMPT2->AddProperty("RINDEX", PhotonEnergy, RefractiveIndex2, nEntries);
162 
163  Air->SetMaterialPropertiesTable(myMPT2);
164
165//
166//      ------------- Volumes --------------
167
168// The experimental Hall
169//
170  G4Box* expHall_box = new G4Box("World",expHall_x,expHall_y,expHall_z);
171
172  G4LogicalVolume* expHall_log
173    = new G4LogicalVolume(expHall_box,Air,"World",0,0,0);
174
175  G4VPhysicalVolume* expHall_phys
176    = new G4PVPlacement(0,G4ThreeVector(),expHall_log,"World",0,false,0);
177
178// The Water Tank
179//     
180  G4Box* waterTank_box = new G4Box("Tank",tank_x,tank_y,tank_z);
181
182  G4LogicalVolume* waterTank_log
183    = new G4LogicalVolume(waterTank_box,Water,"Tank",0,0,0);
184
185  G4VPhysicalVolume* waterTank_phys
186    = new G4PVPlacement(0,G4ThreeVector(),waterTank_log,"Tank",
187                        expHall_log,false,0);
188
189// The Air Bubble
190//   
191  G4Box* bubbleAir_box = new G4Box("Bubble",bubble_x,bubble_y,bubble_z);
192
193  G4LogicalVolume* bubbleAir_log
194    = new G4LogicalVolume(bubbleAir_box,Air,"Bubble",0,0,0);
195
196//G4VPhysicalVolume* bubbleAir_phys =
197      new G4PVPlacement(0,G4ThreeVector(0,2.5*m,0),bubbleAir_log,"Bubble",
198                        waterTank_log,false,0);
199
200//      ------------- Surfaces --------------
201//
202// Water Tank
203//
204  G4OpticalSurface* OpWaterSurface = new G4OpticalSurface("WaterSurface");
205  OpWaterSurface->SetType(dielectric_dielectric);
206  OpWaterSurface->SetFinish(ground);
207  OpWaterSurface->SetModel(unified);
208
209  G4LogicalBorderSurface* WaterSurface = 
210                                 new G4LogicalBorderSurface("WaterSurface",
211                                 waterTank_phys,expHall_phys,OpWaterSurface);
212
213  if(WaterSurface->GetVolume1() == waterTank_phys) G4cout << "Equal" << G4endl;
214  if(WaterSurface->GetVolume2() == expHall_phys  ) G4cout << "Equal" << G4endl;
215
216// Air Bubble
217//
218  G4OpticalSurface* OpAirSurface = new G4OpticalSurface("AirSurface");
219  OpAirSurface->SetType(dielectric_dielectric);
220  OpAirSurface->SetFinish(polished);
221  OpAirSurface->SetModel(glisur);
222
223  G4LogicalSkinSurface* AirSurface = 
224          new G4LogicalSkinSurface("AirSurface", bubbleAir_log, OpAirSurface);
225
226  if(AirSurface->GetLogicalVolume() == bubbleAir_log) G4cout << "Equal" << G4endl;
227  ((G4OpticalSurface*)
228  (AirSurface->GetSurface(bubbleAir_log)->GetSurfaceProperty()))->DumpInfo();
229
230//
231// Generate & Add Material Properties Table attached to the optical surfaces
232//
233  const G4int num = 2;
234  G4double Ephoton[num] = {2.038*eV, 4.144*eV};
235
236  //OpticalWaterSurface
237  G4double RefractiveIndex[num] = {1.35, 1.40};
238  G4double SpecularLobe[num]    = {0.3, 0.3};
239  G4double SpecularSpike[num]   = {0.2, 0.2};
240  G4double Backscatter[num]     = {0.2, 0.2};
241
242  G4MaterialPropertiesTable* myST1 = new G4MaterialPropertiesTable();
243 
244  myST1->AddProperty("RINDEX",                Ephoton, RefractiveIndex, num);
245  myST1->AddProperty("SPECULARLOBECONSTANT",  Ephoton, SpecularLobe,    num);
246  myST1->AddProperty("SPECULARSPIKECONSTANT", Ephoton, SpecularSpike,   num);
247  myST1->AddProperty("BACKSCATTERCONSTANT",   Ephoton, Backscatter,     num);
248
249  OpWaterSurface->SetMaterialPropertiesTable(myST1);
250
251  //OpticalAirSurface
252  G4double Reflectivity[num] = {0.3, 0.5};
253  G4double Efficiency[num]   = {0.8, 1.0};
254
255  G4MaterialPropertiesTable *myST2 = new G4MaterialPropertiesTable();
256
257  myST2->AddProperty("REFLECTIVITY", Ephoton, Reflectivity, num);
258  myST2->AddProperty("EFFICIENCY",   Ephoton, Efficiency,   num);
259
260  OpAirSurface->SetMaterialPropertiesTable(myST2);
261
262//always return the physical World
263  return expHall_phys;
264}
265
266//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Note: See TracBrowser for help on using the repository browser.