ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
HadrontherapyDetectorConstruction.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file HadrontherapyDetectorConstruction.hh
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 // Hadrontherapy advanced example for Geant4
27 // See more at: https://twiki.cern.ch/twiki/bin/view/Geant4/AdvancedExamplesHadrontherapy
28 
29 #ifndef HadrontherapyDetectorConstruction_H
30 #define HadrontherapyDetectorConstruction_H 1
31 
32 #include "G4Box.hh"
33 #include "globals.hh"
34 #include "G4VisAttributes.hh"
35 #include "G4LogicalVolume.hh"
36 #include "G4UnitsTable.hh"
38 
39 class G4VPhysicalVolume;
40 class G4LogicalVolume;
41 class G4PVPlacement;
46 class HadrontherapyLet;
47 
49 {
50 public:
51 
53 
55 
56 public:
59  G4ThreeVector detectorToWorldPosition);
61  HadrontherapyDetectorSD* detectorSD; // Pointer to sensitive detector
62 
64  void VirtualLayer(G4bool Varbool);
68 
70 private:
71 
72  void ConstructPhantom();
73  void ConstructDetector();
74  void ParametersCheck();
75  void CheckOverlaps();
76 
77 public:
78  // Get detector position relative to WORLD
80  {
82  }
84  // Get displacement between phantom and detector by detector position (center of), phantom (center of) and detector sizes
86  {
90  );
91  }
92 
94  // Calculate (and set) detector position by displacement, phantom and detector sizes
95  inline void SetDetectorPosition()
96  {
97  // Adjust detector position
101 
102 
103  }
105  // Check whether detector is inside phantom
106  inline bool IsInside(G4double detectorX,
107  G4double detectorY,
108  G4double detectorZ,
109  G4double phantomX,
110  G4double phantomY,
111  G4double phantomZ,
113  {
114  // Dimensions check... X Y and Z
115  // Firstly check what dimension we are modifying
116  {
117  if (detectorX > phantomX)
118  {
119  G4cout << "Error: Detector X dimension must be smaller or equal to the correspondent of the phantom" << G4endl;
120  return false;
121  }
122  if ( (phantomX - detectorX) < pos.getX())
123  {
124  G4cout << "Error: X dimension doesn't fit with detector to phantom relative position" << G4endl;
125  return false;
126  }
127  }
128 
129  {
130  if (detectorY > phantomY)
131  {
132  G4cout << "Error: Detector Y dimension must be smaller or equal to the correspondent of the phantom" << G4endl;
133  return false;
134  }
135  if ( (phantomY - detectorY) < pos.getY())
136  {
137  G4cout << "Error: Y dimension doesn't fit with detector to phantom relative position" << G4endl;
138  return false;
139  }
140  }
141 
142  {
143  if (detectorZ > phantomZ)
144  {
145  G4cout << "Error: Detector Z dimension must be smaller or equal to the correspondent of the phantom" << G4endl;
146  return false;
147  }
148  if ( (phantomZ - detectorZ) < pos.getZ())
149  {
150  G4cout << "Error: Z dimension doesn't fit with detector to phantom relative position" << G4endl;
151  return false;
152  }
153  }
154 
155  return true;
156  }
158 
160  void SetVoxelSize(G4double sizeX, G4double sizeY, G4double sizeZ);
161  void SetDetectorSize(G4double sizeX, G4double sizeY, G4double sizeZ);
162  void SetPhantomSize(G4double sizeX, G4double sizeY, G4double sizeZ);
164  void SetDetectorToPhantomPosition(G4ThreeVector DetectorToPhantomPosition);
165  void UpdateGeometry();
166  void PrintParameters();
168 
169 private:
172 
175 
179 
183 
187 
191 
195 
196  G4ThreeVector phantomPosition, detectorPosition, detectorToPhantomPosition; // phantom center, detector center, detector to phantom relative position
197 
201 
205 
207 
210 };
211 #endif