Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/hadronic/Hadr03/src/DetectorConstruction.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  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 /// \file DetectorConstruction.cc
 27 /// \brief Implementation of the DetectorConstruction class
 28 //
 29 //
 30 
 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 33 
 34 #include "DetectorConstruction.hh"
 35 
 36 #include "DetectorMessenger.hh"
 37 
 38 #include "G4Box.hh"
 39 #include "G4GeometryManager.hh"
 40 #include "G4LogicalVolume.hh"
 41 #include "G4LogicalVolumeStore.hh"
 42 #include "G4Material.hh"
 43 #include "G4NistManager.hh"
 44 #include "G4PVPlacement.hh"
 45 #include "G4PhysicalVolumeStore.hh"
 46 #include "G4RunManager.hh"
 47 #include "G4SolidStore.hh"
 48 #include "G4SystemOfUnits.hh"
 49 #include "G4UnitsTable.hh"
 50 
 51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 52 
 53 DetectorConstruction::DetectorConstruction()
 54 {
 55   fBoxSize = 10 * m;
 56   DefineMaterials();
 57   SetMaterial("Molybdenum98");
 58   fDetectorMessenger = new DetectorMessenger(this);
 59 }
 60 
 61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 62 
 63 DetectorConstruction::~DetectorConstruction()
 64 {
 65   delete fDetectorMessenger;
 66 }
 67 
 68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 69 
 70 G4VPhysicalVolume* DetectorConstruction::Construct()
 71 {
 72   return ConstructVolumes();
 73 }
 74 
 75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 76 
 77 void DetectorConstruction::DefineMaterials()
 78 {
 79   // define a Material from isotopes
 80   //
 81   MaterialWithSingleIsotope("Molybdenum98", "Mo98", 10.28 * g / cm3, 42, 98);
 82 
 83   // NE213
 84   G4Element* H = new G4Element("Hydrogen", "H", 1., 1.01 * g / mole);
 85   G4Element* C = new G4Element("Carbon", "C", 6., 12.00 * g / mole);
 86   G4Material* ne213 = new G4Material("NE213", 0.874 * g / cm3, 2);
 87   ne213->AddElement(H, 9.2 * perCent);
 88   ne213->AddElement(C, 90.8 * perCent);
 89 
 90   G4Material* hydrogen = new G4Material("hydrogen", 1.0 * g / cm3, 1);
 91   hydrogen->AddElement(H, 1);
 92 
 93   G4Material* carbon = new G4Material("carbon", 1.0 * g / cm3, 1);
 94   carbon->AddElement(C, 1);
 95 
 96   G4Material* plastic = new G4Material("plastic", 1.0 * g / cm3, 2);
 97   plastic->AddElement(H, 1);
 98   plastic->AddElement(C, 1);
 99 
100   // or use G4-NIST materials data base
101   //
102   G4NistManager* man = G4NistManager::Instance();
103   man->FindOrBuildMaterial("G4_B");
104 
105   G4Element* O = man->FindOrBuildElement("O");
106   G4Element* Hf = man->FindOrBuildElement("Hf");
107 
108   G4Material* HfO2 = new G4Material("HfO2", 9.68 * g / cm3, 2);
109   HfO2->AddElement(Hf, 1);
110   HfO2->AddElement(O, 2);
111 
112   /// G4cout << *(G4Material::GetMaterialTable()) << G4endl;
113 }
114 
115 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
116 
117 G4Material* DetectorConstruction::MaterialWithSingleIsotope(G4String name, G4String symbol,
118                                                             G4double density, G4int Z, G4int A)
119 {
120   // define a material from an isotope
121   //
122   G4int ncomponents;
123   G4double abundance, massfraction;
124 
125   G4Isotope* isotope = new G4Isotope(symbol, Z, A);
126 
127   G4Element* element = new G4Element(name, symbol, ncomponents = 1);
128   element->AddIsotope(isotope, abundance = 100. * perCent);
129 
130   G4Material* material = new G4Material(name, density, ncomponents = 1);
131   material->AddElement(element, massfraction = 100. * perCent);
132 
133   return material;
134 }
135 
136 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
137 
138 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
139 {
140   // Cleanup old geometry
141   G4GeometryManager::GetInstance()->OpenGeometry();
142   G4PhysicalVolumeStore::GetInstance()->Clean();
143   G4LogicalVolumeStore::GetInstance()->Clean();
144   G4SolidStore::GetInstance()->Clean();
145 
146   G4Box* sBox = new G4Box("Container",  // its name
147                           fBoxSize / 2, fBoxSize / 2, fBoxSize / 2);  // its dimensions
148 
149   fLBox = new G4LogicalVolume(sBox,  // its shape
150                               fMaterial,  // its material
151                               fMaterial->GetName());  // its name
152 
153   fPBox = new G4PVPlacement(0,  // no rotation
154                             G4ThreeVector(),  // at (0,0,0)
155                             fLBox,  // its logical volume
156                             fMaterial->GetName(),  // its name
157                             0,  // its mother  volume
158                             false,  // no boolean operation
159                             0);  // copy number
160 
161   PrintParameters();
162 
163   // always return the root volume
164   //
165   return fPBox;
166 }
167 
168 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
169 
170 void DetectorConstruction::PrintParameters()
171 {
172   G4cout << "\n The Box is " << G4BestUnit(fBoxSize, "Length") << " of " << fMaterial->GetName()
173          << "\n \n"
174          << fMaterial << G4endl;
175 }
176 
177 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
178 
179 void DetectorConstruction::SetMaterial(G4String materialChoice)
180 {
181   // search the material by its name
182   G4Material* pttoMaterial = G4NistManager::Instance()->FindOrBuildMaterial(materialChoice);
183 
184   if (pttoMaterial) {
185     if (fMaterial != pttoMaterial) {
186       fMaterial = pttoMaterial;
187       if (fLBox) {
188         fLBox->SetMaterial(pttoMaterial);
189       }
190       G4RunManager::GetRunManager()->PhysicsHasBeenModified();
191     }
192   }
193   else {
194     G4cout << "\n--> warning from DetectorConstruction::SetMaterial : " << materialChoice
195            << " not found" << G4endl;
196   }
197 }
198 
199 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
200 
201 void DetectorConstruction::SetSize(G4double value)
202 {
203   fBoxSize = value;
204   G4RunManager::GetRunManager()->ReinitializeGeometry();
205 }
206 
207 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
208