Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/geometry/volumes/src/G4AssemblyVolume.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 /geometry/volumes/src/G4AssemblyVolume.cc (Version 11.3.0) and /geometry/volumes/src/G4AssemblyVolume.cc (Version 10.2.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 //
                                                   >>  26 //
                                                   >>  27 // $Id: G4AssemblyVolume.cc 66872 2013-01-15 01:25:57Z japost $
                                                   >>  28 //
 25 //                                                 29 // 
 26 // Class G4AssemblyVolume - implementation         30 // Class G4AssemblyVolume - implementation
 27 //                                                 31 //
 28 // Author: Radovan Chytracek, CERN - November  <<  32 // ----------------------------------------------------------------------
 29 // ------------------------------------------- << 
 30                                                    33 
 31 #include "G4AssemblyVolume.hh"                     34 #include "G4AssemblyVolume.hh"
 32 #include "G4AssemblyStore.hh"                  << 
 33 #include "G4PVPlacement.hh"                        35 #include "G4PVPlacement.hh"
 34 #include "G4RotationMatrix.hh"                     36 #include "G4RotationMatrix.hh"
 35 #include "G4AffineTransform.hh"                    37 #include "G4AffineTransform.hh"
 36 #include "G4LogicalVolume.hh"                      38 #include "G4LogicalVolume.hh"
 37 #include "G4VPhysicalVolume.hh"                    39 #include "G4VPhysicalVolume.hh"
 38 #include "G4ReflectionFactory.hh"                  40 #include "G4ReflectionFactory.hh"
 39                                                    41 
 40 #include <sstream>                                 42 #include <sstream>
 41                                                    43 
 42 G4ThreadLocal unsigned int G4AssemblyVolume::f     44 G4ThreadLocal unsigned int G4AssemblyVolume::fsInstanceCounter = 0;
 43                                                    45 
 44 // ------------------------------------------- << 
 45 // Default constructor                             46 // Default constructor
 46 //                                                 47 //
 47 G4AssemblyVolume::G4AssemblyVolume()               48 G4AssemblyVolume::G4AssemblyVolume()
                                                   >>  49   : fAssemblyID( 0 )
 48 {                                                  50 {
 49   InstanceCountPlus();                             51   InstanceCountPlus();
 50   SetAssemblyID( GetInstanceCount() );             52   SetAssemblyID( GetInstanceCount() );
 51   SetImprintsCount( 0 );                           53   SetImprintsCount( 0 );
 52   G4AssemblyStore* aStore = G4AssemblyStore::G << 
 53   if (aStore->GetAssembly(fAssemblyID,false) ! << 
 54   {                                            << 
 55     std::ostringstream message;                << 
 56     message << "The assembly has NOT been regi << 
 57             << "          Assembly " << fAssem << 
 58             << " already existing in store !"  << 
 59     G4Exception("G4AssemblyVolume::G4AssemblyV << 
 60                 JustWarning, message);         << 
 61   }                                            << 
 62   else                                         << 
 63   {                                            << 
 64     aStore->Register(this);                    << 
 65   }                                            << 
 66 }                                                  54 }
 67                                                    55 
 68 // ------------------------------------------- << 
 69 // Composing constructor                           56 // Composing constructor
 70 //                                                 57 //
 71 G4AssemblyVolume::G4AssemblyVolume( G4LogicalV     58 G4AssemblyVolume::G4AssemblyVolume( G4LogicalVolume* volume,
 72                                     G4ThreeVec     59                                     G4ThreeVector& translation,
 73                                     G4Rotation     60                                     G4RotationMatrix* rotation )
                                                   >>  61   : fAssemblyID( 0 )
 74 {                                                  62 {
 75   InstanceCountPlus();                             63   InstanceCountPlus();
 76   SetAssemblyID( GetInstanceCount() );             64   SetAssemblyID( GetInstanceCount() );
 77   SetImprintsCount( 0 );                           65   SetImprintsCount( 0 );
 78   AddPlacedVolume(volume, translation, rotatio     66   AddPlacedVolume(volume, translation, rotation);
 79   G4AssemblyStore* aStore = G4AssemblyStore::G << 
 80   if (aStore->GetAssembly(fAssemblyID,false) ! << 
 81   {                                            << 
 82     std::ostringstream message;                << 
 83     message << "The assembly has NOT been regi << 
 84             << "          Assembly " << fAssem << 
 85             << " already existing in store !"  << 
 86     G4Exception("G4Assembly::G4Assembly()", "G << 
 87                 JustWarning, message);         << 
 88   }                                            << 
 89   else                                         << 
 90   {                                            << 
 91     aStore->Register(this);                    << 
 92   }                                            << 
 93 }                                                  67 }
 94                                                    68 
 95 // ------------------------------------------- << 
 96 // Destructor                                      69 // Destructor
 97 //                                                 70 //
 98 G4AssemblyVolume::~G4AssemblyVolume()              71 G4AssemblyVolume::~G4AssemblyVolume()
 99 {                                                  72 {
100   std::size_t howmany = fTriplets.size();      <<  73   unsigned int howmany = fTriplets.size();
101   if( howmany != 0 )                               74   if( howmany != 0 )
102   {                                                75   {
103     for( std::size_t i = 0; i < howmany; ++i ) <<  76     for( unsigned int i = 0; i < howmany; i++ )
104     {                                              77     {
105       G4RotationMatrix* pRotToClean = fTriplet     78       G4RotationMatrix* pRotToClean = fTriplets[i].GetRotation();
106       delete pRotToClean;                      <<  79       if( pRotToClean != 0 )
                                                   >>  80       {
                                                   >>  81         delete pRotToClean;
                                                   >>  82       }
107     }                                              83     }
108   }                                                84   }
109   fTriplets.clear();                               85   fTriplets.clear();
110                                                    86   
111   // No need to delete physical volumes, as th <<  87   howmany = fPVStore.size();
112   // by physical-volumes store.                <<  88   if( howmany != 0 )
113   //                                           <<  89   {
                                                   >>  90     for( unsigned int j = 0; j < howmany; j++ )
                                                   >>  91     {
                                                   >>  92       delete fPVStore[j];
                                                   >>  93     }
                                                   >>  94   }
114   fPVStore.clear();                                95   fPVStore.clear();
115   InstanceCountMinus();                            96   InstanceCountMinus();
116   G4AssemblyStore::GetInstance()->DeRegister(t << 
117 }                                                  97 }
118                                                    98 
119 // ------------------------------------------- << 
120 // Add and place the given volume according to     99 // Add and place the given volume according to the specified
121 // translation and rotation.                      100 // translation and rotation.
122 //                                                101 //
123 // The rotation matrix passed in can be 0 = id    102 // The rotation matrix passed in can be 0 = identity or an address even of an
124 // object on the upper stack frame. During ass    103 // object on the upper stack frame. During assembly imprint, it creates anyway
125 // a new matrix and keeps track of it so it ca    104 // a new matrix and keeps track of it so it can delete it later at destruction
126 // time.                                          105 // time.
127 // This policy has been adopted since user has    106 // This policy has been adopted since user has no control on the way the
128 // rotations are combined.                        107 // rotations are combined.
129 //                                                108 //
130 void G4AssemblyVolume::AddPlacedVolume( G4Logi    109 void G4AssemblyVolume::AddPlacedVolume( G4LogicalVolume*  pVolume,
131                                         G4Thre    110                                         G4ThreeVector&    translation,
132                                         G4Rota    111                                         G4RotationMatrix* pRotation )
133 {                                                 112 {
134   auto toStore = new G4RotationMatrix;         << 113   G4RotationMatrix*  toStore  = new G4RotationMatrix;
135                                                   114   
136   if( pRotation != nullptr )  { *toStore = *pR << 115   if( pRotation != 0 )  { *toStore = *pRotation; }
137                                                   116   
138   G4AssemblyTriplet toAdd( pVolume, translatio    117   G4AssemblyTriplet toAdd( pVolume, translation, toStore );
139   fTriplets.push_back( toAdd );                   118   fTriplets.push_back( toAdd );
140 }                                                 119 }
141                                                   120 
142 // ------------------------------------------- << 
143 // Add and place the given volume according to    121 // Add and place the given volume according to the specified transformation
144 //                                                122 //
145 void G4AssemblyVolume::AddPlacedVolume( G4Logi    123 void G4AssemblyVolume::AddPlacedVolume( G4LogicalVolume*  pVolume,
146                                         G4Tran    124                                         G4Transform3D&    transformation )
147 {                                                 125 {
148   // Decompose transformation                     126   // Decompose transformation
149   G4Scale3D     scale;                            127   G4Scale3D     scale;
150   G4Rotate3D    rotation;                         128   G4Rotate3D    rotation;
151   G4Translate3D translation;                      129   G4Translate3D translation;
152   transformation.getDecomposition(scale, rotat    130   transformation.getDecomposition(scale, rotation, translation);
153                                                   131 
154   G4ThreeVector v = translation.getTranslation << 132   G4ThreeVector      v = translation.getTranslation();
155   auto r = new G4RotationMatrix;               << 133   G4RotationMatrix*  r = new G4RotationMatrix;
156       *r = rotation.getRotation();             << 134                     *r = rotation.getRotation();
157                                                   135   
158   G4bool isReflection = false;                    136   G4bool isReflection = false;
159   if (scale(0,0)*scale(1,1)*scale(2,2) < 0.)      137   if (scale(0,0)*scale(1,1)*scale(2,2) < 0.)  { isReflection = true; }
160                                                   138 
161   G4AssemblyTriplet toAdd( pVolume, v, r, isRe    139   G4AssemblyTriplet toAdd( pVolume, v, r, isReflection );
162   fTriplets.push_back( toAdd );                   140   fTriplets.push_back( toAdd );
163 }                                                 141 }
164                                                   142 
165 // ------------------------------------------- << 
166 // Add and place the given assembly volume acc    143 // Add and place the given assembly volume according to the specified
167 // translation and rotation.                      144 // translation and rotation.
168 //                                                145 //
169 void G4AssemblyVolume::AddPlacedAssembly( G4As    146 void G4AssemblyVolume::AddPlacedAssembly( G4AssemblyVolume* pAssembly,
170                                           G4Th    147                                           G4ThreeVector&    translation,
171                                           G4Ro    148                                           G4RotationMatrix* pRotation )
172 {                                                 149 {
173   auto toStore = new G4RotationMatrix;         << 150   G4RotationMatrix*  toStore  = new G4RotationMatrix;
174                                                   151   
175   if( pRotation != nullptr )  { *toStore = *pR << 152   if( pRotation != 0 )  { *toStore = *pRotation; }
176                                                   153   
177   G4AssemblyTriplet toAdd( pAssembly, translat    154   G4AssemblyTriplet toAdd( pAssembly, translation, toStore );
178   fTriplets.push_back( toAdd );                   155   fTriplets.push_back( toAdd );
179 }                                                 156 }
180                                                   157 
181 // ------------------------------------------- << 
182 // Add and place the given assembly volume acc    158 // Add and place the given assembly volume according to the specified 
183 // transformation                                 159 // transformation
184 //                                                160 //
185 void G4AssemblyVolume::AddPlacedAssembly( G4As    161 void G4AssemblyVolume::AddPlacedAssembly( G4AssemblyVolume* pAssembly,
186                                           G4Tr    162                                           G4Transform3D&    transformation )
187 {                                                 163 {
188   // Decompose transformation                     164   // Decompose transformation
189   //                                              165   //
190   G4Scale3D     scale;                            166   G4Scale3D     scale;
191   G4Rotate3D    rotation;                         167   G4Rotate3D    rotation;
192   G4Translate3D translation;                      168   G4Translate3D translation;
193   transformation.getDecomposition(scale, rotat    169   transformation.getDecomposition(scale, rotation, translation);
194                                                   170 
195   G4ThreeVector      v = translation.getTransl    171   G4ThreeVector      v = translation.getTranslation();
196   auto r = new G4RotationMatrix;               << 172   G4RotationMatrix*  r = new G4RotationMatrix;
197       *r = rotation.getRotation();             << 173                     *r = rotation.getRotation();
198                                                   174   
199   G4bool isReflection = false;                    175   G4bool isReflection = false;
200   if (scale(0,0)*scale(1,1)*scale(2,2) < 0.)      176   if (scale(0,0)*scale(1,1)*scale(2,2) < 0.)  { isReflection = true; }
201                                                   177   
202   G4AssemblyTriplet toAdd( pAssembly, v, r, is    178   G4AssemblyTriplet toAdd( pAssembly, v, r, isReflection );
203   fTriplets.push_back( toAdd );                   179   fTriplets.push_back( toAdd );
204 }                                                 180 }
205                                                   181 
206 // ------------------------------------------- << 
207 // Create an instance of an assembly volume in    182 // Create an instance of an assembly volume inside of the specified
208 // mother volume. This works analogically to m    183 // mother volume. This works analogically to making stamp imprints.
209 // This method makes use of the Geant4 affine     184 // This method makes use of the Geant4 affine transformation class.
210 // The algorithm is defined as follows:           185 // The algorithm is defined as follows:
211 //                                                186 //  
212 // Having rotation matrix Rm and translation v    187 // Having rotation matrix Rm and translation vector Tm to be applied
213 // inside the mother and rotation matrix Ra an    188 // inside the mother and rotation matrix Ra and translation vector Ta
214 // to be applied inside the assembly itself fo    189 // to be applied inside the assembly itself for each of the participating
215 // volumes the resulting transformation is        190 // volumes the resulting transformation is
216 //                                                191 //  
217 // Tfinal = Ta * Tm                               192 // Tfinal = Ta * Tm
218 //                                                193 //  
219 // where Ta and Tm are constructed as             194 // where Ta and Tm are constructed as
220 //                                                195 //  
221 //        -1                                      196 //        -1                                     -1
222 // Ta = Ra  * Ta           and            Tm =    197 // Ta = Ra  * Ta           and            Tm = Rm  * Tm
223 //                                                198 //  
224 // which in words means that we create first t    199 // which in words means that we create first the affine transformations
225 // by inverse rotation matrices and translatio    200 // by inverse rotation matrices and translations for mother and assembly.
226 // The resulting final transformation to be ap    201 // The resulting final transformation to be applied to each of the
227 // participating volumes is their product.        202 // participating volumes is their product.
228 //                                                203 //
229 // IMPORTANT NOTE!                                204 // IMPORTANT NOTE!
230 // The order of multiplication is reversed whe    205 // The order of multiplication is reversed when comparing to CLHEP 3D
231 // transformation matrix(G4Transform3D class).    206 // transformation matrix(G4Transform3D class).
232 //                                                207 //  
233 // The rotation matrix passed in can be 0 = id    208 // The rotation matrix passed in can be 0 = identity or an address even of an
234 // object on the upper stack frame. During ass    209 // object on the upper stack frame. During assembly imprint, it creates anyway
235 // a new matrix and keeps track of it so it ca    210 // a new matrix and keeps track of it so it can delete it later at destruction
236 // time.                                          211 // time.
237 // This policy has been adopted since user has    212 // This policy has been adopted since user has no control on the way the
238 // rotations are combined.                        213 // rotations are combined.
239 //                                                214 // 
240 // If the assembly volume contains assembly (a    215 // If the assembly volume contains assembly (a'), the function is called
241 // recursively with composed transformation:      216 // recursively with composed transformation:
242 //                                                217 //
243 // Tanew =  Ta * Ta'                              218 // Tanew =  Ta * Ta'
244 //                                                219 //
245 void G4AssemblyVolume::MakeImprint( G4Assembly    220 void G4AssemblyVolume::MakeImprint( G4AssemblyVolume* pAssembly,
246                                     G4LogicalV    221                                     G4LogicalVolume*  pMotherLV,
247                                     G4Transfor    222                                     G4Transform3D&    transformation,
248                                     G4int copy    223                                     G4int copyNumBase,
249                                     G4bool sur    224                                     G4bool surfCheck )
250 {                                                 225 {
251   std::size_t numberOfDaughters;               << 226   unsigned int  numberOfDaughters;
252                                                << 227   
253   if( copyNumBase == 0 )                          228   if( copyNumBase == 0 )
254   {                                               229   {
255     numberOfDaughters = pMotherLV->GetNoDaught    230     numberOfDaughters = pMotherLV->GetNoDaughters();
256   }                                               231   }
257   else                                            232   else
258   {                                               233   {
259     numberOfDaughters = copyNumBase;              234     numberOfDaughters = copyNumBase;
260   }                                               235   }
261                                                   236 
262   // We start from the first available index      237   // We start from the first available index
263   //                                              238   //
264   ++numberOfDaughters;                         << 239   numberOfDaughters++;
265                                                   240 
266   ImprintsCountPlus();                            241   ImprintsCountPlus();
267                                                   242   
268   auto triplets = pAssembly->fTriplets;        << 243   std::vector<G4AssemblyTriplet> triplets = pAssembly->fTriplets;
269                                                << 
270   // store the transformation in a container ( << 
271   fImprintsTransf[GetImprintsCount()] = transf << 
272                                                   244 
273   for( std::size_t i = 0; i < triplets.size(); << 245   for( unsigned int   i = 0; i < triplets.size(); i++ )
274   {                                               246   {
275     G4Transform3D Ta( *(triplets[i].GetRotatio    247     G4Transform3D Ta( *(triplets[i].GetRotation()),
276                       triplets[i].GetTranslati    248                       triplets[i].GetTranslation() );
277     if ( triplets[i].IsReflection() )  { Ta =     249     if ( triplets[i].IsReflection() )  { Ta = Ta * G4ReflectZ3D(); }
278                                                   250 
279     G4Transform3D Tfinal = transformation * Ta    251     G4Transform3D Tfinal = transformation * Ta;
280                                                   252     
281     if ( triplets[i].GetVolume() != nullptr )  << 253     if ( triplets[i].GetVolume() )
282     {                                             254     {
283       // Generate the unique name for the next    255       // Generate the unique name for the next PV instance
284       // The name has format:                     256       // The name has format:
285       //                                          257       //
286       // av_WWW_impr_XXX_YYY_ZZZ                  258       // av_WWW_impr_XXX_YYY_ZZZ
287       // where the fields mean:                   259       // where the fields mean:
288       // WWW - assembly volume instance number    260       // WWW - assembly volume instance number
289       // XXX - assembly volume imprint number     261       // XXX - assembly volume imprint number
290       // YYY - the name of a log. volume we wa    262       // YYY - the name of a log. volume we want to make a placement of
291       // ZZZ - the log. volume index inside th    263       // ZZZ - the log. volume index inside the assembly volume
292       //                                          264       //
293       std::stringstream pvName;                   265       std::stringstream pvName;
294       pvName << "av_"                             266       pvName << "av_"
295              << GetAssemblyID()                   267              << GetAssemblyID()
296              << "_impr_"                          268              << "_impr_"
297              << GetImprintsCount()                269              << GetImprintsCount()
298              << "_"                               270              << "_"
299              << triplets[i].GetVolume()->GetNa    271              << triplets[i].GetVolume()->GetName().c_str()
300              << "_pv_"                            272              << "_pv_"
301              << i                                 273              << i
302              << std::ends;                        274              << std::ends;
303                                                   275 
304       // Generate a new physical volume instan    276       // Generate a new physical volume instance inside a mother
305       // (as we allow 3D transformation use G4    277       // (as we allow 3D transformation use G4ReflectionFactory to 
306       //  take into account eventual reflectio    278       //  take into account eventual reflection)
307       //                                          279       //
308       G4PhysicalVolumesPair pvPlaced              280       G4PhysicalVolumesPair pvPlaced
309         = G4ReflectionFactory::Instance()->Pla    281         = G4ReflectionFactory::Instance()->Place( Tfinal,
310                                                   282                                                   pvName.str().c_str(),
311                                                   283                                                   triplets[i].GetVolume(),
312                                                   284                                                   pMotherLV,
313                                                   285                                                   false,
314                                                << 286                                                   numberOfDaughters + i,
315                                                   287                                                   surfCheck );
316                                                   288 
317       // Register the physical volume created     289       // Register the physical volume created by us so we can delete it later
318       //                                          290       //
319       fPVStore.push_back( pvPlaced.first );       291       fPVStore.push_back( pvPlaced.first );
320       if ( pvPlaced.second != nullptr )  { fPV << 292       if ( pvPlaced.second )  { fPVStore.push_back( pvPlaced.second ); }
321                                                << 
322       // Here I want to save the imprint trans << 
323       // imprintTrans[GetImprintsCount()] = tr << 
324                                                << 
325     }                                             293     }
326     else if ( triplets[i].GetAssembly() != nul << 294     else if ( triplets[i].GetAssembly() )
327     {                                             295     {
328       // Place volumes in this assembly with c    296       // Place volumes in this assembly with composed transformation
329       //                                          297       //
330       MakeImprint( triplets[i].GetAssembly(),     298       MakeImprint( triplets[i].GetAssembly(), pMotherLV,
331                    Tfinal, (G4int)i*100+copyNu << 299                    Tfinal, i*100+copyNumBase, surfCheck ); 
332     }                                             300     }
333     else                                          301     else
334     {                                             302     {
335       G4Exception("G4AssemblyVolume::MakeImpri    303       G4Exception("G4AssemblyVolume::MakeImprint(..)",
336                   "GeomVol0003", FatalExceptio    304                   "GeomVol0003", FatalException,
337                   "Triplet has no volume and n    305                   "Triplet has no volume and no assembly");
338     }                                             306     }  
339   }                                               307   }  
340 }                                                 308 }    
341                                                   309 
342 // ------------------------------------------- << 
343 void G4AssemblyVolume::MakeImprint( G4LogicalV    310 void G4AssemblyVolume::MakeImprint( G4LogicalVolume*  pMotherLV,
344                                     G4ThreeVec    311                                     G4ThreeVector&    translationInMother,
345                                     G4Rotation    312                                     G4RotationMatrix* pRotationInMother,
346                                     G4int copy    313                                     G4int copyNumBase,
347                                     G4bool sur    314                                     G4bool surfCheck )
348 {                                                 315 {
349   // If needed user can specify explicitely th    316   // If needed user can specify explicitely the base count from which to start
350   // off for the generation of phys. vol. copy    317   // off for the generation of phys. vol. copy numbers.
351   // The old behaviour is preserved when copyN    318   // The old behaviour is preserved when copyNumBase == 0, e.g. the generated
352   // copy numbers start from the count equal t    319   // copy numbers start from the count equal to current number of daughter
353   // volumes before an imprint is made            320   // volumes before an imprint is made
354                                                   321 
355   // Compose transformation                       322   // Compose transformation
356   //                                              323   //
357   if( pRotationInMother == nullptr )           << 324   if( pRotationInMother == 0 )
358   {                                               325   {
359     // Make it by default an indentity matrix     326     // Make it by default an indentity matrix
360     //                                            327     //
361     pRotationInMother =                           328     pRotationInMother =
362       const_cast<G4RotationMatrix*>( &G4Rotati    329       const_cast<G4RotationMatrix*>( &G4RotationMatrix::IDENTITY );
363   }                                               330   }
364                                                   331 
365   G4Transform3D transform( *pRotationInMother,    332   G4Transform3D transform( *pRotationInMother,
366                             translationInMothe    333                             translationInMother );
367   MakeImprint(this, pMotherLV, transform, copy    334   MakeImprint(this, pMotherLV, transform, copyNumBase, surfCheck);
368 }                                                 335 }
369                                                   336 
370 // ------------------------------------------- << 
371 void G4AssemblyVolume::MakeImprint( G4LogicalV    337 void G4AssemblyVolume::MakeImprint( G4LogicalVolume*  pMotherLV,
372                                     G4Transfor    338                                     G4Transform3D&    transformation,
373                                     G4int copy    339                                     G4int copyNumBase,
374                                     G4bool sur    340                                     G4bool surfCheck )
375 {                                                 341 {
376   // If needed user can specify explicitely th    342   // If needed user can specify explicitely the base count from which to start
377   // off for the generation of phys. vol. copy    343   // off for the generation of phys. vol. copy numbers.
378   // The old behaviour is preserved when copyN    344   // The old behaviour is preserved when copyNumBase == 0, e.g. the generated
379   // copy numbers start from the count equal t    345   // copy numbers start from the count equal to current number of daughter
380   // volumes before a imprint is made             346   // volumes before a imprint is made
381                                                   347 
382   MakeImprint(this, pMotherLV, transformation,    348   MakeImprint(this, pMotherLV, transformation, copyNumBase, surfCheck);
383 }                                                 349 }
384                                                   350 
385 // ------------------------------------------- << 
386 unsigned int G4AssemblyVolume::GetInstanceCoun    351 unsigned int G4AssemblyVolume::GetInstanceCount() const
387 {                                                 352 {
388   return G4AssemblyVolume::fsInstanceCounter;     353   return G4AssemblyVolume::fsInstanceCounter;
389 }                                                 354 }
390                                                   355 
391 // ------------------------------------------- << 356 void         G4AssemblyVolume::SetInstanceCount( unsigned int value )
392 void G4AssemblyVolume::SetInstanceCount( unsig << 
393 {                                                 357 {
394   G4AssemblyVolume::fsInstanceCounter = value;    358   G4AssemblyVolume::fsInstanceCounter = value;
395 }                                                 359 }
396                                                   360 
397 // ------------------------------------------- << 361 void         G4AssemblyVolume::InstanceCountPlus()
398 void G4AssemblyVolume::InstanceCountPlus()     << 
399 {                                                 362 {
400   G4AssemblyVolume::fsInstanceCounter++;          363   G4AssemblyVolume::fsInstanceCounter++;
401 }                                                 364 }
402                                                   365 
403 // ------------------------------------------- << 366 void         G4AssemblyVolume::InstanceCountMinus()
404 void G4AssemblyVolume::InstanceCountMinus()    << 
405 {                                                 367 {
406   G4AssemblyVolume::fsInstanceCounter--;          368   G4AssemblyVolume::fsInstanceCounter--;
407 }                                                 369 }
408                                                   370