Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/parameterisations/Par03/src/Par03RunAction.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.                             *
 10 // *                                                                  *
 11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                                                  *
 18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // ********************************************************************
 25 //
 26 #include "Par03RunAction.hh"
 27 
 28 #include "Par03DetectorConstruction.hh"
 29 
 30 #include "G4AnalysisManager.hh"
 31 
 32 Par03RunAction::Par03RunAction(Par03DetectorConstruction* aDetector)
 33   : G4UserRunAction(), fDetector(aDetector)
 34 {
 35   // Create analysis manager
 36   G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
 37   analysisManager->SetDefaultFileType("root");
 38 
 39   // Default filename, can be overriden with /analysis/setFileName
 40   analysisManager->SetFileName("Par03Output");
 41 }
 42 
 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 44 
 45 Par03RunAction::~Par03RunAction() = default;
 46 
 47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 48 
 49 void Par03RunAction::BeginOfRunAction(const G4Run*)
 50 {
 51   // Get detector dimensions
 52   G4int cellNumZ = fDetector->GetNbOfLayers();
 53   G4int cellNumRho = fDetector->GetNbOfRhoCells();
 54   G4double cellSizeZ = fDetector->GetLength() / cellNumZ;
 55   G4double cellSizeRho = fDetector->GetRadius() / cellNumRho;
 56   // Default max value of energy stored in histogram (in GeV)
 57   G4double maxEnergy = 100;
 58 
 59   // Get analysis manager
 60   auto analysisManager = G4AnalysisManager::Instance();
 61 
 62   // Creating control histograms
 63   analysisManager->CreateH1("energyParticle", "Primary energy;E_{MC} (GeV);Entries", 256, 0,
 64                             1.1 * maxEnergy);
 65   analysisManager->CreateH1("energyDeposited", "Deposited energy;E_{MC} (GeV);Entries", 256, 0,
 66                             1.1 * maxEnergy);
 67   analysisManager->CreateH1(
 68     "energyRatio", "Ratio of energy deposited to primary;E_{dep} /  E_{MC};Entries", 1024, 0, 1);
 69   analysisManager->CreateH1("time", "Simulation time; time (s);Entries", 2048, 0, 30);
 70   analysisManager->CreateH1("longProfile", "Longitudinal profile;t (mm);#LTE#GT (MeV)", cellNumZ,
 71                             -0.5 * cellSizeZ, (cellNumZ - 0.5) * cellSizeZ);
 72   analysisManager->CreateH1("transProfile", "Transverse profile;r (mm);#LTE#GT (MeV)", cellNumRho,
 73                             -0.5 * cellSizeRho, (cellNumRho - 0.5) * cellSizeRho);
 74   analysisManager->CreateH1("longFirstMoment",
 75                             "First moment of longitudinal distribution;#LT#lambda#GT (mm);Entries",
 76                             1024, -0.5 * cellSizeZ,
 77                             cellNumZ * cellSizeZ / 2);  // arbitrary scaling of max value on axis
 78   analysisManager->CreateH1("transFirstMoment",
 79                             "First moment of transverse distribution;#LTr#GT "
 80                             "(mm);Entries",
 81                             1024, -0.5 * cellSizeRho,
 82                             cellNumRho * cellSizeRho
 83                               / 10);  // arbitrary scaling of max value on axis
 84   analysisManager->CreateH1(
 85     "longSecondMoment",
 86     "Second moment of longitudinal distribution;#LT#lambda^{2}#GT "
 87     "(mm^{2});Entries",
 88     1024, 0, std::pow(cellNumZ * cellSizeZ, 2) / 25);  // arbitrary scaling of max value on axis
 89   analysisManager->CreateH1(
 90     "transSecondMoment", "Second moment of transverse distribution;#LTr^{2}#GT (mm^{2});Entries",
 91     1024, 0, std::pow(cellNumRho * cellSizeRho, 2) / 25);  // arbitrary scaling of max value on axis
 92   analysisManager->CreateH1("hitType", "hit type;type (0=full, 1= fast);Entries", 2, -0.5, 1.5);
 93 
 94   // Open an output file
 95   analysisManager->OpenFile();
 96 }
 97 
 98 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 99 
100 void Par03RunAction::EndOfRunAction(const G4Run*)
101 {
102   auto analysisManager = G4AnalysisManager::Instance();
103   analysisManager->Write();
104   analysisManager->CloseFile();
105   analysisManager->Clear();
106 }