ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
DetectorConstruction.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file DetectorConstruction.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 // 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 // and papers
31 // M. Batmunkh et al. J Radiat Res Appl Sci 8 (2015) 498-507
32 // O. Belov et al. Physica Medica 32 (2016) 1510-1520
33 // The Geant4-DNA web site is available at http://geant4-dna.org
34 //
35 // -------------------------------------------------------------------
36 // November 2016
37 // -------------------------------------------------------------------
38 //
41 
42 #ifndef DetectorConstruction_h
43 #define DetectorConstruction_h 1
44 
45 #include "G4VUserDetectorConstruction.hh"
46 #include "G4VSolid.hh"
47 #include "G4UnionSolid.hh"
48 #include "G4VPhysicalVolume.hh"
49 #include "G4LogicalVolume.hh"
50 #include "G4Box.hh"
51 #include "G4Sphere.hh"
52 #include "G4Material.hh"
53 #include "G4NistManager.hh"
54 #include "G4PVPlacement.hh"
55 #include "G4UserLimits.hh"
56 #include "G4VisAttributes.hh"
57 //
58 #include "G4PVParameterised.hh"
59 #include "NeuronLoadDataFile.hh"
60 #include "G4Tubs.hh"
61 #include "G4Orb.hh"
62 #include "G4Ellipsoid.hh"
63 #include "G4LogicalVolumeStore.hh"
64 #include "G4PhysicalVolumeStore.hh"
65 #include "G4SolidStore.hh"
66 
67 const G4int kMT = 100000;
68 
69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
70 
71 class G4Region;
72 
73 class DetectorConstruction : public G4VUserDetectorConstruction
74 {
75 public:
76 
77  DetectorConstruction();
78  virtual ~DetectorConstruction();
79  virtual G4VPhysicalVolume* Construct();
80  G4Region* GetTargetRegion() {return fpRegion;}
81 
82 // For neuron constructions!
83  NeuronLoadDataFile * GetNeuronLoadParamz() {return fNeuronLoadParamz;};
84  const G4VPhysicalVolume* GetsomaPV(G4int i) {return fsomaPV[i];};
85  const G4VPhysicalVolume* GetdendritePV(G4int i) {return fdendritePV[i];};
86  const G4VPhysicalVolume* GetaxonPV(G4int i) {return faxonPV[i];};
87  G4Material* GetTargetMaterial() {return fpWaterMaterial;};
88 
89 private:
90  G4Material* fpDefaultMaterial;
91  G4Material* fpWaterMaterial;
92  G4Region* fpRegion;
93 
94  void DefineMaterials();
95  G4VPhysicalVolume* ConstructDetector();
96 
97  // For neuron constructions!
98  NeuronLoadDataFile * fNeuronLoadParamz;
99  G4bool fCheckOverlaps; // option to activate checking of volumes overlaps
100  G4VSolid * fsomaS[kMT], *fdendriteS[kMT], *faxonS[kMT];
101  G4LogicalVolume * fsomaLV[kMT], *fdendriteLV[kMT], *faxonLV[kMT];
102  G4VPhysicalVolume * fsomaPV[kMT], *fdendritePV[kMT], *faxonPV[kMT];
103 
104  G4VisAttributes * fSomaColour;
105  G4VisAttributes * fDendColour;
106  G4VisAttributes * fAxonColour;
107  G4VisAttributes * fSpineColour;
108  G4VisAttributes * fNeuronColour;
109 
110 public:
111 
112  G4int GetnbSomacomp() {return fnbSomacomp;} //
113  G4double GetMassSomacomp (G4int i) {return fMassSomacomp[i];} //
114  G4double GetMassSomaTot () {return fMassSomaTot;} //
115  G4ThreeVector GetPosSomacomp(G4int i) {return fPosSomacomp[i];} //
116 
117  G4int GetnbDendritecomp() {return fnbDendritecomp;} //
118  G4double GetMassDendcomp (G4int i) {return fMassDendcomp[i];} //
119  G4double GetMassDendTot () {return fMassDendTot;} //
120  G4ThreeVector GetPosDendcomp(G4int i) {return fPosDendcomp[i];} //
121  G4double GetDistADendSoma (G4int i) {return fDistADendSoma[i];} //
122  G4double GetDistBDendSoma (G4int i) {return fDistBDendSoma[i];} //
123 
124  G4int GetnbAxoncomp() {return fnbAxoncomp;} //
125  G4double GetMassAxoncomp (G4int i) {return fMassAxoncomp[i];} //
126  G4double GetMassAxonTot () {return fMassAxonTot;} //
127  G4ThreeVector GetPosAxoncomp(G4int i) {return fPosAxoncomp[i];} //
128  G4double GetDistAxonsoma (G4int i) {return fDistAxonsoma[i];} //
129 
130  G4double GetTotVolNeuron() {return fTotVolNeuron;}
131  G4double GetTotSurfNeuron () {return fTotSurfNeuron;}
132  G4double GetTotMassNeuron () {return fTotMassNeuron;} //
133  G4double GetTotVolSlice() {return fTotVolSlice;}
134  G4double GetTotSurfSlice () {return fTotSurfSlice;}
135  G4double GetTotMassSlice() {return fTotMassSlice;} //
136  G4double GetTotVolMedium() {return fTotVolMedium;}
137  G4double GetTotSurfMedium() {return fTotSurfMedium;} //
138  G4double GetTotMassMedium() {return fTotMassMedium;} //
139 
140 private:
141 
142  G4int fnbSomacomp;
143  G4int fnbDendritecomp;
144  G4int fnbAxoncomp;
145 
146  G4ThreeVector * fPosSomacomp ;
147  G4double * fMassSomacomp ;
148  G4double fMassSomaTot ;
149 
150  G4double * fDistADendSoma ;
151  G4double * fDistBDendSoma ;
152  G4double * fMassDendcomp ;
153  G4double fMassDendTot ;
154  G4ThreeVector * fPosDendcomp ;
155 
156  G4double * fDistAxonsoma ;
157  G4double * fMassAxoncomp ;
158  G4double fMassAxonTot ;
159  G4ThreeVector * fPosAxoncomp ;
160 
161  G4double fTotVolNeuron ;
162  G4double fTotSurfNeuron ;
163  G4double fTotMassNeuron ;
164  G4double fTotVolSlice;
165  G4double fTotSurfSlice ;
166  G4double fTotMassSlice;
167  G4double fTotVolMedium ;
168  G4double fTotSurfMedium ;
169  G4double fTotMassMedium ;
170 
171 };
172 #endif