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 /// \file DicomPhantomParameterisationColour.c 28 /// \brief Implementation of the DicomPhantomP 29 30 #include "DicomPhantomParameterisationColour.h 31 32 #include "DicomHandler.hh" 33 34 #include "G4LogicalVolume.hh" 35 #include "G4Material.hh" 36 #include "G4VPhysicalVolume.hh" 37 #include "G4VVisManager.hh" 38 #include "G4VisAttributes.hh" 39 #include "globals.hh" 40 41 //....oooOO0OOooo........oooOO0OOooo........oo 42 43 G4String DicomPhantomParameterisationColour::f 44 DicomHandler::GetDicomDataPath() + "/ColourM 45 46 //....oooOO0OOooo........oooOO0OOooo........oo 47 DicomPhantomParameterisationColour::DicomPhant 48 : G4PhantomParameterisation() 49 { 50 ReadColourData(colourFile); 51 SetSkipEqualMaterials(false); 52 } 53 54 //....oooOO0OOooo........oooOO0OOooo........oo 55 DicomPhantomParameterisationColour::~DicomPhan 56 57 //....oooOO0OOooo........oooOO0OOooo........oo 58 void DicomPhantomParameterisationColour::ReadC 59 { 60 //----- Add a G4VisAttributes for materials 61 G4VisAttributes* blankAtt = new G4VisAttribu 62 blankAtt->SetVisibility(FALSE); 63 fColours["Default"] = blankAtt; 64 65 //----- Read file 66 std::ifstream fin(colourFile.c_str()); 67 G4int nMate; 68 G4String mateName; 69 G4double cred, cgreen, cblue, copacity; 70 fin >> nMate; 71 for (G4int ii = 0; ii < nMate; ii++) { 72 fin >> mateName; 73 if (fin.eof()) break; 74 fin >> cred >> cgreen >> cblue >> copacity 75 G4Colour colour(cred, cgreen, cblue, copac 76 G4VisAttributes* visAtt = new G4VisAttribu 77 visAtt->SetVisibility(true); 78 fColours[mateName] = visAtt; 79 fColours2[ii] = new G4VisAttributes(*visAt 80 } 81 } 82 83 //....oooOO0OOooo........oooOO0OOooo........oo 84 G4Material* DicomPhantomParameterisationColour 85 86 87 { 88 G4Material* mate = G4PhantomParameterisation 89 90 if (G4VVisManager::GetConcreteInstance() && 91 G4String mateName = mate->GetName(); 92 std::string::size_type iuu = mateName.find 93 if (iuu != std::string::npos) mateName = m 94 95 if (0 < fColours.count(mateName)) 96 physVol->GetLogicalVolume()->SetVisAttri 97 else { 98 bool found = false; 99 for (const auto& itr : fColours) { 100 G4String mat_color = itr.first; 101 auto len = mat_color.length(); 102 if (mateName.find(mat_color) == 0 && m 103 physVol->GetLogicalVolume()->SetVisA 104 found = true; 105 } 106 if (found) break; 107 } 108 if (!found) { 109 G4int matIndex = G4int(GetMaterialInde 110 static uintmax_t n = 0; 111 if (n++ < 100) 112 G4cout << "Unknown material name " < 113 if (fColours2.find(matIndex) != fColou 114 physVol->GetLogicalVolume()->SetVisA 115 else 116 physVol->GetLogicalVolume()->SetVisA 117 } 118 } 119 physVol->GetLogicalVolume()->SetMaterial(m 120 } 121 122 return mate; 123 } 124