| 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 publication:
29 // Med. Phys. 37 (2010) 4692-4708 29 // Med. Phys. 37 (2010) 4692-4708
30 // The Geant4-DNA web site is available at htt 30 // The Geant4-DNA web site is available at http://geant4-dna.org
31 // 31 //
32 // If you use this example, please cite the fo 32 // If you use this example, please cite the following publication:
33 // Rad. Prot. Dos. 133 (2009) 2-11 33 // Rad. Prot. Dos. 133 (2009) 2-11
34 34
35 #ifndef DetectorConstruction_h 35 #ifndef DetectorConstruction_h
36 #define DetectorConstruction_h 1 36 #define DetectorConstruction_h 1
37 37
38 #include "CellParameterisation.hh" <<
39 #include "EMField.hh" <<
40 #include "G4VUserDetectorConstruction.hh" 38 #include "G4VUserDetectorConstruction.hh"
41 #include "G4Box.hh" 39 #include "G4Box.hh"
42 #include "G4Cons.hh" 40 #include "G4Cons.hh"
43 #include "G4Material.hh" 41 #include "G4Material.hh"
44 #include "G4PVPlacement.hh" 42 #include "G4PVPlacement.hh"
45 #include "G4UserLimits.hh" 43 #include "G4UserLimits.hh"
46 #include "G4PVParameterised.hh" 44 #include "G4PVParameterised.hh"
>> 45 #include "PhantomConfiguration.hh"
>> 46 #include "CellParameterisation.hh"
>> 47
>> 48 #include "EMField.hh"
47 #include "G4EqMagElectricField.hh" 49 #include "G4EqMagElectricField.hh"
48 #include "G4PropagatorInField.hh" 50 #include "G4PropagatorInField.hh"
49 #include "G4TransportationManager.hh" 51 #include "G4TransportationManager.hh"
50 #include "G4ChordFinder.hh" 52 #include "G4ChordFinder.hh"
51 #include "G4ClassicalRK4.hh" 53 #include "G4ClassicalRK4.hh"
52 54
>> 55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 56
53 class DetectorConstruction : public G4VUserDet 57 class DetectorConstruction : public G4VUserDetectorConstruction
54 { 58 {
55 public: 59 public:
56 60
57 explicit DetectorConstruction(); << 61 DetectorConstruction();
58 ~DetectorConstruction() override; << 62 ~DetectorConstruction();
59 63
60 G4VPhysicalVolume* Construct() override; << 64 G4VPhysicalVolume* Construct();
61 <<
62 void DefineMaterials(); <<
63 65
64 void ConstructSDandField() override; << 66 void ConstructSDandField();
65 67
66 void SetMassNucleus(G4double mN) {fMassNucle << 68 void SetMassNucleus(G4float mN){ fMassNucleus = mN;}
67 G4double GetMassNucleus() const {return fMa << 69 G4float GetMassNucleus(){return fMassNucleus;}
68 70
69 void SetMassCytoplasm(G4double mC) {fMassCyt << 71 void SetMassCytoplasm(G4float mC){ fMassCytoplasm = mC;}
70 G4double GetMassCytoplasm() const {return f << 72 G4float GetMassCytoplasm(){return fMassCytoplasm;}
71 73
72 void SetNbOfPixelsInPhantom(G4int nP) {fNbOf << 74 void SetNbOfPixelsInPhantom(G4int nP){ fNbOfPixelsInPhantom = nP;}
73 G4int GetNbOfPixelsInPhantom() const {return << 75 G4int GetNbOfPixelsInPhantom(){return fNbOfPixelsInPhantom;}
74 76
75 const G4LogicalVolume* GetLogicalCollDetYoke << 77 // Returns the logical volumes
76 const G4LogicalVolume* GetLogicalIsobutane() << 78
77 const G4LogicalVolume* GetLogicalCollDetGap4 << 79 G4LogicalVolume* GetLogicalCollDetYoke() {return fLogicYoke2;};
78 const G4LogicalVolume* GetLogicalPolyprop() << 80 G4LogicalVolume* GetLogicalIsobutane() {return fLogicBoiteIso;};
79 const G4LogicalVolume* GetLogicalKgm() << 81 G4LogicalVolume* GetLogicalCollDetGap4() {return fLogic4Gap;};
80 << 82 G4LogicalVolume* GetLogicalPolyprop() {return fLogicBoite3;};
81 const G4Material * GetNucleusMaterial1() c << 83 G4LogicalVolume* GetLogicalKgm() {return fLogicKgm;};
82 const G4Material * GetNucleusMaterial2() c <<
83 const G4Material * GetNucleusMaterial3() c <<
84 const G4Material * GetCytoplasmMaterial1() c <<
85 const G4Material * GetCytoplasmMaterial2() c <<
86 const G4Material * GetCytoplasmMaterial3() c <<
87 <<
88 const CellParameterisation * GetCellParamete <<
89 {return fMyCellParameterisation;}; <<
90 <<
91 G4VPhysicalVolume* ConstructLine(); <<
92 84
93 private: 85 private:
94 << 86
95 CellParameterisation * fMyCellParameterisati << 87 G4float fMassPhantom;
96 << 88 G4float fMassNucleus;
97 // Materials << 89 G4float fMassCytoplasm;
98 G4Material* fDefaultMaterial; << 90
99 G4Material* fCollimatorMaterial; << 91 G4double fDensityPhantom;
100 G4Material* fBoiteMaterial; << 92 G4double fDensityNucleus;
101 G4Material* fCathodeMaterial; << 93 G4double fDensityCytoplasm;
102 G4Material* fVerreMaterial; << 94 G4int fNbOfPixelsInPhantom;
103 G4Material* fVerre2Material; << 95
104 G4Material* fKgmMaterial; << 96 G4double fWorldSizeXY;
105 G4Material* fBoite2Material; << 97 G4double fWorldSizeZ;
106 G4Material* fBoite3Material; << 98 G4double fCollObjSizeXY;
107 G4Material* fNucleusMaterial1; << 99 G4double fCollObjSizeZ;
108 G4Material* fCytoplasmMaterial1; << 100
109 G4Material* fNucleusMaterial2; << 101 G4double fCiblePositionX;
110 G4Material* fCytoplasmMaterial2; << 102 G4double fCiblePositionY;
111 G4Material* fNucleusMaterial3; << 103 G4double fCiblePositionZ;
112 G4Material* fCytoplasmMaterial3; << 104
113 << 105 G4double fLineAngle;
>> 106
>> 107 // Materials
>> 108
>> 109 G4Material* fDefaultMaterial;
>> 110 G4Material* fCollimatorMaterial;
>> 111 G4Material* fBoiteMaterial;
>> 112 G4Material* fCathodeMaterial;
>> 113 G4Material* fVerreMaterial;
>> 114 G4Material* fVerre2Material;
>> 115 G4Material* fKgmMaterial;
>> 116 G4Material* fBoite2Material;
>> 117 G4Material* fBoite3Material;
>> 118 G4Material* fNucleusMaterial1;
>> 119 G4Material* fCytoplasmMaterial1;
>> 120 G4Material* fNucleusMaterial2;
>> 121 G4Material* fCytoplasmMaterial2;
>> 122 G4Material* fNucleusMaterial3;
>> 123 G4Material* fCytoplasmMaterial3;
>> 124
114 // Volumes 125 // Volumes
115 126
116 G4VPhysicalVolume* fPhysiWorld; 127 G4VPhysicalVolume* fPhysiWorld;
117 G4LogicalVolume* fLogicWorld; 128 G4LogicalVolume* fLogicWorld;
118 G4Box* fSolidWorld; 129 G4Box* fSolidWorld;
119 130
120 G4VPhysicalVolume* fPhysiVol; 131 G4VPhysicalVolume* fPhysiVol;
121 G4LogicalVolume* fLogicVol; 132 G4LogicalVolume* fLogicVol;
122 G4Box* fSolidVol; 133 G4Box* fSolidVol;
123 134
124 G4VPhysicalVolume* fPhysiBoite; 135 G4VPhysicalVolume* fPhysiBoite;
125 G4LogicalVolume* fLogicBoite; 136 G4LogicalVolume* fLogicBoite;
126 G4Box* fSolidBoite; 137 G4Box* fSolidBoite;
127 138
128 G4VPhysicalVolume* fPhysiYoke1; 139 G4VPhysicalVolume* fPhysiYoke1;
129 G4LogicalVolume* fLogicYoke1; 140 G4LogicalVolume* fLogicYoke1;
130 G4Box* fSolidYoke1; 141 G4Box* fSolidYoke1;
131 142
132 G4VPhysicalVolume* fPhysi1Gap; 143 G4VPhysicalVolume* fPhysi1Gap;
133 G4LogicalVolume* fLogic1Gap; 144 G4LogicalVolume* fLogic1Gap;
134 G4Cons* fSolid1Gap; 145 G4Cons* fSolid1Gap;
135 146
136 G4VPhysicalVolume* fPhysi2Gap; 147 G4VPhysicalVolume* fPhysi2Gap;
137 G4LogicalVolume* fLogic2Gap; 148 G4LogicalVolume* fLogic2Gap;
138 G4Cons* fSolid2Gap; 149 G4Cons* fSolid2Gap;
139 150
140 G4VPhysicalVolume* fPhysi3Gap; 151 G4VPhysicalVolume* fPhysi3Gap;
141 G4LogicalVolume* fLogic3Gap; 152 G4LogicalVolume* fLogic3Gap;
142 G4Cons* fSolid3Gap; 153 G4Cons* fSolid3Gap;
143 154
144 G4VPhysicalVolume* fPhysiYoke2; 155 G4VPhysicalVolume* fPhysiYoke2;
145 G4LogicalVolume* fLogicYoke2; 156 G4LogicalVolume* fLogicYoke2;
146 G4Box* fSolidYoke2; 157 G4Box* fSolidYoke2;
147 158
148 G4VPhysicalVolume* fPhysi4Gap; 159 G4VPhysicalVolume* fPhysi4Gap;
149 G4LogicalVolume* fLogic4Gap; 160 G4LogicalVolume* fLogic4Gap;
150 G4Cons* fSolid4Gap; 161 G4Cons* fSolid4Gap;
151 162
152 G4VPhysicalVolume* fPhysi5Gap; 163 G4VPhysicalVolume* fPhysi5Gap;
153 G4LogicalVolume* fLogic5Gap; 164 G4LogicalVolume* fLogic5Gap;
154 G4Cons* fSolid5Gap; 165 G4Cons* fSolid5Gap;
155 166
156 G4VPhysicalVolume* fPhysiBoiteIso; 167 G4VPhysicalVolume* fPhysiBoiteIso;
157 G4LogicalVolume* fLogicBoiteIso; 168 G4LogicalVolume* fLogicBoiteIso;
158 G4Box* fSolidBoiteIso; 169 G4Box* fSolidBoiteIso;
159 170
160 G4VPhysicalVolume* fPhysiCathode; 171 G4VPhysicalVolume* fPhysiCathode;
161 G4LogicalVolume* fLogicCathode; 172 G4LogicalVolume* fLogicCathode;
162 G4Box* fSolidCathode; 173 G4Box* fSolidCathode;
163 174
164 G4VPhysicalVolume* fPhysiIso; 175 G4VPhysicalVolume* fPhysiIso;
165 G4LogicalVolume* fLogicIso; 176 G4LogicalVolume* fLogicIso;
166 G4Box* fSolidIso; 177 G4Box* fSolidIso;
167 178
168 G4VPhysicalVolume* fPhysiVerre; 179 G4VPhysicalVolume* fPhysiVerre;
169 G4LogicalVolume* fLogicVerre; 180 G4LogicalVolume* fLogicVerre;
170 G4Box* fSolidVerre; 181 G4Box* fSolidVerre;
171 182
172 G4VPhysicalVolume* fPhysiBoite2; 183 G4VPhysicalVolume* fPhysiBoite2;
173 G4LogicalVolume* fLogicBoite2; 184 G4LogicalVolume* fLogicBoite2;
174 G4Box* fSolidBoite2; 185 G4Box* fSolidBoite2;
175 186
176 G4VPhysicalVolume* fPhysiBoite3; 187 G4VPhysicalVolume* fPhysiBoite3;
177 G4LogicalVolume* fLogicBoite3; 188 G4LogicalVolume* fLogicBoite3;
178 G4Box* fSolidBoite3; 189 G4Box* fSolidBoite3;
179 190
180 G4VPhysicalVolume* fPhysiKgm; 191 G4VPhysicalVolume* fPhysiKgm;
181 G4LogicalVolume* fLogicKgm; 192 G4LogicalVolume* fLogicKgm;
182 G4Box* fSolidKgm; 193 G4Box* fSolidKgm;
183 194
184 G4VPhysicalVolume* fPhysiVerre2; 195 G4VPhysicalVolume* fPhysiVerre2;
185 G4LogicalVolume* fLogicVerre2; 196 G4LogicalVolume* fLogicVerre2;
186 G4Box* fSolidVerre2; 197 G4Box* fSolidVerre2;
187 198
188 // CELL PHANTOM 199 // CELL PHANTOM
189 200
190 G4VPhysicalVolume* fPhysiPhantom; 201 G4VPhysicalVolume* fPhysiPhantom;
191 G4LogicalVolume* fLogicPhantom; 202 G4LogicalVolume* fLogicPhantom;
192 G4Box* fSolidPhantom; 203 G4Box* fSolidPhantom;
193 204
194 G4double fMassNucleus; << 205 //
195 G4double fMassCytoplasm; <<
196 206
197 G4double fDensityPhantom; << 207 PhantomConfiguration fMyPhantomConfiguration;
198 G4double fDensityNucleus; << 208 CellParameterisation* fPhantomParam ;
199 G4double fDensityCytoplasm; <<
200 <<
201 G4double fWorldSizeXY; <<
202 G4double fWorldSizeZ; <<
203 G4double fCollObjSizeXY; <<
204 G4double fCollObjSizeZ; <<
205 <<
206 G4double fCiblePositionX; <<
207 G4double fCiblePositionY; <<
208 G4double fCiblePositionZ; <<
209 209
210 G4double fLineAngle; << 210 static G4ThreadLocal EMField * fField;
>> 211
>> 212 G4FieldManager * fFieldMgr;
>> 213 G4MagIntegratorStepper * fStepper;
>> 214 G4EqMagElectricField * fEquation;
>> 215 G4MagInt_Driver * fIntgrDriver;
>> 216 G4ChordFinder * fChordFinder ;
>> 217 G4PropagatorInField * fPropInField;
>> 218
>> 219 //
211 220
212 G4int fNbOfPixelsInPhantom; << 221 void DefineMaterials();
>> 222 G4VPhysicalVolume* ConstructLine();
>> 223
213 }; 224 };
214 225
215 #endif 226 #endif
216 227