VecGeomG4PropagatingNavigator.h

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
#ifndef O2_SIMSETUP_VECGEOMG4PROPAGATINGNAVIGATOR_H_
#define O2_SIMSETUP_VECGEOMG4PROPAGATINGNAVIGATOR_H_
 
#include "VecGeomG4NavigatorBase.h"
 
namespace o2::simsetup
{
 
class VecGeomG4PropagatingNavigator : public VecGeomG4NavigatorBase
{
 public:
  VecGeomG4PropagatingNavigator(VecGeomG4Map const& map, bool zeroSafety);
  ~VecGeomG4PropagatingNavigator() override;
 
  G4double ComputeStep(const G4ThreeVector& globalPoint, const G4ThreeVector& direction,
                       const G4double proposedStepLength, G4double& newSafety) override;
 
  G4VPhysicalVolume* ResetHierarchyAndLocate(const G4ThreeVector& point, const G4ThreeVector& direction,
                                             const G4TouchableHistory& history) override;
 
  G4VPhysicalVolume* LocateGlobalPointAndSetup(const G4ThreeVector& point, const G4ThreeVector* direction = nullptr,
                                               const G4bool relativeSearch = true,
                                               const G4bool ignoreDirection = true) override;
 
  void LocateGlobalPointWithinVolume(const G4ThreeVector& position) override;
 
  G4double ComputeSafety(const G4ThreeVector& globalPoint, const G4double proposedMaxLength = DBL_MAX,
                         const G4bool keepState = true) override;
 
  G4ThreeVector GetLocalExitNormal(G4bool* valid) override;
  G4ThreeVector GetGlobalExitNormal(const G4ThreeVector& point, G4bool* valid) override;
 
 private:
  vecgeom::NavigationState mCurState;   
  vecgeom::NavigationState mNextState;  
  vecgeom::NavigationState mPrevState;  
  vecgeom::NavigationState mEmptyState; 
 
  bool mZeroSafety = false;
  bool mHaveNextState = false; 
 
  G4ThreeVector mNextPoint{-1e8, -1e8, -1e8}; 
  G4ThreeVector mLastDirection{0, 0, 1};      
 
  bool mWouldEnter = false;                      
  bool mWouldExit = false;                       
  bool mOnBoundary = false;                      
  bool mForceReInit = false;                     
  bool mCrossed = false;                         
  bool mExitBlockPending = false;                
  G4ThreeVector mLocatedPoint{-1e8, -1e8, -1e8}; 
 
  int mZeroSteps = 0;
  long mNudgedSteps = 0;
  long mGlobalRelocates = 0;
};
 
} // namespace o2::simsetup
 
#endif