| 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 // G4ErrorPlaneSurfaceTarget class implementat 26 // G4ErrorPlaneSurfaceTarget class implementation
27 // 27 //
28 // Created: P.Arce, September 2004 28 // Created: P.Arce, September 2004
29 // ------------------------------------------- 29 // --------------------------------------------------------------------
30 30
31 #include "G4ErrorPlaneSurfaceTarget.hh" 31 #include "G4ErrorPlaneSurfaceTarget.hh"
32 32
33 #ifdef G4VERBOSE 33 #ifdef G4VERBOSE
34 #include "G4ErrorPropagatorData.hh" // for ver 34 #include "G4ErrorPropagatorData.hh" // for verbosity checking
35 #endif 35 #endif
36 36
37 #include "G4Point3D.hh" 37 #include "G4Point3D.hh"
38 #include "G4ThreeVector.hh" 38 #include "G4ThreeVector.hh"
39 39
40 //-------------------------------------------- 40 //---------------------------------------------------------------------
41 41
42 G4ErrorPlaneSurfaceTarget:: 42 G4ErrorPlaneSurfaceTarget::
43 G4ErrorPlaneSurfaceTarget(G4double aa, G4doubl 43 G4ErrorPlaneSurfaceTarget(G4double aa, G4double ab, G4double ac, G4double ad)
44 : G4Plane3D( aa, ab, ac, ad ) 44 : G4Plane3D( aa, ab, ac, ad )
45 { 45 {
46 theType = G4ErrorTarget_PlaneSurface; 46 theType = G4ErrorTarget_PlaneSurface;
47 47
48 #ifdef G4VERBOSE 48 #ifdef G4VERBOSE
49 if(G4ErrorPropagatorData::verbose() >= 2 ) 49 if(G4ErrorPropagatorData::verbose() >= 2 )
50 { 50 {
51 Dump( " $$$ creating G4ErrorPlaneSurfaceTa 51 Dump( " $$$ creating G4ErrorPlaneSurfaceTarget from parameters");
52 } 52 }
53 #endif 53 #endif
54 } 54 }
55 55
56 //-------------------------------------------- 56 //---------------------------------------------------------------------
57 57
58 G4ErrorPlaneSurfaceTarget:: 58 G4ErrorPlaneSurfaceTarget::
59 G4ErrorPlaneSurfaceTarget(const G4Normal3D &no 59 G4ErrorPlaneSurfaceTarget(const G4Normal3D &norm, const G4Point3D &pt)
60 : G4Plane3D( norm, pt ) 60 : G4Plane3D( norm, pt )
61 { 61 {
62 theType = G4ErrorTarget_PlaneSurface; 62 theType = G4ErrorTarget_PlaneSurface;
63 63
64 #ifdef G4VERBOSE 64 #ifdef G4VERBOSE
65 if(G4ErrorPropagatorData::verbose() >= 2 ) 65 if(G4ErrorPropagatorData::verbose() >= 2 )
66 { 66 {
67 Dump( " $$$ creating G4ErrorPlaneSurfaceTa 67 Dump( " $$$ creating G4ErrorPlaneSurfaceTarget from point and normal");
68 } 68 }
69 #endif 69 #endif
70 } 70 }
71 71
72 //-------------------------------------------- 72 //---------------------------------------------------------------------
73 73
74 G4ErrorPlaneSurfaceTarget:: 74 G4ErrorPlaneSurfaceTarget::
75 G4ErrorPlaneSurfaceTarget(const G4Point3D &p1, 75 G4ErrorPlaneSurfaceTarget(const G4Point3D &p1,
76 const G4Point3D &p2, 76 const G4Point3D &p2,
77 const G4Point3D &p3) 77 const G4Point3D &p3)
78 : G4Plane3D( p1, p2, p3 ) 78 : G4Plane3D( p1, p2, p3 )
79 { 79 {
80 theType = G4ErrorTarget_PlaneSurface; 80 theType = G4ErrorTarget_PlaneSurface;
81 81
82 #ifdef G4VERBOSE 82 #ifdef G4VERBOSE
83 if(G4ErrorPropagatorData::verbose() >= 2 ) 83 if(G4ErrorPropagatorData::verbose() >= 2 )
84 { 84 {
85 Dump( " $$$ creating G4ErrorPlaneSurfaceTa 85 Dump( " $$$ creating G4ErrorPlaneSurfaceTarget from three points");
86 } 86 }
87 #endif 87 #endif
88 } 88 }
89 89
90 //-------------------------------------------- 90 //---------------------------------------------------------------------
91 91
92 G4ThreeVector G4ErrorPlaneSurfaceTarget:: 92 G4ThreeVector G4ErrorPlaneSurfaceTarget::
93 Intersect( const G4ThreeVector& pt, const G4Th 93 Intersect( const G4ThreeVector& pt, const G4ThreeVector& dir ) const
94 { 94 {
95 G4double lam = GetDistanceFromPoint( pt, dir 95 G4double lam = GetDistanceFromPoint( pt, dir );
96 G4Point3D inters = pt + lam * dir; 96 G4Point3D inters = pt + lam * dir;
97 97
98 #ifdef G4VERBOSE 98 #ifdef G4VERBOSE
99 if(G4ErrorPropagatorData::verbose() >= 4 ) 99 if(G4ErrorPropagatorData::verbose() >= 4 )
100 { 100 {
101 G4cout << " $$$ creating G4ErrorPlaneSurfa 101 G4cout << " $$$ creating G4ErrorPlaneSurfaceTarget::Intersect "
102 << inters << G4endl; 102 << inters << G4endl;
103 } 103 }
104 #endif 104 #endif
105 105
106 return inters; 106 return inters;
107 } 107 }
108 108
109 //-------------------------------------------- 109 //---------------------------------------------------------------------
110 110
111 G4double G4ErrorPlaneSurfaceTarget:: 111 G4double G4ErrorPlaneSurfaceTarget::
112 GetDistanceFromPoint( const G4ThreeVector& pt, 112 GetDistanceFromPoint( const G4ThreeVector& pt, const G4ThreeVector& dir ) const
113 { 113 {
114 if( std::fabs( dir.mag() -1. ) > 1.E-6 ) 114 if( std::fabs( dir.mag() -1. ) > 1.E-6 )
115 { 115 {
116 std::ostringstream message; 116 std::ostringstream message;
117 message << "Direction is not a unit vector 117 message << "Direction is not a unit vector: " << dir << " !";
118 G4Exception("G4ErrorPlaneSurfaceTarget::Ge 118 G4Exception("G4ErrorPlaneSurfaceTarget::GetDistanceFromPoint()",
119 "GeomMgt1002", JustWarning, me 119 "GeomMgt1002", JustWarning, message);
120 } 120 }
121 G4double dist = -(a_ * pt.x() + b_ * pt.y() 121 G4double dist = -(a_ * pt.x() + b_ * pt.y() + c_ * pt.z() + d_)
122 / (a_ * dir.x() + b_ * dir.y( 122 / (a_ * dir.x() + b_ * dir.y() + c_ * dir.z() );
123 123
124 #ifdef G4VERBOSE 124 #ifdef G4VERBOSE
125 if(G4ErrorPropagatorData::verbose() >= 3 ) 125 if(G4ErrorPropagatorData::verbose() >= 3 )
126 { 126 {
127 G4cout << " G4ErrorPlaneSurfaceTarget::Get 127 G4cout << " G4ErrorPlaneSurfaceTarget::GetDistanceFromPoint()" << G4endl
128 << " Point: " << pt << ", Directi 128 << " Point: " << pt << ", Direction: " << dir << G4endl
129 << " Distance: " << dist << G4end 129 << " Distance: " << dist << G4endl;
130 } 130 }
131 #endif 131 #endif
132 132
133 return dist; 133 return dist;
134 } 134 }
135 135
136 //-------------------------------------------- 136 //---------------------------------------------------------------------
137 137
138 G4double G4ErrorPlaneSurfaceTarget:: 138 G4double G4ErrorPlaneSurfaceTarget::
139 GetDistanceFromPoint( const G4ThreeVector& pt 139 GetDistanceFromPoint( const G4ThreeVector& pt ) const
140 { 140 {
141 G4ThreeVector vec = point() - pt; 141 G4ThreeVector vec = point() - pt;
142 G4double dist = std::fabs(vec * normal() / n 142 G4double dist = std::fabs(vec * normal() / normal().mag());
143 143
144 #ifdef G4VERBOSE 144 #ifdef G4VERBOSE
145 if(G4ErrorPropagatorData::verbose() >= 3 ) 145 if(G4ErrorPropagatorData::verbose() >= 3 )
146 { 146 {
147 G4cout << " G4ErrorPlaneSurfaceTarget::Get 147 G4cout << " G4ErrorPlaneSurfaceTarget::GetDistanceFromPoint()" << G4endl
148 << " Point: " << pt << G4endl 148 << " Point: " << pt << G4endl
149 << " Distance: " << dist << G4end 149 << " Distance: " << dist << G4endl;
150 } 150 }
151 #endif 151 #endif
152 152
153 return dist; 153 return dist;
154 } 154 }
155 155
156 //-------------------------------------------- 156 //---------------------------------------------------------------------
157 157
158 G4Plane3D G4ErrorPlaneSurfaceTarget:: 158 G4Plane3D G4ErrorPlaneSurfaceTarget::
159 GetTangentPlane( const G4ThreeVector& ) const 159 GetTangentPlane( const G4ThreeVector& ) const
160 { 160 {
161 return *this; 161 return *this;
162 } 162 }
163 163
164 //-------------------------------------------- 164 //---------------------------------------------------------------------
165 165
166 void G4ErrorPlaneSurfaceTarget::Dump( const G4 166 void G4ErrorPlaneSurfaceTarget::Dump( const G4String& msg ) const
167 { 167 {
168 G4cout << msg << " point = " << point() 168 G4cout << msg << " point = " << point()
169 << " normal = " << normal() << 169 << " normal = " << normal() << G4endl;
170 } 170 }
171 171