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 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 
33 #include "HistoManager.hh"
34 #include "G4UnitsTable.hh"
35 
36 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
37 
39  : fFileName("Activation")
40 {
41  Book();
42 }
43 
44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
45 
47 {
49 }
50 
51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
52 
53 void HistoManager::Book()
54 {
55  // Create or get analysis manager
56  // The choice of analysis technology is done via selection of a namespace
57  // in HistoManager.hh
59  analysisManager->SetFileName(fFileName);
60  analysisManager->SetVerboseLevel(1);
61  analysisManager->SetActivation(true); //enable inactivation of histograms
62 
63  // Define histograms start values
65  const G4String id[] = {"0","1","2","3","4","5","6","7","8","9",
66  "10","11","12","13",
67  "14","15","16","17","18","19","20","21","22","23",
68  "24","25","26","27","28","29","30","31","32","33",
69  "34","35","36","37","38","39","40","41","42","43" };
70 
71  const G4String title[] =
72  { "dummy", //0
73  "total energy deposit", //1
74  "Edep (MeV/mm) along beam directiom", //2
75  "total kinetic energy emerging", //3
76  "energy spectrum of emerging gamma", //4
77  "energy spectrum of emerging e+-", //5
78  "energy spectrum of emerging neutrons", //6
79  "energy spectrum of emerging protons", //7
80  "energy spectrum of emerging deuterons", //8
81  "energy spectrum of emerging alphas", //9
82  "energy spectrum of all others emerging ions", //10
83  "energy spectrum of all others emerging baryons", //11
84  "energy spectrum of all others emerging mesons", //12
85  "energy spectrum of all others emerging leptons (neutrinos)", //13
86  "dN/dt (becquerel) of emerging gamma", //14
87  "dN/dt (becquerel) of emerging e+- ", //15
88  "dN/dt (becquerel) of emerging neutrons", //16
89  "dN/dt (becquerel) of emerging protons", //17
90  "dN/dt (becquerel) of emerging deuterons", //18
91  "dN/dt (becquerel) of emerging alphas", //19
92  "dN/dt (becquerel) of all others emerging ions", //20
93  "dN/dt (becquerel) of all others emerging baryons", //21
94  "dN/dt (becquerel) of all others emerging mesons", //22
95  "dN/dt (becquerel) of all others emerging leptons (neutrinos)" //23
96  };
97 
98  // Default values (to be reset via /analysis/h1/set command)
99  G4int nbins = 100;
100  G4double vmin = 0.;
101  G4double vmax = 100.;
102 
103  // Create all histograms as inactivated
104  // as we have not yet set nbins, vmin, vmax
105  for (G4int k=0; k<kMaxHisto1; k++) {
106  G4int ih = analysisManager->CreateH1(id[k], title[k], nbins, vmin, vmax);
107  analysisManager->SetH1Activation(ih, false);
108  }
109 
110  // population of species
111  for (G4int k=kMaxHisto1; k<kMaxHisto2; k++) {
112  G4int ih = analysisManager->CreateH1(id[k], " population",nbins,vmin,vmax);
113  analysisManager->SetH1Activation(ih, false);
114  }
115 }
116 
117 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......