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

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

tag geant4.9.4 beta 1 + modifs locales
File size: 10.2 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.17 2008/07/17 00:32:34 gum Exp $
28// GEANT4 tag $Name: geant4-09-04-beta-01 $
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  // Set the Birks Constant for the Water scintillator
151
152  Water->GetIonisation()->SetBirksConstant(0.126*mm/MeV);
153
154//
155// Air
156//
157  G4double RefractiveIndex2[nEntries] =
158            { 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
159              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
160              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
161              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
162              1.00, 1.00, 1.00, 1.00 };
163
164  G4MaterialPropertiesTable* myMPT2 = new G4MaterialPropertiesTable();
165  myMPT2->AddProperty("RINDEX", PhotonEnergy, RefractiveIndex2, nEntries);
166 
167  Air->SetMaterialPropertiesTable(myMPT2);
168
169//
170//      ------------- Volumes --------------
171
172// The experimental Hall
173//
174  G4Box* expHall_box = new G4Box("World",expHall_x,expHall_y,expHall_z);
175
176  G4LogicalVolume* expHall_log
177    = new G4LogicalVolume(expHall_box,Air,"World",0,0,0);
178
179  G4VPhysicalVolume* expHall_phys
180    = new G4PVPlacement(0,G4ThreeVector(),expHall_log,"World",0,false,0);
181
182// The Water Tank
183//     
184  G4Box* waterTank_box = new G4Box("Tank",tank_x,tank_y,tank_z);
185
186  G4LogicalVolume* waterTank_log
187    = new G4LogicalVolume(waterTank_box,Water,"Tank",0,0,0);
188
189  G4VPhysicalVolume* waterTank_phys
190    = new G4PVPlacement(0,G4ThreeVector(),waterTank_log,"Tank",
191                        expHall_log,false,0);
192
193// The Air Bubble
194//   
195  G4Box* bubbleAir_box = new G4Box("Bubble",bubble_x,bubble_y,bubble_z);
196
197  G4LogicalVolume* bubbleAir_log
198    = new G4LogicalVolume(bubbleAir_box,Air,"Bubble",0,0,0);
199
200//G4VPhysicalVolume* bubbleAir_phys =
201      new G4PVPlacement(0,G4ThreeVector(0,2.5*m,0),bubbleAir_log,"Bubble",
202                        waterTank_log,false,0);
203
204//      ------------- Surfaces --------------
205//
206// Water Tank
207//
208  G4OpticalSurface* OpWaterSurface = new G4OpticalSurface("WaterSurface");
209  OpWaterSurface->SetType(dielectric_dielectric);
210  OpWaterSurface->SetFinish(ground);
211  OpWaterSurface->SetModel(unified);
212
213  G4LogicalBorderSurface* WaterSurface = 
214                                 new G4LogicalBorderSurface("WaterSurface",
215                                 waterTank_phys,expHall_phys,OpWaterSurface);
216
217  if(WaterSurface->GetVolume1() == waterTank_phys) G4cout << "Equal" << G4endl;
218  if(WaterSurface->GetVolume2() == expHall_phys  ) G4cout << "Equal" << G4endl;
219
220// Air Bubble
221//
222  G4OpticalSurface* OpAirSurface = new G4OpticalSurface("AirSurface");
223  OpAirSurface->SetType(dielectric_dielectric);
224  OpAirSurface->SetFinish(polished);
225  OpAirSurface->SetModel(glisur);
226
227  G4LogicalSkinSurface* AirSurface = 
228          new G4LogicalSkinSurface("AirSurface", bubbleAir_log, OpAirSurface);
229
230  if(AirSurface->GetLogicalVolume() == bubbleAir_log) G4cout << "Equal" << G4endl;
231  G4OpticalSurface* opticalSurface = dynamic_cast <G4OpticalSurface*>
232        (AirSurface->GetSurface(bubbleAir_log)->GetSurfaceProperty());
233
234  if (opticalSurface) opticalSurface->DumpInfo();
235
236//
237// Generate & Add Material Properties Table attached to the optical surfaces
238//
239  const G4int num = 2;
240  G4double Ephoton[num] = {2.038*eV, 4.144*eV};
241
242  //OpticalWaterSurface
243  G4double RefractiveIndex[num] = {1.35, 1.40};
244  G4double SpecularLobe[num]    = {0.3, 0.3};
245  G4double SpecularSpike[num]   = {0.2, 0.2};
246  G4double Backscatter[num]     = {0.2, 0.2};
247
248  G4MaterialPropertiesTable* myST1 = new G4MaterialPropertiesTable();
249 
250  myST1->AddProperty("RINDEX",                Ephoton, RefractiveIndex, num);
251  myST1->AddProperty("SPECULARLOBECONSTANT",  Ephoton, SpecularLobe,    num);
252  myST1->AddProperty("SPECULARSPIKECONSTANT", Ephoton, SpecularSpike,   num);
253  myST1->AddProperty("BACKSCATTERCONSTANT",   Ephoton, Backscatter,     num);
254
255  OpWaterSurface->SetMaterialPropertiesTable(myST1);
256
257  //OpticalAirSurface
258  G4double Reflectivity[num] = {0.3, 0.5};
259  G4double Efficiency[num]   = {0.8, 1.0};
260
261  G4MaterialPropertiesTable *myST2 = new G4MaterialPropertiesTable();
262
263  myST2->AddProperty("REFLECTIVITY", Ephoton, Reflectivity, num);
264  myST2->AddProperty("EFFICIENCY",   Ephoton, Efficiency,   num);
265
266  OpAirSurface->SetMaterialPropertiesTable(myST2);
267
268//always return the physical World
269  return expHall_phys;
270}
271
272//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Note: See TracBrowser for help on using the repository browser.