Activation/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 kMaxHisto1 = 24, kMaxHisto2 = 44;
0053   const G4String id[] = {"0",  "1",  "2",  "3",  "4",  "5",  "6",  "7",  "8",  "9",  "10",
0054                          "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21",
0055                          "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32",
0056                          "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43"};
0057
0058   const G4String title[] = {
0059     "dummy",  // 0
0060     "total energy deposit",  // 1
0061     "Edep (MeV/mm) along beam directiom",  // 2
0062     "total kinetic energy emerging",  // 3
0063     "energy spectrum of emerging gamma",  // 4
0064     "energy spectrum of emerging e+-",  // 5
0065     "energy spectrum of emerging neutrons",  // 6
0066     "energy spectrum of emerging protons",  // 7
0067     "energy spectrum of emerging deuterons",  // 8
0068     "energy spectrum of emerging alphas",  // 9
0069     "energy spectrum of all others emerging ions",  // 10
0070     "energy spectrum of all others emerging baryons",  // 11
0071     "energy spectrum of all others emerging mesons",  // 12
0072     "energy spectrum of all others emerging leptons (neutrinos)",  // 13
0073     "dN/dt (becquerel) of emerging gamma",  // 14
0074     "dN/dt (becquerel) of emerging e+- ",  // 15
0075     "dN/dt (becquerel) of emerging neutrons",  // 16
0076     "dN/dt (becquerel) of emerging protons",  // 17
0077     "dN/dt (becquerel) of emerging deuterons",  // 18
0078     "dN/dt (becquerel) of emerging alphas",  // 19
0079     "dN/dt (becquerel) of all others emerging ions",  // 20
0080     "dN/dt (becquerel) of all others emerging baryons",  // 21
0081     "dN/dt (becquerel) of all others emerging mesons",  // 22
0082     "dN/dt (becquerel) of all others emerging leptons (neutrinos)"  // 23
0083   };
0084
0085   // Default values (to be reset via /analysis/h1/set command)
0086   G4int nbins = 100;
0087   G4double vmin = 0.;
0088   G4double vmax = 100.;
0089
0090   // Create all histograms as inactivated
0091   // as we have not yet set nbins, vmin, vmax
0092   for (G4int k = 0; k < kMaxHisto1; k++) {
0093     G4int ih = analysisManager->CreateH1(id[k], title[k], nbins, vmin, vmax);
0094     analysisManager->SetH1Activation(ih, false);
0095   }
0096
0097   // population of species
0098   for (G4int k = kMaxHisto1; k < kMaxHisto2; k++) {
0099     G4int ih = analysisManager->CreateH1(id[k], " population", nbins, vmin, vmax);
0100     analysisManager->SetH1Activation(ih, false);
0101   }
0102 }
0103
0104 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......