B4b::RunAction Class Reference

Run action class. More...

#include <Doxymodules_basic.h>

Inheritance diagram for B4b::RunAction:

Public Member Functions
	RunAction ()
	~RunAction () override=default
G4Run *	GenerateRun () override
void	BeginOfRunAction (const G4Run *) override
void	EndOfRunAction (const G4Run *) override

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 according to a specified file extension.

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

Definition at line 142 of file Doxymodules_basic.h.

Constructor & Destructor Documentation

RunAction()

B4b::RunAction::RunAction

(

)

Definition at line 44 of file RunAction.cc.

{
  // 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 Analysis.hh
  auto analysisManager = G4AnalysisManager::Instance();
 
  // 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", 110, 0., 330*MeV);
  analysisManager->CreateH1("Egap" ,"Edep in gap", 100, 0., 30*MeV);
  analysisManager->CreateH1("Labs" ,"trackL in absorber", 100, 0., 50*cm);
  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();
}

~RunAction()

B4b::RunAction::~RunAction ( )

overridedefault

Member Function Documentation

GenerateRun()

G4Run * B4b::RunAction::GenerateRun ( )

override

Definition at line 82 of file RunAction.cc.

{
  return (new RunData);
}

BeginOfRunAction()

void B4b::RunAction::BeginOfRunAction ( const G4Run *  run )

override

Definition at line 89 of file RunAction.cc.

{
  G4cout << "### Run " << run->GetRunID() << " start." << G4endl;
 
  //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.root";
  // Other supported output types:
  // G4String fileName = "B4.csv";
  // G4String fileName = "B4.hdf5";
  // G4String fileName = "B4.xml";
  analysisManager->OpenFile(fileName);
  G4cout << "Using " << analysisManager->GetType() << G4endl;
}

EndOfRunAction()

void B4b::RunAction::EndOfRunAction ( const G4Run *  )

override

Definition at line 112 of file RunAction.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();
 
}

---

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

• .doxygen/Doxymodules_basic.h
• basic/B4/B4b/include/RunAction.hh
• basic/B4/B4b/src/RunAction.cc