Geant4  10.00.p01
UltraRunAction.cc
Go to the documentation of this file.
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 //
27 // --------------------------------------------------------------
28 // GEANT 4 - ULTRA experiment example
29 // --------------------------------------------------------------
30 //
31 // Code developed by:
32 // B. Tome, M.C. Espirito-Santo, A. Trindade, P. Rodrigues
33 //
34 // ****************************************************
35 // * UltraRunAction.cc
36 // ****************************************************
37 //
38 // RunAction class for Ultra; it has also a Messenger Class
39 //
40 #include "UltraRunAction.hh"
41 #include "G4Run.hh"
42 #include "G4RunManager.hh"
43 #include "Randomize.hh"
44 #include "UltraAnalysisManager.hh"
45 #include <ctime>
46 #include "Randomize.hh"
47 
48 
49 // //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
50 
52 {
53  seed = -1; // RANLUX seed
54  luxury = 3; // RANLUX luxury level (3 is default)
55  saveRndm = 1;
56 }
57 
58 // //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
59 
61 {;}
62 
63 // //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
64 
66 {
67  // Get/create analysis manager: need to do that in the master and in the workers
68  G4AnalysisManager* man = G4AnalysisManager::Instance();
69 
70  // Open an output file
71  man->OpenFile("ultra");
72  man->SetFirstHistoId(1);
73 
74  // Create histogram(s)
75  man->CreateH1("1","Optical photons energy (eV)", //histoID,histo name
76  500,0.,5.); //bins' number, xmin, xmax
77  man->CreateH1("2","Number of Detected Photons",
78  10,0.,10.); //bins' number, xmin, xmax
79 
80  if (!IsMaster()) //it is a slave, do nothing else
81  {
82  G4cout << "ooo Run " << aRun->GetRunID() << " starts on slave." << G4endl;
83  return;
84  }
85 
86  //Master or sequential
87  G4cout << "ooo Run " << aRun->GetRunID() << " starts (global)." << G4endl;
88  if (seed<0) //not initialized by anybody else
89  {
90  seed=time(0);
91  G4Random::setTheSeed(seed,luxury);
92  G4Random::showEngineStatus();
93  }
94 
95  // save Rndm status
96  if (saveRndm > 0)
97  G4Random::saveEngineStatus("beginOfRun.rndm");
98 
99  return;
100 }
101 
102 // //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
103 
105 {
106  // Write histograms to file
107 
108 
109  // Save histograms
110  G4AnalysisManager* man = G4AnalysisManager::Instance();
111  man->Write();
112  man->CloseFile();
113  // Complete clean-up
114  delete G4AnalysisManager::Instance();
115 
116  if (!IsMaster())
117  {
118  G4cout << "### Run " << aRun->GetRunID() << " (slave) ended." << G4endl;
119  return;
120  }
121 
122  G4cout << "### Run " << aRun->GetRunID() << " (global) ended." << G4endl;
123  // save Rndm status
124  if (saveRndm == 1)
125  G4Random::saveEngineStatus("endOfRun.rndm");
126  return;
127 }
128 
129 
G4bool SetFirstHistoId(G4int firstId)
G4int CreateH1(const G4String &name, const G4String &title, G4int nbins, G4double xmin, G4double xmax, const G4String &unitName="none", const G4String &fcnName="none", const G4String &binSchemeName="linear")
void BeginOfRunAction(const G4Run *aRun)
G4GLOB_DLL std::ostream G4cout
G4bool IsMaster() const
G4int GetRunID() const
Definition: G4Run.hh:76
Definition: G4Run.hh:46
void EndOfRunAction(const G4Run *aRun)
#define G4endl
Definition: G4ios.hh:61