| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 /// \file ExGflashHistoManager.cc
28 /// \brief Implementation of the ExGflasHistoM
29 //
30 //....oooOO0OOooo........oooOO0OOooo........oo
31 //....oooOO0OOooo........oooOO0OOooo........oo
32
33 #include "ExGflashHistoManager.hh"
34
35 #include "ExGflashDetectorConstruction.hh"
36
37 #include "G4SystemOfUnits.hh"
38 #include "G4UnitsTable.hh"
39
40 //....oooOO0OOooo........oooOO0OOooo........oo
41
42 ExGflashHistoManager::ExGflashHistoManager(ExG
43 : fFileName("gflash01"), fDet(myDet)
44 {
45 Book();
46 }
47
48 //....oooOO0OOooo........oooOO0OOooo........oo
49
50 ExGflashHistoManager::~ExGflashHistoManager()
51
52 //....oooOO0OOooo........oooOO0OOooo........oo
53
54 void ExGflashHistoManager::Book()
55 {
56 G4AnalysisManager* fAnalysisManager = G4Anal
57 fAnalysisManager->SetDefaultFileType("root")
58 fAnalysisManager->SetFileName(fFileName);
59 fAnalysisManager->SetVerboseLevel(1);
60 fAnalysisManager->SetActivation(true); // e
61 //
62 // Creating histograms
63 //
64 G4int nLbin = fDet->GetnLtot();
65 G4int nRbin = fDet->GetnRtot();
66 G4double dLradl = fDet->GetdLradl();
67 G4double dRradl = fDet->GetdRradl();
68 fVerbose = fDet->GetVerbose();
69
70 fAnalysisManager->CreateH1("h0", "total ener
71
72 fAnalysisManager->CreateH1("h1", "The number
73
74 fAnalysisManager->CreateH1("h2", "The energy
75
76 fAnalysisManager->CreateP1("p0", "longit ene
77 0., 2000.);
78
79 fAnalysisManager->CreateP1("p1", "radial ene
80 0., 2000.);
81
82 fAnalysisManager->CreateP1("p2", "Comul long
83 nLbin * dLradl, 0
84
85 fAnalysisManager->CreateP1("p3", "Cuml radia
86 nRbin * dRradl, 0
87
88 if (fVerbose > 1) G4cout << "\n----> Histogr
89 }
90
91 //....oooOO0OOooo........oooOO0OOooo........oo
92
93 void ExGflashHistoManager::SetBinning()
94 {
95 G4int nLbin = fDet->GetnLtot();
96 G4int nRbin = fDet->GetnRtot();
97 G4double dLradl = fDet->GetdLradl();
98 G4double dRradl = fDet->GetdRradl();
99
100 // Get analysis manager
101 auto analysisManager = G4AnalysisManager::In
102
103 // Redifine profiles
104 analysisManager->SetP1(0, nLbin, 0., nLbin *
105 analysisManager->SetP1(1, nRbin, 0., nRbin *
106 analysisManager->SetP1(2, nLbin, 0., nLbin *
107 analysisManager->SetP1(3, nRbin, 0., nRbin *
108 }
109