B4RunAction Class Reference

#include <B4RunAction.hh>

Inheritance diagram for B4RunAction:

Collaboration diagram for B4RunAction:

Public Member Functions
	B4RunAction ()
virtual	~B4RunAction ()
virtual void	BeginOfRunAction (const G4Run *)
virtual void	EndOfRunAction (const G4Run *)
	B4RunAction ()
virtual	~B4RunAction ()
virtual void	BeginOfRunAction (const G4Run *)
virtual void	EndOfRunAction (const G4Run *)
	B4RunAction ()
virtual	~B4RunAction ()
virtual void	BeginOfRunAction (const G4Run *)
virtual void	EndOfRunAction (const G4Run *)
Public Member Functions inherited from G4UserRunAction
	G4UserRunAction ()
virtual	~G4UserRunAction ()
virtual G4Run *	GenerateRun ()
virtual void	SetMaster (G4bool val=true)
G4bool	IsMaster () const

Additional Inherited Members
Protected Attributes inherited from G4UserRunAction
G4bool	isMaster

Detailed Description

Run action class

It accumulates statistic and computes dispersion of the energy deposit and track lengths of charged particles with use of analysis tools: H1D histograms are created in BeginOfRunAction() for the following physics quantities:

• Edep in absorber
• Edep in gap
• Track length in absorber
• Track length in gap The same values are also saved in the ntuple. The histograms and ntuple are saved in the output file in a format accoring to a selected technology in B4Analysis.hh.

In EndOfRunAction(), the accumulated statistic and computed dispersion is printed.

Definition at line 57 of file B4RunAction.hh.

Constructor & Destructor Documentation

B4RunAction::B4RunAction

(

)

Definition at line 41 of file B4RunAction.cc.

 : G4UserRunAction()
{ 
  // set printing event number per each event
  G4RunManager::GetRunManager()->SetPrintProgress(1);     

  // Create analysis manager
  // The choice of analysis technology is done via selectin of a namespace
  // in B4Analysis.hh
  auto analysisManager = G4AnalysisManager::Instance();
  G4cout << "Using " << analysisManager->GetType() << G4endl;

  // Create directories 
  //analysisManager->SetHistoDirectoryName("histograms");
  //analysisManager->SetNtupleDirectoryName("ntuple");
  analysisManager->SetVerboseLevel(1);
  analysisManager->SetNtupleMerging(true);
    // Note: merging ntuples is available only with Root output

  // Book histograms, ntuple
  //
  
  // Creating histograms
  analysisManager->CreateH1("Eabs","Edep in absorber", 100, 0., 800*MeV);
  analysisManager->CreateH1("Egap","Edep in gap", 100, 0., 100*MeV);
  analysisManager->CreateH1("Labs","trackL in absorber", 100, 0., 1*m);
  analysisManager->CreateH1("Lgap","trackL in gap", 100, 0., 50*cm);

  // Creating ntuple
  //
  analysisManager->CreateNtuple("B4", "Edep and TrackL");
  analysisManager->CreateNtupleDColumn("Eabs");
  analysisManager->CreateNtupleDColumn("Egap");
  analysisManager->CreateNtupleDColumn("Labs");
  analysisManager->CreateNtupleDColumn("Lgap");
  analysisManager->FinishNtuple();
}

Here is the call graph for this function:

B4RunAction::~B4RunAction ( )

virtual

Definition at line 81 of file B4RunAction.cc.

{
  delete G4AnalysisManager::Instance();  
}

B4RunAction::B4RunAction

(

)

virtual B4RunAction::~B4RunAction ( )

virtual

B4RunAction::B4RunAction

(

)

virtual B4RunAction::~B4RunAction ( )

virtual

Member Function Documentation

void B4RunAction::BeginOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

Definition at line 88 of file B4RunAction.cc.

{ 
  //inform the runManager to save random number seed
  //G4RunManager::GetRunManager()->SetRandomNumberStore(true);
  
  // Get analysis manager
  auto analysisManager = G4AnalysisManager::Instance();

  // Open an output file
  //
  G4String fileName = "B4";
  analysisManager->OpenFile(fileName);
}

virtual void B4RunAction::BeginOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

virtual void B4RunAction::BeginOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

virtual void B4RunAction::EndOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

void B4RunAction::EndOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

Definition at line 104 of file B4RunAction.cc.

{
  // print histogram statistics
  //
  auto analysisManager = G4AnalysisManager::Instance();
  if ( analysisManager->GetH1(1) ) {
    G4cout << G4endl << " ----> print histograms statistic ";
    if(isMaster) {
      G4cout << "for the entire run " << G4endl << G4endl; 
    }
    else {
      G4cout << "for the local thread " << G4endl << G4endl; 
    }
    
    G4cout << " EAbs : mean = " 
       << G4BestUnit(analysisManager->GetH1(0)->mean(), "Energy") 
       << " rms = " 
       << G4BestUnit(analysisManager->GetH1(0)->rms(),  "Energy") << G4endl;
    
    G4cout << " EGap : mean = " 
       << G4BestUnit(analysisManager->GetH1(1)->mean(), "Energy") 
       << " rms = " 
       << G4BestUnit(analysisManager->GetH1(1)->rms(),  "Energy") << G4endl;
    
    G4cout << " LAbs : mean = " 
      << G4BestUnit(analysisManager->GetH1(2)->mean(), "Length") 
      << " rms = " 
      << G4BestUnit(analysisManager->GetH1(2)->rms(),  "Length") << G4endl;

    G4cout << " LGap : mean = " 
      << G4BestUnit(analysisManager->GetH1(3)->mean(), "Length") 
      << " rms = " 
      << G4BestUnit(analysisManager->GetH1(3)->rms(),  "Length") << G4endl;
  }

  // save histograms & ntuple
  //
  analysisManager->Write();
  analysisManager->CloseFile();
}

virtual void B4RunAction::EndOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

---

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

• source/geant4.10.03.p03/examples/basic/B4/B4a/include/B4RunAction.hh
• source/geant4.10.03.p03/examples/basic/B4/B4c/include/B4RunAction.hh
• source/geant4.10.03.p03/examples/basic/B4/B4d/include/B4RunAction.hh
• source/geant4.10.03.p03/examples/basic/B4/B4a/src/B4RunAction.cc
• source/geant4.10.03.p03/examples/basic/B4/B4c/src/B4RunAction.cc
• source/geant4.10.03.p03/examples/basic/B4/B4d/src/B4RunAction.cc