TestEm5/src/HistoManager.cc

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 /// \file HistoManager.cc
0027 /// \brief Implementation of the HistoManager class
0028
0029 #include "HistoManager.hh"
0030
0031 #include "G4UnitsTable.hh"
0032
0033 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0034
0035 HistoManager::HistoManager()
0036 {
0037   Book();
0038 }
0039
0040 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0041
0042 void HistoManager::Book()
0043 {
0044   // Create or get analysis manager
0045   G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
0046   analysisManager->SetDefaultFileType("root");
0047   analysisManager->SetFileName(fFileName);
0048   analysisManager->SetVerboseLevel(1);
0049   analysisManager->SetActivation(true);  // enable inactivation of histograms
0050
0051   // Define histograms start values
0052   const G4int kMaxHisto = 62;
0053   const G4String id[] = {
0054     "0",  "1",  "2",  "3",  "4",  "5",  "6",  "7",  "8",  "9",  "10", "11", "12", "13", "14", "15",
0055     "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31",
0056     "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47",
0057     "48", "49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60", "61"};
0058
0059   const G4String title[] = {
0060     "dummy",  // 0
0061     "energy deposit in absorber: dN/dE",  // 1
0062     "energy of charged secondaries at creation",  // 2
0063     "energy of neutral secondaries at creation",  // 3
0064     "energy of charged at creation (log scale)",  // 4
0065     "energy of neutral at creation (log scale)",  // 5
0066     "x_vertex of charged secondaries (all)",  // 6
0067     "x_vertex of charged secondaries (not absorbed)",  // 7
0068     "dummy",
0069     "dummy",  // 8-9
0070     "(transmit, charged) : kinetic energy at exit: dN/dE",  // 10
0071     "(transmit, charged) : ener fluence: dE(MeV)/dOmega",  // 11
0072     "(transmit, charged) : space angle: dN/dOmega",  // 12
0073     "(transmit, charged) : projected angle at exit",  // 13
0074     "(transmit, charged) : projected position at exit",  // 14
0075     "(transmit, charged) : radius at exit",  // 15
0076     "dummy",  // 16
0077     "dummy",  // 17
0078     "dummy",  // 18
0079     "dummy",  // 19
0080     "(transmit, neutral) : kinetic energy at exit",  // 20
0081     "(transmit, neutral) : ener fluence: dE(MeV)/dOmega",  // 21
0082     "(transmit, neutral) : space angle: dN/dOmega",  // 22
0083     "(transmit, neutral) : projected angle at exit",  // 23
0084     "dummy",
0085     "dummy",
0086     "dummy",
0087     "dummy",
0088     "dummy",
0089     "dummy",  // 24-29
0090     "(reflect , charged) : kinetic energy at exit",  // 30
0091     "(reflect , charged) : ener fluence: dE(MeV)/dOmega",  // 31
0092     "(reflect , charged) : space angle: dN/dOmega",  // 32
0093     "(reflect , charged) : projected angle at exit",  // 33
0094     "dummy",
0095     "dummy",
0096     "dummy",
0097     "dummy",
0098     "dummy",
0099     "dummy",  // 34-39
0100     "(reflect , neutral) : kinetic energy at exit",  // 40
0101     "(reflect , neutral) : ener fluence: dE(MeV)/dOmega",  // 41
0102     "(reflect , neutral) : space angle: dN/dOmega",  // 42
0103     "(reflect , neutral) : projected angle at exit",  // 43
0104     "dummy",  // 44
0105     "dummy",  // 45
0106     "dummy",  // 46
0107     "dummy",  // 47
0108     "dummy",  // 48
0109     "dummy",  // 49
0110     "energy of Auger e- at creation",  // 50
0111     "energy of fluorescence gamma at creation",  // 51
0112     "energy of Auger e- at creation (log scale)",  // 52
0113     "energy of fluorescence gamma at creation (log scale)",  // 53
0114     "energy of PIXE Auger e- at creation",  // 54
0115     "energy of PIXE gamma at creation",  // 55
0116     "energy of PIXE Auger e- at creation (log scale)",  // 56
0117     "energy of PIXE gamma at creation (log scale)",  // 57
0118     "energy of DNA Auger e- at creation",  // 58
0119     "energy of DNA gamma at creation",  // 59
0120     "energy of DNA Auger e- at creation (log scale)",  // 60
0121     "energy of DNA gamma at creation (log scale)"  // 61
0122   };
0123
0124   // Default values (to be reset via /analysis/h1/set command)
0125   G4int nbins = 100;
0126   G4double vmin = 0.;
0127   G4double vmax = 100.;
0128
0129   // Create all histograms as inactivated
0130   // as we have not yet set nbins, vmin, vmax
0131   for (G4int k = 0; k < kMaxHisto; k++) {
0132     G4int ih = analysisManager->CreateH1("h" + id[k], title[k], nbins, vmin, vmax);
0133     analysisManager->SetH1Activation(ih, false);
0134   }
0135 }