| 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 and of QinetiQ Ltd, *
20 // * subject to DEFCON 705 IPR conditions. << 20 // * subject DEFCON 705 IPR conditions. *
21 // * By using, copying, modifying or distri 21 // * By using, copying, modifying or distributing the software (or *
22 // * any work based on the software) you ag 22 // * any work based on the software) you agree to acknowledge its *
23 // * use in resulting scientific publicati 23 // * use in resulting scientific publications, and indicate your *
24 // * acceptance of all terms of the Geant4 Sof 24 // * acceptance of all terms of the Geant4 Software license. *
25 // ******************************************* 25 // ********************************************************************
26 // 26 //
27 // G4VFacet implementation. << 27 // $Id: G4VFacet.cc,v 1.3 2006/06/29 18:49:31 gunter Exp $
>> 28 // GEANT4 tag $Name: geant4-08-01-patch-01 $
28 // 29 //
29 // 31 October 2004, P R Truscott, QinetiQ Ltd, << 30 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
30 // 12 October 2012, M Gayer, CERN, - Reviewed <<
31 // ------------------------------------------- <<
32 <<
33 #include "G4VFacet.hh" <<
34 #include "globals.hh" <<
35 #include "geomdefs.hh" <<
36 #include "G4GeometryTolerance.hh" <<
37 <<
38 using namespace std; <<
39 <<
40 const G4double G4VFacet::dirTolerance = 1.0E-1 <<
41 <<
42 ////////////////////////////////////////////// <<
43 // 31 //
44 G4VFacet::G4VFacet() << 32 // MODULE: G4VFacet.hh
45 { <<
46 kCarTolerance = G4GeometryTolerance::GetInst <<
47 } <<
48 <<
49 ////////////////////////////////////////////// <<
50 // 33 //
51 G4VFacet::~G4VFacet() = default; << 34 // Date: 15/06/2005
52 << 35 // Author: P R Truscott
53 ////////////////////////////////////////////// << 36 // Organisation: QinetiQ Ltd, UK
>> 37 // Customer: UK Ministry of Defence : RAO CRP TD Electronic Systems
>> 38 // Contract: C/MAT/N03517
54 // 39 //
55 G4bool G4VFacet::operator== (const G4VFacet& r << 40 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
56 { <<
57 G4double tolerance = kCarTolerance*kCarToler <<
58 <<
59 if (GetNumberOfVertices() != right.GetNumber <<
60 return false; <<
61 else if ((GetCircumcentre()-right.GetCircumc <<
62 return false; <<
63 else if (std::fabs((right.GetSurfaceNormal() <<
64 return false; <<
65 <<
66 G4bool coincident = true; <<
67 G4int i = 0; <<
68 do // Loop checking, 13.08.2015, G.Cosmo <<
69 { <<
70 coincident = false; <<
71 G4int j = 0; <<
72 do // Loop checking, 13.08.2015, G.Cosm <<
73 { <<
74 coincident = (GetVertex(i)-right.GetVert <<
75 } while (!coincident && ++j < GetNumberOfV <<
76 } while (coincident && ++i < GetNumberOfVert <<
77 <<
78 return coincident; <<
79 } <<
80 <<
81 ////////////////////////////////////////////// <<
82 // 41 //
83 void G4VFacet::ApplyTranslation(const G4ThreeV << 42 // CHANGE HISTORY
84 { << 43 // --------------
85 G4int n = GetNumberOfVertices(); << 44 //
86 for (G4int i = 0; i < n; ++i) << 45 // 31 October 2004, P R Truscott, QinetiQ Ltd, UK - Created.
87 { << 46 //
88 SetVertex(i, GetVertex(i) + v); << 47 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
89 } << 48
90 } << 49 #include "G4VFacet.hh"
>> 50 #include "globals.hh"
91 51
92 ////////////////////////////////////////////// 52 ///////////////////////////////////////////////////////////////////////////////
93 // 53 //
94 std::ostream& G4VFacet::StreamInfo(std::ostrea << 54 std::ostream &G4VFacet::StreamInfo(std::ostream &os) const
95 { 55 {
96 os << G4endl; << 56 os <<G4endl;
97 os << "************************************* << 57 os <<"***********************************************************************"
98 << G4endl; << 58 <<G4endl;
99 os << "FACET TYPE = " << GetEntityType << 59 os <<"FACET TYPE = " <<geometryType <<G4endl;
100 os << "ABSOLUTE VECTORS = " << G4endl; << 60 os <<"ABSOLUTE VECTORS = " <<G4endl;
101 G4int n = GetNumberOfVertices(); << 61 os <<"P0 = " <<P0 <<G4endl;
102 for (G4int i = 0; i < n; ++i) << 62 for (G4ThreeVectorList::const_iterator it=P.begin(); it!=P.end(); it++)
103 os << "P[" << i << "] = " << GetVerte << 63 os <<"P[" <<it-P.begin()+1 <<"] = " <<*it <<G4endl;
104 os << "************************************* << 64
105 << G4endl; << 65 os <<"RELATIVE VECTORS = " <<G4endl;
106 << 66 for (G4ThreeVectorList::const_iterator it=E.begin(); it!=E.end(); it++)
>> 67 os <<"E[" <<it-E.begin()+1 <<"] = " <<*it <<G4endl;
>> 68
>> 69 os <<"***********************************************************************"
>> 70 <<G4endl;
>> 71
107 return os; 72 return os;
108 } <<
109 <<
110 G4bool G4VFacet::IsInside (const G4ThreeVector <<
111 { <<
112 G4ThreeVector d = p-GetVertex(0); <<
113 G4double displacement = d.dot(GetSurfaceNorm <<
114 return displacement <= 0.0; <<
115 } 73 }
116 74