Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm11/src/HistoManager.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/extended/electromagnetic/TestEm11/src/HistoManager.cc (Version 11.3.0) and /examples/extended/electromagnetic/TestEm11/src/HistoManager.cc (Version 10.3.p2)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 /// \file electromagnetic/TestEm11/src/HistoMa     26 /// \file electromagnetic/TestEm11/src/HistoManager.cc
 27 /// \brief Implementation of the HistoManager      27 /// \brief Implementation of the HistoManager class
 28 //                                                 28 //
 29 //                                                 29 //
 30 //                                             <<  30 // $Id: HistoManager.cc 98749 2016-08-09 13:43:36Z gcosmo $
                                                   >>  31 // 
 31 //....oooOO0OOooo........oooOO0OOooo........oo     32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32 //....oooOO0OOooo........oooOO0OOooo........oo     33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 33                                                    34 
 34 #include "HistoManager.hh"                         35 #include "HistoManager.hh"
 35                                                << 
 36 #include "DetectorConstruction.hh"                 36 #include "DetectorConstruction.hh"
 37                                                << 
 38 #include "G4UnitsTable.hh"                         37 #include "G4UnitsTable.hh"
 39                                                    38 
 40 //....oooOO0OOooo........oooOO0OOooo........oo     39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 41                                                    40 
 42 HistoManager::HistoManager()                       41 HistoManager::HistoManager()
                                                   >>  42   : fFileName("testem11")
 43 {                                                  43 {
 44   Book();                                          44   Book();
 45 }                                                  45 }
 46                                                    46 
 47 //....oooOO0OOooo........oooOO0OOooo........oo     47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 48                                                    48 
                                                   >>  49 HistoManager::~HistoManager()
                                                   >>  50 {
                                                   >>  51   delete G4AnalysisManager::Instance();
                                                   >>  52 }
                                                   >>  53 
                                                   >>  54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >>  55 
 49 void HistoManager::Book()                          56 void HistoManager::Book()
 50 {                                                  57 {
 51   // Create or get analysis manager                58   // Create or get analysis manager
 52   // The choice of analysis technology is done     59   // The choice of analysis technology is done via selection of a namespace
 53   // in HistoManager.hh                            60   // in HistoManager.hh
 54   G4AnalysisManager* analysisManager = G4Analy     61   G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
 55   analysisManager->SetDefaultFileType("root"); << 
 56   analysisManager->SetFileName(fFileName);         62   analysisManager->SetFileName(fFileName);
 57   analysisManager->SetVerboseLevel(1);             63   analysisManager->SetVerboseLevel(1);
 58   analysisManager->SetActivation(true);            64   analysisManager->SetActivation(true);
 59                                                    65 
 60   // Define histograms start values                66   // Define histograms start values
 61   const G4int kMaxHisto = 10;                      67   const G4int kMaxHisto = 10;
 62   const G4String id[] = {"0",  "1",  "2",  "3" <<  68   const G4String id[] = { "0", "1", "2", "3", "4", "5", "6", "7", "8", "9",
 63                          "12", "13", "14", "15 <<  69                          "10","11","12","13","14","15","16","17","18","19",
 64                                                <<  70                          "20","21","22"};
 65   const G4String title[] = {                   <<  71 
 66     "dummy",  // 0                             <<  72   const G4String title[] = 
 67     "Edep (MeV/mm) along absorbers",  // 1     <<  73                 { "dummy",                                        //0
 68     "total Energy deposited in absorbers",  // <<  74                   "Edep (MeV/mm) along absorbers",                //1
 69     "true track length of the primary particle <<  75                   "total Energy deposited in absorbers",          //2
 70     "true step size of the primary particle",  <<  76                   "true track length of the primary particle",    //3
 71     "projected range of the primary particle", <<  77                   "true step size of the primary particle",       //4
 72     "true track length of charged secondaries" <<  78                   "projected range of the primary particle",      //5
 73     "true step size of charged secondaries",   <<  79                   "true track length of charged secondaries",     //6
 74     "Edep (MeV.cm2/g) along x/r0",  // 8       <<  80                   "true step size of charged secondaries",        //7
 75     "total Energy leakage",  // 9              <<  81                   "Edep (MeV.cm2/g) along x/r0",                  //8
 76     "total Energy : Edep + Eleak"  // 10       <<  82                   "dummy",                                        //9
 77   };                                           <<  83                   "dummy"                                         //10
                                                   >>  84                  };
 78                                                    85 
 79   // Default values (to be reset via /analysis     86   // Default values (to be reset via /analysis/h1/set command)
 80   G4int nbins = 100;                               87   G4int nbins = 100;
 81   G4double vmin = 0.;                              88   G4double vmin = 0.;
 82   G4double vmax = 100.;                            89   G4double vmax = 100.;
 83                                                    90 
 84   // Create all histograms as inactivated      <<  91   // Create all histograms as inactivated 
 85   // as we have not yet set nbins, vmin, vmax      92   // as we have not yet set nbins, vmin, vmax
 86   for (G4int k = 0; k <= kMaxHisto; k++) {     <<  93   for (G4int k=0; k<=kMaxHisto; k++) {
 87     G4int ih = analysisManager->CreateH1(id[k]     94     G4int ih = analysisManager->CreateH1(id[k], title[k], nbins, vmin, vmax);
 88     analysisManager->SetH1Activation(ih, false     95     analysisManager->SetH1Activation(ih, false);
 89   }                                                96   }
 90                                                <<  97   
 91   G4String title2;                                 98   G4String title2;
 92   for (G4int k = 1; k < kMaxAbsor; k++) {      <<  99   for (G4int k=1; k<kMaxAbsor; k++) {
 93     title2 = "Edep in absorber " + id[k];         100     title2 = "Edep in absorber " + id[k];
 94     G4int ih = analysisManager->CreateH1(id[kM << 101     G4int ih 
                                                   >> 102       = analysisManager->CreateH1(id[kMaxHisto+k], title2, nbins, vmin, vmax);
 95     analysisManager->SetH1Activation(ih, false    103     analysisManager->SetH1Activation(ih, false);
 96   }                                               104   }
 97 }                                                 105 }
 98                                                   106 
 99 //....oooOO0OOooo........oooOO0OOooo........oo    107 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
100                                                   108