Par03RunAction Class Reference

Run action. More...

#include <Doxymodules_parameterisations.h>

Inheritance diagram for Par03RunAction:

Public Member Functions
	Par03RunAction (Par03DetectorConstruction *aDetector)
	Constructor. Defines the histograms.
virtual	~Par03RunAction ()
virtual void	BeginOfRunAction (const G4Run *) final
	Open the file for the analysis.
virtual void	EndOfRunAction (const G4Run *) final
	Write and close the file.

Private Attributes
Par03DetectorConstruction *	fDetector
	Pointer to detector construction to retrieve the detector dimensions to setup the histograms.

Detailed Description

Run action.

Create analysis file and define control histograms for showers in detectors. Histograms are configured taking into account the dimensions of the detector (and number of readout cells).

Definition at line 67 of file Doxymodules_parameterisations.h.

Constructor & Destructor Documentation

Par03RunAction()

Par03RunAction::Par03RunAction

(

Par03DetectorConstruction *

aDetector

)

Constructor. Defines the histograms.

Definition at line 31 of file Par03RunAction.cc.

  : G4UserRunAction()
  , fDetector(aDetector)
{
  // Create analysis manager
  G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
  analysisManager->SetDefaultFileType("root");
 
  // Default filename, can be overriden with /analysis/setFileName
  analysisManager->SetFileName("Par03Output");
}

~Par03RunAction()

Par03RunAction::~Par03RunAction ( )

virtualdefault

Member Function Documentation

BeginOfRunAction()

void Par03RunAction::BeginOfRunAction ( const G4Run *  )

finalvirtual

Open the file for the analysis.

Definition at line 49 of file Par03RunAction.cc.

{
  // Get detector dimensions
  G4int cellNumZ       = fDetector->GetNbOfLayers();
  G4int cellNumRho     = fDetector->GetNbOfRhoCells();
  G4double cellSizeZ   = fDetector->GetLength() / cellNumZ;
  G4double cellSizeRho = fDetector->GetRadius() / cellNumRho;
  // Default max value of energy stored in histogram (in GeV)
  G4double maxEnergy = 100;
 
  // Get analysis manager
  auto analysisManager = G4AnalysisManager::Instance();
 
  // Creating control histograms
  analysisManager->CreateH1("energyParticle",
                            "Primary energy;E_{MC} (GeV);Entries", 256, 0,
                            1.1 * maxEnergy);
  analysisManager->CreateH1("energyDeposited",
                            "Deposited energy;E_{MC} (GeV);Entries", 256, 0,
                            1.1 * maxEnergy);
  analysisManager->CreateH1(
    "energyRatio",
    "Ratio of energy deposited to primary;E_{dep} /  E_{MC};Entries", 1024, 0,
    1);
  analysisManager->CreateH1("time", "Simulation time; time (s);Entries", 2048,
                            0, 30);
  analysisManager->CreateH1(
    "longProfile", "Longitudinal profile;t (mm);#LTE#GT (MeV)", cellNumZ,
    -0.5 * cellSizeZ, (cellNumZ - 0.5) * cellSizeZ);
  analysisManager->CreateH1(
    "transProfile", "Transverse profile;r (mm);#LTE#GT (MeV)", cellNumRho,
    -0.5 * cellSizeRho, (cellNumRho - 0.5) * cellSizeRho);
  analysisManager->CreateH1(
    "longFirstMoment",
    "First moment of longitudinal distribution;#LT#lambda#GT (mm);Entries",
    1024, -0.5 * cellSizeZ,
    cellNumZ * cellSizeZ / 2);  // arbitrary scaling of max value on axis
  analysisManager->CreateH1("transFirstMoment",
                            "First moment of transverse distribution;#LTr#GT "
                            "(mm);Entries",
                            1024, -0.5 * cellSizeRho,
                            cellNumRho * cellSizeRho /
                              10);  // arbitrary scaling of max value on axis
  analysisManager->CreateH1(
    "longSecondMoment",
    "Second moment of longitudinal distribution;#LT#lambda^{2}#GT "
    "(mm^{2});Entries",
    1024, 0,
    std::pow(cellNumZ * cellSizeZ, 2) /
      25);  // arbitrary scaling of max value on axis
  analysisManager->CreateH1(
    "transSecondMoment",
    "Second moment of transverse distribution;#LTr^{2}#GT (mm^{2});Entries",
    1024, 0,
    std::pow(cellNumRho * cellSizeRho, 2) /
      25);  // arbitrary scaling of max value on axis
  analysisManager->CreateH1(
    "hitType", "hit type;type (0=full, 1= fast);Entries", 2, -0.5, 1.5);
 
  // Open an output file
  analysisManager->OpenFile();
}

EndOfRunAction()

void Par03RunAction::EndOfRunAction ( const G4Run *  )

finalvirtual

Write and close the file.

Definition at line 114 of file Par03RunAction.cc.

{
  auto analysisManager = G4AnalysisManager::Instance();
  analysisManager->Write();
  analysisManager->CloseFile();
  analysisManager->Clear();
}

Member Data Documentation

fDetector

Par03DetectorConstruction* Par03RunAction::fDetector

private

Pointer to detector construction to retrieve the detector dimensions to setup the histograms.

Definition at line 57 of file Par03RunAction.hh.

---

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

• .doxygen/Doxymodules_parameterisations.h
• extended/parameterisations/Par03/include/Par03RunAction.hh
• extended/parameterisations/Par03/src/Par03RunAction.cc