ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
DetectorConstruction.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file DetectorConstruction.cc
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 //
28 
29 #include "DetectorConstruction.hh"
30 #include "DetectorMessenger.hh"
31 
32 #include "G4NistManager.hh"
33 #include "G4Sphere.hh"
34 #include "G4LogicalVolume.hh"
35 #include "G4VPhysicalVolume.hh"
36 #include "G4PVPlacement.hh"
37 #include "G4PVReplica.hh"
38 
39 #include "G4GeometryManager.hh"
40 #include "G4PhysicalVolumeStore.hh"
41 #include "G4LogicalVolumeStore.hh"
42 #include "G4SolidStore.hh"
43 #include "G4RunManager.hh"
44 
45 #include "G4UnitsTable.hh"
46 #include "G4PhysicalConstants.hh"
47 #include "G4SystemOfUnits.hh"
48 #include "G4UserLimits.hh"
49 
50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
51 
52 DetectorConstruction::DetectorConstruction()
53  : G4VUserDetectorConstruction(),
54  fAbsorMaterial(0),
55  fAbsor(0),
56  fDetectorMessenger(0)
57 {
58  //default tracking cut
59  fTrackingCut = 7.4*eV;
60 
61  // default parameter values
62  fAbsorRadius = 1*m;
63 
64  DefineMaterials();
65  SetMaterial("G4_WATER");
66 
67  // create commands for interactive definition of the detector
68  fDetectorMessenger = new DetectorMessenger(this);
69 }
70 
71 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
72 
73 DetectorConstruction::~DetectorConstruction()
74 { delete fDetectorMessenger;}
75 
76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
77 
78 G4VPhysicalVolume* DetectorConstruction::Construct()
79 {
80  return ConstructVolumes();
81 }
82 
83 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
84 
85 void DetectorConstruction::DefineMaterials()
86 {
87  G4NistManager* man = G4NistManager::Instance();
88 
89  man->FindOrBuildMaterial("G4_WATER");
90 
91 }
92 
93 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
94 
95 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
96 {
97  G4GeometryManager::GetInstance()->OpenGeometry();
98  G4PhysicalVolumeStore::GetInstance()->Clean();
99  G4LogicalVolumeStore::GetInstance()->Clean();
100  G4SolidStore::GetInstance()->Clean();
101 
102  // Spherical absorber
103  //
104  G4Sphere*
105  sAbsor = new G4Sphere("Absorber", //name
106  0., fAbsorRadius, 0., twopi, 0., pi); //size
107 
108  fLogicalAbsor = new G4LogicalVolume(sAbsor, //solid
109  fAbsorMaterial, //material
110  "Absorber"); //name
111 
112  fAbsor = new G4PVPlacement(0, //no rotation
113  G4ThreeVector(), //at (0,0,0)
114  fLogicalAbsor, //logical volume
115  "Absorber", //name
116  0, //mother volume
117  false, //no boolean operation
118  0); //copy number
119 
120  PrintParameters();
121 
122  fLogicalAbsor->SetUserLimits(new G4UserLimits(DBL_MAX,DBL_MAX,DBL_MAX,
123  fTrackingCut));
124 
125  //
126  //always return the root volume
127  //
128  return fAbsor;
129 }
130 
131 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
132 
133 void DetectorConstruction::PrintParameters() const
134 {
135  G4cout << "\n---------------------------------------------------------\n";
136  G4cout << "---> The tracking cut is set to "
137  << G4BestUnit(fTrackingCut,"Energy") << G4endl;
138  G4cout << "---> The Absorber is a sphere of "
139  << G4BestUnit(fAbsorRadius,"Length") << " radius of "
140  << fAbsorMaterial->GetName() << " made of"
141  << "\n \n" << fAbsorMaterial << G4endl;
142  G4cout << "\n---------------------------------------------------------\n";
143 }
144 
145 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
146 
147 void DetectorConstruction::SetTrackingCut(G4double value)
148 {
149  fTrackingCut = value;
150  G4RunManager::GetRunManager()->ReinitializeGeometry();
151 }
152 
153 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
154 
155 void DetectorConstruction::SetRadius(G4double value)
156 {
157  fAbsorRadius = value;
158  G4RunManager::GetRunManager()->ReinitializeGeometry();
159 }
160 
161 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
162 
163 void DetectorConstruction::SetMaterial(G4String materialChoice)
164 {
165  // search the material by its name
166  G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);
167  if (pttoMaterial) fAbsorMaterial = pttoMaterial;
168  G4RunManager::GetRunManager()->PhysicsHasBeenModified();
169 }