ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4MIRDUpperSpine.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4MIRDUpperSpine.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 "G4MIRDUpperSpine.hh"
31 
32 #include "globals.hh"
33 #include "G4SystemOfUnits.hh"
34 #include "G4SDManager.hh"
35 #include "G4VisAttributes.hh"
36 #include "G4VisAttributes.hh"
38 #include "G4EllipticalTube.hh"
39 #include "G4RotationMatrix.hh"
40 #include "G4ThreeVector.hh"
41 #include "G4VPhysicalVolume.hh"
42 #include "G4PVPlacement.hh"
43 #include "G4Box.hh"
44 #include "G4SubtractionSolid.hh"
45 #include "G4HumanPhantomColour.hh"
46 
48 {
49 }
50 
52 {
53 
54 }
55 
56 
58  G4VPhysicalVolume* mother,
59  const G4String& colourName
60  , G4bool wireFrame,G4bool)
61 {
63 
64  G4cout<<"Construct "<<volumeName<<" with mother volume "<<mother->GetName()<<G4endl;
65 
66 
67  G4Material* skeleton = material -> GetMaterial("skeleton");
68 
69  delete material;
70 
71  G4double dx = 2. *cm;
72  G4double dy = 2.5 *cm;
73  G4double dz = 4.25*cm;
74 
75  G4EllipticalTube* upperSpine = new G4EllipticalTube("UpperSpine",dx, dy, dz);
76 
77  G4double xx = 20. * cm;
78  G4double yy = 10. * cm;
79  G4double zz = 5. * cm;
80 
81  G4Box* subtraction = new G4Box("box", xx/2., yy/2., zz/2.);
82 
83  G4RotationMatrix* matrix = new G4RotationMatrix();
84  matrix -> rotateX(-25.* deg);
85 
86  G4SubtractionSolid* upper_spine = new G4SubtractionSolid("upperspine",upperSpine, subtraction,
87  matrix, G4ThreeVector(0., -2.5 * cm, 5.5* cm));
88 
89  G4LogicalVolume* logicUpperSpine = new G4LogicalVolume(upper_spine, skeleton,
90  "logical" + volumeName,
91  0, 0, 0);
92  // Define rotation and position here!
93  G4VPhysicalVolume* physUpperSpine = new G4PVPlacement(0,
94  G4ThreeVector(0.0, 5.5 *cm, -3.5 *cm),
95  "physicalUpperSpine",
96  logicUpperSpine,
97  mother,
98  false,
99  0, true);
100 
101  // Visualization Attributes
102  //G4VisAttributes* UpperSpineVisAtt = new G4VisAttributes(G4Colour(0.46,0.53,0.6));
103  G4HumanPhantomColour* colourPointer = new G4HumanPhantomColour();
104  G4Colour colour = colourPointer -> GetColour(colourName);
105  G4VisAttributes* UpperSpineVisAtt = new G4VisAttributes(colour);
106 
107  UpperSpineVisAtt->SetForceSolid(wireFrame);
108  logicUpperSpine->SetVisAttributes(UpperSpineVisAtt);
109 
110  G4cout << "UpperSpine created !!!!!!" << G4endl;
111 
112  // Testing UpperSpine Volume
113  G4double UpperSpineVol = logicUpperSpine->GetSolid()->GetCubicVolume();
114  G4cout << "Volume of UpperSpine = " << UpperSpineVol/cm3 << " cm^3" << G4endl;
115 
116  // Testing UpperSpine Material
117  G4String UpperSpineMat = logicUpperSpine->GetMaterial()->GetName();
118  G4cout << "Material of UpperSpine = " << UpperSpineMat << G4endl;
119 
120  // Testing Density
121  G4double UpperSpineDensity = logicUpperSpine->GetMaterial()->GetDensity();
122  G4cout << "Density of Material = " << UpperSpineDensity*cm3/g << " g/cm^3" << G4endl;
123 
124  // Testing Mass
125  G4double UpperSpineMass = (UpperSpineVol)*UpperSpineDensity;
126  G4cout << "Mass of UpperSpine = " << UpperSpineMass/gram << " g" << G4endl;
127 
128 
129  return physUpperSpine;
130 }