TG4RootNavMgr.cxx

Go to the documentation of this file.

// @(#)root/g4root:$Id$
// Author: Andrei Gheata   07/08/06
 
/*************************************************************************
 * Copyright (C) 1995-2000, Rene Brun and Fons Rademakers.               *
 * All rights reserved.                                                  *
 *                                                                       *
 * For the licensing terms see $ROOTSYS/LICENSE.                         *
 * For the list of contributors see $ROOTSYS/README/CREDITS.             *
 *************************************************************************/

 
#include "TG4RootNavMgr.h"
#include "TG4RootDetectorConstruction.h"
#include "TG4RootNavigator.h"
 
#include "G4PropagatorInField.hh"
#include "G4RunManager.hh"
#include "G4TransportationManager.hh"
 
#include "TGeoManager.h"

// ClassImp(TG4RootNavMgr)
 
G4ThreadLocal TG4RootNavMgr* TG4RootNavMgr::fRootNavMgr = 0;
TG4RootNavMgr* TG4RootNavMgr::fgMasterInstance = 0;
 
//______________________________________________________________________________
TG4RootNavMgr::TG4RootNavMgr()
  : TObject(),
    fGeometry(0),
    fNavigator(0),
    fDetConstruction(0),
    fConnected(kFALSE)
{
}
 
//______________________________________________________________________________
TG4RootNavMgr::TG4RootNavMgr(
  TGeoManager* geom, TG4RootDetectorConstruction* detConstruction)
  : TObject(),
    fGeometry(geom),
    fNavigator(0),
    fDetConstruction(detConstruction),
    fConnected(kFALSE)
{
  if (!detConstruction) {
    fDetConstruction = new TG4RootDetectorConstruction(geom);
    SetNavigator(new TG4RootNavigator());
  }
  else {
    SetNavigator(new TG4RootNavigator(detConstruction));
  }
}
 
//______________________________________________________________________________
TG4RootNavMgr::~TG4RootNavMgr()
{
  //   if (fNavigator) delete fNavigator;
  if (fDetConstruction) delete fDetConstruction;
  fRootNavMgr = 0;
 
  G4bool isMaster = !G4Threading::IsWorkerThread();
  if (isMaster) {
    fgMasterInstance = 0;
  }
}
 
//______________________________________________________________________________
TG4RootNavMgr* TG4RootNavMgr::GetInstance(TGeoManager* geom)
{
  if (fRootNavMgr) return fRootNavMgr;
  // Check if we have to create one.
  if (!geom) return NULL;
  fRootNavMgr = new TG4RootNavMgr(geom);
  G4bool isMaster = !G4Threading::IsWorkerThread();
  Printf("isMaster=%d", isMaster);
  if (isMaster) {
    fgMasterInstance = fRootNavMgr;
  }
  return fRootNavMgr;
}
 
//______________________________________________________________________________
TG4RootNavMgr* TG4RootNavMgr::GetInstance(const TG4RootNavMgr& navMgr)
{
  if (fRootNavMgr) return fRootNavMgr;
  // Check if we have to create one.
  fRootNavMgr = new TG4RootNavMgr(navMgr.fGeometry, navMgr.fDetConstruction);
  G4bool isMaster = !G4Threading::IsWorkerThread();
  if (isMaster) {
    fgMasterInstance = fRootNavMgr;
  }
  return fRootNavMgr;
}
 
//_____________________________________________________________________________
TG4RootNavMgr* TG4RootNavMgr::GetMasterInstance()
{
  return fgMasterInstance;
}
 
//______________________________________________________________________________
Bool_t TG4RootNavMgr::ConnectToG4()
{
  if (fConnected) {
    Info("ConnectToG4", "Already connected");
    return kTRUE;
  }
  if (!fDetConstruction) {
    Error("ConnectToG4", "No detector construction set !");
    return kFALSE;
  }
 
  if (!fNavigator) {
    Error(
      "ConnectToG4", "Navigator has to be created befor connecting to G4 !!!");
    return kFALSE;
  }
 
  G4RunManager* runManager = G4RunManager::GetRunManager();
  if (!runManager) {
    Error("ConnectToG4", "Unable to connect: G4RunManager not instantiated");
    return kFALSE;
  }
  G4bool isMaster = !G4Threading::IsWorkerThread();
  if (isMaster) runManager->SetUserInitialization(fDetConstruction);
  Info("ConnectToG4",
    "ROOT detector construction class connected to G4RunManager");
  fConnected = kTRUE;
  return kTRUE;
}
 
//______________________________________________________________________________
void TG4RootNavMgr::SetNavigator(TG4RootNavigator* nav)
{
  if (fConnected) {
    Error("SetNavigator",
      "Navigator set after instantiation of G4RunManager. Won't set!!!");
    return;
  }
  G4TransportationManager* trMgr =
    G4TransportationManager::GetTransportationManager();
  //   G4Navigator *oldNav = trMgr->GetNavigatorForTracking();
  trMgr->SetNavigatorForTracking(nav);
  G4FieldManager* fieldMgr =
    trMgr->GetPropagatorInField()->GetCurrentFieldManager();
  delete trMgr->GetPropagatorInField();
  trMgr->SetPropagatorInField(new G4PropagatorInField(nav, fieldMgr));
  trMgr->ActivateNavigator(nav);
  G4EventManager* evtMgr = G4EventManager::GetEventManager();
  if (evtMgr)
    evtMgr->GetTrackingManager()->GetSteppingManager()->SetNavigator(nav);
  fNavigator = nav;
  //   trMgr->DeRegisterNavigator(oldNav);
  Info("SetNavigator",
    "TG4RootNavigator created and registered to G4TransportationManager");
}
 
//______________________________________________________________________________
void TG4RootNavMgr::Initialize(
  TVirtualUserPostDetConstruction* sdinit, Int_t nthreads)
{
  Info("Initialize", "Creating G4 hierarchy ...");
  if (fDetConstruction) fDetConstruction->Initialize(sdinit);
  if (nthreads > 1) gGeoManager->SetMaxThreads(nthreads);
}
 
//______________________________________________________________________________
void TG4RootNavMgr::LocateGlobalPointAndSetup(Double_t* pt, Double_t* dir)
{
  G4ThreeVector point(pt[0], pt[1], pt[2]);
  if (dir) {
    G4ThreeVector direction(dir[0], dir[1], dir[2]);
    fNavigator->LocateGlobalPointAndSetup(point, &direction);
  }
  else {
    fNavigator->LocateGlobalPointAndSetup(point);
  }
  fNavigator->PrintState();
}
 
//______________________________________________________________________________
void TG4RootNavMgr::SetVerboseLevel(Int_t level)
{
  fNavigator->SetVerboseLevel(level);
}
 
//______________________________________________________________________________
void TG4RootNavMgr::PrintG4State() const
{
  G4NavigationHistory* history = fNavigator->GetHistory();
  G4cout << *history << G4endl;
}
 
//______________________________________________________________________________
void TG4RootNavMgr::SetG4TrackingManager(G4TrackingManager* trackingManager)
{
  fNavigator->SetG4TrackingManager(trackingManager);
}
 
//______________________________________________________________________________
void TG4RootNavMgr::SetGeometryRestoreFunction(
  std::function<Bool_t(Int_t)> restoreGeoStateFunction)
{
  fNavigator->SetGeometryRestoreFunction(restoreGeoStateFunction);
}