| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 /// \file electromagnetic/TestEm8/src/RunActio 26 /// \file electromagnetic/TestEm8/src/RunAction.cc
27 /// \brief Implementation of the RunAction cla 27 /// \brief Implementation of the RunAction class
28 // 28 //
>> 29 // $Id: RunAction.cc 67268 2013-02-13 11:38:40Z ihrivnac $
29 // 30 //
30 //-------------------------------------------- 31 //---------------------------------------------------------------------------
31 // 32 //
32 // ClassName: RunAction 33 // ClassName: RunAction
33 // 34 //
34 // Author: V.Ivanchenko 01.09.2010 35 // Author: V.Ivanchenko 01.09.2010
35 // 36 //
36 //-------------------------------------------- 37 //----------------------------------------------------------------------------
37 // 38 //
38 39
39 //....oooOO0OOooo........oooOO0OOooo........oo << 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
40 //....oooOO0OOooo........oooOO0OOooo........oo << 41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
41 42
42 #include "RunAction.hh" 43 #include "RunAction.hh"
>> 44 #include "HistoManager.hh"
>> 45 #include "G4UImanager.hh"
>> 46 #include "G4VVisManager.hh"
43 47
44 #include "Run.hh" << 48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
45 #include "TestParameters.hh" <<
46 <<
47 #include "G4Run.hh" <<
48 #include "G4RunManager.hh" <<
49 #include "G4SystemOfUnits.hh" <<
50 #include "Randomize.hh" <<
51 <<
52 //....oooOO0OOooo........oooOO0OOooo........oo <<
53 49
54 RunAction::RunAction() 50 RunAction::RunAction()
55 { << 51 : G4UserRunAction()
56 TestParameters::GetPointer(); << 52 {}
57 fAnalysisManager = G4AnalysisManager::Instan <<
58 fAnalysisManager->SetDefaultFileType("root") <<
59 fAnalysisManager->SetFileName("testem8"); <<
60 fAnalysisManager->SetVerboseLevel(1); <<
61 fAnalysisManager->SetActivation(true); <<
62 } <<
63 <<
64 //....oooOO0OOooo........oooOO0OOooo........oo <<
65 <<
66 void RunAction::Book() <<
67 { <<
68 // Always creating analysis manager <<
69 TestParameters* param = TestParameters::GetP <<
70 G4int nBinsE = param->GetNumberBins(); <<
71 G4int nBinsCluster = param->GetNumberBinsClu <<
72 G4int nMaxCluster = param->GetMaxCluster(); <<
73 G4double maxEnergy = param->GetMaxEnergy(); <<
74 G4double factorALICE = param->GetFactorALICE <<
75 <<
76 // G4cout << "### maxenergy(keV)= " << maxEn <<
77 // << " factorALICE= " << factorAL <<
78 <<
79 // Creating an 1-dimensional histograms in t <<
80 fAnalysisManager->SetFirstHistoId(1); <<
81 fAnalysisManager->CreateH1("h1", "Energy dep <<
82 maxEnergy / keV); <<
83 fAnalysisManager->CreateH1("h2", "Number of <<
84 G4double(nMaxClus <<
85 fAnalysisManager->CreateH1("h3", "Energy dep <<
86 maxEnergy * facto <<
87 fAnalysisManager->OpenFile(); <<
88 } <<
89 53
90 //....oooOO0OOooo........oooOO0OOooo........oo << 54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
91 55
92 G4Run* RunAction::GenerateRun() << 56 RunAction::~RunAction()
93 { << 57 {}
94 fRun = new Run(); <<
95 return fRun; <<
96 } <<
97 58
98 //....oooOO0OOooo........oooOO0OOooo........oo << 59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
99 60
100 void RunAction::BeginOfRunAction(const G4Run* 61 void RunAction::BeginOfRunAction(const G4Run* aRun)
101 { 62 {
102 G4int id = aRun->GetRunID(); 63 G4int id = aRun->GetRunID();
103 G4cout << "### Run " << id << " start analys << 64 G4cout << "### Run " << id << " start" << G4endl;
>> 65 (HistoManager::GetPointer())->BeginOfRun();
>> 66
>> 67 #ifdef G4VIS_USE
>> 68 G4UImanager* UI = G4UImanager::GetUIpointer();
104 69
105 fRun->BeginOfRun(); << 70 G4VVisManager* pVVisManager = G4VVisManager::GetConcreteInstance();
>> 71
>> 72 if(pVVisManager)
>> 73 {
>> 74 UI->ApplyCommand("/vis/scene/notifyHandlers");
>> 75 }
>> 76 #endif
106 77
107 // histograms <<
108 Book(); <<
109 } 78 }
110 79
111 //....oooOO0OOooo........oooOO0OOooo........oo << 80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
112 81
113 void RunAction::EndOfRunAction(const G4Run*) 82 void RunAction::EndOfRunAction(const G4Run*)
114 { 83 {
115 // print Run summary << 84 G4cout << "RunAction: End of run actions are started" << G4endl;
116 G4cout << "RunAction: End of run actions are << 85 (HistoManager::GetPointer())->EndOfRun();
117 << " Nevt= " << fRun->GetNumberOfEv << 86
118 << G4endl; << 87 #ifdef G4VIS_USE
119 if (isMaster) { << 88 if (G4VVisManager::GetConcreteInstance()) {
120 fRun->EndOfRun(); << 89 G4UImanager::GetUIpointer()->ApplyCommand("/vis/viewer/update");
121 } <<
122 // save histos and close analysis <<
123 if (fAnalysisManager->IsActive()) { <<
124 fAnalysisManager->Write(); <<
125 fAnalysisManager->CloseFile(); <<
126 fAnalysisManager->Clear(); <<
127 } 90 }
>> 91 #endif
128 } 92 }
129 93
130 //....oooOO0OOooo........oooOO0OOooo........oo << 94 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
131 95