Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/dna/microyz/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 // This example is provided by the Geant4-DNA collaboration
 27 // Any report or published results obtained using the Geant4-DNA software
 28 // shall cite the following Geant4-DNA collaboration publications:
 29 // Med. Phys. 45 (2018) e722-e739
 30 // Phys. Med. 31 (2015) 861-874
 31 // Med. Phys. 37 (2010) 4692-4708
 32 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 157–178
 33 // The Geant4-DNA web site is available at http://geant4-dna.org
 34 //
 35 /// \file DetectorConstruction.cc
 36 /// \brief Implementation of the DetectorConstruction class
 37 
 38 #include "DetectorConstruction.hh"
 39 
 40 #include "DetectorMessenger.hh"
 41 #include "TrackerSD.hh"
 42 
 43 #include "G4Box.hh"
 44 #include "G4GeometryManager.hh"
 45 #include "G4GeometryTolerance.hh"
 46 #include "G4NistManager.hh"
 47 #include "G4PVPlacement.hh"
 48 #include "G4SDManager.hh"
 49 #include "G4SystemOfUnits.hh"
 50 #include "G4UnitsTable.hh"
 51 #include "G4UserLimits.hh"
 52 
 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 54 
 55 DetectorConstruction::DetectorConstruction() : G4VUserDetectorConstruction()
 56 {
 57   // Create commands for interactive definition of the detector
 58   fDetectorMessenger = std::make_unique<DetectorMessenger>(this);
 59 }
 60 
 61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 62 
 63 DetectorConstruction::~DetectorConstruction() = default;
 64 
 65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 66 
 67 G4VPhysicalVolume* DetectorConstruction::Construct()
 68 {
 69   // Define materials
 70   DefineMaterials();
 71 
 72   // Define volumes
 73   return DefineVolumes();
 74 }
 75 
 76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 77 
 78 void DetectorConstruction::DefineMaterials()
 79 {
 80   // Water is defined from NIST material database
 81   G4NistManager* man = G4NistManager::Instance();
 82   G4Material* H2O = man->FindOrBuildMaterial("G4_WATER");
 83   fpWaterMaterial = H2O;
 84 
 85   // Print materials
 86   G4cout << *(G4Material::GetMaterialTable()) << G4endl;
 87 }
 88 
 89 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 90 
 91 G4VPhysicalVolume* DetectorConstruction::DefineVolumes()
 92 {
 93   G4double worldLength = 10 * m;
 94 
 95   // World
 96 
 97   G4GeometryManager::GetInstance()->SetWorldMaximumExtent(worldLength);
 98 
 99   G4cout << "Computed tolerance = "
100          << G4GeometryTolerance::GetInstance()->GetSurfaceTolerance() / nm << " nm" << G4endl;
101 
102   auto* worldS = new G4Box("world",  // its name
103                            worldLength / 2, worldLength / 2, worldLength / 2);  // its size
104 
105   auto* worldLV = new G4LogicalVolume(worldS,  // its solid
106                                       fpWaterMaterial,  // its material
107                                       "World_LV");  // its name
108 
109   G4VPhysicalVolume* worldPV = new G4PVPlacement(nullptr,  // no rotation
110                                                  G4ThreeVector(),  // at (0,0,0)
111                                                  worldLV,  // its logical volume
112                                                  "World",  // its name
113                                                  nullptr,  // its mother  volume
114                                                  false,  // no boolean operations
115                                                  0,  // copy number
116                                                  false);  // checking overlaps
117 
118   worldLV->SetUserLimits(new G4UserLimits(fpMaxStepSize, DBL_MAX, DBL_MAX, fpTrackingCut));
119 
120   PrintParameters();
121 
122   return worldPV;
123 }
124 
125 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
126 
127 void DetectorConstruction::ConstructSDandField()
128 {
129   // Sensitive detector
130 
131   G4String trackerChamberSDname = "TrackerChamberSD";
132 
133   auto* aTrackerSD = new TrackerSD(trackerChamberSDname, "TrackerHitsCollection");
134   aTrackerSD->SetRadius(fRadius);
135 
136   G4SDManager::GetSDMpointer()->AddNewDetector(aTrackerSD);
137 
138   SetSensitiveDetector("World_LV", aTrackerSD, true);
139 }
140 
141 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
142 
143 void DetectorConstruction::SetTrackingCut(const G4double& value)
144 {
145   fpTrackingCut = value;
146 }
147 
148 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
149 
150 void DetectorConstruction::SetMaxStepSize(const G4double& value)
151 {
152   fpMaxStepSize = value;
153 }
154 
155 void DetectorConstruction::SetTrackerSDRadius(const G4double& value)
156 {
157   fRadius = value;
158 }
159 
160 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
161 
162 void DetectorConstruction::PrintParameters() const
163 {
164   G4cout << "\n---------------------------------------------------------\n";
165   G4cout << "---> The tracking cut is set to " << G4BestUnit(fpTrackingCut, "Energy") << G4endl;
166   G4cout << "---> The maximum step size is set to " << G4BestUnit(fpMaxStepSize, "Length")
167          << G4endl;
168   G4cout << "---> The TrackerSD radius is set to " << G4BestUnit(fRadius, "Length") << G4endl;
169   G4cout << "\n---------------------------------------------------------\n";
170 }
171 
172 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
173