source: trunk/examples/extended/electromagnetic/TestEm7/src/DetectorConstruction.cc @ 953

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

update
File size: 10.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// $Id: DetectorConstruction.cc,v 1.8 2007/01/11 15:41:46 vnivanch Exp $
27// GEANT4 tag $Name: geant4-09-01-patch-02 $
28//
29//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
30//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31
32#include "DetectorConstruction.hh"
33#include "DetectorMessenger.hh"
34
35#include "G4Material.hh"
36#include "G4Box.hh"
37#include "G4LogicalVolume.hh"
38#include "G4PVPlacement.hh"
39#include "G4UniformMagField.hh"
40
41#include "G4GeometryManager.hh"
42#include "G4PhysicalVolumeStore.hh"
43#include "G4LogicalVolumeStore.hh"
44#include "G4SolidStore.hh"
45
46#include "G4NistManager.hh"
47
48#include "G4UnitsTable.hh"
49
50//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
51
52DetectorConstruction::DetectorConstruction()
53{
54  // default parameter values
55  absorSizeX = absorSizeYZ = 20*cm;
56  worldSizeX = worldSizeYZ = 1.2*absorSizeX;
57 
58  worldMaterial = absorMaterial = 0;
59  magField = 0;
60  lAbsor   = 0;
61 
62  tallySize     = G4ThreeVector();
63  tallyMaterial = 0;
64  tallyMass     = 0.; 
65  tallyNumber   = 0;
66  tallyPosition = new G4ThreeVector[MaxTally];
67  lTally        = 0;
68 
69  DefineMaterials();
70  SetMaterial("Water");
71  SetTallyMaterial("Water");
72
73  // create commands for interactive definition of the detector 
74  detectorMessenger = new DetectorMessenger(this);
75}
76
77//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
78
79DetectorConstruction::~DetectorConstruction()
80{ delete [] tallyPosition; delete detectorMessenger;}
81
82//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
83
84G4VPhysicalVolume* DetectorConstruction::Construct()
85{
86  return ConstructVolumes();
87}
88
89//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
90
91void DetectorConstruction::DefineMaterials()
92{ 
93  //
94  // define Elements
95  //
96  G4double z, a;
97
98  G4Element* H = new G4Element("Hydrogen", "H", z= 1, a= 1.008*g/mole);
99  G4Element* N = new G4Element("Nitrogen", "N", z= 7, a= 14.01*g/mole);
100  G4Element* O = new G4Element("Oxygen"  , "O", z= 8, a= 16.00*g/mole);
101
102  //
103  // define Materials.
104  //
105  G4double density, temperature, pressure;
106  G4int    ncomponents, natoms;
107  G4double fractionmass;
108 
109  G4Material* H2O = 
110    new G4Material("Water", density= 1.0*g/cm3, ncomponents=2);
111  H2O->AddElement(H, natoms=2);
112  H2O->AddElement(O, natoms=1);
113  H2O->GetIonisation()->SetMeanExcitationEnergy(75.0*eV);
114
115  G4Material* Air = 
116    new G4Material("Air"  , density= 1.290*mg/cm3, ncomponents=2);
117  Air->AddElement(N, fractionmass=0.7);
118  Air->AddElement(O, fractionmass=0.3);
119
120  density     = universe_mean_density;    //from PhysicalConstants.h
121  pressure    = 3.e-18*pascal;
122  temperature = 2.73*kelvin;
123  G4Material* vacuum = 
124    new G4Material("Galactic",z= 1,a= 1.008*g/mole,density,
125                   kStateGas,temperature,pressure);
126
127  //default materials
128  worldMaterial = vacuum;
129}
130
131//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
132 
133G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
134{
135  G4GeometryManager::GetInstance()->OpenGeometry();
136  G4PhysicalVolumeStore::GetInstance()->Clean();
137  G4LogicalVolumeStore::GetInstance()->Clean();
138  G4SolidStore::GetInstance()->Clean();
139
140  // World
141  //
142  G4Box*
143  sWorld = new G4Box("World",                                   //name
144                   worldSizeX/2,worldSizeYZ/2,worldSizeYZ/2);   //dimensions
145
146  G4LogicalVolume*                                                           
147  lWorld = new G4LogicalVolume(sWorld,                  //shape
148                               worldMaterial,           //material
149                              "World");                 //name
150
151  G4VPhysicalVolume*                                   
152  pWorld = new G4PVPlacement(0,                         //no rotation
153                           G4ThreeVector(),             //at (0,0,0)
154                           lWorld,                      //logical volume
155                           "World",                     //name
156                           0,                           //mother  volume
157                           false,                       //no boolean operation
158                           0);                          //copy number
159  //                       
160  // Absorber
161  //                       
162  G4Box*
163  sAbsor = new G4Box("Absorber",                                //name
164                   absorSizeX/2,absorSizeYZ/2,absorSizeYZ/2);   //dimensions
165                                                             
166  lAbsor = new G4LogicalVolume(sAbsor,                  //shape
167                               absorMaterial,           //material
168                              "Absorber");              //name
169 
170                             
171           new G4PVPlacement(0,                         //no rotation
172                           G4ThreeVector(),             //at (0,0,0)
173                           lAbsor,                      //logical volume
174                           "Absorber",                  //name
175                           lWorld,                      //mother  volume
176                           false,                       //no boolean operation
177                           0);                          //copy number
178  //
179  // Tallies (optional)
180  //
181  if (tallyNumber > 0) {     
182    G4Box*
183    sTally = new G4Box("Tally",tallySize.x()/2,tallySize.y()/2,tallySize.z()/2);
184    lTally = new G4LogicalVolume(sTally,tallyMaterial,"Tally");
185   
186    for (G4int j=0; j<tallyNumber; j++)
187       {
188        new G4PVPlacement(0,                            //no rotation
189                          tallyPosition[j],             //position
190                          lTally,                       //logical volume
191                          "Tally",                      //name
192                          lAbsor,                       //mother  volume
193                          false,                        //no boolean operation
194                          j);                           //copy number
195       }
196       
197    tallyMass = tallySize.x()*tallySize.y()*tallySize.z()
198               *(tallyMaterial->GetDensity());
199  } 
200
201  PrintParameters();
202   
203  //
204  //always return the World volume
205  // 
206  return pWorld;
207}
208
209//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
210
211void DetectorConstruction::PrintParameters()
212{
213  G4cout << *(G4Material::GetMaterialTable()) << G4endl;
214  G4cout << "\n---------------------------------------------------------\n";
215  G4cout << "---> The Absorber is " << G4BestUnit(absorSizeX,"Length")
216         << " of " << absorMaterial->GetName() << G4endl;
217  G4cout << "\n---------------------------------------------------------\n";
218 
219  if (tallyNumber > 0) {
220    G4cout << "---> There are " << tallyNumber << " tallies : "
221           << G4BestUnit(tallySize,"Length")
222           << " of " << tallyMaterial->GetName()
223           << "  (mass : " << G4BestUnit(tallyMass,"Mass") << ")" << G4endl;
224           
225    for (G4int j=0; j<tallyNumber; j++)
226     G4cout << "tally " << j << ": "
227            << "position = " << G4BestUnit(tallyPosition[j],"Length") << G4endl;
228    G4cout << "\n---------------------------------------------------------\n";
229  }         
230}
231
232//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
233
234void DetectorConstruction::SetSizeX(G4double value)
235{
236  absorSizeX = value; worldSizeX = 1.2*absorSizeX;
237}
238 
239//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
240
241void DetectorConstruction::SetSizeYZ(G4double value)
242{
243  absorSizeYZ = value; 
244  worldSizeYZ = 1.2*absorSizeYZ;
245} 
246
247//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
248
249void DetectorConstruction::SetMaterial(G4String materialChoice)
250{
251  // search the material by its name   
252  G4Material* pttoMaterial =
253    G4NistManager::Instance()->FindOrBuildMaterial(materialChoice);
254  if (pttoMaterial) absorMaterial = pttoMaterial;
255}
256
257//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
258
259#include "G4FieldManager.hh"
260#include "G4TransportationManager.hh"
261
262void DetectorConstruction::SetMagField(G4double fieldValue)
263{
264  //apply a global uniform magnetic field along Z axis
265  G4FieldManager* fieldMgr
266   = G4TransportationManager::GetTransportationManager()->GetFieldManager();
267   
268  if (magField) delete magField;        //delete the existing magn field
269 
270  if (fieldValue!=0.)                   // create a new one if non nul
271    {
272      magField = new G4UniformMagField(G4ThreeVector(0.,0.,fieldValue));       
273      fieldMgr->SetDetectorField(magField);
274      fieldMgr->CreateChordFinder(magField);
275    }
276   else
277    {
278      magField = 0;
279      fieldMgr->SetDetectorField(magField);
280    }
281}
282//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
283
284void DetectorConstruction::SetTallySize(G4ThreeVector value)
285{
286  tallySize = value;
287} 
288
289//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
290
291void DetectorConstruction::SetTallyMaterial(G4String materialChoice)
292{
293  // search the material by its name   
294  G4Material* pttoMaterial =
295    G4NistManager::Instance()->FindOrBuildMaterial(materialChoice);
296  if (pttoMaterial) tallyMaterial = pttoMaterial;
297}
298
299//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
300
301void DetectorConstruction::SetTallyPosition(G4ThreeVector value)
302{
303  if (tallyNumber < MaxTally) {
304    tallyPosition[tallyNumber] = value; 
305    tallyNumber++;
306  }
307} 
308
309//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
310
311#include "G4RunManager.hh"
312 
313void DetectorConstruction::UpdateGeometry()
314{
315G4RunManager::GetRunManager()->DefineWorldVolume(ConstructVolumes());
316}
317
318//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Note: See TracBrowser for help on using the repository browser.