Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/analysis/AnaEx02/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/analysis/AnaEx02/src/HistoManager.cc (Version 11.3.0) and /examples/extended/analysis/AnaEx02/src/HistoManager.cc (Version 11.1.2)


  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 analysis/AnaEx02/src/HistoManager.cc     26 /// \file analysis/AnaEx02/src/HistoManager.cc
 27 /// \brief Implementation of the HistoManager      27 /// \brief Implementation of the HistoManager class
 28 //                                                 28 //
 29 //                                                 29 //
 30 //....oooOO0OOooo........oooOO0OOooo........oo     30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 31 //....oooOO0OOooo........oooOO0OOooo........oo     31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32                                                    32 
 33 #include "HistoManager.hh"                     <<  33 #include <TH1D.h>
                                                   >>  34 #include <TFile.h>
                                                   >>  35 #include <TTree.h>
                                                   >>  36 #include <CLHEP/Units/SystemOfUnits.h>
 34                                                    37 
                                                   >>  38 #include "HistoManager.hh"
 35 #include "G4UnitsTable.hh"                         39 #include "G4UnitsTable.hh"
 36                                                    40 
 37 #include <CLHEP/Units/SystemOfUnits.h>         <<  41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 38 #include <TFile.h>                             <<  42 
 39 #include <TH1D.h>                              <<  43 HistoManager::HistoManager()
 40 #include <TTree.h>                             <<  44 {
                                                   >>  45 
                                                   >>  46   // // histograms
                                                   >>  47   // for (G4int k=0; k<kMaxHisto; k++) fHisto[k] = 0;
                                                   >>  48 }
 41                                                    49 
 42 //....oooOO0OOooo........oooOO0OOooo........oo     50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 43                                                    51 
 44 HistoManager::~HistoManager()                      52 HistoManager::~HistoManager()
 45 {                                                  53 {
 46   delete fRootFile;                                54   delete fRootFile;
 47 }                                                  55 }
 48                                                    56 
 49 //....oooOO0OOooo........oooOO0OOooo........oo     57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 50                                                    58 
 51 void HistoManager::Book()                          59 void HistoManager::Book()
 52 {                                                  60 {
 53   // Creating a tree container to handle histo     61   // Creating a tree container to handle histograms and ntuples.
 54   // This tree is associated to an output file     62   // This tree is associated to an output file.
 55   //                                               63   //
 56   G4String fileName = "AnaEx02.root";              64   G4String fileName = "AnaEx02.root";
 57   fRootFile = new TFile(fileName, "RECREATE"); <<  65   fRootFile = new TFile(fileName,"RECREATE");
 58   if (fRootFile == nullptr) {                      66   if (fRootFile == nullptr) {
 59     G4cout << " HistoManager::Book :"              67     G4cout << " HistoManager::Book :"
 60            << " problem creating the ROOT TFil <<  68            << " problem creating the ROOT TFile "
                                                   >>  69            << G4endl;
 61     return;                                        70     return;
 62   }                                                71   }
 63                                                    72 
 64   // id = 0                                        73   // id = 0
 65   fHisto[0] = new TH1D("EAbs", "Edep in absorb <<  74   fHisto[0] = new TH1D("EAbs", "Edep in absorber (MeV)", 100, 0., 800*CLHEP::MeV);
 66   // id = 1                                        75   // id = 1
 67   fHisto[1] = new TH1D("EGap", "Edep in gap (M <<  76   fHisto[1] = new TH1D("EGap", "Edep in gap (MeV)", 100, 0., 100*CLHEP::MeV);
 68   // id = 2                                        77   // id = 2
 69   fHisto[2] = new TH1D("LAbs", "trackL in abso <<  78   fHisto[2] = new TH1D("LAbs", "trackL in absorber (mm)", 100, 0., 1*CLHEP::m);
 70   // id = 3                                        79   // id = 3
 71   fHisto[3] = new TH1D("LGap", "trackL in gap  <<  80   fHisto[3] = new TH1D("LGap", "trackL in gap (mm)", 100, 0., 50*CLHEP::cm);
 72                                                    81 
 73   for (G4int i = 0; i < kMaxHisto; ++i) {      <<  82   for ( G4int i=0; i<kMaxHisto; ++i ) {
 74     if (fHisto[i] == nullptr) {                    83     if (fHisto[i] == nullptr) {
 75       G4cout << "\n can't create histo " << i      84       G4cout << "\n can't create histo " << i << G4endl;
 76     }                                              85     }
 77   }                                                86   }
 78                                                    87 
 79   // create 1st ntuple                             88   // create 1st ntuple
 80   fNtuple1 = new TTree("Ntuple1", "Edep");         89   fNtuple1 = new TTree("Ntuple1", "Edep");
 81   fNtuple1->Branch("Eabs", &fEabs, "Eabs/D");      90   fNtuple1->Branch("Eabs", &fEabs, "Eabs/D");
 82   fNtuple1->Branch("Egap", &fEgap, "Egap/D");      91   fNtuple1->Branch("Egap", &fEgap, "Egap/D");
 83                                                    92 
 84   // create 2nd ntuple                             93   // create 2nd ntuple
 85   fNtuple2 = new TTree("Ntuple2", "TrackL");       94   fNtuple2 = new TTree("Ntuple2", "TrackL");
 86   fNtuple2->Branch("Labs", &fLabs, "Labs/D");      95   fNtuple2->Branch("Labs", &fLabs, "Labs/D");
 87   fNtuple2->Branch("Lgap", &fLgap, "Lgap/D");      96   fNtuple2->Branch("Lgap", &fLgap, "Lgap/D");
 88                                                    97 
 89   G4cout << "\n----> Output file is open in "      98   G4cout << "\n----> Output file is open in " << fileName << G4endl;
 90 }                                                  99 }
 91                                                   100 
 92 //....oooOO0OOooo........oooOO0OOooo........oo    101 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 93                                                   102 
 94 void HistoManager::Save()                         103 void HistoManager::Save()
 95 {                                                 104 {
 96   if (fRootFile == nullptr) {                  << 105   if (fRootFile == nullptr) { return;
 97     return;                                    << 106 }
 98   }                                            << 
 99                                                   107 
100   fRootFile->Write();  // Writing the histogra << 108   fRootFile->Write();       // Writing the histograms to the file
101   fRootFile->Close();  // and closing the tree << 109   fRootFile->Close();       // and closing the tree (and the file)
102                                                   110 
103   G4cout << "\n----> Histograms and ntuples ar    111   G4cout << "\n----> Histograms and ntuples are saved\n" << G4endl;
104 }                                                 112 }
105                                                   113 
106 //....oooOO0OOooo........oooOO0OOooo........oo    114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
107                                                   115 
108 void HistoManager::FillHisto(G4int ih, G4doubl    116 void HistoManager::FillHisto(G4int ih, G4double xbin, G4double weight)
109 {                                                 117 {
110   if (ih >= kMaxHisto) {                          118   if (ih >= kMaxHisto) {
111     G4cerr << "---> warning from HistoManager:    119     G4cerr << "---> warning from HistoManager::FillHisto() : histo " << ih
112            << " does not exist. (xbin=" << xbi << 120            << " does not exist. (xbin=" << xbin << " weight=" << weight << ")"
                                                   >> 121            << G4endl;
113     return;                                       122     return;
114   }                                               123   }
115   if (fHisto[ih] != nullptr) {                 << 124   if  (fHisto[ih] != nullptr) { fHisto[ih]->Fill(xbin, weight); }
116     fHisto[ih]->Fill(xbin, weight);            << 
117   }                                            << 
118 }                                                 125 }
119                                                   126 
120 //....oooOO0OOooo........oooOO0OOooo........oo    127 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
121                                                   128 
122 void HistoManager::Normalize(G4int ih, G4doubl    129 void HistoManager::Normalize(G4int ih, G4double fac)
123 {                                                 130 {
124   if (ih >= kMaxHisto) {                          131   if (ih >= kMaxHisto) {
125     G4cout << "---> warning from HistoManager:    132     G4cout << "---> warning from HistoManager::Normalize() : histo " << ih
126            << " does not exist. (fac=" << fac     133            << " does not exist. (fac=" << fac << ")" << G4endl;
127     return;                                       134     return;
128   }                                               135   }
129   if (fHisto[ih] != nullptr) {                    136   if (fHisto[ih] != nullptr) {
130     fHisto[ih]->Scale(fac);                       137     fHisto[ih]->Scale(fac);
131   }                                               138   }
132 }                                                 139 }
133                                                   140 
134 //....oooOO0OOooo........oooOO0OOooo........oo    141 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
135                                                   142 
136 void HistoManager::FillNtuple(G4double energyA << 143 void HistoManager::FillNtuple(G4double energyAbs, G4double energyGap,
137                               G4double trackLG << 144                               G4double trackLAbs , G4double trackLGap )
138 {                                                 145 {
139   fEabs = energyAbs;                              146   fEabs = energyAbs;
140   fEgap = energyGap;                              147   fEgap = energyGap;
141   fLabs = trackLAbs;                              148   fLabs = trackLAbs;
142   fLgap = trackLGap;                              149   fLgap = trackLGap;
143                                                   150 
144   if (fNtuple1 != nullptr) {                   << 151   if (fNtuple1 != nullptr) { fNtuple1->Fill(); }
145     fNtuple1->Fill();                          << 152   if (fNtuple2 != nullptr) { fNtuple2->Fill(); }
146   }                                            << 
147   if (fNtuple2 != nullptr) {                   << 
148     fNtuple2->Fill();                          << 
149   }                                            << 
150 }                                                 153 }
151                                                   154 
152 //....oooOO0OOooo........oooOO0OOooo........oo    155 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
153                                                   156 
154 void HistoManager::PrintStatistic()               157 void HistoManager::PrintStatistic()
155 {                                                 158 {
156   G4cout << "\n ----> print histograms statist    159   G4cout << "\n ----> print histograms statistic \n" << G4endl;
157   for (auto h1 : fHisto) {                        160   for (auto h1 : fHisto) {
158     const G4String name = h1->GetName();          161     const G4String name = h1->GetName();
159                                                   162 
160     G4String unitCategory;                        163     G4String unitCategory;
161     if (name[0] == 'E') {                      << 164     if (name[0] == 'E' ) { unitCategory = "Energy"; }
162       unitCategory = "Energy";                 << 165     if (name[0] == 'L' ) { unitCategory = "Length"; }
163     }                                          << 
164     if (name[0] == 'L') {                      << 
165       unitCategory = "Length";                 << 
166     }                                          << 
167                                                   166 
168     G4cout << name << ": mean = " << G4BestUni << 167     G4cout << name
169            << " rms = " << G4BestUnit(h1->GetR << 168            << ": mean = " << G4BestUnit(h1->GetMean(), unitCategory)
                                                   >> 169            << " rms = " << G4BestUnit(h1->GetRMS(), unitCategory )
                                                   >> 170            << G4endl;
170   }                                               171   }
171 }                                                 172 }
172                                                   173 
173 //....oooOO0OOooo........oooOO0OOooo........oo    174 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
174                                                   175