source: trunk/examples/advanced/gammaray_telescope/include/GammaRayTelDetectorConstruction.hh@ 850

Last change on this file since 850 was 807, checked in by garnier, 17 years ago

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