| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // This example is provided by the Geant4-DNA 26 // This example is provided by the Geant4-DNA collaboration
27 // Any report or published results obtained us << 27 // Any report or published results obtained using the Geant4-DNA software
28 // shall cite the following Geant4-DNA collabo 28 // shall cite the following Geant4-DNA collaboration publications:
29 // Med. Phys. 45 (2018) e722-e739 << 29 // Med. Phys. 45 (2018) e722-e739
30 // Phys. Med. 31 (2015) 861-874 << 30 // Phys. Med. 31 (2015) 861-874
31 // Med. Phys. 37 (2010) 4692-4708 << 31 // Med. Phys. 37 (2010) 4692-4708
32 // Int. J. Model. Simul. Sci. Comput. 1 (2010) << 32 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 157\u2013178
33 // 33 //
34 // The Geant4-DNA web site is available at htt 34 // The Geant4-DNA web site is available at http://geant4-dna.org
35 // 35 //
36 /// \file DetectorConstruction.cc 36 /// \file DetectorConstruction.cc
37 /// \brief Implementation of the DetectorConst 37 /// \brief Implementation of the DetectorConstruction class
38 << 38
39 #include "DetectorConstruction.hh" 39 #include "DetectorConstruction.hh"
40 <<
41 #include "DetectorMessenger.hh" 40 #include "DetectorMessenger.hh"
42 #include "TrackerSD.hh" 41 #include "TrackerSD.hh"
43 42
44 #include "G4Box.hh" << 43 #include "G4Material.hh"
45 #include "G4GeometryManager.hh" <<
46 #include "G4GeometryTolerance.hh" <<
47 #include "G4NistManager.hh" 44 #include "G4NistManager.hh"
>> 45
>> 46 #include "G4Box.hh"
48 #include "G4PVPlacement.hh" 47 #include "G4PVPlacement.hh"
49 #include "G4SDManager.hh" << 48
>> 49 #include "G4GeometryTolerance.hh"
>> 50 #include "G4GeometryManager.hh"
>> 51
50 #include "G4SystemOfUnits.hh" 52 #include "G4SystemOfUnits.hh"
51 #include "G4UnitsTable.hh" <<
52 #include "G4UserLimits.hh" 53 #include "G4UserLimits.hh"
>> 54 #include "G4UnitsTable.hh"
53 55
54 //....oooOO0OOooo........oooOO0OOooo........oo << 56 #include "G4SDManager.hh"
55 57
56 DetectorConstruction::DetectorConstruction() : << 58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 59
>> 60 DetectorConstruction::DetectorConstruction()
>> 61 :G4VUserDetectorConstruction(),
>> 62 fDetectorMessenger(0)
57 { 63 {
58 // Default tracking cut 64 // Default tracking cut
59 fTrackingCut = 9 * eV; << 65 fTrackingCut = 9*eV;
60 66
61 // Create commands for interactive definitio << 67 // Create commands for interactive definition of the detector
62 fDetectorMessenger = new DetectorMessenger(t 68 fDetectorMessenger = new DetectorMessenger(this);
>> 69
63 } 70 }
64 71
65 //....oooOO0OOooo........oooOO0OOooo........oo 72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
66 << 73
67 DetectorConstruction::~DetectorConstruction() 74 DetectorConstruction::~DetectorConstruction()
68 { 75 {
69 delete fDetectorMessenger; 76 delete fDetectorMessenger;
70 } 77 }
71 78
72 //....oooOO0OOooo........oooOO0OOooo........oo 79 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
73 << 80
74 G4VPhysicalVolume* DetectorConstruction::Const 81 G4VPhysicalVolume* DetectorConstruction::Construct()
75 { 82 {
76 // Define materials 83 // Define materials
77 DefineMaterials(); 84 DefineMaterials();
78 85
79 // Define volumes 86 // Define volumes
80 return DefineVolumes(); 87 return DefineVolumes();
81 } 88 }
82 89
83 //....oooOO0OOooo........oooOO0OOooo........oo 90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
84 91
85 void DetectorConstruction::DefineMaterials() 92 void DetectorConstruction::DefineMaterials()
86 { 93 {
87 // Water is defined from NIST material datab 94 // Water is defined from NIST material database
88 G4NistManager* man = G4NistManager::Instance << 95 G4NistManager * man = G4NistManager::Instance();
89 <<
90 G4Material* H2O = man->FindOrBuildMaterial(" <<
91 96
>> 97 G4Material * H2O = man->FindOrBuildMaterial("G4_WATER");
>> 98
92 fWaterMaterial = H2O; 99 fWaterMaterial = H2O;
93 100
94 // Print materials 101 // Print materials
95 G4cout << *(G4Material::GetMaterialTable()) 102 G4cout << *(G4Material::GetMaterialTable()) << G4endl;
96 } 103 }
97 104
98 //....oooOO0OOooo........oooOO0OOooo........oo 105 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
99 106
100 G4VPhysicalVolume* DetectorConstruction::Defin 107 G4VPhysicalVolume* DetectorConstruction::DefineVolumes()
101 { 108 {
102 G4double worldLength = 1 * km; << 109 G4double worldLength = 1*km;
103 110
104 // World 111 // World
105 112
106 G4GeometryManager::GetInstance()->SetWorldMa 113 G4GeometryManager::GetInstance()->SetWorldMaximumExtent(worldLength);
107 114
108 G4cout << "Computed tolerance = " 115 G4cout << "Computed tolerance = "
109 << G4GeometryTolerance::GetInstance() << 116 << G4GeometryTolerance::GetInstance()->GetSurfaceTolerance()/nm
110 << 117 << " nm" << G4endl;
111 G4Box* worldS = new G4Box("world", // its n <<
112 worldLength / 2, w <<
113 <<
114 G4LogicalVolume* worldLV = new G4LogicalVolu <<
115 <<
116 <<
117 <<
118 G4VPhysicalVolume* worldPV = new G4PVPlaceme <<
119 <<
120 <<
121 <<
122 <<
123 <<
124 <<
125 <<
126 118
127 worldLV->SetUserLimits(new G4UserLimits(DBL_ << 119 G4Box* worldS
>> 120 = new G4Box("world", //its name
>> 121 worldLength/2,worldLength/2,worldLength/2); //its size
>> 122
>> 123 G4LogicalVolume* worldLV
>> 124 = new G4LogicalVolume(
>> 125 worldS, //its solid
>> 126 fWaterMaterial, //its material
>> 127 "World_LV"); //its name
>> 128
>> 129 G4VPhysicalVolume* worldPV
>> 130 = new G4PVPlacement(
>> 131 0, // no rotation
>> 132 G4ThreeVector(), // at (0,0,0)
>> 133 worldLV, // its logical volume
>> 134 "World", // its name
>> 135 0, // its mother volume
>> 136 false, // no boolean operations
>> 137 0, // copy number
>> 138 false); // checking overlaps
>> 139
>> 140 worldLV->SetUserLimits(new G4UserLimits(DBL_MAX,DBL_MAX,DBL_MAX,
>> 141 fTrackingCut));
128 142
129 PrintParameters(); 143 PrintParameters();
130 << 144
131 return worldPV; 145 return worldPV;
132 } 146 }
133 147
134 //....oooOO0OOooo........oooOO0OOooo........oo 148 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
135 << 149
136 void DetectorConstruction::ConstructSDandField 150 void DetectorConstruction::ConstructSDandField()
137 { 151 {
138 // Sensitive detectors 152 // Sensitive detectors
139 153
140 G4String trackerChamberSDname = "TrackerCham 154 G4String trackerChamberSDname = "TrackerChamberSD";
141 << 155
142 TrackerSD* aTrackerSD = new TrackerSD(tracke << 156 TrackerSD* aTrackerSD = new TrackerSD(trackerChamberSDname,
143 << 157 "TrackerHitsCollection");
>> 158
144 G4SDManager::GetSDMpointer()->AddNewDetector 159 G4SDManager::GetSDMpointer()->AddNewDetector(aTrackerSD);
145 160
>> 161 // Setting aTrackerSD to all logical volumes with the same name
>> 162 // of "Chamber_LV".
146 SetSensitiveDetector("World_LV", aTrackerSD, 163 SetSensitiveDetector("World_LV", aTrackerSD, true);
147 } << 164
>> 165 }
148 166
149 //....oooOO0OOooo........oooOO0OOooo........oo 167 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
150 168
151 void DetectorConstruction::SetTrackingCut(G4do 169 void DetectorConstruction::SetTrackingCut(G4double value)
152 { 170 {
153 fTrackingCut = value; 171 fTrackingCut = value;
154 } 172 }
155 173
156 //....oooOO0OOooo........oooOO0OOooo........oo 174 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
157 175
158 void DetectorConstruction::PrintParameters() c 176 void DetectorConstruction::PrintParameters() const
159 { 177 {
160 G4cout << "\n------------------------------- 178 G4cout << "\n---------------------------------------------------------\n";
161 G4cout << "---> The tracking cut is set to " << 179 G4cout << "---> The tracking cut is set to "
>> 180 << G4BestUnit(fTrackingCut,"Energy") << G4endl;
162 G4cout << "\n------------------------------- 181 G4cout << "\n---------------------------------------------------------\n";
163 } 182 }
164 183
165 //....oooOO0OOooo........oooOO0OOooo........oo 184 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 185
166 186