2 // ********************************************************************
 
    3 // * License and Disclaimer                                           *
 
    5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
 
    6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
 
    7 // * conditions of the Geant4 Software License,  included in the file *
 
    8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
 
    9 // * include a list of copyright holders.                             *
 
   11 // * Neither the authors of this software system, nor their employing *
 
   12 // * institutes,nor the agencies providing financial support for this *
 
   13 // * work  make  any representation or  warranty, express or implied, *
 
   14 // * regarding  this  software system or assume any liability for its *
 
   15 // * use.  Please see the license in the file  LICENSE  and URL above *
 
   16 // * for the full disclaimer and the limitation of liability.         *
 
   18 // * This  code  implementation is the result of  the  scientific and *
 
   19 // * technical work of the GEANT4 collaboration.                      *
 
   20 // * By using,  copying,  modifying or  distributing the software (or *
 
   21 // * any work based  on the software)  you  agree  to acknowledge its *
 
   22 // * use  in  resulting  scientific  publications,  and indicate your *
 
   23 // * acceptance of all terms of the Geant4 Software license.          *
 
   24 // ********************************************************************
 
   27 // $Id: G4TouchableHistory.icc 67974 2013-03-13 10:17:37Z gcosmo $
 
   30 // class G4TouchableHistory inline implementation
 
   31 // ----------------------------------------------------------------------
 
   33 extern G4GEOM_DLL G4ThreadLocal G4Allocator<G4TouchableHistory> *aTouchableHistoryAllocator;
 
   36 void  G4TouchableHistory::UpdateYourself( G4VPhysicalVolume*   pPhysVol,
 
   37                                     const G4NavigationHistory* pHistory ) 
 
   40   G4AffineTransform tf(fhistory.GetTopTransform().Inverse());
 
   43     // This means that the track has left the World Volume.
 
   44     // Since the Navigation History does not already reflect this,
 
   45     // we must correct this problem here.
 
   47     fhistory.SetFirstEntry(pPhysVol);
 
   49   ftlate = tf.NetTranslation();
 
   50   frot = tf.NetRotation();
 
   54 G4int G4TouchableHistory::CalculateHistoryIndex( G4int stackDepth ) const
 
   56   return (fhistory.GetDepth()-stackDepth); // was -1
 
   60 G4VPhysicalVolume* G4TouchableHistory::GetVolume( G4int depth ) const
 
   62   return fhistory.GetVolume(CalculateHistoryIndex(depth));
 
   66 G4VSolid* G4TouchableHistory::GetSolid( G4int depth ) const
 
   68   return fhistory.GetVolume(CalculateHistoryIndex(depth))
 
   69                             ->GetLogicalVolume()->GetSolid();
 
   73 G4int G4TouchableHistory::GetReplicaNumber( G4int depth ) const
 
   75   return fhistory.GetReplicaNo(CalculateHistoryIndex(depth));
 
   79 G4int G4TouchableHistory::GetHistoryDepth()  const
 
   81   return  fhistory.GetDepth();
 
   85 G4int G4TouchableHistory::MoveUpHistory( G4int num_levels )
 
   87   G4int maxLevelsMove = fhistory.GetDepth();
 
   88   G4int minLevelsMove = 0;              // Cannot redescend today!
 
   89                                         // Soon it will be possible
 
   90                                         // by adding a data member here
 
   92   if( num_levels > maxLevelsMove )
 
   94     num_levels = maxLevelsMove;
 
   96   else if( num_levels < minLevelsMove )
 
   98     num_levels = minLevelsMove;
 
  100   fhistory.BackLevel( num_levels ); 
 
  106 const G4NavigationHistory* G4TouchableHistory::GetHistory() const
 
  111 // There is no provision in case this class is subclassed.
 
  112 // If it is subclassed, this will fail and may not give errors!
 
  115 void* G4TouchableHistory::operator new(size_t)
 
  117   // Once the Navigator calls InitializeAllocator,
 
  118   // the if block below can be removed.
 
  120   if (!aTouchableHistoryAllocator)
 
  122      aTouchableHistoryAllocator = new G4Allocator<G4TouchableHistory>;
 
  124   return (void *) aTouchableHistoryAllocator->MallocSingle();
 
  128 void G4TouchableHistory::operator delete(void *aTH)
 
  130   aTouchableHistoryAllocator->FreeSingle((G4TouchableHistory *) aTH);