| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // Created by Stewart Boogert on 05/03/2023. 26 // Created by Stewart Boogert on 05/03/2023.
27 // 27 //
28 28
29 #include "G4VVtkPipeline.hh" // Move basic has 29 #include "G4VVtkPipeline.hh" // Move basic hashing functionality to G4VtkUtility.(hh/cc
30 #include "G4VtkUtility.hh" 30 #include "G4VtkUtility.hh"
31 31
32 #include "G4Transform3D.hh" 32 #include "G4Transform3D.hh"
33 #include "G4Colour.hh" 33 #include "G4Colour.hh"
34 34
35 #include <vtkMatrix4x4.h> << 35 #include "vtkMatrix4x4.h"
36 #include <vtkPlane.h> << 36 #include "vtkPlane.h"
37 37
38 vtkSmartPointer<vtkMatrix4x4> G4Transform3DToV 38 vtkSmartPointer<vtkMatrix4x4> G4Transform3DToVtkMatrix4x4(const G4Transform3D& g4Trans)
39 { 39 {
40 auto transform = vtkSmartPointer<vtkMatrix4x 40 auto transform = vtkSmartPointer<vtkMatrix4x4>();
41 41
42 double transformArray[16] = {g4Trans.xx(), g 42 double transformArray[16] = {g4Trans.xx(), g4Trans.xy(), g4Trans.xy(), g4Trans.dx(),
43 g4Trans.yx(), g 43 g4Trans.yx(), g4Trans.yy(), g4Trans.yy(), g4Trans.dy(),
44 g4Trans.zx(), g 44 g4Trans.zx(), g4Trans.zy(), g4Trans.zy(), g4Trans.dz(),
45 0., 0 45 0., 0., 0., 1.};
46 transform->DeepCopy(transformArray); 46 transform->DeepCopy(transformArray);
47 return transform; 47 return transform;
48 } 48 }
49 49
50 vtkSmartPointer<vtkPlane> G4Plane3DToVtkPlane( 50 vtkSmartPointer<vtkPlane> G4Plane3DToVtkPlane(const G4Plane3D& g4plane)
51 { 51 {
52 auto plane = vtkSmartPointer<vtkPlane>::New( 52 auto plane = vtkSmartPointer<vtkPlane>::New();
53 53
54 auto g4normal = g4plane.normal(); 54 auto g4normal = g4plane.normal();
55 auto g4point = g4plane.point(); 55 auto g4point = g4plane.point();
56 56
57 plane->SetNormal(g4normal.x(), g4normal.y(), 57 plane->SetNormal(g4normal.x(), g4normal.y(), g4normal.z());
58 plane->SetOrigin(g4point.x(), g4point.y(), g 58 plane->SetOrigin(g4point.x(), g4point.y(), g4point.z());
59 59
60 return plane; 60 return plane;
61 } 61 }
62 62
63 G4Plane3D VtkPlaneToG4Plane3D(vtkPlane* vtkPla 63 G4Plane3D VtkPlaneToG4Plane3D(vtkPlane* vtkPlane)
64 { 64 {
65 auto point = vtkPlane->GetOrigin(); 65 auto point = vtkPlane->GetOrigin();
66 auto normal = vtkPlane->GetNormal(); 66 auto normal = vtkPlane->GetNormal();
67 67
68 G4Plane3D plane = 68 G4Plane3D plane =
69 G4Plane3D(G4Normal3D(normal[0], normal[1], 69 G4Plane3D(G4Normal3D(normal[0], normal[1], normal[2]), G4Point3D(point[0], point[1], point[2]));
70 return plane; 70 return plane;
71 } 71 }
72 72
73 void MaxBounds(G4double* maxBound, G4double* b 73 void MaxBounds(G4double* maxBound, G4double* boundToCheck)
74 { 74 {
75 if (boundToCheck[0] < maxBound[0]) maxBound[ 75 if (boundToCheck[0] < maxBound[0]) maxBound[0] = boundToCheck[0];
76 if (boundToCheck[1] > maxBound[1]) maxBound[ 76 if (boundToCheck[1] > maxBound[1]) maxBound[1] = boundToCheck[1];
77 77
78 if (boundToCheck[2] < maxBound[2]) maxBound[ 78 if (boundToCheck[2] < maxBound[2]) maxBound[2] = boundToCheck[2];
79 if (boundToCheck[3] > maxBound[3]) maxBound[ 79 if (boundToCheck[3] > maxBound[3]) maxBound[3] = boundToCheck[3];
80 80
81 if (boundToCheck[4] < maxBound[4]) maxBound[ 81 if (boundToCheck[4] < maxBound[4]) maxBound[4] = boundToCheck[4];
82 if (boundToCheck[5] > maxBound[5]) maxBound[ 82 if (boundToCheck[5] > maxBound[5]) maxBound[5] = boundToCheck[5];
83 } 83 }
84 84
85 85
86 std::size_t MakeHash(const G4ThreeVector &v) 86 std::size_t MakeHash(const G4ThreeVector &v)
87 { 87 {
88 88
89 std::size_t xhash = std::hash<G4double>{}(v. 89 std::size_t xhash = std::hash<G4double>{}(v.x());
90 std::size_t yhash = std::hash<G4double>{}(v. 90 std::size_t yhash = std::hash<G4double>{}(v.y());
91 std::size_t zhash = std::hash<G4double>{}(v. 91 std::size_t zhash = std::hash<G4double>{}(v.z());
92 92
93 std::hash_combine(xhash, yhash); 93 std::hash_combine(xhash, yhash);
94 std::hash_combine(yhash, zhash); 94 std::hash_combine(yhash, zhash);
95 95
96 return xhash; 96 return xhash;
97 } 97 }
98 98
99 std::size_t MakeHash(const G4Colour &c) 99 std::size_t MakeHash(const G4Colour &c)
100 { 100 {
101 101
102 std::size_t rhash = std::hash<G4double>{}(c. 102 std::size_t rhash = std::hash<G4double>{}(c.GetRed());
103 std::size_t ghash = std::hash<G4double>{}(c. 103 std::size_t ghash = std::hash<G4double>{}(c.GetGreen());
104 std::size_t bhash = std::hash<G4double>{}(c. 104 std::size_t bhash = std::hash<G4double>{}(c.GetBlue());
105 std::size_t ahash = std::hash<G4double>{}(c. 105 std::size_t ahash = std::hash<G4double>{}(c.GetAlpha());
106 106
107 std::hash_combine(rhash, ghash); 107 std::hash_combine(rhash, ghash);
108 std::hash_combine(rhash, bhash); 108 std::hash_combine(rhash, bhash);
109 std::hash_combine(rhash, ahash); 109 std::hash_combine(rhash, ahash);
110 110
111 return rhash; 111 return rhash;
112 } 112 }