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 virtual ~DetectorConstruction(); 59 62 60 G4VPhysicalVolume* Construct() override; << 63 virtual G4VPhysicalVolume* Construct(); 61 << 62 void DefineMaterials(); << 63 64 64 void ConstructSDandField() override; << 65 virtual 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() const {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() const {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() const {return fNbOfPixelsInPhantom;} 74 75 75 const G4LogicalVolume* GetLogicalCollDetYoke 76 const G4LogicalVolume* GetLogicalCollDetYoke() const {return fLogicYoke2;}; 76 const G4LogicalVolume* GetLogicalIsobutane() 77 const G4LogicalVolume* GetLogicalIsobutane() const {return fLogicBoiteIso;}; 77 const G4LogicalVolume* GetLogicalCollDetGap4 78 const G4LogicalVolume* GetLogicalCollDetGap4() const {return fLogic4Gap;}; 78 const G4LogicalVolume* GetLogicalPolyprop() 79 const G4LogicalVolume* GetLogicalPolyprop() const {return fLogicBoite3;}; 79 const G4LogicalVolume* GetLogicalKgm() 80 const G4LogicalVolume* GetLogicalKgm() const {return fLogicKgm;}; 80 81 81 const G4Material * GetNucleusMaterial1() c 82 const G4Material * GetNucleusMaterial1() const {return fNucleusMaterial1;}; 82 const G4Material * GetNucleusMaterial2() c 83 const G4Material * GetNucleusMaterial2() const {return fNucleusMaterial2;}; 83 const G4Material * GetNucleusMaterial3() c 84 const G4Material * GetNucleusMaterial3() const {return fNucleusMaterial3;}; 84 const G4Material * GetCytoplasmMaterial1() c 85 const G4Material * GetCytoplasmMaterial1() const {return fCytoplasmMaterial1;}; 85 const G4Material * GetCytoplasmMaterial2() c 86 const G4Material * GetCytoplasmMaterial2() const {return fCytoplasmMaterial2;}; 86 const G4Material * GetCytoplasmMaterial3() c 87 const G4Material * GetCytoplasmMaterial3() const {return fCytoplasmMaterial3;}; 87 88 88 const CellParameterisation * GetCellParamete 89 const CellParameterisation * GetCellParameterisation() const 89 {return fMyCellParameterisation;}; 90 {return fMyCellParameterisation;}; 90 << 91 G4VPhysicalVolume* ConstructLine(); << 92 91 93 private: 92 private: >> 93 >> 94 G4double fMassNucleus; >> 95 G4double fMassCytoplasm; >> 96 >> 97 G4double fDensityPhantom; >> 98 G4double fDensityNucleus; >> 99 G4double fDensityCytoplasm; >> 100 G4int fNbOfPixelsInPhantom; >> 101 >> 102 G4double fWorldSizeXY; >> 103 G4double fWorldSizeZ; >> 104 G4double fCollObjSizeXY; >> 105 G4double fCollObjSizeZ; >> 106 >> 107 G4double fCiblePositionX; >> 108 G4double fCiblePositionY; >> 109 G4double fCiblePositionZ; 94 110 95 CellParameterisation * fMyCellParameterisati << 111 G4double fLineAngle; 96 << 112 97 // Materials << 113 // Materials >> 114 98 G4Material* fDefaultMaterial; 115 G4Material* fDefaultMaterial; 99 G4Material* fCollimatorMaterial; 116 G4Material* fCollimatorMaterial; 100 G4Material* fBoiteMaterial; 117 G4Material* fBoiteMaterial; 101 G4Material* fCathodeMaterial; 118 G4Material* fCathodeMaterial; 102 G4Material* fVerreMaterial; 119 G4Material* fVerreMaterial; 103 G4Material* fVerre2Material; 120 G4Material* fVerre2Material; 104 G4Material* fKgmMaterial; 121 G4Material* fKgmMaterial; 105 G4Material* fBoite2Material; 122 G4Material* fBoite2Material; 106 G4Material* fBoite3Material; 123 G4Material* fBoite3Material; 107 G4Material* fNucleusMaterial1; 124 G4Material* fNucleusMaterial1; 108 G4Material* fCytoplasmMaterial1; 125 G4Material* fCytoplasmMaterial1; 109 G4Material* fNucleusMaterial2; 126 G4Material* fNucleusMaterial2; 110 G4Material* fCytoplasmMaterial2; 127 G4Material* fCytoplasmMaterial2; 111 G4Material* fNucleusMaterial3; 128 G4Material* fNucleusMaterial3; 112 G4Material* fCytoplasmMaterial3; 129 G4Material* fCytoplasmMaterial3; 113 << 130 114 // Volumes 131 // Volumes 115 132 116 G4VPhysicalVolume* fPhysiWorld; 133 G4VPhysicalVolume* fPhysiWorld; 117 G4LogicalVolume* fLogicWorld; 134 G4LogicalVolume* fLogicWorld; 118 G4Box* fSolidWorld; 135 G4Box* fSolidWorld; 119 136 120 G4VPhysicalVolume* fPhysiVol; 137 G4VPhysicalVolume* fPhysiVol; 121 G4LogicalVolume* fLogicVol; 138 G4LogicalVolume* fLogicVol; 122 G4Box* fSolidVol; 139 G4Box* fSolidVol; 123 140 124 G4VPhysicalVolume* fPhysiBoite; 141 G4VPhysicalVolume* fPhysiBoite; 125 G4LogicalVolume* fLogicBoite; 142 G4LogicalVolume* fLogicBoite; 126 G4Box* fSolidBoite; 143 G4Box* fSolidBoite; 127 144 128 G4VPhysicalVolume* fPhysiYoke1; 145 G4VPhysicalVolume* fPhysiYoke1; 129 G4LogicalVolume* fLogicYoke1; 146 G4LogicalVolume* fLogicYoke1; 130 G4Box* fSolidYoke1; 147 G4Box* fSolidYoke1; 131 148 132 G4VPhysicalVolume* fPhysi1Gap; 149 G4VPhysicalVolume* fPhysi1Gap; 133 G4LogicalVolume* fLogic1Gap; 150 G4LogicalVolume* fLogic1Gap; 134 G4Cons* fSolid1Gap; 151 G4Cons* fSolid1Gap; 135 152 136 G4VPhysicalVolume* fPhysi2Gap; 153 G4VPhysicalVolume* fPhysi2Gap; 137 G4LogicalVolume* fLogic2Gap; 154 G4LogicalVolume* fLogic2Gap; 138 G4Cons* fSolid2Gap; 155 G4Cons* fSolid2Gap; 139 156 140 G4VPhysicalVolume* fPhysi3Gap; 157 G4VPhysicalVolume* fPhysi3Gap; 141 G4LogicalVolume* fLogic3Gap; 158 G4LogicalVolume* fLogic3Gap; 142 G4Cons* fSolid3Gap; 159 G4Cons* fSolid3Gap; 143 160 144 G4VPhysicalVolume* fPhysiYoke2; 161 G4VPhysicalVolume* fPhysiYoke2; 145 G4LogicalVolume* fLogicYoke2; 162 G4LogicalVolume* fLogicYoke2; 146 G4Box* fSolidYoke2; 163 G4Box* fSolidYoke2; 147 164 148 G4VPhysicalVolume* fPhysi4Gap; 165 G4VPhysicalVolume* fPhysi4Gap; 149 G4LogicalVolume* fLogic4Gap; 166 G4LogicalVolume* fLogic4Gap; 150 G4Cons* fSolid4Gap; 167 G4Cons* fSolid4Gap; 151 168 152 G4VPhysicalVolume* fPhysi5Gap; 169 G4VPhysicalVolume* fPhysi5Gap; 153 G4LogicalVolume* fLogic5Gap; 170 G4LogicalVolume* fLogic5Gap; 154 G4Cons* fSolid5Gap; 171 G4Cons* fSolid5Gap; 155 172 156 G4VPhysicalVolume* fPhysiBoiteIso; 173 G4VPhysicalVolume* fPhysiBoiteIso; 157 G4LogicalVolume* fLogicBoiteIso; 174 G4LogicalVolume* fLogicBoiteIso; 158 G4Box* fSolidBoiteIso; 175 G4Box* fSolidBoiteIso; 159 176 160 G4VPhysicalVolume* fPhysiCathode; 177 G4VPhysicalVolume* fPhysiCathode; 161 G4LogicalVolume* fLogicCathode; 178 G4LogicalVolume* fLogicCathode; 162 G4Box* fSolidCathode; 179 G4Box* fSolidCathode; 163 180 164 G4VPhysicalVolume* fPhysiIso; 181 G4VPhysicalVolume* fPhysiIso; 165 G4LogicalVolume* fLogicIso; 182 G4LogicalVolume* fLogicIso; 166 G4Box* fSolidIso; 183 G4Box* fSolidIso; 167 184 168 G4VPhysicalVolume* fPhysiVerre; 185 G4VPhysicalVolume* fPhysiVerre; 169 G4LogicalVolume* fLogicVerre; 186 G4LogicalVolume* fLogicVerre; 170 G4Box* fSolidVerre; 187 G4Box* fSolidVerre; 171 188 172 G4VPhysicalVolume* fPhysiBoite2; 189 G4VPhysicalVolume* fPhysiBoite2; 173 G4LogicalVolume* fLogicBoite2; 190 G4LogicalVolume* fLogicBoite2; 174 G4Box* fSolidBoite2; 191 G4Box* fSolidBoite2; 175 192 176 G4VPhysicalVolume* fPhysiBoite3; 193 G4VPhysicalVolume* fPhysiBoite3; 177 G4LogicalVolume* fLogicBoite3; 194 G4LogicalVolume* fLogicBoite3; 178 G4Box* fSolidBoite3; 195 G4Box* fSolidBoite3; 179 196 180 G4VPhysicalVolume* fPhysiKgm; 197 G4VPhysicalVolume* fPhysiKgm; 181 G4LogicalVolume* fLogicKgm; 198 G4LogicalVolume* fLogicKgm; 182 G4Box* fSolidKgm; 199 G4Box* fSolidKgm; 183 200 184 G4VPhysicalVolume* fPhysiVerre2; 201 G4VPhysicalVolume* fPhysiVerre2; 185 G4LogicalVolume* fLogicVerre2; 202 G4LogicalVolume* fLogicVerre2; 186 G4Box* fSolidVerre2; 203 G4Box* fSolidVerre2; 187 204 188 // CELL PHANTOM 205 // CELL PHANTOM 189 206 190 G4VPhysicalVolume* fPhysiPhantom; 207 G4VPhysicalVolume* fPhysiPhantom; 191 G4LogicalVolume* fLogicPhantom; 208 G4LogicalVolume* fLogicPhantom; 192 G4Box* fSolidPhantom; 209 G4Box* fSolidPhantom; 193 210 194 G4double fMassNucleus; << 211 CellParameterisation * fMyCellParameterisation; 195 G4double fMassCytoplasm; << 196 << 197 G4double fDensityPhantom; << 198 G4double fDensityNucleus; << 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 212 210 G4double fLineAngle; << 213 // EM FIELD 211 214 212 G4int fNbOfPixelsInPhantom; << 215 static G4ThreadLocal EMField * fField; >> 216 >> 217 void DefineMaterials(); >> 218 G4VPhysicalVolume* ConstructLine(); >> 219 213 }; 220 }; 214 221 215 #endif 222 #endif 216 223