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 #include "G4MIRDRightArmBone.hh" 31 #include "globals.hh" 32 #include "G4SystemOfUnits.hh" 33 #include "G4SDManager.hh" 34 #include "G4VisAttributes.hh" 35 #include "G4EllipticalTube.hh" 36 #include "G4RotationMatrix.hh" 37 #include "G4ThreeVector.hh" 38 #include "G4VPhysicalVolume.hh" 39 #include "G4PVPlacement.hh" 40 #include "G4UnionSolid.hh" 41 #include "G4EllipticalCone.hh" 42 #include "G4HumanPhantomColour.hh" 43 #include "G4HumanPhantomMaterial.hh" 44 45 G4VPhysicalVolume* G4MIRDRightArmBone::Constru 46 const G4String& colourName,G4bool 47 { 48 auto* material = new G4HumanPhantomMaterial( 49 50 G4cout<<"Construct "<<volumeName<<" with mot 51 52 auto* skeleton = material -> GetMaterial("sk 53 54 delete material; 55 56 G4double dx = 1.4 * cm;//a 57 G4double dy = 2.7 * cm;//b 58 // G4double dz= 46. * cm;//z0 59 60 auto* rightArm = new G4EllipticalTube("OneRi 61 62 auto* logicRightArmBone = new G4LogicalVolum 63 skeleton, 64 "logical" + volumeName, 65 nullptr, nullptr, nullptr); 66 67 auto* rm = new G4RotationMatrix(); 68 rm->rotateX(180.*degree); 69 G4VPhysicalVolume* physRightArmBone = new G4 70 G4ThreeVector(-18.4 * cm, 0.0, 71 //-x0 72 "physicalRightArmBone", 73 logicRightArmBone, 74 mother, 75 false,0,true); 76 77 // Visualization Attributes 78 //G4VisAttributes* RightArmBoneVisAtt = new 79 auto* colourPointer = new G4HumanPhantomColo 80 G4Colour colour = colourPointer -> GetColour 81 auto* RightArmBoneVisAtt = new G4VisAttribut 82 83 RightArmBoneVisAtt->SetForceSolid(wireFrame) 84 logicRightArmBone->SetVisAttributes(RightArm 85 86 G4cout << "RightArmBone created !!!!!!" << G 87 88 // Testing RightArmBone Volume 89 G4double RightArmBoneVol = logicRightArmBone 90 G4cout << "Volume of RightArmBone = " << Rig 91 92 // Testing RightArmBone Material 93 G4String RightArmBoneMat = logicRightArmBone 94 G4cout << "Material of RightArmBone = " << R 95 96 // Testing Density 97 G4double RightArmBoneDensity = logicRightArm 98 G4cout << "Density of Material = " << RightA 99 100 // Testing Mass 101 G4double RightArmBoneMass = (RightArmBoneVol 102 G4cout << "Mass of RightArmBone = " << Right 103 104 return physRightArmBone; 105 } 106