Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm5/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/TestEm5/src/HistoManager.cc (Version 11.3.0) and /examples/extended/electromagnetic/TestEm5/src/HistoManager.cc (Version 10.4.p1)


  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/TestEm5/src/HistoMan     26 /// \file electromagnetic/TestEm5/src/HistoManager.cc
 27 /// \brief Implementation of the HistoManager      27 /// \brief Implementation of the HistoManager class
 28 //                                                 28 //
 29 //                                                 29 //
                                                   >>  30 // $Id: HistoManager.cc 104417 2017-05-30 08:30:48Z gcosmo $
                                                   >>  31 // 
 30 //....oooOO0OOooo........oooOO0OOooo........oo     32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 31 //....oooOO0OOooo........oooOO0OOooo........oo     33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32                                                    34 
 33 #include "HistoManager.hh"                         35 #include "HistoManager.hh"
 34                                                << 
 35 #include "G4UnitsTable.hh"                         36 #include "G4UnitsTable.hh"
 36                                                    37 
 37 //....oooOO0OOooo........oooOO0OOooo........oo     38 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 38                                                    39 
 39 HistoManager::HistoManager()                       40 HistoManager::HistoManager()
                                                   >>  41   : fFileName("testem5")
 40 {                                                  42 {
 41   Book();                                          43   Book();
 42 }                                                  44 }
 43                                                    45 
 44 //....oooOO0OOooo........oooOO0OOooo........oo     46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 45                                                    47 
                                                   >>  48 HistoManager::~HistoManager()
                                                   >>  49 {
                                                   >>  50   delete G4AnalysisManager::Instance();
                                                   >>  51 }
                                                   >>  52 
                                                   >>  53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >>  54 
 46 void HistoManager::Book()                          55 void HistoManager::Book()
 47 {                                                  56 {
 48   // Create or get analysis manager                57   // Create or get analysis manager
 49   // The choice of analysis technology is done     58   // The choice of analysis technology is done via selection of a namespace
 50   // in HistoManager.hh                            59   // in HistoManager.hh
 51   G4AnalysisManager* analysisManager = G4Analy     60   G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
 52   analysisManager->SetDefaultFileType("root"); << 
 53   analysisManager->SetFileName(fFileName);         61   analysisManager->SetFileName(fFileName);
 54   analysisManager->SetVerboseLevel(1);             62   analysisManager->SetVerboseLevel(1);
 55   analysisManager->SetActivation(true);  // en <<  63   analysisManager->SetActivation(true);    // enable inactivation of histograms
 56                                                    64 
 57   // Define histograms start values                65   // Define histograms start values
 58   const G4int kMaxHisto = 62;                      66   const G4int kMaxHisto = 62;
 59   const G4String id[] = {                      <<  67   const G4String id[] = { "0", "1", "2", "3", "4", "5", "6", "7", "8", "9",
 60     "0",  "1",  "2",  "3",  "4",  "5",  "6",   <<  68                          "10","11","12","13","14","15","16","17","18","19",
 61     "16", "17", "18", "19", "20", "21", "22",  <<  69                          "20","21","22","23","24","25","26","27","28","29",
 62     "32", "33", "34", "35", "36", "37", "38",  <<  70                          "30","31","32","33","34","35","36","37","38","39",
 63     "48", "49", "50", "51", "52", "53", "54",  <<  71                          "40","41","42","43","44","45","46","47","48","49",
 64                                                <<  72                          "50","51","52","53","54","55","56","57","58","59",
 65   const G4String title[] = {                   <<  73                          "60","61"
 66     "dummy",  // 0                             <<  74                         };
 67     "energy deposit in absorber: dN/dE",  // 1 <<  75                         
 68     "energy of charged secondaries at creation <<  76   const G4String title[] =
 69     "energy of neutral secondaries at creation <<  77                 { "dummy",                                                //0
 70     "energy of charged at creation (log scale) <<  78                   "energy deposit in absorber: dN/dE",                    //1
 71     "energy of neutral at creation (log scale) <<  79                   "energy of charged secondaries at creation",            //2
 72     "x_vertex of charged secondaries (all)",   <<  80                   "energy of neutral secondaries at creation",            //3
 73     "x_vertex of charged secondaries (not abso <<  81                   "energy of charged at creation (log scale)",            //4
 74     "dummy",                                   <<  82                   "energy of neutral at creation (log scale)",            //5
 75     "dummy",  // 8-9                           <<  83                   "x_vertex of charged secondaries (all)",                //6
 76     "(transmit, charged) : kinetic energy at e <<  84                   "x_vertex of charged secondaries (not absorbed)",       //7
 77     "(transmit, charged) : ener fluence: dE(Me <<  85                   "dummy","dummy",                                        //8-9
 78     "(transmit, charged) : space angle: dN/dOm <<  86                   "(transmit, charged) : kinetic energy at exit: dN/dE",  //10
 79     "(transmit, charged) : projected angle at  <<  87                   "(transmit, charged) : ener fluence: dE(MeV)/dOmega",   //11
 80     "(transmit, charged) : projected position  <<  88                   "(transmit, charged) : space angle: dN/dOmega",         //12
 81     "(transmit, charged) : radius at exit",  / <<  89                   "(transmit, charged) : projected angle at exit",        //13
 82     "dummy",  // 16                            <<  90                   "(transmit, charged) : projected position at exit",     //14
 83     "dummy",  // 17                            <<  91                   "(transmit, charged) : radius at exit",                 //15
 84     "dummy",  // 18                            <<  92                   "dummy",                                                //16
 85     "dummy",  // 19                            <<  93                   "dummy",                                                //17
 86     "(transmit, neutral) : kinetic energy at e <<  94                   "dummy",                                                //18
 87     "(transmit, neutral) : ener fluence: dE(Me <<  95                   "dummy",                                                //19
 88     "(transmit, neutral) : space angle: dN/dOm <<  96                   "(transmit, neutral) : kinetic energy at exit",         //20
 89     "(transmit, neutral) : projected angle at  <<  97                   "(transmit, neutral) : ener fluence: dE(MeV)/dOmega",   //21
 90     "dummy",                                   <<  98                   "(transmit, neutral) : space angle: dN/dOmega",         //22
 91     "dummy",                                   <<  99                   "(transmit, neutral) : projected angle at exit",        //23
 92     "dummy",                                   << 100                   "dummy","dummy","dummy","dummy","dummy","dummy",       //24-29
 93     "dummy",                                   << 101                   "(reflect , charged) : kinetic energy at exit",         //30
 94     "dummy",                                   << 102                   "(reflect , charged) : ener fluence: dE(MeV)/dOmega",   //31
 95     "dummy",  // 24-29                         << 103                   "(reflect , charged) : space angle: dN/dOmega",         //32
 96     "(reflect , charged) : kinetic energy at e << 104                   "(reflect , charged) : projected angle at exit",        //33
 97     "(reflect , charged) : ener fluence: dE(Me << 105                   "dummy","dummy","dummy","dummy","dummy","dummy",       //34-39
 98     "(reflect , charged) : space angle: dN/dOm << 106                   "(reflect , neutral) : kinetic energy at exit",         //40
 99     "(reflect , charged) : projected angle at  << 107                   "(reflect , neutral) : ener fluence: dE(MeV)/dOmega",   //41
100     "dummy",                                   << 108                   "(reflect , neutral) : space angle: dN/dOmega",         //42
101     "dummy",                                   << 109                   "(reflect , neutral) : projected angle at exit",        //43
102     "dummy",                                   << 110                   "dummy",                                                //44
103     "dummy",                                   << 111                   "dummy",                                                //45
104     "dummy",                                   << 112                   "dummy",                                                //46
105     "dummy",  // 34-39                         << 113                   "dummy",                                                //47
106     "(reflect , neutral) : kinetic energy at e << 114                   "dummy",                                                //48
107     "(reflect , neutral) : ener fluence: dE(Me << 115                   "dummy",                                                //49
108     "(reflect , neutral) : space angle: dN/dOm << 116                   "energy of Auger e- at creation",                       //50
109     "(reflect , neutral) : projected angle at  << 117                   "energy of fluorescence gamma at creation",             //51
110     "dummy",  // 44                            << 118                   "energy of Auger e- at creation (log scale)",           //52
111     "dummy",  // 45                            << 119                   "energy of fluorescence gamma at creation (log scale)", //53
112     "dummy",  // 46                            << 120                   "energy of PIXE Auger e- at creation",                  //54
113     "dummy",  // 47                            << 121                   "energy of PIXE gamma at creation",                     //55
114     "dummy",  // 48                            << 122                   "energy of PIXE Auger e- at creation (log scale)",      //56
115     "dummy",  // 49                            << 123                   "energy of PIXE gamma at creation (log scale)",         //57
116     "energy of Auger e- at creation",  // 50   << 124                   "energy of DNA Auger e- at creation",                   //58
117     "energy of fluorescence gamma at creation" << 125                   "energy of DNA gamma at creation",                      //59
118     "energy of Auger e- at creation (log scale << 126                   "energy of DNA Auger e- at creation (log scale)",       //60
119     "energy of fluorescence gamma at creation  << 127                   "energy of DNA gamma at creation (log scale)"           //61
120     "energy of PIXE Auger e- at creation",  // << 128                  };
121     "energy of PIXE gamma at creation",  // 55 << 
122     "energy of PIXE Auger e- at creation (log  << 
123     "energy of PIXE gamma at creation (log sca << 
124     "energy of DNA Auger e- at creation",  //  << 
125     "energy of DNA gamma at creation",  // 59  << 
126     "energy of DNA Auger e- at creation (log s << 
127     "energy of DNA gamma at creation (log scal << 
128   };                                           << 
129                                                   129 
130   // Default values (to be reset via /analysis << 130   // Default values (to be reset via /analysis/h1/set command)               
131   G4int nbins = 100;                              131   G4int nbins = 100;
132   G4double vmin = 0.;                             132   G4double vmin = 0.;
133   G4double vmax = 100.;                           133   G4double vmax = 100.;
134                                                   134 
135   // Create all histograms as inactivated      << 135   // Create all histograms as inactivated 
136   // as we have not yet set nbins, vmin, vmax     136   // as we have not yet set nbins, vmin, vmax
137   for (G4int k = 0; k < kMaxHisto; k++) {      << 137   for (G4int k=0; k<kMaxHisto; k++) {
138     G4int ih = analysisManager->CreateH1("h" + << 138     G4int ih = analysisManager->CreateH1("h"+id[k], title[k], nbins,vmin,vmax);
139     analysisManager->SetH1Activation(ih, false    139     analysisManager->SetH1Activation(ih, false);
140   }                                               140   }
141 }                                                 141 }
142                                                   142