source: trunk/source/visualization/test/src/MyDetectorConstruction.cc@ 858

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

geant4.8.2 beta
File size: 29.3 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: MyDetectorConstruction.cc,v 1.39 2006/11/14 16:03:55 allison Exp $
28// GEANT4 tag $Name: HEAD $
29//
30//
31#include "MyDetectorConstruction.hh"
32
33#include "MyDetectorMessenger.hh"
34//#include "MyCalorimeterSD.hh"
35//#include "MyTrackerSD.hh"
36//#include "MyCalorimeterHit.hh"
37//#include "MyTrackerHit.hh"
38
39#include "G4Material.hh"
40#include "G4MaterialTable.hh"
41#include "G4Element.hh"
42#include "G4ElementTable.hh"
43#include "G4Box.hh"
44#include "G4Tubs.hh"
45#include "G4Cons.hh"
46#include "G4Sphere.hh"
47#include "G4Orb.hh"
48#include "G4Trap.hh"
49#include "G4EllipticalTube.hh"
50#include "G4EllipticalCone.hh"
51#include "G4Polyhedra.hh"
52#include "G4Polycone.hh"
53#include "G4Tet.hh"
54#include "G4TwistedTubs.hh"
55#include "G4TwistedBox.hh"
56#include "G4TwistedTrd.hh"
57#include "G4TwistedTrap.hh"
58#include "G4IntersectionSolid.hh"
59#include "G4SubtractionSolid.hh"
60#include "G4UnionSolid.hh"
61#include "G4DisplacedSolid.hh"
62#include "G4LogicalVolume.hh"
63#include "G4RotationMatrix.hh"
64#include "G4ThreeVector.hh"
65#include "G4Transform3D.hh"
66#include "G4PVPlacement.hh"
67#include "G4PVReplica.hh"
68#include "G4SDManager.hh"
69#include "G4VisAttributes.hh"
70#include "G4Colour.hh"
71#include "G4VisExtent.hh"
72
73#include "G4Polyhedron.hh"
74
75MyDetectorConstruction::MyDetectorConstruction()
76{
77 new MyDetectorMessenger(this);
78
79 expHall_x = 10. * m;
80 expHall_y = 10. * m;
81 expHall_z = 10. * m;
82
83 calBox_x = 100.*cm;
84 calBox_y = 50.*cm;
85 calBox_z = 50.*cm;
86 rotAngle = 30.*deg;
87 calPos = 200.*cm;
88 calMaterialName = "Pb";
89
90 trackerRadius = 50.*cm;
91 trackerHight = 100.*cm;
92 trackerPos = -200.*cm;
93}
94
95MyDetectorConstruction::~MyDetectorConstruction()
96{
97 for (size_t i = 0; i < materialPointerStore.size(); i++) {
98 delete materialPointerStore[i];
99 }
100}
101
102G4VPhysicalVolume* MyDetectorConstruction::Construct()
103{
104
105 //------------------------------------------------------ materials
106
107 G4double a, iz, z, density;
108 G4String name, symbol;
109 G4int nel;
110
111 a = 14.01*g/mole;
112 G4Element* elN = new G4Element(name="Nitrogen", symbol="N", iz=7., a);
113 a = 16.00*g/mole;
114 G4Element* elO = new G4Element(name="Oxygen", symbol="O", iz=8., a);
115
116 density = 1.29e-03*g/cm3;
117 G4Material* Air = new G4Material(name="Air", density, nel=2);
118 materialPointerStore.push_back(Air);
119 Air->AddElement(elN, .7);
120 Air->AddElement(elO, .3);
121
122 a = 207.19*g/mole;
123 density = 11.35*g/cm3;
124 G4Material* Pb = new G4Material(name="Lead", z=82., a, density);
125 materialPointerStore.push_back(Pb);
126
127 a = 39.95*g/mole;
128 density = 1.782e-03*g/cm3;
129 G4Material* Ar = new G4Material(name="ArgonGas", z=18., a, density);
130 materialPointerStore.push_back(Ar);
131
132 a = 26.98*g/mole;
133 density = 2.7*g/cm3;
134 G4Material* Al = new G4Material(name="Aluminum", z=13., a, density);
135 materialPointerStore.push_back(Al);
136
137 a = 55.85*g/mole;
138 density = 7.87*g/cm3;
139 G4Material* Fe = new G4Material(name="Iron", z=26., a, density);
140 materialPointerStore.push_back(Fe);
141
142 //------------------------------------------------------ volumes
143
144 //------------------------------ experimental hall
145 G4Box * experimentalHall_box
146 = new G4Box("expHall_b",expHall_x,expHall_y,expHall_z);
147 G4LogicalVolume * experimentalHall_log
148 = new G4LogicalVolume(experimentalHall_box,Air,"expHall_L",0,0,0);
149 // G4VisAttributes * experimentalHallVisAtt
150 // = new G4VisAttributes(G4Colour(1.0,1.0,1.0));
151 // experimentalHallVisAtt->SetForceWireframe(true);
152 // experimentalHall_log->SetVisAttributes(experimentalHallVisAtt);
153 experimentalHall_log -> SetVisAttributes (&G4VisAttributes::GetInvisible());
154 G4VPhysicalVolume * experimentalHall_phys
155 = new G4PVPlacement(0,G4ThreeVector(),"expHall_P",
156 experimentalHall_log,0,false,0);
157
158 //------------------------------ calorimeter boxes
159 G4Box * calorimeter_box
160 = new G4Box("calorimeter_b",calBox_x,calBox_y,calBox_z);
161 G4Material* calMat;
162 if(calMaterialName=="Pb")
163 { calMat = Pb; }
164 else if(calMaterialName=="Al")
165 { calMat = Al; }
166 else if(calMaterialName=="Fe")
167 { calMat = Fe; }
168 else
169 { calMat = Air; }
170 G4LogicalVolume * calorimeter_log
171 = new G4LogicalVolume(calorimeter_box,calMat,"calo_L",0,0,0);
172 G4VisAttributes * calorimeterVisAtt
173 = new G4VisAttributes(G4Colour(0.,0.,1.));
174 // calorimeterVisAtt->SetForceWireframe(true);
175 calorimeter_log->SetVisAttributes(calorimeterVisAtt);
176 G4VPhysicalVolume* calo_phys;
177 for(G4int i=0;i<3;i++)
178 {
179 G4RotationMatrix rm;
180 rm.rotateZ(i*rotAngle);
181 rm.print(G4cout);
182 calo_phys =
183 new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(0.*cm,i*calPos,0.*cm)),
184 "calo_phys",calorimeter_log,experimentalHall_phys,
185 false,i);
186 calo_phys->GetObjectRotationValue().print(G4cout);
187 }
188
189 //------------------------------ tracker tube
190 G4Tubs * tracker_tube
191 = new G4Tubs("tracker_tube",0.*cm,trackerRadius,trackerHight,
192 0.*deg,360.*deg);
193 /*
194 G4cout << "Tracker tube polyhedron:\n"
195 << (HepPolyhedron)(*(tracker_tube->GetPolyhedron())) << G4endl;
196 */
197 G4LogicalVolume * tracker_log
198 = new G4LogicalVolume(tracker_tube,Ar,"tracker_L",0,0,0);
199 G4VisAttributes * trackerVisAtt
200 = new G4VisAttributes(G4Colour(0.,0.,1.));
201 // trackerVisAtt->SetForceWireframe(true);
202 tracker_log->SetVisAttributes(trackerVisAtt);
203 //////////////............
204 G4RotationMatrix* tracker_rm = new G4RotationMatrix;
205 tracker_rm->rotateY(-30.*deg);
206 tracker_rm->print(G4cout);
207 G4VPhysicalVolume* tracker_phys =
208 //new G4PVPlacement(tracker_rm,G4ThreeVector(0.*cm,trackerPos,200.*cm),
209 new G4PVPlacement(G4Transform3D(*tracker_rm,G4ThreeVector(0.*cm,trackerPos,200.*cm)),
210 //////////////............
211 //new G4PVPlacement(0,G4ThreeVector(0.*cm,trackerPos,0.*cm),
212 "tracker_phys",tracker_log,experimentalHall_phys,
213 false,0);
214 tracker_phys->GetObjectRotationValue().print(G4cout);
215
216 //------------------------------ displaced box
217 G4Box * undisplaced_box = new G4Box("undisplaced_box",30.*cm,50.*cm,70.*cm);
218 G4RotationMatrix rm_db;
219 rm_db.rotateZ(20.*deg);
220 G4DisplacedSolid* displaced_box = new G4DisplacedSolid
221 ("displaced_box",undisplaced_box,
222 G4Transform3D(rm_db,G4ThreeVector(200.*cm,0.,0.)));
223 G4cout << "Displaced box extent:\n" << displaced_box->GetExtent() << G4endl;
224 G4LogicalVolume * displaced_box_log
225 = new G4LogicalVolume(displaced_box,Ar,"displaced_box_L",0,0,0);
226 new G4PVPlacement(0,G4ThreeVector(0.*cm,-200.*cm,0.*cm),
227 "displaced_box_phys",displaced_box_log,
228 experimentalHall_phys,
229 false,0);
230
231 //-------------------------------------------- Boolean solids
232
233 G4Tubs* cylinder1 = new G4Tubs("Cylinder #1",20*cm,50*cm,30*cm,0,2*pi);
234 G4Box* box1 = new G4Box("Box #1",20*cm,30*cm,40*cm);
235 G4Box* box2 = new G4Box("Box #2",10*cm,20*cm,35*cm);
236 G4RotationMatrix* rm1 = new G4RotationMatrix;
237 rm1->rotateZ(20*deg);
238 G4RotationMatrix* rm2 = new G4RotationMatrix;
239 rm2->rotateZ(60*deg);
240
241 G4IntersectionSolid* cyl1Ibox1 =
242 new G4IntersectionSolid("cylinder1-intersection-box1", cylinder1, box1,
243 rm1,G4ThreeVector(30.01*cm,30.01*cm,0.01*cm));
244 G4cout << "cylinder1-intersection-box1 extent:\n"
245 << cyl1Ibox1->GetExtent() << G4endl;
246 G4IntersectionSolid* cyl1Ibox1Ibox2 =
247 new G4IntersectionSolid("cylinder1-intersection-box1-intersection-box2", cyl1Ibox1, box2,
248 rm2,G4ThreeVector(0.01,40.01*cm,0.01*cm));
249 G4cout << "cylinder1-intersection-box1-intersection-box2 extent:\n"
250 << cyl1Ibox1Ibox2->GetExtent() << G4endl;
251 G4LogicalVolume * intersection_log
252 = new G4LogicalVolume(cyl1Ibox1Ibox2,Ar,"intersection_L",0,0,0);
253 const G4VisAttributes* bool_red =
254 new G4VisAttributes(G4Colour(1.,0.,0.));
255 intersection_log->SetVisAttributes(bool_red);
256 new G4PVPlacement(0,G4ThreeVector(100.*cm,-50*cm,0.*cm),
257 "intersection_phys",intersection_log,experimentalHall_phys,
258 true,0);
259
260 G4SubtractionSolid* cyl1Sbox1 =
261 new G4SubtractionSolid("cylinder1-subtraction-box1", cylinder1, box1,
262 rm1,G4ThreeVector(30.*cm,30.*cm,1.*cm));
263 G4SubtractionSolid* cyl1Sbox1Sbox2 =
264 new G4SubtractionSolid("cylinder1-subtraction-box1-subtraction-box2", cyl1Sbox1, box2,
265 rm2,G4ThreeVector(0.,40*cm,2.*cm));
266 G4LogicalVolume * subtraction_log
267 = new G4LogicalVolume(cyl1Sbox1Sbox2,Ar,"subtraction_L",0,0,0);
268 const G4VisAttributes* bool_green =
269 new G4VisAttributes(G4Colour(0.,1.,0.));
270 subtraction_log->SetVisAttributes(bool_green);
271 new G4PVPlacement(0,G4ThreeVector(200.*cm,-50*cm,0.*cm),
272 "subtraction_phys",subtraction_log,experimentalHall_phys,
273 true,0);
274
275 G4UnionSolid* cyl1Ubox1 =
276 new G4UnionSolid("cylinder1-union-box1", cylinder1, box1,
277 rm1,G4ThreeVector(30.*cm,30.*cm,1.*cm));
278 G4UnionSolid* cyl1Ubox1Ubox2 =
279 new G4UnionSolid("cylinder1-union-box1-union-box2", cyl1Ubox1, box2,
280 rm2,G4ThreeVector(0.,40*cm,2.*cm));
281 G4LogicalVolume * union_log
282 = new G4LogicalVolume(cyl1Ubox1Ubox2,Ar,"union_L",0,0,0);
283 const G4VisAttributes* bool_blue =
284 new G4VisAttributes(G4Colour(0.,0.,1.));
285 union_log->SetVisAttributes(bool_blue);
286 new G4PVPlacement(0,G4ThreeVector(350.*cm,-50*cm,0.*cm),
287 "union_phys",union_log,experimentalHall_phys,
288 true,0);
289
290 //----------- Tubes, replicas(!?) and daughter boxes
291
292 const G4double eps = 10 * mm;
293 const G4double alp = 10 * mrad;
294
295 G4VisAttributes * grey
296 = new G4VisAttributes(G4Colour(0.5,0.5,0.5,0.1));
297
298
299 G4double tube_dPhi = twopi;
300 G4Tubs* tube = new G4Tubs("tube",20*cm,50*cm,30*cm,0.,tube_dPhi);
301 G4LogicalVolume * tube_log
302 = new G4LogicalVolume(tube,Ar,"tube_L",0,0,0);
303 G4VisAttributes * tube_VisAtt
304 = new G4VisAttributes(G4Colour(0.,1.,0.,0.1));
305 // tube_log->SetVisAttributes(tube_VisAtt);
306 tube_log->SetVisAttributes(G4VisAttributes::GetInvisible());
307 new G4PVPlacement(0,G4ThreeVector(-200.*cm,0.,0.*cm),
308 "tube_phys",tube_log,experimentalHall_phys,
309 false,0);
310
311 G4double divided_tube_dPhi = tube_dPhi / 6.;
312 G4Tubs* divided_tube = new G4Tubs
313 ("divided_tube",
314 20*cm,50*cm,30*cm,-divided_tube_dPhi/2.,divided_tube_dPhi);
315 G4LogicalVolume * divided_tube_log
316 = new G4LogicalVolume(divided_tube,Ar,"divided_tube_L",0,0,0);
317 divided_tube_log->SetVisAttributes(grey);
318 new G4PVReplica("divided_tube_phys",divided_tube_log,tube_log,
319 kPhi,6,divided_tube_dPhi);
320 /************
321 G4int iCopy;
322 for (iCopy = 0; iCopy < 6; iCopy++) {
323 new G4PVPlacement
324 (G4Transform3D
325 (G4RotationMatrix().rotateZ(divided_tube_dPhi/2.+iCopy*pi/3.),
326 G4ThreeVector()),
327 divided_tube_log,"divided_tube_phys",tube_log,
328 false,iCopy);
329 }
330 *********/
331
332 G4double divided_tube_inset_dPhi = divided_tube_dPhi - 2. * alp;
333 G4Tubs* divided_tube_inset = new G4Tubs
334 ("divided_tube_inset",
335 20*cm+eps,50*cm-eps,30*cm-eps,
336 -divided_tube_inset_dPhi/2.,
337 divided_tube_inset_dPhi);
338 G4LogicalVolume * divided_tube_inset_log
339 = new G4LogicalVolume(divided_tube_inset,Ar,"divided_tube_inset_L",0,0,0);
340 G4VisAttributes * divided_tube_inset_VisAtt
341 = new G4VisAttributes(G4Colour(1.,0.,0.,0.2));
342 // divided_tube_inset_log->SetVisAttributes(divided_tube_inset_VisAtt);
343 divided_tube_inset_log->SetVisAttributes(G4VisAttributes::Invisible);
344 new G4PVPlacement(0,G4ThreeVector(),
345 divided_tube_inset_log,"divided_tube_inset_phys",
346 divided_tube_log,false,0);
347
348 G4double sub_divided_tube_dPhi = divided_tube_inset_dPhi / 4.;
349 G4Tubs* sub_divided_tube = new G4Tubs
350 ("sub_divided_tube",
351 20*cm+eps,50*cm-eps,30*cm-eps,
352 -sub_divided_tube_dPhi/2.,sub_divided_tube_dPhi);
353 G4LogicalVolume * sub_divided_tube_log
354 = new G4LogicalVolume(sub_divided_tube,Ar,"sub_divided_tube_L",0,0,0);
355 sub_divided_tube_log->SetVisAttributes(grey);
356 new G4PVReplica("sub_divided_tube_phys",
357 sub_divided_tube_log,divided_tube_inset_log,
358 kPhi,4,sub_divided_tube_dPhi,-divided_tube_inset_dPhi/2.);
359 /************
360 for (iCopy = 0; iCopy < 4; iCopy++) {
361 new G4PVPlacement
362 (G4Transform3D
363 (G4RotationMatrix().rotateZ
364 (-divided_tube_inset_dPhi/2.
365 +(iCopy+0.5)*sub_divided_tube_dPhi),
366 G4ThreeVector()),
367 sub_divided_tube_log,"sub_divided_tube_phys",divided_tube_inset_log,
368 false,iCopy);
369 }
370 ************/
371
372 G4Box* daughter_box = new G4Box("daughter_box",4.*cm,3.*cm,25.*cm);
373 G4LogicalVolume * daughter_box_log
374 = new G4LogicalVolume(daughter_box,Ar,"daughter_box_L",0,0,0);
375 G4VisAttributes * daughter_box_VisAtt
376 = new G4VisAttributes(G4Colour(0.,0.,1.,0.3));
377 daughter_box_log->SetVisAttributes(daughter_box_VisAtt);
378 G4Box* grand_daughter_box = new G4Box("grand_daughter_box",1*cm,2*cm,5*cm);
379 G4LogicalVolume * grand_daughter_box_log
380 = new G4LogicalVolume(grand_daughter_box,Ar,"grand_daughter_box_L",0,0,0);
381 G4VisAttributes * grand_daughter_box_VisAtt
382 = new G4VisAttributes(G4Colour(1.,1.,0.));
383 grand_daughter_box_log->SetVisAttributes(grand_daughter_box_VisAtt);
384 G4Box* grand_daughter_box2 = new G4Box("grand_daughter_box2",1*cm,2*cm,5*cm);
385 G4LogicalVolume * grand_daughter_box2_log
386 = new G4LogicalVolume
387 (grand_daughter_box2,Ar,"grand_daughter_box2_L",0,0,0);
388 G4VisAttributes * grand_daughter_box2_VisAtt
389 = new G4VisAttributes(G4Colour(1.,0.,1.));
390 grand_daughter_box2_log->SetVisAttributes(grand_daughter_box2_VisAtt);
391 new G4PVPlacement(0,G4ThreeVector(-2*cm,0.,0.),
392 grand_daughter_box_log,"grand_daughter_box_phys",
393 daughter_box_log,false,0);
394 new G4PVPlacement(0,G4ThreeVector(2*cm,0.,0.),
395 grand_daughter_box2_log,"grand_daughter_box2_phys",
396 daughter_box_log,false,0);
397 new G4PVPlacement(0,G4ThreeVector(40*cm,0.,0.),
398 daughter_box_log,"daughter_box_phys",
399 sub_divided_tube_log,false,0);
400
401 //------------------------------------------------ Extra placements
402 // For good measure, as a test of drawn volume path, place in one of
403 // the earlier volumes...
404 new G4PVPlacement(0,G4ThreeVector(),
405 "divided_tube_inset_phys",divided_tube_inset_log,
406 calo_phys,false,0); // Place in PV.
407 new G4PVPlacement(0,G4ThreeVector(),
408 grand_daughter_box2_log,"grand_daughter_box2_phys",
409 calorimeter_log,false,0); // Place in LV. Same effect.
410
411 //------------------------------------------------ sensitive detectors
412 /****************
413 G4SDManager* SDman = G4SDManager::GetSDMpointer();
414
415 G4String calorimeterSDname = "example2/calorimeter";
416 MyCalorimeterSD * myCalorimeterSD = new MyCalorimeterSD( calorimeterSDname );
417 SDman->AddNewDetector( myCalorimeterSD );
418 calorimeter_log->SetSensitiveDetector( myCalorimeterSD );
419
420 G4String trackerSDname = "example2/tracker";
421 MyTrackerSD * myTrackerSD = new MyTrackerSD( trackerSDname );
422 SDman->AddNewDetector( myTrackerSD );
423 tracker_log->SetSensitiveDetector( myTrackerSD );
424 ******************/
425
426 //-------------------------------------------- Sphere
427
428 G4Sphere* PD_vol_crystal
429 = new G4Sphere("Test_Sphere",
430 100.*cm, // inner radius
431 200.*cm, // outer radius
432 0., // start phi
433 90.*degree, // delta phi
434 0., // start theta
435 90.*degree // delta theta
436 );
437
438
439 G4LogicalVolume* PD_log_crystal
440 = new G4LogicalVolume(PD_vol_crystal,Ar,"Test_Sphere");
441
442 G4VisAttributes * PD_att_crystal
443 = new G4VisAttributes(G4Colour(1.,0.,1.));
444 PD_att_crystal->SetForceAuxEdgeVisible(true);
445 PD_att_crystal->SetForceLineSegmentsPerCircle(100);
446 PD_log_crystal->SetVisAttributes(PD_att_crystal);
447
448 G4RotationMatrix rm;
449
450 new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(200.*cm,300.*cm,0)),
451 "PD_physical", PD_log_crystal,
452 experimentalHall_phys,false,0);
453
454//-------------------------------------------- Polyhedra and Polycone
455
456 const G4int numRZ = 10;
457 G4double polyhedra_r[] = {0,5,3,4,9,9,3,3,2,0};
458 G4double polyhedra_z[] = {0,1,2,3,0,5,4,3,2,1};
459 for (int i = 0; i < numRZ; ++i) {
460 polyhedra_r[i] *= 10*cm;
461 polyhedra_z[i] *= 10.*cm;
462 }
463
464 G4Polyhedra* polyhedra_solid
465 = new G4Polyhedra("Polyhedra_Test",
466 0.,270.*deg,6,numRZ,polyhedra_r,polyhedra_z);
467 //0.,twopi,6,numRZ,polyhedra_r,polyhedra_z);
468 G4cout << polyhedra_solid->StreamInfo(G4cout) << G4endl;
469
470 G4LogicalVolume* polyhedra_log
471 = new G4LogicalVolume(polyhedra_solid,Ar,"Polyhedra_Test");
472
473 G4VisAttributes * polyhedra_atts
474 = new G4VisAttributes(G4Colour(0.,1.,1.));
475 polyhedra_log->SetVisAttributes(polyhedra_atts);
476
477 new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(200.*cm,100.*cm,0)),
478 "Polyhedra_Test", polyhedra_log,
479 experimentalHall_phys,false,0);
480
481 const G4int numRZ1 = 10;
482 G4double polycone_r[] = {1,5,3,4,9,9,3,3,2,1};
483 G4double polycone_z[] = {0,1,2,3,0,5,4,3,2,1};
484 for (int i = 0; i < numRZ1; ++i) {
485 polycone_r[i] *= 10*cm;
486 polycone_z[i] *= 10.*cm;
487 }
488
489 G4Polycone* polycone_solid
490 = new G4Polycone("Polycone_Test",
491 0.*deg,270.*deg,numRZ1,polycone_r,polycone_z);
492 //0.,twopi,numRZ,polycone_r,polycone_z);
493 G4cout << polycone_solid->StreamInfo(G4cout) << G4endl;
494
495 G4LogicalVolume* polycone_log
496 = new G4LogicalVolume(polycone_solid,Ar,"Polycone_Test");
497
498 G4VisAttributes * polycone_atts
499 = new G4VisAttributes(G4Colour(1.,0.5,0.5));
500 polycone_log->SetVisAttributes(polycone_atts);
501
502 new G4PVPlacement(G4Transform3D(rm,G4ThreeVector(400.*cm,120.*cm,0)),
503 "Polycone_Test", polycone_log,
504 experimentalHall_phys,false,0);
505
506 //-------------------------------------------- Orb
507
508 G4Orb* orb
509 = new G4Orb("Test orb", 100.*cm);
510
511 G4LogicalVolume* orb_log
512 = new G4LogicalVolume(orb, Ar,"Test_orb");
513
514 G4VisAttributes * orb_att
515 = new G4VisAttributes(G4Colour(1.,0.,1.));
516 orb_att->SetForceAuxEdgeVisible(true);
517 orb_log->SetVisAttributes(orb_att);
518
519 new G4PVPlacement(G4Translate3D(-300.*cm,200.*cm,0),
520 "Test orb", orb_log,
521 experimentalHall_phys,false,0);
522
523 //-------------------------------------------- Trapezoid
524
525 G4Trap* trap1_solid = new G4Trap
526 ("trap1_solid",
527 50.*cm, // pDz Half-length along the z-axis
528 0., // pTheta Polar angle of the line joining the centres of the faces
529 // at -/+pDz
530 0., // pPhi Azimuthal angle of the line joing the centre of the face at
531 // -pDz to the centre of the face at +pDz
532 40.*cm, // pDy1 Half-length along y of the face at -pDz
533 30.*cm, // pDx1 Half-length along x of the side at y=-pDy1 of the face at -pDz
534 30.*cm, // pDx2 Half-length along x of the side at y=+pDy1 of the face at -pDz
535 20.*deg, // pAlp1 Angle with respect to the y axis from the centre of the side
536 // at y=-pDy1 to the centre at y=+pDy1 of the face at -pDz
537 40.*cm, // pDy2 Half-length along y of the face at +pDz
538 30.*cm, // pDx3 Half-length along x of the side at y=-pDy2 of the face at +pDz
539 30.*cm, // pDx4 Half-length along x of the side at y=+pDy2 of the face at +pDz
540 20.*deg // pAlp2 Angle with respect to the y axis from the centre of the side
541 // at y=-pDy2 to the centre at y=+pDy2 of the face at +pDz
542 );
543 G4LogicalVolume* trap1_log = new G4LogicalVolume (trap1_solid,Ar,"trap1_log");
544 new G4PVPlacement
545 (G4Transform3D(rm,G4ThreeVector(-200.*cm,200.*cm,-200.*cm)),
546 "trap1_phys", trap1_log,
547 experimentalHall_phys,false,0);
548
549 G4Trap* trap2_solid = new G4Trap
550 ("trap2_solid",
551 50.*cm, // pDz Half-length along the z-axis
552 20.*deg, // pTheta Polar angle of the line joining the centres of the faces
553 // at -/+pDz
554 90.*deg, // pPhi Azimuthal angle of the line joing the centre of the face at
555 // -pDz to the centre of the face at +pDz
556 40.*cm, // pDy1 Half-length along y of the face at -pDz
557 30.*cm, // pDx1 Half-length along x of the side at y=-pDy1 of the face at -pDz
558 30.*cm, // pDx2 Half-length along x of the side at y=+pDy1 of the face at -pDz
559 0., // pAlp1 Angle with respect to the y axis from the centre of the side
560 // at y=-pDy1 to the centre at y=+pDy1 of the face at -pDz
561 40.*cm, // pDy2 Half-length along y of the face at +pDz
562 30.*cm, // pDx3 Half-length along x of the side at y=-pDy2 of the face at +pDz
563 30.*cm, // pDx4 Half-length along x of the side at y=+pDy2 of the face at +pDz
564 0. // pAlp2 Angle with respect to the y axis from the centre of the side
565 // at y=-pDy2 to the centre at y=+pDy2 of the face at +pDz
566 );
567 G4LogicalVolume* trap2_log = new G4LogicalVolume (trap2_solid,Ar,"trap2_log");
568 new G4PVPlacement
569 (G4Transform3D(rm,G4ThreeVector(-200.*cm,400.*cm,-200.*cm)),
570 "trap2_phys", trap2_log,
571 experimentalHall_phys,false,0);
572
573 //-------------------------------------------- Elliptical Tube
574 G4VSolid* eTube = new G4EllipticalTube("e-tube",100.*cm,50.*cm,100.*cm);
575 G4LogicalVolume* eTubeLog = new G4LogicalVolume(eTube,Ar,"e-tube-log");
576 new G4PVPlacement(G4Translate3D(G4ThreeVector(-400.*cm,0.,0)),
577 "e-tube-phys", eTubeLog,
578 experimentalHall_phys,false,0);
579
580 //-------------------------------------------- Elliptical Cone
581 G4VSolid* eCone = new G4EllipticalCone("e-cone",50.*cm,100.*cm,100.*cm,50.*cm);
582 //G4VSolid* eCone = new G4EllipticalCone("e-cone",1.*mm,0.5*mm,40.*mm,20.*mm);
583 G4LogicalVolume* eConeLog = new G4LogicalVolume(eCone,Ar,"e-cone-log");
584 new G4PVPlacement(G4Translate3D(G4ThreeVector(-500.*cm,400.*cm,0)),
585 "e-cone-phys", eConeLog,
586 experimentalHall_phys,false,0);
587
588 //-------------------------------------------- G4Cons
589 G4VSolid* eCons = new G4Cons("e-cons",50.*cm,70.*cm,100.*cm,140.*cm,200.*cm,0,twopi);
590 G4LogicalVolume* eConsLog = new G4LogicalVolume(eCons,Ar,"e-cons-log");
591 new G4PVPlacement(G4Translate3D(G4ThreeVector(-300.*cm,-200.*cm,60*cm)),
592 "e-cons-phys", eConsLog,
593 experimentalHall_phys,false,0);
594
595 /*
596 //--------------------------- Simple shared logical volume tree
597
598 G4VSolid* boxTwo = new G4Box("Box2", 100.*cm, 100.*cm, 100.*cm);
599 G4LogicalVolume* log2 = new G4LogicalVolume(boxTwo,Ar,"Log2");
600 new G4PVPlacement(G4Translate3D(G4ThreeVector(600.*cm,0.,0.)),
601 "B",log2,
602 experimentalHall_phys,false,0);
603 new G4PVPlacement(G4Translate3D(G4ThreeVector(900.*cm,0.,0.)),
604 "C",log2,
605 experimentalHall_phys,false,0);
606 G4VSolid* boxThree = new G4Box("Box3", 50.*cm, 50.*cm, 50.*cm);
607 G4LogicalVolume* log3 = new G4LogicalVolume(boxThree,Ar,"Log3");
608 new G4PVPlacement(G4Translate3D(),
609 log3,"D",
610 log2,false,0);
611 G4VSolid* boxFour = new G4Box("Box3", 30.*cm, 30.*cm, 30.*cm);
612 G4LogicalVolume* log4 = new G4LogicalVolume(boxFour,Ar,"Log4");
613 new G4PVPlacement(G4Translate3D(),
614 log4,"E",
615 log3,false,0);
616 */
617
618 /*
619 //--------------------------- Boolean for logical volume test
620
621 G4double alSLayer_x = 30.*cm;
622 G4double alSLayer_y = 50.*cm;
623 G4double alSLayer_z = 2.5*cm;
624
625 G4Box* box1B = new
626 G4Box("alSLayer_box",alSLayer_x/2.,alSLayer_y/2.,alSLayer_z/2.);
627
628 G4Tubs* Cylinder1 = new G4Tubs("Cylinder#1",0.*mm,alSLayer_x/2.,alSLayer_x/2.,0.,2.*M_PI);
629
630 G4UnionSolid* b1UnionC1 = new G4UnionSolid("Box+Cylinder", box1B,
631 Cylinder1);
632
633 G4LogicalVolume* alSpaceCraft_log = new
634 G4LogicalVolume(b1UnionC1,Al,"alLayer_log",0,0,0);
635
636 new G4PVPlacement(0,G4ThreeVector(900.*cm, 200.*cm, 0.),alSpaceCraft_log,"alLayer_phys",experimentalHall_log,false,0);
637 // new G4PVPlacement(G4Translate3D(G4ThreeVector(900.*cm, 200.*cm, 0.)),"alLayer_phys",alSpaceCraft_log,experimentalHall_log,false,0);
638 */
639
640 //----------- Radially replicated tube sector
641
642 G4double rMin = 50.*cm;
643 G4double DeltaR = 50.*cm;
644 G4VSolid* rrTubs = new G4Tubs
645 ("rrTubs",rMin,rMin + DeltaR,200*cm,180*deg,90*deg);
646 G4LogicalVolume* rrTubsLog = new G4LogicalVolume
647 (rrTubs,Ar,"rrTubs-log");
648 rrTubsLog->SetVisAttributes(G4VisAttributes::Invisible);
649 new G4PVPlacement(G4Translate3D(G4ThreeVector(400.*cm,-200.*cm,0)),
650 "rrTubs-phys", rrTubsLog,
651 experimentalHall_phys,false,0);
652 G4double deltaR = DeltaR / 6.;
653 G4Tubs* drrTubs = new G4Tubs
654 ("drrTubs",rMin,rMin + deltaR,200*cm,180*deg,90*deg);
655 G4LogicalVolume * drrTubsLog = new G4LogicalVolume
656 (drrTubs,Ar,"drrTubs-log");
657 drrTubsLog->SetVisAttributes(G4Colour::Red());
658 new G4PVReplica("drrTubs-phys",drrTubsLog,rrTubsLog,
659 kRho,6,deltaR,rMin);
660
661
662 //----------- tetrahedron
663
664 G4VSolid* tet = new G4Tet
665 ("tet",
666 G4ThreeVector(),
667 G4ThreeVector(0.,100.*cm,0.),
668 G4ThreeVector(0.,0.,100.*cm),
669 G4ThreeVector(100*cm,0.,0.));
670 G4LogicalVolume* tet_log = new G4LogicalVolume
671 (tet,Ar,"tet-log");
672 tet_log->SetVisAttributes(G4VisAttributes(G4Colour(0.,1.,0.)));
673 new G4PVPlacement
674 (G4Translate3D(G4ThreeVector(300.*cm,-400.*cm,0.)),
675 "tet-phys", tet_log,
676 experimentalHall_phys,false,0);
677
678 //----------- Twisted solids
679
680 G4VSolid* aVolume;
681 G4LogicalVolume* aLog;
682 G4double fTrackerR1 ; // r1
683 G4double fTrackerR2 ; // r2
684 G4double fTrackerpDz ; // Full length of Tracker (pDz)
685 G4double fTrackerpDx1 ; // twisted Trapezoid
686 G4double fTrackerpDx2 ;
687 G4double fTrackerpDy1 ;
688 G4double fTrackerpDy2 ;
689 G4double fTwistAngle ;
690 G4double fPhi ;
691
692 G4double myScale = 5.;
693
694 fTwistAngle = 90*deg ;
695 fTrackerpDz = myScale*20*cm ;
696 fTrackerR1 = myScale*7*cm ;
697 fTrackerR2 = myScale*10*cm ;
698 fPhi = 180*deg ;
699
700 aVolume = new G4TwistedTubs
701 ("aTwistedTubs", fTwistAngle, fTrackerR1, fTrackerR2, fTrackerpDz, fPhi ) ;
702 //G4cout << "aTwistedTubs: volume: " << aVolume->GetCubicVolume() << G4endl;
703 aLog = new G4LogicalVolume
704 (aVolume,Ar,"aTwistedTubs-log");
705 aLog->SetVisAttributes(G4VisAttributes(G4Colour(1.,1.,0.)));
706 // Comment out this G4PVPlacement to get voxelisation crash...
707 new G4PVPlacement
708 (G4Translate3D(G4ThreeVector(200.*cm,-400.*cm,0.)),
709 "aTwistedTubs-phys", aLog, experimentalHall_phys,false,0);
710
711 fTwistAngle = 50*deg ;
712 fTrackerpDx1 = myScale*4*cm ;
713 fTrackerpDy1 = myScale*6*cm ;
714 fTrackerpDz = myScale*15*cm ;
715
716 aVolume = new G4TwistedBox
717 ("aTwistedBox",fTwistAngle,fTrackerpDx1,fTrackerpDy1,fTrackerpDz) ;
718 aLog = new G4LogicalVolume
719 (aVolume,Ar,"aTwistedBox-log");
720 aLog->SetVisAttributes(G4VisAttributes(G4Colour(1.,1.,0.)));
721 // Comment out this G4PVPlacement to get voxelisation crash...
722 new G4PVPlacement
723 (G4Translate3D(G4ThreeVector(100.*cm,-400.*cm,0.)),
724 "aTwistedBox-phys", aLog, experimentalHall_phys,false,0);
725
726 fTrackerpDx1 = myScale*4*cm ;
727 fTrackerpDx2 = myScale*7*cm ;
728 fTrackerpDy1 = myScale*2*cm ;
729 fTrackerpDy2 = myScale*4*cm ;
730 fTrackerpDz = myScale*15*cm ;
731 fTwistAngle = 50*deg ;
732
733 aVolume = new G4TwistedTrd
734 ("aTwistedTrd",fTrackerpDx1,fTrackerpDx2,fTrackerpDy1,fTrackerpDy2,
735 fTrackerpDz,fTwistAngle);
736 aLog = new G4LogicalVolume
737 (aVolume,Ar,"aTwistedTrd-log");
738 aLog->SetVisAttributes(G4VisAttributes(G4Colour(1.,1.,0.)));
739 new G4PVPlacement
740 (G4Translate3D(G4ThreeVector(000.*cm,-400.*cm,0.)),
741 "aTwistedTrd-phys", aLog, experimentalHall_phys,false,0);
742
743 fTrackerpDx1 = myScale*4*cm ;
744 fTrackerpDx2 = myScale*7*cm ;
745 fTrackerpDy1 = myScale*2*cm ;
746 fTrackerpDy2 = myScale*4*cm ;
747 fTrackerpDz = myScale*15*cm ;
748 fTwistAngle = 50*deg ;
749
750 aVolume = new G4TwistedTrd
751 ("anotherTwistedTrd",fTrackerpDx1,fTrackerpDx2,fTrackerpDy1,fTrackerpDy2,
752 fTrackerpDz,fTwistAngle);
753 aLog = new G4LogicalVolume
754 (aVolume,Ar,"anotherTwistedTrd-log");
755 aLog->SetVisAttributes(G4VisAttributes(G4Colour(1.,1.,0.)));
756 new G4PVPlacement
757 (G4Translate3D(G4ThreeVector(-100.*cm,-400.*cm,0.)),
758 "anotherTwistedTrd-phys", aLog, experimentalHall_phys,false,0);
759
760 //-------------------------------------------- return
761 return experimentalHall_phys;
762}
Note: See TracBrowser for help on using the repository browser.