G4MaterialScanner Class Reference

#include <G4MaterialScanner.hh>

Public Member Functions
	G4MaterialScanner ()
G4ThreeVector	GetEyePosition () const
G4int	GetNPhi () const
G4int	GetNTheta () const
G4double	GetPhiMin () const
G4double	GetPhiSpan () const
const G4String &	GetRegionName () const
G4bool	GetRegionSensitive () const
G4double	GetThetaMin () const
G4double	GetThetaSpan () const
void	Scan ()
void	SetEyePosition (const G4ThreeVector &val)
void	SetNPhi (G4int val)
void	SetNTheta (G4int val)
void	SetPhiMin (G4double val)
void	SetPhiSpan (G4double val)
G4bool	SetRegionName (const G4String &val)
void	SetRegionSensitive (G4bool val=true)
void	SetThetaMin (G4double val)
void	SetThetaSpan (G4double val)
	~G4MaterialScanner ()

Private Member Functions
void	DoScan ()
void	RestoreUserActions ()
void	StoreUserActions ()

Private Attributes
G4ThreeVector	eyeDirection
G4ThreeVector	eyePosition
G4int	nPhi = 37
G4int	nTheta = 91
G4double	phiMin = 0.0
G4double	phiSpan = 0.0
G4String	regionName = "notDefined"
G4bool	regionSensitive = false
G4EventManager *	theEventManager = nullptr
G4UserEventAction *	theMatScannerEventAction = nullptr
G4UserStackingAction *	theMatScannerStackingAction = nullptr
G4MSSteppingAction *	theMatScannerSteppingAction = nullptr
G4UserTrackingAction *	theMatScannerTrackingAction = nullptr
G4MatScanMessenger *	theMessenger = nullptr
G4RayShooter *	theRayShooter = nullptr
G4Region *	theRegion = nullptr
G4double	thetaMin = 0.0
G4double	thetaSpan = 0.0
G4UserEventAction *	theUserEventAction = nullptr
G4UserStackingAction *	theUserStackingAction = nullptr
G4UserSteppingAction *	theUserSteppingAction = nullptr
G4UserTrackingAction *	theUserTrackingAction = nullptr

Detailed Description

Definition at line 52 of file G4MaterialScanner.hh.

Constructor & Destructor Documentation

G4MaterialScanner()

G4MaterialScanner::G4MaterialScanner

(

)

Definition at line 50 of file G4MaterialScanner.cc.

{
  theRayShooter   = new G4RayShooter();
  theMessenger    = new G4MatScanMessenger(this);
  theEventManager = G4EventManager::GetEventManager();
 
  eyePosition = G4ThreeVector(0., 0., 0.);
  thetaSpan   = 90. * deg;
  phiSpan     = 360. * deg;
}

References deg, eyePosition, G4EventManager::GetEventManager(), phiSpan, theEventManager, theMessenger, theRayShooter, and thetaSpan.

~G4MaterialScanner()

G4MaterialScanner::~G4MaterialScanner

(

)

Definition at line 62 of file G4MaterialScanner.cc.

{
  delete theRayShooter;
  delete theMatScannerSteppingAction;
  delete theMessenger;
}

References theMatScannerSteppingAction, theMessenger, and theRayShooter.

Member Function Documentation

DoScan()

void G4MaterialScanner::DoScan ( )

private

= G4Geantino::GeantinoDefinition()->GetProcessManager()->GetProcessList();

Definition at line 129 of file G4MaterialScanner.cc.

{
  // Confirm material table is updated
  G4RunManagerKernel::GetRunManagerKernel()->UpdateRegion();
 
 
  // Close geometry and set the application state
  G4GeometryManager* geomManager = G4GeometryManager::GetInstance();
  geomManager->OpenGeometry();
  geomManager->CloseGeometry(1, 0);
 
  G4ThreeVector center(0, 0, 0);
  G4Navigator* navigator = G4TransportationManager::GetTransportationManager()
                             ->GetNavigatorForTracking();
  navigator->LocateGlobalPointAndSetup(center, 0, false);
 
  G4StateManager* theStateMan = G4StateManager::GetStateManager();
  theStateMan->SetNewState(G4State_GeomClosed);
 
  // Event loop
  G4int iEvent = 0;
  for(G4int iTheta = 0; iTheta < nTheta; ++iTheta)
  {
    G4double theta = thetaMin;
    if(iTheta > 0)
      theta += G4double(iTheta) * thetaSpan / G4double(nTheta - 1);
    G4double aveLength = 0.;
    G4double aveX0     = 0.;
    G4double aveLambda = 0.;
    G4cout << G4endl;
    G4cout
      << "         Theta(deg)    Phi(deg)  Length(mm)          x0     lambda0"
      << G4endl;
    G4cout << G4endl;
    for(G4int iPhi = 0; iPhi < nPhi; ++iPhi)
    {
      G4Event* anEvent = new G4Event(iEvent++);
      G4double phi     = phiMin;
      if(iPhi > 0)
        phi += G4double(iPhi) * phiSpan / G4double(nPhi - 1);
      eyeDirection =
        G4ThreeVector(std::cos(theta) * std::cos(phi),
                      std::cos(theta) * std::sin(phi), std::sin(theta));
      theRayShooter->Shoot(anEvent, eyePosition, eyeDirection);
      theMatScannerSteppingAction->Initialize(regionSensitive, theRegion);
      theEventManager->ProcessOneEvent(anEvent);
      G4double length = theMatScannerSteppingAction->GetTotalStepLength();
      G4double x0     = theMatScannerSteppingAction->GetX0();
      G4double lambda = theMatScannerSteppingAction->GetLambda0();
 
      G4cout << "        " << std::setw(11) << theta / deg << " "
             << std::setw(11) << phi / deg << " " << std::setw(11)
             << length / mm << " " << std::setw(11) << x0 << " "
             << std::setw(11) << lambda << G4endl;
      aveLength += length / mm;
      aveX0 += x0;
      aveLambda += lambda;
    }
    if(nPhi > 1)
    {
      G4cout << G4endl;
      G4cout << " ave. for theta = " << std::setw(11) << theta / deg << " : "
             << std::setw(11) << aveLength / nPhi << " " << std::setw(11)
             << aveX0 / nPhi << " " << std::setw(11) << aveLambda / nPhi
             << G4endl;
    }
  }
 
  theStateMan->SetNewState(G4State_Idle);
  return;
}

References G4GeometryManager::CloseGeometry(), deg, eyeDirection, eyePosition, G4cout, G4endl, G4State_GeomClosed, G4State_Idle, G4GeometryManager::GetInstance(), G4MSSteppingAction::GetLambda0(), G4TransportationManager::GetNavigatorForTracking(), G4RunManagerKernel::GetRunManagerKernel(), G4StateManager::GetStateManager(), G4MSSteppingAction::GetTotalStepLength(), G4TransportationManager::GetTransportationManager(), G4MSSteppingAction::GetX0(), G4MSSteppingAction::Initialize(), G4InuclParticleNames::lambda, mm, write_gdml::navigator, nPhi, nTheta, G4GeometryManager::OpenGeometry(), phiMin, phiSpan, G4EventManager::ProcessOneEvent(), regionSensitive, G4StateManager::SetNewState(), G4RayShooter::Shoot(), theEventManager, theMatScannerSteppingAction, theRayShooter, theRegion, thetaMin, thetaSpan, and G4RunManagerKernel::UpdateRegion().

Referenced by Scan().

GetEyePosition()

G4ThreeVector G4MaterialScanner::GetEyePosition ( ) const

inline

Definition at line 64 of file G4MaterialScanner.hh.

{ return eyePosition; }

References eyePosition.

Referenced by G4MatScanMessenger::GetCurrentValue().

GetNPhi()

G4int G4MaterialScanner::GetNPhi ( ) const

inline

Definition at line 72 of file G4MaterialScanner.hh.

{ return nPhi; }

References nPhi.

Referenced by G4MatScanMessenger::GetCurrentValue(), and G4MatScanMessenger::SetNewValue().

GetNTheta()

G4int G4MaterialScanner::GetNTheta ( ) const

inline

Definition at line 66 of file G4MaterialScanner.hh.

{ return nTheta; }

References nTheta.

Referenced by G4MatScanMessenger::GetCurrentValue(), and G4MatScanMessenger::SetNewValue().

GetPhiMin()

G4double G4MaterialScanner::GetPhiMin ( ) const

inline

Definition at line 74 of file G4MaterialScanner.hh.

{ return phiMin; }

References phiMin.

Referenced by G4MatScanMessenger::GetCurrentValue(), and G4MatScanMessenger::SetNewValue().

GetPhiSpan()

G4double G4MaterialScanner::GetPhiSpan ( ) const

inline

Definition at line 76 of file G4MaterialScanner.hh.

{ return phiSpan; }

References phiSpan.

Referenced by G4MatScanMessenger::GetCurrentValue(), and G4MatScanMessenger::SetNewValue().

GetRegionName()

const G4String & G4MaterialScanner::GetRegionName ( ) const

inline

Definition at line 80 of file G4MaterialScanner.hh.

{ return regionName; }

References regionName.

Referenced by G4MatScanMessenger::GetCurrentValue().

GetRegionSensitive()

G4bool G4MaterialScanner::GetRegionSensitive ( ) const

inline

Definition at line 78 of file G4MaterialScanner.hh.

{ return regionSensitive; }

References regionSensitive.

Referenced by G4MatScanMessenger::GetCurrentValue().

GetThetaMin()

G4double G4MaterialScanner::GetThetaMin ( ) const

inline

Definition at line 68 of file G4MaterialScanner.hh.

{ return thetaMin; }

References thetaMin.

Referenced by G4MatScanMessenger::GetCurrentValue(), and G4MatScanMessenger::SetNewValue().

GetThetaSpan()

G4double G4MaterialScanner::GetThetaSpan ( ) const

inline

Definition at line 70 of file G4MaterialScanner.hh.

{ return thetaSpan; }

References thetaSpan.

Referenced by G4MatScanMessenger::GetCurrentValue(), and G4MatScanMessenger::SetNewValue().

RestoreUserActions()

void G4MaterialScanner::RestoreUserActions ( )

private

Definition at line 114 of file G4MaterialScanner.cc.

{
  theEventManager->SetUserAction(theUserEventAction);
  theEventManager->SetUserAction(theUserStackingAction);
  theEventManager->SetUserAction(theUserTrackingAction);
  theEventManager->SetUserAction(theUserSteppingAction);
 
  G4SDManager* theSDMan = G4SDManager::GetSDMpointerIfExist();
  if(theSDMan != nullptr)
  {
    theSDMan->Activate("/", true);
  }
}

References G4SDManager::Activate(), G4SDManager::GetSDMpointerIfExist(), G4EventManager::SetUserAction(), theEventManager, theUserEventAction, theUserStackingAction, theUserSteppingAction, and theUserTrackingAction.

Referenced by Scan().

Scan()

void G4MaterialScanner::Scan

(

)

Definition at line 70 of file G4MaterialScanner.cc.

{
  G4StateManager* theStateMan     = G4StateManager::GetStateManager();
  G4ApplicationState currentState = theStateMan->GetCurrentState();
  if(currentState != G4State_Idle)
  {
    G4cerr << "Illegal application state - Scan() ignored." << G4endl;
    return;
  }
 
  if(theMatScannerSteppingAction == nullptr)
  {
    theMatScannerSteppingAction = new G4MSSteppingAction();
  }
  StoreUserActions();
  DoScan();
  RestoreUserActions();
}

References DoScan(), G4cerr, G4endl, G4State_Idle, G4StateManager::GetCurrentState(), G4StateManager::GetStateManager(), RestoreUserActions(), StoreUserActions(), and theMatScannerSteppingAction.

Referenced by G4MatScanMessenger::SetNewValue().

SetEyePosition()

void G4MaterialScanner::SetEyePosition ( const G4ThreeVector &  val )

inline

Definition at line 63 of file G4MaterialScanner.hh.

{ eyePosition = val; }

References eyePosition.

Referenced by G4MatScanMessenger::SetNewValue().

SetNPhi()

void G4MaterialScanner::SetNPhi ( G4int  val )

inline

Definition at line 71 of file G4MaterialScanner.hh.

{ nPhi = val; }

References nPhi.

Referenced by G4MatScanMessenger::SetNewValue().

SetNTheta()

void G4MaterialScanner::SetNTheta ( G4int  val )

inline

Definition at line 65 of file G4MaterialScanner.hh.

{ nTheta = val; }

References nTheta.

Referenced by G4MatScanMessenger::SetNewValue().

SetPhiMin()

void G4MaterialScanner::SetPhiMin ( G4double  val )

inline

Definition at line 73 of file G4MaterialScanner.hh.

{ phiMin = val; }

References phiMin.

Referenced by G4MatScanMessenger::SetNewValue().

SetPhiSpan()

void G4MaterialScanner::SetPhiSpan ( G4double  val )

inline

Definition at line 75 of file G4MaterialScanner.hh.

{ phiSpan = val; }

References phiSpan.

Referenced by G4MatScanMessenger::SetNewValue().

SetRegionName()

G4bool G4MaterialScanner::SetRegionName

(

const G4String &

val

)

Definition at line 209 of file G4MaterialScanner.cc.

{
  G4Region* aRegion = G4RegionStore::GetInstance()->GetRegion(val);
  if(aRegion != nullptr)
  {
    theRegion  = aRegion;
    regionName = val;
    return true;
  }
  else
  {
    G4cerr << "Region <" << val << "> not found. Command ignored." << G4endl;
    G4cerr << "Defined regions are : " << G4endl;
    for(std::size_t i = 0; i < G4RegionStore::GetInstance()->size(); ++i)
    {
      G4cerr << " " << (*(G4RegionStore::GetInstance()))[i]->GetName();
    }
    G4cerr << G4endl;
    return false;
  }
}

References G4cerr, G4endl, G4RegionStore::GetInstance(), G4RegionStore::GetRegion(), regionName, and theRegion.

Referenced by G4MatScanMessenger::SetNewValue().

SetRegionSensitive()

void G4MaterialScanner::SetRegionSensitive ( G4bool  val = true )

inline

Definition at line 77 of file G4MaterialScanner.hh.

{ regionSensitive = val; }

References regionSensitive.

Referenced by G4MatScanMessenger::SetNewValue().

SetThetaMin()

void G4MaterialScanner::SetThetaMin ( G4double  val )

inline

Definition at line 67 of file G4MaterialScanner.hh.

{ thetaMin = val; }

References thetaMin.

Referenced by G4MatScanMessenger::SetNewValue().

SetThetaSpan()

void G4MaterialScanner::SetThetaSpan ( G4double  val )

inline

Definition at line 69 of file G4MaterialScanner.hh.

{ thetaSpan = val; }

References thetaSpan.

Referenced by G4MatScanMessenger::SetNewValue().

StoreUserActions()

void G4MaterialScanner::StoreUserActions ( )

private

Definition at line 90 of file G4MaterialScanner.cc.

{
  theUserEventAction    = theEventManager->GetUserEventAction();
  theUserStackingAction = theEventManager->GetUserStackingAction();
  theUserTrackingAction = theEventManager->GetUserTrackingAction();
  theUserSteppingAction = theEventManager->GetUserSteppingAction();
 
  theEventManager->SetUserAction(theMatScannerEventAction);
  theEventManager->SetUserAction(theMatScannerStackingAction);
  theEventManager->SetUserAction(theMatScannerTrackingAction);
  theEventManager->SetUserAction(theMatScannerSteppingAction);
 
  G4SDManager* theSDMan = G4SDManager::GetSDMpointerIfExist();
  if(theSDMan != nullptr)
  {
    theSDMan->Activate("/", false);
  }
 
  G4GeometryManager* theGeomMan = G4GeometryManager::GetInstance();
  theGeomMan->OpenGeometry();
  theGeomMan->CloseGeometry(true);
}

References G4SDManager::Activate(), G4GeometryManager::CloseGeometry(), G4GeometryManager::GetInstance(), G4SDManager::GetSDMpointerIfExist(), G4EventManager::GetUserEventAction(), G4EventManager::GetUserStackingAction(), G4EventManager::GetUserSteppingAction(), G4EventManager::GetUserTrackingAction(), G4GeometryManager::OpenGeometry(), G4EventManager::SetUserAction(), theEventManager, theMatScannerEventAction, theMatScannerStackingAction, theMatScannerSteppingAction, theMatScannerTrackingAction, theUserEventAction, theUserStackingAction, theUserSteppingAction, and theUserTrackingAction.

Referenced by Scan().

Field Documentation

eyeDirection

G4ThreeVector G4MaterialScanner::eyeDirection

private

Definition at line 116 of file G4MaterialScanner.hh.

Referenced by DoScan().

eyePosition

G4ThreeVector G4MaterialScanner::eyePosition

private

Definition at line 108 of file G4MaterialScanner.hh.

Referenced by DoScan(), G4MaterialScanner(), GetEyePosition(), and SetEyePosition().

nPhi

G4int G4MaterialScanner::nPhi = 37

private

Definition at line 112 of file G4MaterialScanner.hh.

Referenced by DoScan(), GetNPhi(), and SetNPhi().

nTheta

G4int G4MaterialScanner::nTheta = 91

private

Definition at line 109 of file G4MaterialScanner.hh.

Referenced by DoScan(), GetNTheta(), and SetNTheta().

phiMin

G4double G4MaterialScanner::phiMin = 0.0

private

Definition at line 113 of file G4MaterialScanner.hh.

Referenced by DoScan(), GetPhiMin(), and SetPhiMin().

phiSpan

G4double G4MaterialScanner::phiSpan = 0.0

private

Definition at line 114 of file G4MaterialScanner.hh.

Referenced by DoScan(), G4MaterialScanner(), GetPhiSpan(), and SetPhiSpan().

regionName

G4String G4MaterialScanner::regionName = "notDefined"

private

Definition at line 119 of file G4MaterialScanner.hh.

Referenced by GetRegionName(), and SetRegionName().

regionSensitive

G4bool G4MaterialScanner::regionSensitive = false

private

Definition at line 118 of file G4MaterialScanner.hh.

Referenced by DoScan(), GetRegionSensitive(), and SetRegionSensitive().

theEventManager

G4EventManager* G4MaterialScanner::theEventManager = nullptr

private

Definition at line 96 of file G4MaterialScanner.hh.

Referenced by DoScan(), G4MaterialScanner(), RestoreUserActions(), and StoreUserActions().

theMatScannerEventAction

G4UserEventAction* G4MaterialScanner::theMatScannerEventAction = nullptr

private

Definition at line 103 of file G4MaterialScanner.hh.

Referenced by StoreUserActions().

theMatScannerStackingAction

G4UserStackingAction* G4MaterialScanner::theMatScannerStackingAction = nullptr

private

Definition at line 104 of file G4MaterialScanner.hh.

Referenced by StoreUserActions().

theMatScannerSteppingAction

G4MSSteppingAction* G4MaterialScanner::theMatScannerSteppingAction = nullptr

private

Definition at line 106 of file G4MaterialScanner.hh.

Referenced by DoScan(), Scan(), StoreUserActions(), and ~G4MaterialScanner().

theMatScannerTrackingAction

G4UserTrackingAction* G4MaterialScanner::theMatScannerTrackingAction = nullptr

private

Definition at line 105 of file G4MaterialScanner.hh.

Referenced by StoreUserActions().

theMessenger

G4MatScanMessenger* G4MaterialScanner::theMessenger = nullptr

private

Definition at line 94 of file G4MaterialScanner.hh.

Referenced by G4MaterialScanner(), and ~G4MaterialScanner().

theRayShooter

G4RayShooter* G4MaterialScanner::theRayShooter = nullptr

private

Definition at line 93 of file G4MaterialScanner.hh.

Referenced by DoScan(), G4MaterialScanner(), and ~G4MaterialScanner().

theRegion

G4Region* G4MaterialScanner::theRegion = nullptr

private

Definition at line 120 of file G4MaterialScanner.hh.

Referenced by DoScan(), and SetRegionName().

thetaMin

G4double G4MaterialScanner::thetaMin = 0.0

private

Definition at line 110 of file G4MaterialScanner.hh.

Referenced by DoScan(), GetThetaMin(), and SetThetaMin().

thetaSpan

G4double G4MaterialScanner::thetaSpan = 0.0

private

Definition at line 111 of file G4MaterialScanner.hh.

Referenced by DoScan(), G4MaterialScanner(), GetThetaSpan(), and SetThetaSpan().

theUserEventAction

G4UserEventAction* G4MaterialScanner::theUserEventAction = nullptr

private

Definition at line 98 of file G4MaterialScanner.hh.

Referenced by RestoreUserActions(), and StoreUserActions().

theUserStackingAction

G4UserStackingAction* G4MaterialScanner::theUserStackingAction = nullptr

private

Definition at line 99 of file G4MaterialScanner.hh.

Referenced by RestoreUserActions(), and StoreUserActions().

theUserSteppingAction

G4UserSteppingAction* G4MaterialScanner::theUserSteppingAction = nullptr

private

Definition at line 101 of file G4MaterialScanner.hh.

Referenced by RestoreUserActions(), and StoreUserActions().

theUserTrackingAction

G4UserTrackingAction* G4MaterialScanner::theUserTrackingAction = nullptr

private

Definition at line 100 of file G4MaterialScanner.hh.

Referenced by RestoreUserActions(), and StoreUserActions().

---

The documentation for this class was generated from the following files:

• source/run/include/G4MaterialScanner.hh
• source/run/src/G4MaterialScanner.cc