source: trunk/examples/advanced/cosmicray_charging/src/LISAScienceModuleStructures.icc

Last change on this file was 807, checked in by garnier, 16 years ago

update
File size: 25.2 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// *                                                                  *
28// * cosmicray_charging advanced example for Geant4                   *
29// * (adapted simulation of test-mass charging in the LISA mission)   *
30// *                                                                  *
31// * Henrique Araujo (h.araujo@imperial.ac.uk) & Peter Wass           *
32// * Imperial College London                                          *
33// *                                                                  *
34// * LISADetectorConstruction class                                   *
35// *                                                                  *
36// ********************************************************************
37//
38// HISTORY
39// 22/02/2004: migrated from LISA-V04
40// 24/11/2004: migrated to cmath
41// 08/12/2005: removed compilation warnings
42//
43// ********************************************************************
44
45
46//**************************************************************************
47// Science Module Structure (SMS)
48//**************************************************************************
49//    Primary Structure
50//    Lower Deck
51//    Upper Deck
52//    Thermal Shield
53//    Solar Array
54//    Optical Solar Reflectors (OSR)
55//    Radiator Panels
56//**************************************************************************
57
58
59  // Structure Tubes
60  G4double StructTube_dia       =  100.*mm;
61  G4double StructTube_thi       =   11.*mm;
62  G4double StructTube_hei       =  490.*mm;
63  G4double StructTube_sep       = 1480.*mm;
64  G4double StructTubeSpacer_dia =  200.*mm;
65  G4double StructTubeSpacer_thi =   50.*mm;
66  G4double StructTubeSpacer_hei =   16.*mm;
67
68  // Diagonal Panels
69  G4double DiagPanelSkin_thi   = 0.5*mm;
70  G4double DiagPanelSkin_len   = 1420.*mm;
71  G4double DiagPanelSkin_hei   = 490.*mm;
72  G4double DiagPanelCore_thi   = 19.*mm;
73  G4double DiagPanelCore_len   = DiagPanelSkin_len;
74  G4double DiagPanelCore_hei   = DiagPanelSkin_hei;
75  G4double DiagPanelCutout_dia = 466.*mm;
76  G4double DiagPanelCutout_off = 50.*mm;
77
78  // Radial Panels
79  G4double RadialPanelSkin_thi = 0.5*mm;
80  G4double RadialPanelSkin_hei = 490.*mm;
81  G4double RadialPanelSkin_top = 424.*mm;
82  G4double RadialPanelSkin_bot = 159.*mm;
83  G4double RadialPanelCore_thi =  19.*mm;
84  G4double RadialPanel_roff    = 910.*mm;
85
86  // Upper and Lower Decks
87  G4double UpperDeck_dia     = 2700.*mm;
88  G4double UpperDeckCore_thi =  29.0*mm;
89  G4double LowerDeckCore_thi =  29.0*mm;
90  G4double UpperDeck_zoff    =   0.5*mm;
91  G4double UpperDeckSkin_thi =   0.5*mm;
92  G4double LowerDeck_dia     = 2174.*mm;
93  G4double LowerDeckSkin_thi =   0.5*mm;
94  G4double LowerDeck_zoff    =   0.5*mm;
95  G4double MLIBlanket_dia    = LowerDeck_dia;
96  G4double MLIBlanket_thi    =   1.0*mm;
97
98  // Thermal Shield
99  G4double ThermalShieldFacesheet_thi = 0.6*mm;
100  G4double ThermalShieldCore_thi      = 20.*mm;
101  G4double ThermalShieldFoam_thi      = 20.*mm;
102  G4double ThermalShield_zoff         = StructTubeSpacer_hei;
103  G4double ThermalShield_thi          = 3*ThermalShieldFacesheet_thi +
104                                        ThermalShieldCore_thi +
105                                        ThermalShieldFoam_thi;
106  // Solar Array and OSR Assy
107  G4double SolarArray_i   = 0.5*1683*mm;
108  G4double SolarArray_o   = 0.5*UpperDeck_dia;
109  G4double SolarArray_thi = 0.5*mm;
110  G4double OSRAssy_o      = SolarArray_i;
111  G4double OSRAssy_thi    = 0.5*mm;
112
113  // Radiator Panels
114  G4double RadPanel_top  = UpperDeck_dia+12.*mm;
115  G4double RadPanel_bot  = LowerDeck_dia+12.*mm;
116  G4double RadPanel_len  = 510.*mm;
117  G4double RadPanel_thi  = 1.25*mm;
118  // G4double RadPanel_zoff = 20*mm;
119  G4double RadPanelCutout_roff  = 0.5*(RadPanel_top+RadPanel_bot)/2.;
120  G4double RadPanelCutout1_dia  = 466.*mm;
121  G4double RadPanelCutout1_len  = 200.*mm;
122  G4double RadPanelCutout2_hei  = 500.*mm;
123  G4double RadPanelCutout2_wid  = 400.*mm;
124  G4double RadPanelCutout2_len  = 400.*mm;
125  G4double RadPanelCutout2_zoff = 330.*mm;
126
127
128
129
130  // Structure Tubes **********************************************************
131
132  // Struct Tube: 2.16 kg
133  // (ID adjusted to make up for mass of missing top part of tube)
134  G4Tubs* StructTube_sol = new G4Tubs("StructTube_sol",
135    0.5*StructTube_dia-StructTube_thi, 0.5*StructTube_dia,
136    0.5*StructTube_hei, 0., 360.*deg);
137  G4LogicalVolume* StructTube_log = new G4LogicalVolume
138    (StructTube_sol, Al6061, "StructTube_log");
139  G4VPhysicalVolume* StructTube_phys;
140  StructTube_phys = new G4PVPlacement(0, spacecraft_pos +
141    G4ThreeVector(+0.5*StructTube_sep/std::cos(30.*deg),0,0),
142    "StructTube_phys", StructTube_log, wld_phys, false, 0);
143  StructTube_phys = new G4PVPlacement(0, spacecraft_pos +
144    G4ThreeVector(-0.5*StructTube_sep*std::tan(30.*deg),+0.5*StructTube_sep,0),
145    "StructTube_phys", StructTube_log, wld_phys, false, 1);
146  StructTube_phys = new G4PVPlacement(0, spacecraft_pos +
147    G4ThreeVector(-0.5*StructTube_sep*std::tan(30.*deg),-0.5*StructTube_sep,0),
148    "StructTube_phys", StructTube_log, wld_phys, false, 2);
149  StructTube_log->SetVisAttributes(sol_white_vat);
150
151  // Struct Tube Spacer: 1.0 kg
152  // (ID adjusted to include mass of missing TopCap)
153  G4Tubs* StructTubeSpacer_sol = new G4Tubs("StructTubeSpacer_sol",
154    0.5*StructTubeSpacer_dia-StructTubeSpacer_thi,
155    0.5*StructTubeSpacer_dia, 0.5*StructTubeSpacer_hei, 0., 360.*deg);
156  G4LogicalVolume* StructTubeSpacer_log = new G4LogicalVolume
157    (StructTubeSpacer_sol, Al6061, "StructTubeSpacer_log");
158  G4VPhysicalVolume* StructTubeSpacer_phys;
159  StructTubeSpacer_phys = new G4PVPlacement
160    (0, spacecraft_pos + G4ThreeVector(+0.5*StructTube_sep/std::cos(30.*deg),0,
161    +0.5*StructTube_hei+UpperDeckCore_thi+2*UpperDeckSkin_thi
162    +0.5*StructTubeSpacer_hei+UpperDeck_zoff),
163    "StructTubeSpacer_phys", StructTubeSpacer_log, wld_phys, false, 0);
164  StructTubeSpacer_phys = new G4PVPlacement
165    (0, spacecraft_pos + G4ThreeVector(-0.5*StructTube_sep*std::tan(30.*deg),
166    +0.5*StructTube_sep, +0.5*StructTube_hei+UpperDeckCore_thi
167    +2*UpperDeckSkin_thi+0.5*StructTubeSpacer_hei+UpperDeck_zoff),
168    "StructTubeSpacer_phys", StructTubeSpacer_log, wld_phys, false, 1);
169  StructTubeSpacer_phys = new G4PVPlacement
170    (0, spacecraft_pos + G4ThreeVector(-0.5*StructTube_sep*std::tan(30.*deg),
171    -0.5*StructTube_sep,+0.5*StructTube_hei+UpperDeckCore_thi
172    +2*UpperDeckSkin_thi+0.5*StructTubeSpacer_hei+UpperDeck_zoff),
173    "StructTubeSpacer_phys", StructTubeSpacer_log, wld_phys, false, 2);
174  StructTubeSpacer_log->SetVisAttributes(sol_grey_vat);
175
176
177  // Diagonal Panels **********************************************************
178
179  // Diagonal Panel: core (x3)
180  // Honeycomb core is represented by Al6061 at 0.05 g/cm3
181  G4Box* DiagPanelCore_box = new G4Box("DiagPanelCore_box",
182    0.5*DiagPanelCore_len-2.*cm, 0.5*DiagPanelCore_thi, 0.5*DiagPanelCore_hei);
183  G4Tubs* DiagPanelCutout_tub = new G4Tubs("DiagPanelCutout_tub", 0.,
184    0.5*DiagPanelCutout_dia, 2*DiagPanelCore_thi, 0., 360.*deg);
185  G4UnionSolid* DiagPanelCutout_sol = new G4UnionSolid("DiagPanelCutout_sol",
186    DiagPanelCutout_tub, DiagPanelCutout_tub,
187    G4Transform3D(G4RotationMatrix(), G4ThreeVector(DiagPanelCutout_off,0,0)));
188  G4RotationMatrix SMS_rot1; SMS_rot1.rotateX(90.*deg);
189  G4SubtractionSolid* DiagPanelCore_sol =
190    new G4SubtractionSolid("DiagPanelCore_sol", DiagPanelCore_box,
191    DiagPanelCutout_sol, G4Transform3D(SMS_rot1,
192    G4ThreeVector(-0.5*DiagPanelCutout_off,0,0)));
193  G4LogicalVolume* DiagPanelCore_log = new G4LogicalVolume
194    (DiagPanelCore_sol, AlHoneycomb, "DiagPanelCore_log");
195  G4VPhysicalVolume* DiagPanelCore_phys;
196  G4RotationMatrix SMS_rot2; SMS_rot2.rotateZ(90.*deg);
197  DiagPanelCore_phys = new G4PVPlacement(G4Transform3D(SMS_rot2,
198    spacecraft_pos + G4ThreeVector(-0.5*DiagPanelCore_len*std::tan(30.*deg),0,0)),
199    "DiagPanelCore_phys", DiagPanelCore_log, wld_phys, false, 0);
200  G4RotationMatrix SMS_rot3; SMS_rot3.rotateZ(-30.*deg);
201  DiagPanelCore_phys = new G4PVPlacement(G4Transform3D(SMS_rot3,
202    spacecraft_pos + G4ThreeVector
203    (.5*DiagPanelCore_len*std::sin(30.*deg)*std::tan(30.*deg),
204    0.5*DiagPanelCore_len*std::sin(30.*deg),0)), "DiagPanelCore_phys",
205    DiagPanelCore_log, wld_phys, false, 1);
206  G4RotationMatrix SMS_rot4; SMS_rot4.rotateZ(+30.*deg);
207  DiagPanelCore_phys = new G4PVPlacement(G4Transform3D(SMS_rot4,
208    spacecraft_pos + G4ThreeVector
209    (0.5*DiagPanelCore_len*std::sin(30*deg)*std::tan(30*deg),
210    -0.5*DiagPanelCore_len*std::sin(30.*deg),0)),
211    "DiagPanelCore_phys", DiagPanelCore_log, wld_phys, false, 2);
212  DiagPanelCore_log->SetVisAttributes(sol_white_vat);
213
214  // Diagonal Panel: skins (x6)
215  G4Box* DiagPanelSkin_box = new G4Box("DiagPanelSkin_box",
216    0.5*DiagPanelSkin_len-2.*cm, 0.5*DiagPanelSkin_thi, 0.5*DiagPanelSkin_hei);
217  G4SubtractionSolid* DiagPanelSkin_sol =
218    new G4SubtractionSolid("DiagPanelSkin_sol", DiagPanelSkin_box,
219    DiagPanelCutout_sol, G4Transform3D(SMS_rot1,
220    G4ThreeVector(-0.5*DiagPanelCutout_off,0,0)));
221  G4LogicalVolume* DiagPanelSkin_log = new G4LogicalVolume
222    (DiagPanelSkin_sol, CFRP, "DiagPanelSkin_log");
223  G4VPhysicalVolume* DiagPanelSkin_phys;
224  DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot2,
225    spacecraft_pos + G4ThreeVector(-0.5*DiagPanelCore_len*std::tan(30.*deg) +
226    0.5*(DiagPanelCore_thi+DiagPanelSkin_thi),0,0)),
227    "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 0);
228  DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot2,
229    spacecraft_pos + G4ThreeVector(-0.5*DiagPanelCore_len*std::tan(30.*deg) -
230    0.5*(DiagPanelCore_thi+DiagPanelSkin_thi),0,0)),
231    "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 1);
232  DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot3,
233    spacecraft_pos + G4ThreeVector
234    (0.5*DiagPanelCore_len*std::sin(30*deg)*std::tan(30*deg),
235    0.5*DiagPanelCore_len*std::sin(30.*deg),0) +
236    G4ThreeVector(0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::sin(30.*deg),
237    0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::cos(30.*deg),0)),
238    "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 2);
239  DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot3,
240    spacecraft_pos + G4ThreeVector
241    (0.5*DiagPanelCore_len*std::sin(30*deg)*std::tan(30*deg),
242    0.5*DiagPanelCore_len*std::sin(30.*deg),0) +
243    G4ThreeVector(-0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::sin(30.*deg),
244    -0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::cos(30.*deg),0)),
245    "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 3);
246  DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot4,
247    spacecraft_pos + G4ThreeVector
248    (0.5*DiagPanelCore_len*std::sin(30*deg)*std::tan(30*deg),
249    -0.5*DiagPanelCore_len*std::sin(30.*deg),0) +
250    G4ThreeVector(0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::sin(30.*deg),
251    -0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::cos(30.*deg),0)),
252    "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 3);
253  DiagPanelSkin_phys = new G4PVPlacement(G4Transform3D(SMS_rot4,
254    spacecraft_pos + G4ThreeVector
255    (0.5*DiagPanelCore_len*std::sin(30*deg)*std::tan(30*deg),
256    -0.5*DiagPanelCore_len*std::sin(30.*deg),0) +
257    G4ThreeVector(-0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::sin(30.*deg),
258    +0.5*(DiagPanelCore_thi+DiagPanelSkin_thi)*std::cos(30.*deg),0)),
259    "DiagPanelSkin_phys", DiagPanelSkin_log, wld_phys, false, 4);
260  DiagPanelSkin_log->SetVisAttributes(sol_white_vat);
261
262
263
264  // Radial Panels ************************************************************
265
266  // Radial Panels: skins (x3)
267  G4Trap* RadialPanelSkin_sol = new G4Trap("RadialPanelSkin_sol",
268    2*RadialPanelSkin_thi+RadialPanelCore_thi, RadialPanelSkin_hei,
269    RadialPanelSkin_top, RadialPanelSkin_bot);
270  G4LogicalVolume* RadialPanelSkin_log = new G4LogicalVolume
271    (RadialPanelSkin_sol, CFRP, "RadialPanelSkin_log");
272  G4VPhysicalVolume* RadialPanelSkin_phys;
273  G4RotationMatrix SMS_rot5; SMS_rot5.rotateX(-90.*deg);
274  RadialPanelSkin_phys = new G4PVPlacement
275    (G4Transform3D(SMS_rot5, spacecraft_pos +
276    G4ThreeVector(+RadialPanel_roff/std::cos(30.*deg),0,0)),
277    "RadialPanelSkin_phys", RadialPanelSkin_log, wld_phys, false, 0);
278  G4RotationMatrix SMS_rot6;
279  SMS_rot6.rotateX(-90.*deg); SMS_rot6.rotateZ(120.*deg);
280  RadialPanelSkin_phys = new G4PVPlacement
281    (G4Transform3D(SMS_rot6, spacecraft_pos +
282    G4ThreeVector(-RadialPanel_roff*std::tan(30.*deg),+RadialPanel_roff,0)),
283    "RadialPanelSkin_phys", RadialPanelSkin_log, wld_phys, false, 1);
284  G4RotationMatrix SMS_rot7;
285  SMS_rot7.rotateX(-90.*deg); SMS_rot7.rotateZ(-120.*deg);
286  RadialPanelSkin_phys = new G4PVPlacement
287    (G4Transform3D(SMS_rot7, spacecraft_pos +
288    G4ThreeVector(-RadialPanel_roff*std::tan(30.*deg),-RadialPanel_roff,0)),
289    "RadialPanelSkin_phys", RadialPanelSkin_log, wld_phys, false, 2);
290  RadialPanelSkin_log->SetVisAttributes(sol_white_vat);
291  // Radial Panels: core
292  // Honeycomb core is represented by Al6061 at 0.05 g/cm3
293  G4Trap* RadialPanelCore_sol = new G4Trap("RadialPanelCore_sol",
294    RadialPanelCore_thi, RadialPanelSkin_hei,
295    RadialPanelSkin_top, RadialPanelSkin_bot);
296  G4LogicalVolume* RadialPanelCore_log = new G4LogicalVolume
297    (RadialPanelCore_sol, AlHoneycomb, "RadialPanelCore_log");
298  G4VPhysicalVolume* RadialPanelCore_phys;
299  RadialPanelCore_phys = new G4PVPlacement
300    (0, G4ThreeVector(0,0,0), "RadialPanelCore_phys",
301    RadialPanelCore_log, RadialPanelSkin_phys, false, 0);
302  RadialPanelCore_log->SetVisAttributes(sol_white_vat);
303
304
305
306  // Upper and Lower Decks ****************************************************
307
308
309  // Upper Deck - skin: 2*4.56 kg
310  G4Tubs* UpperDeckSkin_sol = new G4Tubs
311    ("UpperDeckSkin_sol", 0., 0.5*UpperDeck_dia,
312    0.5*UpperDeckCore_thi + UpperDeckSkin_thi, 0., 360.*deg);
313  G4LogicalVolume* UpperDeckSkin_log = new G4LogicalVolume
314    (UpperDeckSkin_sol, CFRP, "UpperDeckSkin_log");
315  G4VPhysicalVolume* UpperDeckSkin_phys;
316  UpperDeckSkin_phys = new G4PVPlacement
317    (0, spacecraft_pos + G4ThreeVector(0,0,+0.5*StructTube_hei
318    +0.5*UpperDeckCore_thi+UpperDeckSkin_thi+UpperDeck_zoff),
319    "UpperDeckSkin_phys", UpperDeckSkin_log, wld_phys, false, 0);
320  UpperDeckSkin_log->SetVisAttributes(white_vat);
321  // upper deck - core: 8.3 kg
322  // Honeycomb core is represented by Al6061 at 0.05 g/cm3
323  G4Tubs* UpperDeckCore_sol = new G4Tubs
324    ("UpperDeckCore_sol",0.,.5*UpperDeck_dia,.5*UpperDeckCore_thi,0.,360.*deg);
325  G4LogicalVolume* UpperDeckCore_log = new G4LogicalVolume
326    (UpperDeckCore_sol, AlHoneycomb, "UpperDeckCore_log");
327  G4VPhysicalVolume* UpperDeckCore_phys;
328  UpperDeckCore_phys = new G4PVPlacement(0,G4ThreeVector(),
329    "UpperDeckCore_phys", UpperDeckCore_log, UpperDeckSkin_phys, false, 0);
330  UpperDeckCore_log->SetVisAttributes(white_vat);
331
332  // Lower Deck - skin: 2*2.95 kg
333  G4Tubs* LowerDeckSkin_sol = new G4Tubs
334    ("LowerDeckSkin_sol", 0., 0.5*LowerDeck_dia,
335    0.5*LowerDeckCore_thi + LowerDeckSkin_thi, 0., 360.*deg);
336  G4LogicalVolume* LowerDeckSkin_log = new G4LogicalVolume
337    (LowerDeckSkin_sol, CFRP, "LowerDeckSkin_log");
338  G4VPhysicalVolume* LowerDeckSkin_phys;
339  LowerDeckSkin_phys = new G4PVPlacement
340    (0, spacecraft_pos + G4ThreeVector(0,0,-0.5*StructTube_hei
341    -0.5*LowerDeckCore_thi-LowerDeckSkin_thi-LowerDeck_zoff),
342    "LowerDeckSkin_phys", LowerDeckSkin_log, wld_phys, false, 0);
343  LowerDeckSkin_log->SetVisAttributes(sol_dgrey_vat);
344  // lower deck - core: 5.38 kg
345  // Honeycomb core is represented by Al6061 at 0.05 g/cm3
346  G4Tubs* LowerDeckCore_sol = new G4Tubs
347    ("LowerDeckCore_sol",0.,.5*LowerDeck_dia,.5*LowerDeckCore_thi,0.,360.*deg);
348  G4LogicalVolume* LowerDeckCore_log = new G4LogicalVolume
349    (LowerDeckCore_sol, AlHoneycomb, "LowerDeckCore_log");
350  G4VPhysicalVolume* LowerDeckCore_phys;
351  LowerDeckCore_phys = new G4PVPlacement(0,G4ThreeVector(),
352    "LowerDeckCore_phys", LowerDeckCore_log, LowerDeckSkin_phys, false, 0);
353  LowerDeckCore_log->SetVisAttributes(sol_dgrey_vat);
354
355  // Lower Deck MLI Blanket: 5.3 kg
356  G4Tubs* MLIBlanket_sol = new G4Tubs("MLIBlanket_sol",
357    0., 0.5*MLIBlanket_dia, 0.5*MLIBlanket_thi, 0., 360.*deg);
358  G4LogicalVolume* MLIBlanket_log = new G4LogicalVolume
359    (MLIBlanket_sol, MLImat, "MLIBlanket_log");
360  G4VPhysicalVolume* MLIBlanket_phys;
361  MLIBlanket_phys = new G4PVPlacement
362    (0, spacecraft_pos + G4ThreeVector(0,0,-0.5*StructTube_hei
363    -LowerDeckCore_thi-2*LowerDeckSkin_thi-LowerDeck_zoff-.5*MLIBlanket_thi),
364    "MLIBlanket_phys", MLIBlanket_log, wld_phys, false, 0);
365  MLIBlanket_log->SetVisAttributes(sol_orange_vat);
366
367
368  // Thermal Shield ***********************************************************
369
370  // Thermal Shield (Facesheet+HoneycombCore+Facesheet+Foam+Facesheet)
371  // Thermal Shield - facesheet: 3*5.47 kg
372  G4Tubs* ThermalShieldFacesheet_sol = new G4Tubs("ThermalShieldFacesheet_sol",
373    0., 0.5*UpperDeck_dia, 0.5*ThermalShield_thi, 0., 360.*deg);
374  G4LogicalVolume* ThermalShieldFacesheet_log = new G4LogicalVolume
375    (ThermalShieldFacesheet_sol, CFRP, "ThermalShieldFacesheet_log");
376  G4VPhysicalVolume* ThermalShieldFacesheet_phys;
377  ThermalShieldFacesheet_phys = new G4PVPlacement(0,
378    spacecraft_pos + G4ThreeVector
379    (0,0,+.5*StructTube_hei+UpperDeckCore_thi+2*UpperDeckSkin_thi+
380    UpperDeck_zoff+0.5*ThermalShield_thi+ThermalShield_zoff),
381    "ThermalShieldFacesheet_phys",ThermalShieldFacesheet_log,wld_phys,false,0);
382  ThermalShieldFacesheet_log->SetVisAttributes(red_vat);
383  // Thermal Shield: core
384  // Honeycomb core is represented by Al6061 at 0.05 g/cm3
385  G4Tubs* ThermalShieldCore_sol = new G4Tubs("ThermalShieldCore_sol",
386    0., 0.5*UpperDeck_dia, 0.5*ThermalShieldCore_thi, 0., 360.*deg);
387  G4LogicalVolume* ThermalShieldCore_log = new G4LogicalVolume
388    (ThermalShieldCore_sol, AlHoneycomb, "ThermalShieldCore_log");
389  G4VPhysicalVolume* ThermalShieldCore_phys;
390  ThermalShieldCore_phys = new G4PVPlacement(0,
391    G4ThreeVector(0,0,.5*ThermalShieldFacesheet_thi+.5*ThermalShieldCore_thi),
392    "ThermalShieldCore_phys", ThermalShieldCore_log,
393    ThermalShieldFacesheet_phys, false, 0);
394  ThermalShieldCore_log->SetVisAttributes(red_vat);
395  // Thermal Shield: foam
396  G4Tubs* ThermalShieldFoam_sol = new G4Tubs("ThermalShieldFoam_sol",
397    0., 0.5*UpperDeck_dia, 0.5*ThermalShieldFoam_thi, 0., 360.*deg);
398  G4LogicalVolume* ThermalShieldFoam_log = new G4LogicalVolume
399    (ThermalShieldFoam_sol, foam, "ThermalShieldFoam_log");
400  G4VPhysicalVolume* ThermalShieldFoam_phys;
401  ThermalShieldFoam_phys = new G4PVPlacement(0,
402    G4ThreeVector(0,0,-.5*ThermalShieldFacesheet_thi-.5*ThermalShieldFoam_thi),
403    "ThermalShieldFoam_phys", ThermalShieldFoam_log,
404    ThermalShieldFacesheet_phys, false, 0);
405  ThermalShieldFoam_log->SetVisAttributes(red_vat);
406
407
408
409  // Solar Array and OSR Assy *************************************************
410
411  // Solar Array: 13.5 kg
412  // Thickness for correct mass
413  G4Tubs* SolarArray_sol = new G4Tubs("SolarArray_sol",
414    SolarArray_i, SolarArray_o, 0.5*SolarArray_thi, 0., 360.*deg);
415  G4LogicalVolume* SolarArray_log = new G4LogicalVolume
416    (SolarArray_sol, Scell, "SolarArray_log");
417  G4VPhysicalVolume* SolarArray_phys;
418  SolarArray_phys = new G4PVPlacement(0, spacecraft_pos +
419    G4ThreeVector(0,0,+.5*StructTube_hei+UpperDeckCore_thi+2*UpperDeckSkin_thi+
420    UpperDeck_zoff+ThermalShield_thi+ThermalShield_zoff+0.5*SolarArray_thi),
421    "SolarArray_phys",SolarArray_log, wld_phys, false, 0);
422  SolarArray_log->SetVisAttributes(lblue_vat);
423  // Optical Solar Reflector: NO MASS
424  G4Tubs* OSRAssy_sol = new G4Tubs("OSRAssy_sol",
425    0., OSRAssy_o, 0.5*OSRAssy_thi, 0., 360.*deg);
426  G4LogicalVolume* OSRAssy_log = new G4LogicalVolume
427    (OSRAssy_sol, vacuum, "OSRAssy_log");
428  G4VPhysicalVolume* OSRAssy_phys;
429  OSRAssy_phys = new G4PVPlacement(0, spacecraft_pos +
430    G4ThreeVector(0,0,+.5*StructTube_hei+UpperDeckCore_thi+2*UpperDeckSkin_thi+
431    UpperDeck_zoff+ThermalShield_thi+ThermalShield_zoff+0.5*OSRAssy_thi),
432    "OSRAssy_phys",OSRAssy_log, wld_phys, false, 0);
433  OSRAssy_log->SetVisAttributes(lblue_vat);
434
435
436
437  // Radiator Panels **********************************************************
438
439  // Radiator Panels: 3*3.87 kg
440  // Thickess adjusted for correct mass of FS+Core+FS
441  G4Cons* RadPanel_con = new G4Cons("RadPanel_con",
442     0.5*RadPanel_bot-RadPanel_thi, 0.5*RadPanel_bot,
443     0.5*RadPanel_top-RadPanel_thi, 0.5*RadPanel_top,
444     0.5*RadPanel_len, 0.*deg, 360.*deg);
445  // window cutouts
446  G4Tubs* RadPanelCutout_tub = new G4Tubs("RadPanelCutout_tub", 0.,
447    0.5*RadPanelCutout1_dia, 2*RadPanelCutout1_len, 0., 360.*deg);
448  G4RotationMatrix SMS_rot8; SMS_rot8.rotateY(90.*deg);
449  G4SubtractionSolid* RadPanel1_sol =
450    new G4SubtractionSolid("RadPanel1_sol", RadPanel_con,
451    RadPanelCutout_tub, G4Transform3D(SMS_rot8,
452    G4ThreeVector(-RadPanelCutout_roff,0,0)));
453  G4RotationMatrix SMS_rot9;
454  SMS_rot9.rotateY(90.*deg); SMS_rot9.rotateZ(49.*deg);
455  G4SubtractionSolid* RadPanel2_sol =
456    new G4SubtractionSolid("RadPanel2_sol", RadPanel1_sol,
457    RadPanelCutout_tub, G4Transform3D(SMS_rot9, G4ThreeVector
458    (RadPanelCutout_roff*std::sin(49.*deg),RadPanelCutout_roff*std::cos(49.*deg),0)));
459  G4RotationMatrix SMS_rot10;
460  SMS_rot10.rotateY(90.*deg); SMS_rot10.rotateZ(131.*deg);
461  G4SubtractionSolid* RadPanel3_sol =
462    new G4SubtractionSolid("RadPanel3_sol", RadPanel2_sol,
463    RadPanelCutout_tub, G4Transform3D(SMS_rot10, G4ThreeVector
464    (RadPanelCutout_roff*std::sin(131.*deg),RadPanelCutout_roff*std::cos(131.*deg),0)));
465  // box cutouts
466  G4Box* RadPanelCutout_box = new G4Box("RadPanelCutout_box",
467    0.5*RadPanelCutout2_len, 0.5*RadPanelCutout2_wid, 0.5*RadPanelCutout2_hei);
468  G4SubtractionSolid* RadPanel4_sol =
469    new G4SubtractionSolid("RadPanel4_sol", RadPanel3_sol,
470    RadPanelCutout_box, G4Transform3D(G4RotationMatrix(), G4ThreeVector
471    (RadPanelCutout_roff,0,RadPanelCutout2_zoff)));
472  G4RotationMatrix SMS_rot11; SMS_rot11.rotateZ(+120.*deg);
473  G4SubtractionSolid* RadPanel5_sol =
474    new G4SubtractionSolid("RadPanel5_sol", RadPanel4_sol,
475    RadPanelCutout_box, G4Transform3D(SMS_rot11,
476    G4ThreeVector(RadPanelCutout_roff*std::cos(120.*deg),
477    RadPanelCutout_roff*std::sin(120.*deg), RadPanelCutout2_zoff)));
478  G4RotationMatrix SMS_rot12; SMS_rot12.rotateZ(-120.*deg);
479  G4SubtractionSolid* RadPanel6_sol =
480    new G4SubtractionSolid("RadPanel6_sol", RadPanel5_sol,
481    RadPanelCutout_box, G4Transform3D(SMS_rot12,
482    G4ThreeVector(RadPanelCutout_roff*std::cos(-120.*deg),
483    RadPanelCutout_roff*std::sin(-120.*deg), RadPanelCutout2_zoff)));
484  G4LogicalVolume* RadPanel_log =
485    new G4LogicalVolume(RadPanel6_sol, Al6061, "RadPanel_log");
486  G4VPhysicalVolume* RadPanel_phys;
487  RadPanel_phys =
488    new G4PVPlacement(0, spacecraft_pos + G4ThreeVector(0,0,0),
489    "RadPanel_phys", RadPanel_log, wld_phys, false, 0);
490  RadPanel_log->SetVisAttributes(sol_white_vat);
491
492//*****************************************************************************
493
494
495  // invisible spacecraft top
496  if(0) {
497     UpperDeckSkin_log->SetVisAttributes(G4VisAttributes::Invisible);
498     UpperDeckCore_log->SetVisAttributes(G4VisAttributes::Invisible);
499     ThermalShieldFacesheet_log->SetVisAttributes(G4VisAttributes::Invisible);
500     ThermalShieldCore_log->SetVisAttributes(G4VisAttributes::Invisible);
501     ThermalShieldFoam_log->SetVisAttributes(G4VisAttributes::Invisible);
502     //  SolarArray_log->SetVisAttributes(G4VisAttributes::Invisible);
503     OSRAssy_log->SetVisAttributes(G4VisAttributes::Invisible);
504  }
505 
506 // solid spacecraft top
507  if(0) {
508     UpperDeckSkin_log->SetVisAttributes(sol_white_vat);
509     UpperDeckCore_log->SetVisAttributes(sol_white_vat);
510     ThermalShieldFacesheet_log->SetVisAttributes(sol_red_vat);
511     ThermalShieldCore_log->SetVisAttributes(sol_red_vat);
512     ThermalShieldFoam_log->SetVisAttributes(sol_red_vat);
513     SolarArray_log->SetVisAttributes(sol_lblue_vat);
514     OSRAssy_log->SetVisAttributes(sol_lgreen_vat);
515  }
516
517//*****************************************************************************
Note: See TracBrowser for help on using the repository browser.