Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/persistency/gdml/src/G4GDMLParameterisation.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /persistency/gdml/src/G4GDMLParameterisation.cc (Version 11.3.0) and /persistency/gdml/src/G4GDMLParameterisation.cc (Version 11.0.p2)


  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 // G4GDMLParameterisation implementation           26 // G4GDMLParameterisation implementation
 27 //                                                 27 //
 28 // Author: Zoltan Torzsok, November 2007           28 // Author: Zoltan Torzsok, November 2007
 29 // -------------------------------------------     29 // --------------------------------------------------------------------
 30                                                    30 
 31 #include "G4GDMLParameterisation.hh"               31 #include "G4GDMLParameterisation.hh"
 32 #include "G4PolyconeHistorical.hh"                 32 #include "G4PolyconeHistorical.hh"
 33 #include "G4PolyhedraHistorical.hh"                33 #include "G4PolyhedraHistorical.hh"
 34                                                    34 
 35 // -------------------------------------------     35 // --------------------------------------------------------------------
 36 G4int G4GDMLParameterisation::GetSize() const      36 G4int G4GDMLParameterisation::GetSize() const
 37 {                                                  37 {
 38   return (G4int) parameterList.size();             38   return (G4int) parameterList.size();
 39 }                                                  39 }
 40                                                    40 
 41 // -------------------------------------------     41 // --------------------------------------------------------------------
 42 void G4GDMLParameterisation::AddParameter(cons     42 void G4GDMLParameterisation::AddParameter(const PARAMETER& newParameter)
 43 {                                                  43 {
 44   parameterList.push_back(newParameter);           44   parameterList.push_back(newParameter);
 45 }                                                  45 }
 46                                                    46 
 47 // -------------------------------------------     47 // --------------------------------------------------------------------
 48 void G4GDMLParameterisation::ComputeTransforma     48 void G4GDMLParameterisation::ComputeTransformation(
 49   const G4int index, G4VPhysicalVolume* physvo     49   const G4int index, G4VPhysicalVolume* physvol) const
 50 {                                                  50 {
 51   physvol->SetTranslation(parameterList[index]     51   physvol->SetTranslation(parameterList[index].position);
 52   physvol->SetRotation(parameterList[index].pR     52   physvol->SetRotation(parameterList[index].pRot);
 53 }                                                  53 }
 54                                                    54 
 55 // -------------------------------------------     55 // --------------------------------------------------------------------
 56 void G4GDMLParameterisation::ComputeDimensions     56 void G4GDMLParameterisation::ComputeDimensions(G4Box& box, const G4int index,
 57                                                    57                                                const G4VPhysicalVolume*) const
 58 {                                                  58 {
 59   box.SetXHalfLength(parameterList[index].dime     59   box.SetXHalfLength(parameterList[index].dimension[0]);
 60   box.SetYHalfLength(parameterList[index].dime     60   box.SetYHalfLength(parameterList[index].dimension[1]);
 61   box.SetZHalfLength(parameterList[index].dime     61   box.SetZHalfLength(parameterList[index].dimension[2]);
 62 }                                                  62 }
 63                                                    63 
 64 // -------------------------------------------     64 // --------------------------------------------------------------------
 65 void G4GDMLParameterisation::ComputeDimensions     65 void G4GDMLParameterisation::ComputeDimensions(G4Trd& trd, const G4int index,
 66                                                    66                                                const G4VPhysicalVolume*) const
 67 {                                                  67 {
 68   trd.SetXHalfLength1(parameterList[index].dim     68   trd.SetXHalfLength1(parameterList[index].dimension[0]);
 69   trd.SetXHalfLength2(parameterList[index].dim     69   trd.SetXHalfLength2(parameterList[index].dimension[1]);
 70   trd.SetYHalfLength1(parameterList[index].dim     70   trd.SetYHalfLength1(parameterList[index].dimension[2]);
 71   trd.SetYHalfLength2(parameterList[index].dim     71   trd.SetYHalfLength2(parameterList[index].dimension[3]);
 72   trd.SetZHalfLength(parameterList[index].dime     72   trd.SetZHalfLength(parameterList[index].dimension[4]);
 73 }                                                  73 }
 74                                                    74 
 75 // -------------------------------------------     75 // --------------------------------------------------------------------
 76 void G4GDMLParameterisation::ComputeDimensions     76 void G4GDMLParameterisation::ComputeDimensions(G4Trap& trap, const G4int index,
 77                                                    77                                                const G4VPhysicalVolume*) const
 78 {                                                  78 {
 79   trap.SetAllParameters(parameterList[index].d     79   trap.SetAllParameters(parameterList[index].dimension[0],    // Dz
 80                         parameterList[index].d     80                         parameterList[index].dimension[1],    // Theta
 81                         parameterList[index].d     81                         parameterList[index].dimension[2],    // Phi
 82                         parameterList[index].d     82                         parameterList[index].dimension[3],    // Dy1
 83                         parameterList[index].d     83                         parameterList[index].dimension[4],    // Dx1
 84                         parameterList[index].d     84                         parameterList[index].dimension[5],    // Dx2
 85                         parameterList[index].d     85                         parameterList[index].dimension[6],    // pAlp1,
 86                         parameterList[index].d     86                         parameterList[index].dimension[7],    // pDy2,
 87                         parameterList[index].d     87                         parameterList[index].dimension[8],    // pDx3,
 88                         parameterList[index].d     88                         parameterList[index].dimension[9],    // pDx4,
 89                         parameterList[index].d     89                         parameterList[index].dimension[10]);  // pAlp2
 90 }                                                  90 }
 91                                                    91 
 92 // -------------------------------------------     92 // --------------------------------------------------------------------
 93 void G4GDMLParameterisation::ComputeDimensions     93 void G4GDMLParameterisation::ComputeDimensions(G4Tubs& tubs, const G4int index,
 94                                                    94                                                const G4VPhysicalVolume*) const
 95 {                                                  95 {
 96   tubs.SetInnerRadius(parameterList[index].dim     96   tubs.SetInnerRadius(parameterList[index].dimension[0]);
 97   tubs.SetOuterRadius(parameterList[index].dim     97   tubs.SetOuterRadius(parameterList[index].dimension[1]);
 98   tubs.SetZHalfLength(parameterList[index].dim     98   tubs.SetZHalfLength(parameterList[index].dimension[2]);
 99   tubs.SetStartPhiAngle(parameterList[index].d     99   tubs.SetStartPhiAngle(parameterList[index].dimension[3]);
100   tubs.SetDeltaPhiAngle(parameterList[index].d    100   tubs.SetDeltaPhiAngle(parameterList[index].dimension[4]);
101 }                                                 101 }
102                                                   102 
103 // -------------------------------------------    103 // --------------------------------------------------------------------
104 void G4GDMLParameterisation::ComputeDimensions    104 void G4GDMLParameterisation::ComputeDimensions(G4Cons& cons, const G4int index,
105                                                   105                                                const G4VPhysicalVolume*) const
106 {                                                 106 {
107   cons.SetInnerRadiusMinusZ(parameterList[inde    107   cons.SetInnerRadiusMinusZ(parameterList[index].dimension[0]);
108   cons.SetOuterRadiusMinusZ(parameterList[inde    108   cons.SetOuterRadiusMinusZ(parameterList[index].dimension[1]);
109   cons.SetInnerRadiusPlusZ(parameterList[index    109   cons.SetInnerRadiusPlusZ(parameterList[index].dimension[2]);
110   cons.SetOuterRadiusPlusZ(parameterList[index    110   cons.SetOuterRadiusPlusZ(parameterList[index].dimension[3]);
111   cons.SetZHalfLength(parameterList[index].dim    111   cons.SetZHalfLength(parameterList[index].dimension[4]);
112   cons.SetStartPhiAngle(parameterList[index].d    112   cons.SetStartPhiAngle(parameterList[index].dimension[5]);
113   cons.SetDeltaPhiAngle(parameterList[index].d    113   cons.SetDeltaPhiAngle(parameterList[index].dimension[6]);
114 }                                                 114 }
115                                                   115 
116 // -------------------------------------------    116 // --------------------------------------------------------------------
117 void G4GDMLParameterisation::ComputeDimensions    117 void G4GDMLParameterisation::ComputeDimensions(G4Sphere& sphere,
118                                                   118                                                const G4int index,
119                                                   119                                                const G4VPhysicalVolume*) const
120 {                                                 120 {
121   sphere.SetInnerRadius(parameterList[index].d    121   sphere.SetInnerRadius(parameterList[index].dimension[0]);
122   sphere.SetOuterRadius(parameterList[index].d    122   sphere.SetOuterRadius(parameterList[index].dimension[1]);
123   sphere.SetStartPhiAngle(parameterList[index]    123   sphere.SetStartPhiAngle(parameterList[index].dimension[2]);
124   sphere.SetDeltaPhiAngle(parameterList[index]    124   sphere.SetDeltaPhiAngle(parameterList[index].dimension[3]);
125   sphere.SetStartThetaAngle(parameterList[inde    125   sphere.SetStartThetaAngle(parameterList[index].dimension[4]);
126   sphere.SetDeltaThetaAngle(parameterList[inde    126   sphere.SetDeltaThetaAngle(parameterList[index].dimension[5]);
127 }                                                 127 }
128                                                   128 
129 // -------------------------------------------    129 // --------------------------------------------------------------------
130 void G4GDMLParameterisation::ComputeDimensions    130 void G4GDMLParameterisation::ComputeDimensions(G4Orb& orb, const G4int index,
131                                                   131                                                const G4VPhysicalVolume*) const
132 {                                                 132 {
133   orb.SetRadius(parameterList[index].dimension    133   orb.SetRadius(parameterList[index].dimension[0]);
134 }                                                 134 }
135                                                   135 
136 // -------------------------------------------    136 // --------------------------------------------------------------------
137 void G4GDMLParameterisation::ComputeDimensions    137 void G4GDMLParameterisation::ComputeDimensions(G4Ellipsoid& ellipsoid,
138                                                   138                                                const G4int index,
139                                                   139                                                const G4VPhysicalVolume*) const
140 {                                                 140 {
141   ellipsoid.SetSemiAxis(parameterList[index].d    141   ellipsoid.SetSemiAxis(parameterList[index].dimension[0],
142                         parameterList[index].d    142                         parameterList[index].dimension[1],
143                         parameterList[index].d    143                         parameterList[index].dimension[2]);
144   ellipsoid.SetZCuts(parameterList[index].dime    144   ellipsoid.SetZCuts(parameterList[index].dimension[3],
145                      parameterList[index].dime    145                      parameterList[index].dimension[4]);
146 }                                                 146 }
147                                                   147 
148 // -------------------------------------------    148 // --------------------------------------------------------------------
149 void G4GDMLParameterisation::ComputeDimensions    149 void G4GDMLParameterisation::ComputeDimensions(G4Torus& torus,
150                                                   150                                                const G4int index,
151                                                   151                                                const G4VPhysicalVolume*) const
152 {                                                 152 {
153   torus.SetAllParameters(parameterList[index].    153   torus.SetAllParameters(parameterList[index].dimension[0],   // pRmin
154                          parameterList[index].    154                          parameterList[index].dimension[1],   // pRmax
155                          parameterList[index].    155                          parameterList[index].dimension[2],   // pRtor
156                          parameterList[index].    156                          parameterList[index].dimension[3],   // pSPhi
157                          parameterList[index].    157                          parameterList[index].dimension[4]);  // pDPhi
158 }                                                 158 }
159                                                   159 
160 // -------------------------------------------    160 // --------------------------------------------------------------------
161 void G4GDMLParameterisation::ComputeDimensions    161 void G4GDMLParameterisation::ComputeDimensions(G4Para& para, const G4int index,
162                                                   162                                                const G4VPhysicalVolume*) const
163 {                                                 163 {
164   para.SetXHalfLength(parameterList[index].dim    164   para.SetXHalfLength(parameterList[index].dimension[0]);
165   para.SetYHalfLength(parameterList[index].dim    165   para.SetYHalfLength(parameterList[index].dimension[1]);
166   para.SetZHalfLength(parameterList[index].dim    166   para.SetZHalfLength(parameterList[index].dimension[2]);
167   para.SetAlpha(parameterList[index].dimension    167   para.SetAlpha(parameterList[index].dimension[3]);
168   para.SetTanAlpha(std::tan(parameterList[inde    168   para.SetTanAlpha(std::tan(parameterList[index].dimension[3]));
169   para.SetThetaAndPhi(parameterList[index].dim    169   para.SetThetaAndPhi(parameterList[index].dimension[4],
170                       parameterList[index].dim    170                       parameterList[index].dimension[5]);
171 }                                                 171 }
172                                                   172 
173 // -------------------------------------------    173 // --------------------------------------------------------------------
174 void G4GDMLParameterisation::ComputeDimensions    174 void G4GDMLParameterisation::ComputeDimensions(G4Hype& hype, const G4int index,
175                                                   175                                                const G4VPhysicalVolume*) const
176 {                                                 176 {
177   hype.SetInnerRadius(parameterList[index].dim    177   hype.SetInnerRadius(parameterList[index].dimension[0]);
178   hype.SetOuterRadius(parameterList[index].dim    178   hype.SetOuterRadius(parameterList[index].dimension[1]);
179   hype.SetZHalfLength(parameterList[index].dim    179   hype.SetZHalfLength(parameterList[index].dimension[4]);
180   hype.SetInnerStereo(parameterList[index].dim    180   hype.SetInnerStereo(parameterList[index].dimension[2]);
181   hype.SetOuterStereo(parameterList[index].dim    181   hype.SetOuterStereo(parameterList[index].dimension[3]);
182 }                                                 182 }
183                                                   183 
184 // -------------------------------------------    184 // --------------------------------------------------------------------
185 void G4GDMLParameterisation::ComputeDimensions    185 void G4GDMLParameterisation::ComputeDimensions(G4Polycone& pcone,
186                                                   186                                                const G4int index,
187                                                   187                                                const G4VPhysicalVolume*) const
188 {                                                 188 {
189   G4PolyconeHistorical origparam(*(pcone.GetOr    189   G4PolyconeHistorical origparam(*(pcone.GetOriginalParameters()));
190   origparam.Start_angle   = parameterList[inde    190   origparam.Start_angle   = parameterList[index].dimension[0];
191   origparam.Opening_angle = parameterList[inde    191   origparam.Opening_angle = parameterList[index].dimension[1];
192   origparam.Num_z_planes  = (G4int) parameterL    192   origparam.Num_z_planes  = (G4int) parameterList[index].dimension[2];
193   G4int nZplanes          = origparam.Num_z_pl    193   G4int nZplanes          = origparam.Num_z_planes;
194                                                   194 
195   for(G4int ii = 0; ii < nZplanes; ++ii)          195   for(G4int ii = 0; ii < nZplanes; ++ii)
196   {                                               196   {
197     origparam.Rmin[ii]     = parameterList[ind    197     origparam.Rmin[ii]     = parameterList[index].dimension[3 + ii * 3];
198     origparam.Rmax[ii]     = parameterList[ind    198     origparam.Rmax[ii]     = parameterList[index].dimension[4 + ii * 3];
199     origparam.Z_values[ii] = parameterList[ind    199     origparam.Z_values[ii] = parameterList[index].dimension[5 + ii * 3];
200   }                                               200   }
201                                                   201 
202   pcone.SetOriginalParameters(&origparam);  //    202   pcone.SetOriginalParameters(&origparam);  // copy values & transfer pointers
203   pcone.Reset();                            //    203   pcone.Reset();                            // reset to new solid parameters
204 }                                                 204 }
205                                                   205 
206 // -------------------------------------------    206 // --------------------------------------------------------------------
207 void G4GDMLParameterisation::ComputeDimensions    207 void G4GDMLParameterisation::ComputeDimensions(G4Polyhedra& polyhedra,
208                                                   208                                                const G4int index,
209                                                   209                                                const G4VPhysicalVolume*) const
210 {                                                 210 {
211   G4PolyhedraHistorical origparam(*(polyhedra.    211   G4PolyhedraHistorical origparam(*(polyhedra.GetOriginalParameters()));
212   origparam.Start_angle   = parameterList[inde    212   origparam.Start_angle   = parameterList[index].dimension[0];
213   origparam.Opening_angle = parameterList[inde    213   origparam.Opening_angle = parameterList[index].dimension[1];
214   origparam.Num_z_planes  = (G4int) parameterL    214   origparam.Num_z_planes  = (G4int) parameterList[index].dimension[2];
215   origparam.numSide       = (G4int) parameterL    215   origparam.numSide       = (G4int) parameterList[index].dimension[3];
216                                                   216 
217   G4int nZplanes = origparam.Num_z_planes;        217   G4int nZplanes = origparam.Num_z_planes;
218                                                   218 
219   for(G4int ii = 0; ii < nZplanes; ++ii)          219   for(G4int ii = 0; ii < nZplanes; ++ii)
220   {                                               220   {
221     origparam.Rmin[ii]     = parameterList[ind    221     origparam.Rmin[ii]     = parameterList[index].dimension[4 + ii * 3];
222     origparam.Rmax[ii]     = parameterList[ind    222     origparam.Rmax[ii]     = parameterList[index].dimension[5 + ii * 3];
223     origparam.Z_values[ii] = parameterList[ind    223     origparam.Z_values[ii] = parameterList[index].dimension[6 + ii * 3];
224   }                                               224   }
225   polyhedra.SetOriginalParameters(                225   polyhedra.SetOriginalParameters(
226     &origparam);      // copy values & transfe    226     &origparam);      // copy values & transfer pointers
227   polyhedra.Reset();  // reset to new solid pa    227   polyhedra.Reset();  // reset to new solid parameters
228 }                                                 228 }
229                                                   229