human_phantom/src/G4PhantomHeadBuilder.cc

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 // Previous authors: G. Guerrieri, S. Guatelli and M. G. Pia, INFN Genova, Italy
0027 // Authors (since 2007): S. Guatelli, University of Wollongong, Australia
0028 //
0029 #include "G4PhantomHeadBuilder.hh"
0030 #include "G4VBodyFactory.hh"
0031 #include "G4MIRDBodyFactory.hh"
0032 #include "G4ORNLFemaleBodyFactory.hh"
0033 #include "G4ORNLMaleBodyFactory.hh"
0034 #include "G4PVPlacement.hh"
0035
0036 G4PhantomHeadBuilder::G4PhantomHeadBuilder(): fModel("MIRD")
0037 {
0038   // sex can be "female" or "male"
0039   fBody = nullptr;
0040   fMotherVolume = nullptr;
0041   fHeadVolume = nullptr;
0042 }
0043
0044 G4PhantomHeadBuilder::~G4PhantomHeadBuilder()
0045 {
0046   delete fBody;
0047 }
0048
0049 void G4PhantomHeadBuilder::BuildHead(const G4String& colourName, G4bool solidVis, G4bool sensitivity)
0050 {
0051   if (fMotherVolume == nullptr)
0052     G4Exception("G4PhantomHeadBuilder::BuildHead()", "human_phantom0011", FatalException, "The moder Volume volume is missing !!!!!");
0053   else
0054   {
0055   G4cout <<"MotherVolume: " <<  fMotherVolume -> GetName()<< G4endl;
0056   G4cout << "sensitivity : "<< sensitivity << G4endl;
0057   fHeadVolume = fBody -> CreateOrgan("Head", fMotherVolume, colourName, solidVis, sensitivity);
0058 }
0059 }
0060
0061 void G4PhantomHeadBuilder::BuildSkull(const G4String& colourName, G4bool solidVis, G4bool sensitivity)
0062 {
0063   if (fHeadVolume == nullptr)
0064     G4Exception("G4PhantomHeadBuilder::BuildSkull()", "human_phantom0012", FatalException, "The head volume is missing !!!!!");
0065
0066   G4cout <<"MotherVolume: " << fHeadVolume -> GetName()<< G4endl;
0067   G4cout << "sensitivity : "<< sensitivity << G4endl;
0068   fBody -> CreateOrgan( "Skull", fHeadVolume, colourName, solidVis, sensitivity);
0069 }
0070
0071 void G4PhantomHeadBuilder::BuildBrain(const G4String& colourName, G4bool solidVis, G4bool sensitivity)
0072 {
0073  if (fHeadVolume == nullptr)
0074    G4Exception("G4PhantomHeadBuilder::BuildBrain()", "human_phantom0013", FatalException, "The head volume is missing !!!!!");
0075     else fBody -> CreateOrgan("Brain", fHeadVolume, colourName, solidVis, sensitivity);
0076 }
0077
0078 G4VPhysicalVolume* G4PhantomHeadBuilder::GetPhantom()
0079 {
0080   return fMotherVolume;
0081 }
0082
0083 void G4PhantomHeadBuilder::SetMotherVolume(G4VPhysicalVolume* mother)
0084 {
0085   fMotherVolume = mother;
0086 }
0087
0088
0089 void G4PhantomHeadBuilder::SetModel(G4String modelFlag)
0090 {
0091   fModel = modelFlag;
0092
0093   if(fModel=="MIRD") fBody = new G4MIRDBodyFactory();
0094   else if(fModel=="ORNLFemale")
0095 {
0096 #ifdef G4LIB_USE_GDML
0097 fBody = new G4ORNLFemaleBodyFactory();
0098 #else
0099 G4cout << fModel << " Working with GDML only! set G4LIB_USE_GDML 1" << G4endl;
0100 #endif
0101 }
0102 else if(fModel=="ORNLMale")
0103 {
0104 #ifdef G4LIB_USE_GDML
0105 fBody = new G4ORNLMaleBodyFactory();
0106 #else
0107 G4cout << fModel << " Working with GDML only! set G4LIB_USE_GDML 1" << G4endl;
0108 #endif
0109 }
0110
0111 }
0112