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 //
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 publication:
29 // Med. Phys. 37 (2010) 4692-4708
30 // The Geant4-DNA web site is available at http://geant4-dna.org
31 //
32 //
35 
36 #include "DetectorConstruction.hh"
37 #include "DetectorMessenger.hh"
38 
39 #include "G4SystemOfUnits.hh"
40 #include "G4UserLimits.hh"
41 #include "G4NistManager.hh"
42 #include "G4RunManager.hh"
43 #include "G4LogicalVolumeStore.hh"
44 
45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
46 
48  G4VUserDetectorConstruction(), fpWaterMaterial(nullptr),
49  fLogicWorld(nullptr),fPhysiWorld(nullptr)
50 {
51  // create commands for interactive definition of the detector
53 }
54 
55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
56 
58 {
59  delete fDetectorMessenger;
60 }
61 
62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
63 
65 
66 {
67  if(fPhysiWorld) { return fPhysiWorld; }
69  return ConstructDetector();
70 }
71 
72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
73 
75 {
76 
77  // Water is defined from NIST material database
79 
80  G4Material * H2O = man->FindOrBuildMaterial("G4_WATER");
81 
82  /*
83  If one wishes to test other density value for water material,
84  one should use instead:
85  G4Material * H2O = man->BuildMaterialWithNewDensity("G4_WATER_MODIFIED",
86  "G4_WATER",1.100*g/cm3);
87 
88  Note: any string for "G4_WATER_MODIFIED" parameter is accepted
89  and "G4_WATER" parameter should not be changed
90  Both materials are created and can be selected from dna.mac
91  */
92  fpWaterMaterial = H2O;
93 
94  //G4cout << "-> Density of water material (g/cm3)="
95  // << fpWaterMaterial->GetDensity()/(g/cm/cm/cm) << G4endl;
96 
98 }
99 
100 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
102 {
103  // WORLD VOLUME
105  G4double worldSizeY = worldSizeX;
106  G4double worldSizeZ = worldSizeX;
107 
108  G4Box* solidWorld = new G4Box("World", //its name
109  worldSizeX / 2, worldSizeY / 2, worldSizeZ / 2); //its size
110 
111  fLogicWorld = new G4LogicalVolume(solidWorld, //its solid
112  fpWaterMaterial, //its material
113  "World"); //its name
114 
115  fPhysiWorld = new G4PVPlacement(0, //no rotation
116  G4ThreeVector(), //at (0,0,0)
117  "World", //its name
118  fLogicWorld, //its logical volume
119  0, //its mother volume
120  false, //no boolean operation
121  0); //copy number
122 
123  // Visualization attributes - white
124  G4VisAttributes* worldVisAtt = new G4VisAttributes(G4Colour(1.0, 1.0, 1.0));
125  worldVisAtt->SetVisibility(true);
126  fLogicWorld->SetVisAttributes(worldVisAtt);
127 
128  G4VisAttributes* worldVisAtt1 = new G4VisAttributes(G4Colour(1.0, 0.0, 0.0));
129  worldVisAtt1->SetVisibility(true);
130 
131  //
132  // Shows how to introduce a 20 eV tracking cut
133  //
134  //logicWorld->SetUserLimits(new G4UserLimits(DBL_MAX,DBL_MAX,DBL_MAX,20*eV));
135 
136  return fPhysiWorld;
137 }
138 
139 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
140 
141 void DetectorConstruction::SetMaterial(G4String materialChoice)
142 {
143  // Search the material by its name
145  materialChoice);
146 
147  if (pttoMaterial)
148  {
149  fpWaterMaterial = pttoMaterial;
150  if(fLogicWorld) {
152  }
154  }
155 }