ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
HistoManager.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file HistoManager.cc
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 //
28 //
29 //
30 //
31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 
34 #include "HistoManager.hh"
35 #include "G4UnitsTable.hh"
36 
37 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
38 
39 HistoManager::HistoManager()
40  : fFileName("testem5")
41 {
42  Book();
43 }
44 
45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
46 
47 HistoManager::~HistoManager()
48 {
49  delete G4AnalysisManager::Instance();
50 }
51 
52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
53 
54 void HistoManager::Book()
55 {
56  // Create or get analysis manager
57  // The choice of analysis technology is done via selection of a namespace
58  // in HistoManager.hh
59  G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
60  analysisManager->SetFileName(fFileName);
61  analysisManager->SetVerboseLevel(1);
62  analysisManager->SetActivation(true); // enable inactivation of histograms
63 
64  // Define histograms start values
65  const G4int kMaxHisto = 62;
66  const G4String id[] = { "0", "1", "2", "3", "4", "5", "6", "7", "8", "9",
67  "10","11","12","13","14","15","16","17","18","19",
68  "20","21","22","23","24","25","26","27","28","29",
69  "30","31","32","33","34","35","36","37","38","39",
70  "40","41","42","43","44","45","46","47","48","49",
71  "50","51","52","53","54","55","56","57","58","59",
72  "60","61"
73  };
74 
75  const G4String title[] =
76  { "dummy", //0
77  "energy deposit in absorber: dN/dE", //1
78  "energy of charged secondaries at creation", //2
79  "energy of neutral secondaries at creation", //3
80  "energy of charged at creation (log scale)", //4
81  "energy of neutral at creation (log scale)", //5
82  "x_vertex of charged secondaries (all)", //6
83  "x_vertex of charged secondaries (not absorbed)", //7
84  "dummy","dummy", //8-9
85  "(transmit, charged) : kinetic energy at exit: dN/dE", //10
86  "(transmit, charged) : ener fluence: dE(MeV)/dOmega", //11
87  "(transmit, charged) : space angle: dN/dOmega", //12
88  "(transmit, charged) : projected angle at exit", //13
89  "(transmit, charged) : projected position at exit", //14
90  "(transmit, charged) : radius at exit", //15
91  "dummy", //16
92  "dummy", //17
93  "dummy", //18
94  "dummy", //19
95  "(transmit, neutral) : kinetic energy at exit", //20
96  "(transmit, neutral) : ener fluence: dE(MeV)/dOmega", //21
97  "(transmit, neutral) : space angle: dN/dOmega", //22
98  "(transmit, neutral) : projected angle at exit", //23
99  "dummy","dummy","dummy","dummy","dummy","dummy", //24-29
100  "(reflect , charged) : kinetic energy at exit", //30
101  "(reflect , charged) : ener fluence: dE(MeV)/dOmega", //31
102  "(reflect , charged) : space angle: dN/dOmega", //32
103  "(reflect , charged) : projected angle at exit", //33
104  "dummy","dummy","dummy","dummy","dummy","dummy", //34-39
105  "(reflect , neutral) : kinetic energy at exit", //40
106  "(reflect , neutral) : ener fluence: dE(MeV)/dOmega", //41
107  "(reflect , neutral) : space angle: dN/dOmega", //42
108  "(reflect , neutral) : projected angle at exit", //43
109  "dummy", //44
110  "dummy", //45
111  "dummy", //46
112  "dummy", //47
113  "dummy", //48
114  "dummy", //49
115  "energy of Auger e- at creation", //50
116  "energy of fluorescence gamma at creation", //51
117  "energy of Auger e- at creation (log scale)", //52
118  "energy of fluorescence gamma at creation (log scale)", //53
119  "energy of PIXE Auger e- at creation", //54
120  "energy of PIXE gamma at creation", //55
121  "energy of PIXE Auger e- at creation (log scale)", //56
122  "energy of PIXE gamma at creation (log scale)", //57
123  "energy of DNA Auger e- at creation", //58
124  "energy of DNA gamma at creation", //59
125  "energy of DNA Auger e- at creation (log scale)", //60
126  "energy of DNA gamma at creation (log scale)" //61
127  };
128 
129  // Default values (to be reset via /analysis/h1/set command)
130  G4int nbins = 100;
131  G4double vmin = 0.;
132  G4double vmax = 100.;
133 
134  // Create all histograms as inactivated
135  // as we have not yet set nbins, vmin, vmax
136  for (G4int k=0; k<kMaxHisto; k++) {
137  G4int ih = analysisManager->CreateH1("h"+id[k], title[k], nbins,vmin,vmax);
138  analysisManager->SetH1Activation(ih, false);
139  }
140 }