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 electromagnetic/TestEm5/src/HistoManager.cc
0027 /// \brief Implementation of the HistoManager class
0028 //
0029 //
0030 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0031 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0032
0033 #include "HistoManager.hh"
0034
0035 #include "G4UnitsTable.hh"
0036
0037 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0038
0039 HistoManager::HistoManager()
0040 {
0041   Book();
0042 }
0043
0044 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0045
0046 void HistoManager::Book()
0047 {
0048   // Create or get analysis manager
0049   // The choice of analysis technology is done via selection of a namespace
0050   // in HistoManager.hh
0051   G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
0052   analysisManager->SetDefaultFileType("root");
0053   analysisManager->SetFileName(fFileName);
0054   analysisManager->SetVerboseLevel(1);
0055   analysisManager->SetActivation(true);  // enable inactivation of histograms
0056
0057   // Define histograms start values
0058   const G4int kMaxHisto = 62;
0059   const G4String id[] = {
0060     "0",  "1",  "2",  "3",  "4",  "5",  "6",  "7",  "8",  "9",  "10", "11", "12", "13", "14", "15",
0061     "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31",
0062     "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47",
0063     "48", "49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60", "61"};
0064
0065   const G4String title[] = {
0066     "dummy",  // 0
0067     "energy deposit in absorber: dN/dE",  // 1
0068     "energy of charged secondaries at creation",  // 2
0069     "energy of neutral secondaries at creation",  // 3
0070     "energy of charged at creation (log scale)",  // 4
0071     "energy of neutral at creation (log scale)",  // 5
0072     "x_vertex of charged secondaries (all)",  // 6
0073     "x_vertex of charged secondaries (not absorbed)",  // 7
0074     "dummy",
0075     "dummy",  // 8-9
0076     "(transmit, charged) : kinetic energy at exit: dN/dE",  // 10
0077     "(transmit, charged) : ener fluence: dE(MeV)/dOmega",  // 11
0078     "(transmit, charged) : space angle: dN/dOmega",  // 12
0079     "(transmit, charged) : projected angle at exit",  // 13
0080     "(transmit, charged) : projected position at exit",  // 14
0081     "(transmit, charged) : radius at exit",  // 15
0082     "dummy",  // 16
0083     "dummy",  // 17
0084     "dummy",  // 18
0085     "dummy",  // 19
0086     "(transmit, neutral) : kinetic energy at exit",  // 20
0087     "(transmit, neutral) : ener fluence: dE(MeV)/dOmega",  // 21
0088     "(transmit, neutral) : space angle: dN/dOmega",  // 22
0089     "(transmit, neutral) : projected angle at exit",  // 23
0090     "dummy",
0091     "dummy",
0092     "dummy",
0093     "dummy",
0094     "dummy",
0095     "dummy",  // 24-29
0096     "(reflect , charged) : kinetic energy at exit",  // 30
0097     "(reflect , charged) : ener fluence: dE(MeV)/dOmega",  // 31
0098     "(reflect , charged) : space angle: dN/dOmega",  // 32
0099     "(reflect , charged) : projected angle at exit",  // 33
0100     "dummy",
0101     "dummy",
0102     "dummy",
0103     "dummy",
0104     "dummy",
0105     "dummy",  // 34-39
0106     "(reflect , neutral) : kinetic energy at exit",  // 40
0107     "(reflect , neutral) : ener fluence: dE(MeV)/dOmega",  // 41
0108     "(reflect , neutral) : space angle: dN/dOmega",  // 42
0109     "(reflect , neutral) : projected angle at exit",  // 43
0110     "dummy",  // 44
0111     "dummy",  // 45
0112     "dummy",  // 46
0113     "dummy",  // 47
0114     "dummy",  // 48
0115     "dummy",  // 49
0116     "energy of Auger e- at creation",  // 50
0117     "energy of fluorescence gamma at creation",  // 51
0118     "energy of Auger e- at creation (log scale)",  // 52
0119     "energy of fluorescence gamma at creation (log scale)",  // 53
0120     "energy of PIXE Auger e- at creation",  // 54
0121     "energy of PIXE gamma at creation",  // 55
0122     "energy of PIXE Auger e- at creation (log scale)",  // 56
0123     "energy of PIXE gamma at creation (log scale)",  // 57
0124     "energy of DNA Auger e- at creation",  // 58
0125     "energy of DNA gamma at creation",  // 59
0126     "energy of DNA Auger e- at creation (log scale)",  // 60
0127     "energy of DNA gamma at creation (log scale)"  // 61
0128   };
0129
0130   // Default values (to be reset via /analysis/h1/set command)
0131   G4int nbins = 100;
0132   G4double vmin = 0.;
0133   G4double vmax = 100.;
0134
0135   // Create all histograms as inactivated
0136   // as we have not yet set nbins, vmin, vmax
0137   for (G4int k = 0; k < kMaxHisto; k++) {
0138     G4int ih = analysisManager->CreateH1("h" + id[k], title[k], nbins, vmin, vmax);
0139     analysisManager->SetH1Activation(ih, false);
0140   }
0141 }