Context Navigation

BuildShapes.cc @ 1347

Last change on this file since 1347 was 1347, checked in by garnier, 13 years ago
geant4 tag 9.4
File size: 13.1 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: BuildShapes.cc,v 1.6 2006/06/29 21:34:20 gunter Exp $
28	// GEANT4 tag $Name: geant4-09-04-ref-00 $
29	//
30	//
31
32	#include "BuildShapes.hh"
33
34	#include "G4PVPlacement.hh"
35	#include "G4LogicalVolume.hh"
36
37	#include "G4Material.hh"
38	#include "G4VisAttributes.hh"
39
40	#include "G4Box.hh"
41	#include "G4Cons.hh"
42	#include "G4Tubs.hh"
43	#include "G4Trd.hh"
44	#include "G4Trap.hh"
45	#include "G4Sphere.hh"
46	#include "G4Para.hh"
47
48	// #define USE_BREP_SOLID
49
50	#if defined USE_BREP_SOLID
51
52	#include"G4BREPSolidPCone.h"
53	#include"G4BREPSolidPolyhedra.h"
54
55	#endif
56
57	G4VPhysicalVolume* TheWorld ()
58	{
59	G4double a, iz, density;
60	G4String name, symbol;
61	G4int nel;
62	a = 14.01*g/mole;
63	G4Element* elN = new G4Element(name="Nitrogen", symbol="N", iz=7., a);
64	a = 16.00*g/mole;
65	G4Element* elO = new G4Element(name="Oxigen", symbol="O", iz=8., a);
66	density = 1.29e-03*g/cm3;
67	G4Material* Air = new G4Material(name="Air", density, nel=2);
68	Air->AddElement(elN, .7);
69	Air->AddElement(elO, .3);
70
71	G4Box *myWorldBox= new G4Box("WBox",10.0, 10.0, 10.0); // Big enough!
72	G4LogicalVolume *myWorldLog=new G4LogicalVolume(myWorldBox,Air,
73	"WLog", 0, 0, 0);
74	myWorldLog -> SetVisAttributes (&G4VisAttributes::Invisible);
75	G4PVPlacement *myWorldPhys=new G4PVPlacement(0,G4ThreeVector(),
76	"WPhys",
77	myWorldLog,
78	0,false,0);
79	return myWorldPhys;
80	}
81
82	G4VPhysicalVolume* BuildBox()
83	{
84	G4VPhysicalVolume* theWorld = TheWorld ();
85
86	//----- material
87	G4double z = 13.;
88	G4double a = 26.98*g/mole;
89	G4double density = 2.7*g/cm3;
90	G4Material myMaterial("Aluminium", z, a, density);
91
92	//----- Vis Attributes
93	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(1,0,0));
94
95	//----- Solid
96	G4Box *box = new G4Box ("Box", 1 ,4, 9 );
97
98	//----- Logical Volume
99	G4LogicalVolume *boxLog=new G4LogicalVolume( box, &myMaterial,
100	"BoxLog",0,0,0);
101	boxLog -> SetVisAttributes ( pVA );
102
103	//----- physical volume
104	new G4PVPlacement( 0,G4ThreeVector(),
105	"BoxPhys",
106	boxLog,
107	theWorld,
108	false,0);
109
110	return theWorld;
111
112	} // BuildBox
113
114
115	G4VPhysicalVolume* BuildCylinder()
116	{
117	G4VPhysicalVolume* theWorld = TheWorld ();
118
119	//----- material
120	G4double z = 13.;
121	G4double a = 26.98*g/mole;
122	G4double density = 2.7*g/cm3;
123	G4Material myMaterial("Aluminium", z, a, density);
124
125	//----- Vis Attributes
126	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(0,1,0));
127
128	//----- Solid
129	G4Tubs *tubs = new G4Tubs ("Tubs", 0. , 1.5, 2.0 , 0., twopi);
130
131	//----- Logical Volume
132	G4LogicalVolume *tubsLog=new G4LogicalVolume( tubs, &myMaterial,
133	"tubsLog",0,0,0);
134	tubsLog -> SetVisAttributes ( pVA );
135
136	//----- physical volume
137	new G4PVPlacement( 0,G4ThreeVector(),
138	"tubsPhys",
139	tubsLog,
140	theWorld,
141	false,0);
142
143	return theWorld;
144
145	}// BuildCyllinder
146
147
148
149	G4VPhysicalVolume* BuildTubs()
150	{
151	G4VPhysicalVolume* theWorld = TheWorld ();
152
153	//----- material
154	G4double z = 13.;
155	G4double a = 26.98*g/mole;
156	G4double density = 2.7*g/cm3;
157	G4Material myMaterial("Aluminium", z, a, density);
158
159	//----- Vis Attributes
160	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(0,1,0));
161
162	//----- Solid
163	G4Tubs *tubs = new G4Tubs ("Tubs", 1.0 , 1.5, 2.0 , 1.5, 4.5 );
164
165	//----- Logical Volume
166	G4LogicalVolume *tubsLog=new G4LogicalVolume( tubs, &myMaterial,
167	"tubsLog",0,0,0);
168	tubsLog -> SetVisAttributes ( pVA );
169
170	//----- physical volume
171	new G4PVPlacement( 0,G4ThreeVector(),
172	"tubsPhys",
173	tubsLog,
174	theWorld,
175	false,0);
176
177	return theWorld;
178
179	}// BuildTubs
180
181
182
183	G4VPhysicalVolume* BuildCons()
184	{
185	G4VPhysicalVolume* theWorld = TheWorld ();
186
187	//----- material
188	G4double z = 13.;
189	G4double a = 26.98*g/mole;
190	G4double density = 2.7*g/cm3;
191	G4Material myMaterial("Aluminium", z, a, density);
192
193	//----- Vis Attributes
194	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(0,0,1));
195
196	//----- Solid
197	G4Cons *cons = new G4Cons ("Cons", 1.0 , 1.5, 0.5, 1.0 ,2.0, 1.5, 4.5 );
198
199	//----- Logical Volume
200	G4LogicalVolume *consLog=new G4LogicalVolume( cons, &myMaterial,
201	"consLog",0,0,0);
202	consLog -> SetVisAttributes ( pVA );
203
204	//----- physical volume
205	new G4PVPlacement( 0,G4ThreeVector(),
206	"consPhys",
207	consLog,
208	theWorld,
209	false,0);
210
211	return theWorld;
212
213	}// BuildCons
214
215
216	G4VPhysicalVolume* BuildTrd()
217	{
218	G4VPhysicalVolume* theWorld = TheWorld ();
219
220	//----- material
221	G4double z = 13.;
222	G4double a = 26.98*g/mole;
223	G4double density = 2.7*g/cm3;
224	G4Material myMaterial("Aluminium", z, a, density);
225
226	//----- Vis Attributes
227	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(0,1,1));
228
229	//----- Solid
230	G4Trd *trd = new G4Trd ("Trd", 1.0 , 0.5, 1.0, 0.5, 4.0 );
231
232	//----- Logical Volume
233	G4LogicalVolume *trdLog=new G4LogicalVolume( trd, &myMaterial,
234	"trdLog",0,0,0);
235	trdLog -> SetVisAttributes ( pVA );
236
237	//----- physical volume
238	new G4PVPlacement( 0,G4ThreeVector(),
239	"trdPhys",
240	trdLog,
241	theWorld,
242	false,0);
243
244	return theWorld;
245
246	}// BuildTrd
247
248
249
250	G4VPhysicalVolume* BuildTrap()
251	{
252	G4VPhysicalVolume* theWorld = TheWorld ();
253
254	//----- material
255	G4double z = 13.;
256	G4double a = 26.98*g/mole;
257	G4double density = 2.7*g/cm3;
258	G4Material myMaterial("Aluminium", z, a, density);
259
260	//----- Vis Attributes
261	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(1,0,1));
262
263	//----- Solid
264	G4Trap *trap = new G4Trap ("Trap", 2.0, 0.26, 0.0, 1.0, \
265	0.7, 1.2, 0.0, 0.5, 0.35, 0.6, 0.0 );
266
267	//----- Logical Volume
268	G4LogicalVolume *trapLog=new G4LogicalVolume( trap, &myMaterial,
269	"trapLog",0,0,0);
270	trapLog -> SetVisAttributes ( pVA );
271
272	//----- physical volume
273	new G4PVPlacement( 0,G4ThreeVector(),
274	"trapPhys",
275	trapLog,
276	theWorld,
277	false,0);
278
279	return theWorld;
280
281	}// BuildTrap
282
283
284	G4VPhysicalVolume* BuildSphereFull()
285	{
286	G4VPhysicalVolume* theWorld = TheWorld ();
287
288	//----- material
289	G4double z = 13.;
290	G4double a = 26.98*g/mole;
291	G4double density = 2.7*g/cm3;
292	G4Material myMaterial("Aluminium", z, a, density);
293
294	//----- Vis Attributes
295	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(1,1,0));
296
297	//----- Solid
298	G4Sphere *Sphere = new G4Sphere ("Sphere", 1.6, 1.8, 0.0, (twopi), 0.0, (twopi) );
299
300	//----- Logical Volume
301	G4LogicalVolume *SphereLog=new G4LogicalVolume( Sphere, &myMaterial,
302	"SphereLog",0,0,0);
303	SphereLog -> SetVisAttributes ( pVA );
304
305	//----- physical volume
306	new G4PVPlacement( 0,G4ThreeVector(),
307	"SpherePhys",
308	SphereLog,
309	theWorld,
310	false,0);
311
312	return theWorld;
313
314	} // BuildSphereFull
315
316
317	G4VPhysicalVolume* BuildSphereSeg()
318	{
319	G4VPhysicalVolume* theWorld = TheWorld ();
320
321	//----- material
322	G4double z = 13.;
323	G4double a = 26.98*g/mole;
324	G4double density = 2.7*g/cm3;
325	G4Material myMaterial("Aluminium", z, a, density);
326
327	//----- Vis Attributes
328	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(1,1,0));
329
330	//----- Solid
331	G4Sphere *Sphere = new G4Sphere ("Sphere", 1.6, 1.8, 0.0, 6.0, 0.5, 2.5 );
332
333	//----- Logical Volume
334	G4LogicalVolume *SphereLog=new G4LogicalVolume( Sphere, &myMaterial,
335	"SphereLog",0,0,0);
336	SphereLog -> SetVisAttributes ( pVA );
337
338	//----- physical volume
339	new G4PVPlacement( 0,G4ThreeVector(),
340	"SpherePhys",
341	SphereLog,
342	theWorld,
343	false,0);
344
345	return theWorld;
346
347	} // BuildSphereSeg
348
349
350
351	G4VPhysicalVolume* BuildPara()
352	{
353	G4VPhysicalVolume* theWorld = TheWorld ();
354
355	//----- material
356	G4double z = 13.;
357	G4double a = 26.98*g/mole;
358	G4double density = 2.7*g/cm3;
359	G4Material myMaterial("Aluminium", z, a, density);
360
361	//----- Vis Attributes
362	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(1,0,0));
363
364	//----- Solid
365	// G4Para *para = new G4Para ("Para", 1.0 ,1.0, 3.0, 0.0, 0.0, 0.0 );
366	G4Para *para = new G4Para ("Para", 1.0 ,1.0, 3.0, pi/6.0, pi/12.0, pi/12.0 );
367
368	//----- Logical Volume
369	G4LogicalVolume *paraLog=new G4LogicalVolume( para, &myMaterial,
370	"ParaLog",0,0,0);
371	paraLog -> SetVisAttributes ( pVA );
372
373	//----- physical volume
374	new G4PVPlacement( 0,G4ThreeVector(),
375	"ParaPhys",
376	paraLog,
377	theWorld,
378	false,0);
379
380	return theWorld;
381
382	} // BuildPara
383
384
385
386
387	G4VPhysicalVolume* BuildPCon()
388	{
389	#if defined USE_BREP_SOLID
390	G4VPhysicalVolume* theWorld = TheWorld ();
391
392	//----- material
393	G4double z = 13.;
394	G4double a = 26.98*g/mole;
395	G4double density = 2.7*g/cm3;
396	G4Material myMaterial("Aluminium", z, a, density);
397
398	//----- shape
399	G4double sphi = 1.8 ;
400	G4double dphi = 4.71 ;
401	G4int nz = 4 ;
402	G4double DzArray [] = { -3.0, 1.0, 1.0, 3.0 };
403	G4double RminArray [] = { 1.0, 0.5, 1.5, 1.0 };
404	G4double RmaxArray [] = { 1.5, 1.0, 2.0, 1.5 };
405
406	//----- Vis Attributes
407	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(0,1,0));
408
409	//----- Solid
410	G4BREPSolidPCone *pcon = new G4BREPSolidPCone ("Pcon", sphi , dphi, nz, DzArray[0], \
411	DzArray, RminArray, RmaxArray );
412
413	//----- Logical Volume
414	G4LogicalVolume *pconLog=new G4LogicalVolume( pcon, &myMaterial,
415	"PconLog",0,0,0);
416	pconLog -> SetVisAttributes ( pVA );
417
418	//----- physical volume
419	new G4PVPlacement( 0,G4ThreeVector(),
420	"PconPhys",
421	pconLog,
422	theWorld,
423	false,0);
424
425	return theWorld;
426
427	#else
428	return NULL ;
429	#endif
430
431	} // BuildPCon
432
433
434	G4VPhysicalVolume* BuildPGon()
435	{
436	#if defined USE_BREP_SOLID
437	G4VPhysicalVolume* theWorld = TheWorld ();
438
439	//----- material
440	G4double z = 13.;
441	G4double a = 26.98*g/mole;
442	G4double density = 2.7*g/cm3;
443	G4Material myMaterial("Aluminium", z, a, density);
444
445	//----- shape
446	G4double sphi = 1.8 ;
447	G4double dphi = 4.71 ;
448	G4int nside = 5 ;
449	G4int nz = 4 ;
450	G4double DzArray [] = { -3.0, 1.0, 1.0, 3.0 };
451	G4double RminArray [] = { 1.0, 0.5, 1.5, 1.0 };
452	G4double RmaxArray [] = { 1.5, 1.0, 2.0, 1.5 };
453
454	//----- Vis Attributes
455	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(1,0,0));
456
457	//----- Solid
458	G4BREPSolidPolyhedra *pgon = new G4BREPSolidPolyhedra ("Pcon", sphi , dphi, nside, nz, \
459	DzArray[0], \
460	DzArray, RminArray, RmaxArray );
461
462	//----- Logical Volume
463	G4LogicalVolume *pgonLog=new G4LogicalVolume( pgon, &myMaterial,
464	"PgonLog",0,0,0);
465	pgonLog -> SetVisAttributes ( pVA );
466
467	//----- physical volume
468	new G4PVPlacement( 0,G4ThreeVector(),
469	"PgonPhys",
470	pgonLog,
471	theWorld,
472	false,0);
473
474	return theWorld;
475	#else
476	return NULL ;
477	#endif
478	} // BuildPGon
479
480	G4VPhysicalVolume* BuildForcedWireframeBox()
481	{
482	G4VPhysicalVolume* theWorld = TheWorld ();
483
484	//----- material
485	G4double z = 13.;
486	G4double a = 26.98*g/mole;
487	G4double density = 2.7*g/cm3;
488	G4Material myMaterial("Aluminium", z, a, density);
489
490	//----- Vis Attributes
491	G4VisAttributes* pVA = new G4VisAttributes (G4Colour(1,0,0)); // for box
492	pVA->SetForceWireframe( true );
493	G4VisAttributes* pVA2 = new G4VisAttributes (G4Colour(0,1,1)); // for tubs
494
495	//----- Solid
496	G4Box *box = new G4Box ("Box", 2 ,4, 6 );
497	G4Tubs *tubs = new G4Tubs ("Tubs", 1.0 , 1.5, 2.0 , 1.5, 4.5 );
498
499	//----- Logical Volume
500	G4LogicalVolume *boxLog=new G4LogicalVolume( box, &myMaterial,
501	"BoxLog",0,0,0);
502	boxLog -> SetVisAttributes ( pVA );
503	G4LogicalVolume *tubsLog=new G4LogicalVolume( tubs, &myMaterial,
504	"tubsLog",0,0,0);
505	tubsLog -> SetVisAttributes ( pVA2 );
506
507	//----- physical volume
508	G4PVPlacement *boxPhys=new G4PVPlacement( 0,G4ThreeVector(),
509	"BoxPhys",
510	boxLog,
511	theWorld,
512	false,0);
513	new G4PVPlacement( 0,G4ThreeVector(),
514	"tubsPhys",
515	tubsLog,
516	boxPhys ,false,0);
517	return theWorld;
518
519	} // G4VPhysicalVolume* BuidlForcedWireFrameBox()

Note: See TracBrowser for help on using the repository browser.

Context Navigation

source: trunk/source/visualization/test/src/BuildShapes.cc @ 1347

Download in other formats: