G4MIRDHead Class Reference

#include <G4MIRDHead.hh>

Inheritance diagram for G4MIRDHead:

Public Member Functions
	G4MIRDHead ()
	~G4MIRDHead ()
G4VPhysicalVolume *	Construct (const G4String &, G4VPhysicalVolume *, const G4String &, G4bool, G4bool)
Public Member Functions inherited from G4VOrgan
	G4VOrgan ()
virtual	~G4VOrgan ()

Detailed Description

Definition at line 44 of file G4MIRDHead.hh.

Constructor & Destructor Documentation

G4MIRDHead::G4MIRDHead

(

)

Definition at line 50 of file G4MIRDHead.cc.

{

}

G4MIRDHead::~G4MIRDHead

(

)

Definition at line 55 of file G4MIRDHead.cc.

{

}

Member Function Documentation

G4VPhysicalVolume * G4MIRDHead::Construct ( const G4String &  volumeName, G4VPhysicalVolume *  mother, const G4String &  colourName, G4bool  wireFrame, G4bool    )

virtual

Implements G4VOrgan.

Definition at line 61 of file G4MIRDHead.cc.

References python.hepunit::cm, python.hepunit::cm3, python.hepunit::degree, g(), G4cout, G4endl, G4VSolid::GetCubicVolume(), G4Material::GetDensity(), G4LogicalVolume::GetMaterial(), G4VPhysicalVolume::GetName(), G4Material::GetName(), G4LogicalVolume::GetSolid(), python.hepunit::gram, eplot::material, CLHEP::HepRotation::rotateX(), CLHEP::HepRotation::rotateY(), G4VisAttributes::SetForceSolid(), and G4LogicalVolume::SetVisAttributes().

{
  G4cout << "Construct " << volumeName <<" with mother "<<mother->GetName()<<G4endl;
 
  G4HumanPhantomMaterial * material = new G4HumanPhantomMaterial();
  
  G4Material* soft = material -> GetMaterial("soft_tissue");
  
  // MIRD male model
  // Ellipsoid
  G4double ax = 7.0 * cm;
  G4double by = 10.0 * cm;
  G4double cz = 8.50 * cm;
  G4double zcut1 = 0.0 * cm;
  G4double zcut2 = 8.5 * cm;

  G4Ellipsoid* head1 = new G4Ellipsoid("Head1", ax, by, cz, zcut1, zcut2);

  G4double dx = 7.0 * cm;
  G4double dy = 10.0 * cm;
  G4double dz = 7.75 * cm;
 

  G4EllipticalTube* head2 = new G4EllipticalTube("Head2", dx, dy, dz);

  G4UnionSolid* head = new G4UnionSolid("Head",head2,head1,
                    0, // Rotation 
                    G4ThreeVector(0.* cm, 0.*cm, 7.7500 * cm) );

  G4LogicalVolume* logicHead = new G4LogicalVolume(head, soft,"logical" + volumeName,
                           0, 0,0);
  // Define rotation and position here!
  G4RotationMatrix* rm = new G4RotationMatrix();
  rm->rotateX(180.*degree); 
  rm->rotateY(180.*degree); 
  
  G4VPhysicalVolume* physHead = new G4PVPlacement(rm,
                          //G4ThreeVector(0.* cm,0.*cm, 70.75*cm), //FA
                          G4ThreeVector(0.* cm,0.*cm, 77.75*cm),
                          "physicalHead",
                          logicHead,
                          mother,
                          false,
                          0, true);

 
  // Visualization Attributes

  G4HumanPhantomColour* colourPointer = new G4HumanPhantomColour();
  G4Colour colour = colourPointer -> GetColour(colourName);
  G4VisAttributes* HeadVisAtt = new G4VisAttributes(colour);

  HeadVisAtt->SetForceSolid(wireFrame);
  // HeadVisAtt->SetLineWidth(0.7* mm);
  //HeadVisAtt-> SetForceAuxEdgeVisible(true);
  logicHead->SetVisAttributes(HeadVisAtt);

  // Testing Head Volume
  G4double HeadVol = logicHead->GetSolid()->GetCubicVolume();
  G4cout << "Volume of Head = " << HeadVol/cm3 << " cm^3" << G4endl;
  
  // Testing Head Material
  G4String HeadMat = logicHead->GetMaterial()->GetName();
  G4cout << "Material of Head = " << HeadMat << G4endl;
  
  // Testing Density
  G4double HeadDensity = logicHead->GetMaterial()->GetDensity();
  G4cout << "Density of Material = " << HeadDensity*cm3/g << " g/cm^3" << G4endl;

  // Testing Mass
  G4double HeadMass = (HeadVol)*HeadDensity;
  G4cout << "Mass of Head = " << HeadMass/gram << " g" << G4endl;
  
  return physHead;
}

---

The documentation for this class was generated from the following files:

• G4MIRDHead.hh
• G4MIRDHead.cc