ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4MIRDMaleGenitalia.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4MIRDMaleGenitalia.cc
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 // Previous authors: G. Guerrieri, S. Guatelli and M. G. Pia, INFN Genova, Italy
27 // Authors (since 2007): S. Guatelli, University of Wollongong, Australia
28 //
29 //
30 #include "G4MIRDMaleGenitalia.hh"
31 
32 #include "globals.hh"
33 #include "G4SystemOfUnits.hh"
34 #include "G4SDManager.hh"
35 #include "G4VisAttributes.hh"
36 #include "G4Ellipsoid.hh"
37 #include "G4ThreeVector.hh"
38 #include "G4VPhysicalVolume.hh"
39 #include "G4RotationMatrix.hh"
40 #include "G4Material.hh"
41 #include "G4LogicalVolume.hh"
43 #include "G4VPhysicalVolume.hh"
44 #include "G4PVPlacement.hh"
45 #include "G4UnionSolid.hh"
46 #include "G4HumanPhantomColour.hh"
47 #include "G4Trap.hh"
48 #include "G4Cons.hh"
49 #include "G4SubtractionSolid.hh"
50 
52 {
53 }
54 
56 {
57 
58 }
59 
60 
62  const G4String& colourName, G4bool wireFrame, G4bool)
63 {
64  G4cout<<"Construct "<<volumeName<<" with mother volume "<<mother->GetName()<<G4endl;
65 
67  G4Material* soft = material -> GetMaterial("soft_tissue");
68  delete material;
69 
70 
71  G4double pDz=2.4*cm;
72  G4double pTheta=0*degree;
73  G4double pPhi=0*degree;
74  G4double pDy1=4.76*cm;
75  G4double pDx1=9.52*cm;
76  G4double pDx2=9.52*cm;
77  G4double pAlp1=0*degree;
78  G4double pDy2=5*cm;
79  G4double pDx3=10*cm;
80  G4double pDx4=10*cm;
81  G4double pAlp2=0*degree;
82 
83  G4Trap* genitaliaTrap= new G4Trap("GenitaliaTrap",
84  pDz,pTheta,pPhi,pDy1,
85  pDx1,pDx2,pAlp1,pDy2,
86  pDx3,pDx4,pAlp2);
87 
88 
89  G4double rmin1 = 0.* cm;
90  G4double rmin2 = 0.* cm;
91  G4double dz= 5 * cm;
92  G4double rmax1= 9.51 * cm;
93  G4double rmax2= 10.01 * cm;
94  G4double startphi= 0.* degree;
95  G4double deltaphi= 360. * degree;
96 
97  G4Cons* genitaliaLegL = new G4Cons("GenitaliaLegL",
98  rmin1, rmax1,
99  rmin2, rmax2, dz/2.,
100  startphi, deltaphi);
101 
102  G4Cons* genitaliaLegR = new G4Cons("GenitaliaLegR",
103  rmin1, rmax1,
104  rmin2, rmax2, dz/2.,
105  startphi, deltaphi);
106 
107  G4UnionSolid* genitaliaLegs = new G4UnionSolid("GenitaliaLegs",genitaliaLegL,genitaliaLegR,
108  0,
109  G4ThreeVector(20.* cm, 0.*cm,0* cm) );
110 
111 
112 
113  G4SubtractionSolid* MaleGenitalia = new G4SubtractionSolid("MaleGenitalia",genitaliaTrap,genitaliaLegs,
114  0,//
115  G4ThreeVector(-10.* cm, -5.*cm,0* cm) );
116 
117 
118  G4LogicalVolume* logicMaleGenitalia = new G4LogicalVolume(MaleGenitalia,
119  soft,
120  "logical" + volumeName,
121  0, 0, 0);
122 
123  // Define rotation and position here!
124  G4VPhysicalVolume* physMaleGenitalia = new G4PVPlacement(0,
125  G4ThreeVector(0*cm,5.*cm, -2.4*cm),
126  "physicalMaleGenitalia",
127  logicMaleGenitalia,
128  mother,
129  false,
130  0, true);
131 
132  // Visualization Attributes
133  //G4VisAttributes* MaleGenitaliaVisAtt = new G4VisAttributes(G4Colour(0.85,0.44,0.84));
134  G4HumanPhantomColour* colourPointer = new G4HumanPhantomColour();
135  G4Colour colour = colourPointer -> GetColour(colourName);
136  delete colourPointer;
137 
138  G4VisAttributes* MaleGenitaliaVisAtt = new G4VisAttributes(colour);
139  MaleGenitaliaVisAtt->SetForceSolid(wireFrame);
140  logicMaleGenitalia->SetVisAttributes(MaleGenitaliaVisAtt);
141 
142  G4cout << "MaleGenitalia created !!!!!!" << G4endl;
143 
144  // Testing MaleGenitalia Volume
145  G4double MaleGenitaliaVol = logicMaleGenitalia->GetSolid()->GetCubicVolume();
146  G4cout << "Volume of MaleGenitalia = " << MaleGenitaliaVol/cm3 << " cm^3" << G4endl;
147 
148  // Testing MaleGenitalia Material
149  G4String MaleGenitaliaMat = logicMaleGenitalia->GetMaterial()->GetName();
150  G4cout << "Material of MaleGenitalia = " << MaleGenitaliaMat << G4endl;
151 
152  // Testing Density
153  G4double MaleGenitaliaDensity = logicMaleGenitalia->GetMaterial()->GetDensity();
154  G4cout << "Density of Material = " << MaleGenitaliaDensity*cm3/g << " g/cm^3" << G4endl;
155 
156  // Testing Mass
157  G4double MaleGenitaliaMass = (MaleGenitaliaVol)*MaleGenitaliaDensity;
158  G4cout << "Mass of MaleGenitalia = " << MaleGenitaliaMass/gram << " g" << G4endl;
159 
160  return physMaleGenitalia;
161 }