Changeset 921 for trunk/source/geometry/solids/specific
- Timestamp:
- Feb 16, 2009, 10:14:30 AM (16 years ago)
- Location:
- trunk/source/geometry/solids/specific
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/geometry/solids/specific/History
r850 r921 1 $Id: History,v 1.1 45 2008/08/12 08:57:31gcosmo Exp $1 $Id: History,v 1.150 2008/11/21 09:26:53 gcosmo Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 17 17 * Reverse chronological order (last date on top), please * 18 18 ---------------------------------------------------------- 19 20 21-Nov-2008, G.Cosmo (geom-specific-V09-01-19) 21 - Added missing accessors to G4EllipticalCone, required for detector 22 persistency. 23 24 13-Nov-2008, G.Cosmo (geom-specific-V09-01-18) 25 - Corrected initialisation of algorithm in G4TriangularFacet constructor. 26 27 30-Oct-2008, I.Hrivnacova (geom-specific-V09-01-17) 28 - G4ExtrudedSolid: fixed bug in the decomposition of polygonal sides for 29 quadrangular facets in MakeFacets(). Addresses problem report #1029. 30 31 23-Sep-2008, T.Nikitina (geom-specific-V09-01-16) 32 - Corrected algorithm in G4TriangularFacet::GetPointOnFace() according to 33 suggestion advanced in problem report #1025. Fixes a problem of false 34 overlaps detection related to G4ExtrudedSolid and G4TessellatedSolid. 35 36 12-Sep-2008, G.Cosmo (geom-specific-V09-01-15) 37 - G4VFacet: corrected increment of indeces in operator==(), following 38 problem report #1024. 19 39 20 40 25-Jul-2008, I.Hrivnacova (geom-specific-V09-01-14) -
trunk/source/geometry/solids/specific/include/G4EllipticalCone.hh
r850 r921 25 25 // 26 26 // 27 // $Id: G4EllipticalCone.hh,v 1.1 1 2007/08/20 15:21:40 tnikitinExp $28 // GEANT4 tag $Name: HEAD$27 // $Id: G4EllipticalCone.hh,v 1.12 2008/11/21 09:26:22 gcosmo Exp $ 28 // GEANT4 tag $Name: geant4-09-02-cand-01 $ 29 29 // 30 30 // … … 96 96 // 97 97 inline G4double GetSemiAxisMax () const; 98 inline G4double GetSemiAxisX () const; 99 inline G4double GetSemiAxisY () const; 100 inline G4double GetZMax() const; 98 101 inline G4double GetZTopCut() const; 99 102 inline void SetSemiAxis (G4double x, G4double y, G4double z); -
trunk/source/geometry/solids/specific/include/G4EllipticalCone.icc
r850 r921 25 25 // 26 26 // 27 // $Id: G4EllipticalCone.icc,v 1. 6 2006/10/20 13:45:20gcosmo Exp $28 // GEANT4 tag $Name: HEAD$27 // $Id: G4EllipticalCone.icc,v 1.7 2008/11/21 09:26:22 gcosmo Exp $ 28 // GEANT4 tag $Name: geant4-09-02-cand-01 $ 29 29 // 30 30 // … … 43 43 { 44 44 return ySemiAxis > xSemiAxis ? ySemiAxis : xSemiAxis; 45 } 46 47 inline 48 G4double G4EllipticalCone::GetSemiAxisX () const 49 { 50 return xSemiAxis; 51 } 52 53 inline 54 G4double G4EllipticalCone::GetSemiAxisY () const 55 { 56 return ySemiAxis; 57 } 58 59 inline 60 G4double G4EllipticalCone::GetZMax() const 61 { 62 return zheight; 45 63 } 46 64 -
trunk/source/geometry/solids/specific/src/G4ExtrudedSolid.cc
r850 r921 25 25 // 26 26 // 27 // $Id: G4ExtrudedSolid.cc,v 1.1 7 2008/08/12 08:54:57 gcosmoExp $28 // GEANT4 tag $Name: HEAD$27 // $Id: G4ExtrudedSolid.cc,v 1.18 2008/10/30 11:47:45 ivana Exp $ 28 // GEANT4 tag $Name: geant4-09-02-cand-01 $ 29 29 // 30 30 // … … 562 562 563 563 good = AddFacet( new G4QuadrangularFacet( GetVertex(fNz-1, 3), GetVertex(fNz-1, 2), 564 GetVertex(fNz-1, 1), GetVertex( 1, 0),564 GetVertex(fNz-1, 1), GetVertex(fNz-1, 0), 565 565 ABSOLUTE) ); 566 566 if ( ! good ) { return false; } -
trunk/source/geometry/solids/specific/src/G4TriangularFacet.cc
r850 r921 25 25 // ******************************************************************** 26 26 // 27 // $Id: G4TriangularFacet.cc,v 1.1 0 2007/12/10 16:30:35 gunterExp $28 // GEANT4 tag $Name: HEAD$27 // $Id: G4TriangularFacet.cc,v 1.12 2008/11/13 08:25:07 gcosmo Exp $ 28 // GEANT4 tag $Name: geant4-09-02-cand-01 $ 29 29 // 30 30 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% … … 75 75 : G4VFacet() 76 76 { 77 if (!tGeomAlg) { tGeomAlg = G4TessellatedGeometryAlgorithms::GetInstance(); }77 tGeomAlg = G4TessellatedGeometryAlgorithms::GetInstance(); 78 78 P0 = Pt0; 79 79 nVertices = 3; … … 710 710 G4ThreeVector G4TriangularFacet::GetPointOnFace() const 711 711 { 712 G4double lambda0 = CLHEP::RandFlat::shoot(0.,1.); 713 G4double lambda1 = CLHEP::RandFlat::shoot(0.,lambda0); 714 712 G4double alpha = CLHEP::RandFlat::shoot(0.,1.); 713 G4double beta = CLHEP::RandFlat::shoot(0.,1); 714 G4double lambda1=alpha*beta; 715 G4double lambda0=alpha-lambda1; 716 715 717 return (P0 + lambda0*E[0] + lambda1*E[1]); 716 718 } -
trunk/source/geometry/solids/specific/src/G4VFacet.cc
r850 r921 25 25 // ******************************************************************** 26 26 // 27 // $Id: G4VFacet.cc,v 1. 6 2007/08/23 14:45:03gcosmo Exp $28 // GEANT4 tag $Name: HEAD$27 // $Id: G4VFacet.cc,v 1.7 2008/09/12 07:16:22 gcosmo Exp $ 28 // GEANT4 tag $Name: geant4-09-02-cand-01 $ 29 29 // 30 30 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% … … 96 96 { 97 97 coincident = (GetVertex(i)-right.GetVertex(j)).mag2() < tolerance; 98 } while (!coincident && j++< nVertices);99 } while (coincident && i++< nVertices);98 } while (!coincident && ++j < nVertices); 99 } while (coincident && ++i < nVertices); 100 100 101 101 return coincident; … … 117 117 std::ostream &G4VFacet::StreamInfo(std::ostream &os) const 118 118 { 119 os << G4endl;120 os << "***********************************************************************"121 << G4endl;122 os << "FACET TYPE = " <<geometryType <<G4endl;123 os << "ABSOLUTE VECTORS = " <<G4endl;124 os << "P0 = " <<P0 <<G4endl;119 os << G4endl; 120 os << "*********************************************************************" 121 << G4endl; 122 os << "FACET TYPE = " << geometryType << G4endl; 123 os << "ABSOLUTE VECTORS = " << G4endl; 124 os << "P0 = " << P0 << G4endl; 125 125 for (G4ThreeVectorList::const_iterator it=P.begin(); it!=P.end(); it++) 126 os <<"P[" <<it-P.begin()+1 <<"] = " <<*it <<G4endl;127 128 os << "RELATIVE VECTORS = " <<G4endl;126 { os << "P[" << it-P.begin()+1 << "] = " << *it << G4endl; } 127 128 os << "RELATIVE VECTORS = " << G4endl; 129 129 for (G4ThreeVectorList::const_iterator it=E.begin(); it!=E.end(); it++) 130 os <<"E[" <<it-E.begin()+1 <<"] = " <<*it <<G4endl;131 132 os << "***********************************************************************"133 << G4endl;130 { os << "E[" << it-E.begin()+1 << "] = " << *it << G4endl; } 131 132 os << "*********************************************************************" 133 << G4endl; 134 134 135 135 return os;
Note: See TracChangeset
for help on using the changeset viewer.