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