Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 // Previous authors: G. Guerrieri, S. Guatelli 28 // Authors (since 2007): S. Guatelli, Universi 29 // 30 31 #include "G4MIRDRightScapula.hh" 32 33 #include "globals.hh" 34 #include "G4SystemOfUnits.hh" 35 #include "G4SDManager.hh" 36 #include "G4Cons.hh" 37 #include "G4VisAttributes.hh" 38 #include "G4HumanPhantomMaterial.hh" 39 #include "G4EllipticalTube.hh" 40 #include "G4ThreeVector.hh" 41 #include "G4VPhysicalVolume.hh" 42 #include "G4RotationMatrix.hh" 43 #include "G4LogicalVolume.hh" 44 #include "G4PVPlacement.hh" 45 #include "G4SubtractionSolid.hh" 46 #include "G4HumanPhantomColour.hh" 47 #include "G4Box.hh" 48 49 G4VPhysicalVolume* G4MIRDRightScapula::Constru 50 const G4String& colourName, G4boo 51 { 52 G4cout<<"Construct "<<volumeName<<" with mot 53 54 auto* material = new G4HumanPhantomMaterial( 55 auto* skeleton = material -> GetMaterial("sk 56 57 G4double ax_in = 17.* cm; 58 G4double by_in = 9.8* cm; 59 G4double ax_out = 19.*cm; 60 G4double by_out = 9.8*cm; 61 G4double dz= 16.4* cm; 62 63 auto* inner_scapula = new G4EllipticalTube(" 64 auto* outer_scapula = new G4EllipticalTube(" 65 66 auto* subtraction = new G4Box("subtraction", 67 68 G4double xx = -ax_out * 0.242 ; //(sin 14deg 69 G4double yy = - ax_out * 0.97; // (cos 14 d 70 71 auto* rm = new G4RotationMatrix(); 72 rm -> rotateZ(14.* degree); 73 74 auto* scapula_first = new G4SubtractionSoli 75 outer_scapula, 76 subtraction, 77 rm, 78 G4ThreeVector(xx, yy, 0. * 79 80 G4double xx2 = ax_out * 0.62470 ; //(cos 51. 81 G4double yy2 = ax_out * 0.78087; // (sin 51. 82 83 auto* rm2 = new G4RotationMatrix(); 84 rm2 -> rotateZ(38.6598* degree); 85 86 auto* scapula_bone = new G4SubtractionSolid 87 scapula_first, 88 subtraction, 89 rm2, 90 G4ThreeVector(xx2, yy2, 0. 91 92 auto* scapula = new G4SubtractionSolid("Sca 93 scapula_bone, 94 inner_scapula); 95 96 auto* logicRightScapula = new G4LogicalVolum 97 skeleton, 98 "logical" + volumeName, 99 nullptr, nullptr, nullptr); 100 101 G4VPhysicalVolume* physRightScapula = new G4 102 G4ThreeVector(0. * cm, 0. * cm 103 "physicalRightScapula", 104 logicRightScapula, 105 mother, 106 false, 107 0, true); 108 109 // Visualization Attributes 110 //G4VisAttributes* RightScapulaVisAtt = new 111 auto* colourPointer = new G4HumanPhantomColo 112 G4Colour colour = colourPointer -> GetColour 113 auto* RightScapulaVisAtt = new G4VisAttribut 114 RightScapulaVisAtt->SetForceSolid(wireFrame) 115 logicRightScapula->SetVisAttributes(RightSca 116 G4cout << "RightScapula created !!!!!!" << G 117 118 // Testing RightScapula Volume 119 G4double RightScapulaVol = logicRightScapula 120 G4cout << "Volume of RightScapula = " << Rig 121 122 // Testing RightScapula Material 123 G4String RightScapulaMat = logicRightScapula 124 G4cout << "Material of RightScapula = " << R 125 126 // Testing Density 127 G4double RightScapulaDensity = logicRightSca 128 G4cout << "Density of Material = " << RightS 129 130 // Testing Mass 131 G4double RightScapulaMass = (RightScapulaVol 132 G4cout << "Mass of RightScapula = " << Right 133 134 return physRightScapula; 135 } 136