Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
GammaRayTelDetectorConstruction.hh
Go to the documentation of this file.
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 //
27 // $Id$
28 // ------------------------------------------------------------
29 // GEANT 4 class header file
30 // CERN Geneva Switzerland
31 //
32 //
33 // ------------ GammaRayTelDetectorConstruction ------
34 // by F.Longo, R.Giannitrapani & G.Santin (13 nov 2000)
35 //
36 // ************************************************************
37 
38 #ifndef GammaRayTelDetectorConstruction_h
39 #define GammaRayTelDetectorConstruction_h 1
40 
41 #include "G4VUserDetectorConstruction.hh"
42 #include "globals.hh"
43 
44 class G4Box;
45 class G4LogicalVolume;
46 class G4VPhysicalVolume;
47 class G4Material;
48 class G4UniformMagField;
49 class G4Region;
50 class GammaRayTelDetectorMessenger;
51 class GammaRayTelTrackerSD;
52 class GammaRayTelAnticoincidenceSD;
53 class GammaRayTelCalorimeterSD;
54 class GammaRayTelTrackerROGeometry;
55 
56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
57 
58 class GammaRayTelDetectorConstruction : public G4VUserDetectorConstruction
59 {
60 public:
61 
62  GammaRayTelDetectorConstruction();
63  ~GammaRayTelDetectorConstruction();
64 
65 public:
66 
67  void SetNbOfTKRLayers (G4int); // TKR number of layers, material, detector
68  void SetTKRTileSizeXY (G4double);
69  void SetNbOfTKRTiles (G4int);
70  void SetTKRSiliconThickness(G4double);
71  void SetTKRSiliconPitch(G4double);
72 
73  void SetTKRLayerDistance (G4double);
74  void SetTKRViewsDistance (G4double);
75 
76  void SetConverterMaterial (G4String); // TKR Converter material & thickness
77  void SetConverterThickness(G4double);
78 
79  void SetNbOfCALLayers (G4int); // CAL material, lenght, thickness
80  void SetNbOfCALBars (G4int);
81  void SetCALBarThickness(G4double);
82 
83  void SetACDThickness (G4double); //ACD Thickness
84 
85  void SetMagField(G4double); // Magnetic Field
86 
87 
88  G4VPhysicalVolume* Construct();
89  void UpdateGeometry();
90 
91 public:
92 
93  void PrintPayloadParameters();
94 
95 
96  G4double GetWorldSizeZ() {return WorldSizeZ;};
97  G4double GetWorldSizeXY() {return WorldSizeXY;};
98 
99  G4double GetPayloadSizeZ() {return PayloadSizeZ;};
100  G4double GetPayloadSizeXY() {return PayloadSizeXY;};
101 
102  G4double GetTKRSizeZ() {return TKRSizeZ;};
103  G4double GetTKRSizeXY() {return TKRSizeXY;};
104 
105  G4double GetCALSizeZ() {return CALSizeZ;};
106  G4double GetCALTKRDistance() {return CALTKRDistance;};
107 
108  G4double GetTKRSiliconThickness() {return TKRSiliconThickness;};
109  G4double GetTKRSiliconTileXY() {return TKRSiliconTileXY;};
110  G4double GetTKRSiliconPitch() {return TKRSiliconPitch;};
111  G4int GetNbOfTKRLayers() {return NbOfTKRLayers;};
112  G4int GetNbOfTKRTiles() {return NbOfTKRTiles;};
113  G4int GetNbOfTKRStrips() {return NbOfTKRStrips;};
114  G4double GetTKRLayerDistance() {return TKRLayerDistance;};
115  G4double GetTKRViewsDistance() {return TKRViewsDistance;};
116 
117  G4double GetTKRActiveTileXY() {return TKRActiveTileXY;};
118  G4double GetTKRActiveTileZ() {return TKRActiveTileZ;};
119  G4double GetSiliconGuardRing() {return SiliconGuardRing;}
120  G4double GetTilesSeparation() {return TilesSeparation;};
121 
122  G4Material* GetConverterMaterial() {return ConverterMaterial;};
123  G4double GetConverterThickness() {return ConverterThickness;};
124 
125  G4double GetCALBarThickness() {return CALBarThickness;};
126  G4int GetNbOfCALLayers() {return NbOfCALLayers;};
127  G4int GetNbOfCALBars() {return NbOfCALBars;};
128 
129  G4double GetACDThickness() {return ACDThickness;};
130  G4int GetNbOfACDTopTiles() {return NbOfACDTopTiles;};
131  G4int GetNbOfACDLateralTiles() {return NbOfACDLateralTiles;};
132 
133 private:
134 
135  G4Material* ConverterMaterial;
136  G4double ConverterThickness;
137 
138  G4double TKRSiliconThickness;
139  G4double TKRSiliconTileXY;
140  G4double TKRSiliconPitch;
141 
142  G4double TKRSizeXY;
143  G4double TKRSizeZ;
144 
145  G4double TKRLayerDistance;
146  G4double TKRViewsDistance;
147  G4double TKRSupportThickness;
148 
149  G4int NbOfTKRLayers;
150  G4int NbOfTKRTiles;
151 
152  G4double CALBarThickness;
153  G4int NbOfCALLayers;
154  G4int NbOfCALBars;
155  G4double CALSizeXY;
156  G4double CALSizeZ;
157 
158  G4double CALBarX;
159  G4double CALBarY;
160  G4double CALBarZ;
161 
162  G4double ACDThickness;
163  G4double ACTSizeXY;
164  G4double ACTSizeZ;
165 
166  G4double ACL1SizeX;
167  G4double ACL1SizeY;
168  G4double ACL1SizeZ;
169 
170  G4double ACL2SizeX;
171  G4double ACL2SizeY;
172  G4double ACL2SizeZ;
173 
174 
175  G4int NbOfACDLateralTiles;
176  G4int NbOfACDTopTiles;
177 
178  G4double TilesSeparation;
179  G4double ACDTKRDistance;
180  G4double CALTKRDistance;
181  G4double TKRActiveTileXY;
182  G4double TKRActiveTileZ;
183 
184  G4double SiliconGuardRing;
185  G4int NbOfTKRStrips;
186 
187  G4double TKRXStripX;
188  G4double TKRYStripX;
189  G4double TKRXStripY;
190  G4double TKRYStripY;
191  G4double TKRZStrip;
192 
193  G4double PayloadSizeZ;
194  G4double PayloadSizeXY;
195 
196  G4Material* defaultMaterial;
197  G4Material* CALMaterial;
198  G4Material* TKRMaterial;
199  G4Material* ACDMaterial;
200  G4double WorldSizeXY;
201  G4double WorldSizeZ;
202 
203  G4Box* solidWorld; // World
204  G4LogicalVolume* logicWorld;
205  G4VPhysicalVolume* physiWorld;
206 
207  G4Box* solidPayload; // Payload
208  G4LogicalVolume* logicPayload;
209  G4VPhysicalVolume* physiPayload;
210 
211  G4Box* solidTKR; // Tracker
212  G4LogicalVolume* logicTKR;
213  G4VPhysicalVolume* physiTKR;
214 
215  G4Box* solidCAL; // Calorimeter
216  G4LogicalVolume* logicCAL;
217  G4VPhysicalVolume* physiCAL;
218 
219  G4Box* solidACT; // Top Anticoincidence
220  G4LogicalVolume* logicACT;
221  G4VPhysicalVolume* physiACT;
222 
223  G4Box* solidACL1; // Lateral Anticoincidence
224  G4LogicalVolume* logicACL1;
225  G4VPhysicalVolume* physiACL1;
226 
227  G4Box* solidACL2;
228  G4LogicalVolume* logicACL2;
229  G4VPhysicalVolume* physiACL2;
230 
231  G4Box* solidTKRDetectorX; // Tracker PLANE X
232  G4LogicalVolume* logicTKRDetectorX;
233  G4VPhysicalVolume* physiTKRDetectorX;
234 
235  G4Box* solidTKRDetectorY; // Tracker PLANE Y
236  G4LogicalVolume* logicTKRDetectorY;
237  G4VPhysicalVolume* physiTKRDetectorY;
238 
239  G4Box* solidCALLayerX; // Calorimeter PLANE X
240  G4LogicalVolume* logicCALLayerX;
241  G4VPhysicalVolume* physiCALLayerX;
242 
243  G4Box* solidCALLayerY; // Calorimeter PLANE Y
244  G4LogicalVolume* logicCALLayerY;
245  G4VPhysicalVolume* physiCALLayerY;
246 
247  G4Box* solidCALDetectorX; // Calorimeter DETECTOR X
248  G4LogicalVolume* logicCALDetectorX;
249  G4VPhysicalVolume* physiCALDetectorX;
250 
251  G4Box* solidCALDetectorY; // Calorimeter DETECTOR Y
252  G4LogicalVolume* logicCALDetectorY;
253  G4VPhysicalVolume* physiCALDetectorY;
254 
255  G4Box* solidPlane; // Support Plane
256  G4LogicalVolume* logicPlane;
257  G4VPhysicalVolume* physiPlane;
258 
259  G4Box* solidConverter; // Converter
260  G4LogicalVolume* logicConverter;
261  G4VPhysicalVolume* physiConverter;
262 
263  G4UniformMagField* magField; //pointer to the magnetic field
264 
265  GammaRayTelDetectorMessenger* detectorMessenger; //pointer to the Messenger
266  GammaRayTelTrackerSD* trackerSD; //pointer to the sensitive detector
267  GammaRayTelCalorimeterSD* calorimeterSD; //pointer to the sensitive detector
268  GammaRayTelAnticoincidenceSD* anticoincidenceSD; //pointer to the sensitive detector
269 
270  //G4Region* aTKRRegion; // TKR cut region
271  //G4Region* aCALRegion; // CAL cut region
272 
273 private:
274 
275  void DefineMaterials();
276  void ComputePayloadParameters();
277  G4VPhysicalVolume* ConstructPayload();
278 };
279 
280 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
281 
282 inline void GammaRayTelDetectorConstruction::ComputePayloadParameters()
283 {
284  // Compute derived parameters of the payload
285 
286  TKRSupportThickness =TKRLayerDistance -2.*TKRSiliconThickness
287  - TKRViewsDistance - ConverterThickness;
288  TKRSizeXY = NbOfTKRTiles*TKRSiliconTileXY + (NbOfTKRTiles+1)*TilesSeparation;
289  TKRSizeZ = NbOfTKRLayers*TKRLayerDistance;
290 
291  TKRActiveTileXY = TKRSiliconTileXY - 2*SiliconGuardRing;
292  TKRActiveTileZ = TKRSiliconThickness;
293  NbOfTKRStrips = G4int(TKRActiveTileXY/TKRSiliconPitch);
294 
295  SiliconGuardRing = TKRActiveTileXY - NbOfTKRStrips*TKRSiliconPitch;
296  TKRActiveTileXY = TKRSiliconTileXY - 2*SiliconGuardRing;
297 
298  TKRXStripX = TKRYStripY = TKRSiliconPitch;
299  TKRYStripX = TKRXStripY = TKRActiveTileXY;
300  TKRZStrip = TKRSiliconThickness;
301 
302  CALSizeXY = TKRSizeXY;
303  CALSizeZ = 2.*NbOfCALLayers*CALBarThickness;
304 
305  CALBarX = CALSizeXY;
306  CALBarY = CALSizeXY/(NbOfCALBars);
307  CALBarZ = CALBarThickness;
308 
309  ACTSizeXY = TKRSizeXY + 2*ACDTKRDistance + 2*ACDThickness;
310  ACTSizeZ = ACDThickness;
311 
312  ACL1SizeX = TKRSizeXY + 2*ACDTKRDistance + ACDThickness;
313  ACL1SizeY = ACDThickness;
314  ACL1SizeZ = TKRSizeZ + CALSizeZ + ACDTKRDistance + CALTKRDistance;
315 
316  ACL2SizeX = ACDThickness;
317  ACL2SizeY = TKRSizeXY + 2*ACDTKRDistance + ACDThickness;
318  ACL2SizeZ = TKRSizeZ + CALSizeZ + ACDTKRDistance + CALTKRDistance;
319 
320  PayloadSizeZ = 1.1*(ACL1SizeZ + ACTSizeZ);
321  PayloadSizeXY = (ACTSizeXY);
322 
323  WorldSizeZ = 1.5*PayloadSizeZ; WorldSizeXY = 1.5*PayloadSizeXY;
324 
325 }
326 
327 #endif
328 
329 
330 
331 
332 
333 
334 
335 
336 
Generated on Sat May 25 2013 14:32:09 for Geant4 by   doxygen 1.8.4