Geant4 Cross Reference |
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 hadronic/Hadr01/include/Histo.hh << 26 // $Id: Histo.hh,v 1.4 2006/08/11 14:38:11 vnivanch Exp $ 27 /// \brief Definition of the Histo class << 27 // GEANT4 tag $Name: geant4-08-03-patch-01 $ 28 // << 29 28 30 #ifndef Histo_h 29 #ifndef Histo_h 31 #define Histo_h 1 30 #define Histo_h 1 32 31 33 //-------------------------------------------- 32 //--------------------------------------------------------------------------- 34 // 33 // 35 // ClassName: Histo 34 // ClassName: Histo 36 // 35 // 37 // Description: Utility class to hold and mani 36 // Description: Utility class to hold and manipulate histograms/nTuples 38 // 37 // 39 // Author: V.Ivanchenko 30/10/03 38 // Author: V.Ivanchenko 30/10/03 40 // 39 // 41 //-------------------------------------------- 40 //---------------------------------------------------------------------------- 42 // 41 // 43 42 44 //....oooOO0OOooo........oooOO0OOooo........oo 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 45 44 46 #include "G4DataVector.hh" << 47 #include "globals.hh" 45 #include "globals.hh" 48 << 49 #include <vector> 46 #include <vector> >> 47 #include "G4DynamicParticle.hh" >> 48 #include "G4VPhysicalVolume.hh" >> 49 #include "G4DataVector.hh" >> 50 #include "G4Track.hh" >> 51 50 52 51 //....oooOO0OOooo........oooOO0OOooo........oo 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 52 54 53 class G4RootAnalysisManager; << 54 class HistoMessenger; 55 class HistoMessenger; 55 56 >> 57 namespace AIDA { >> 58 class ITree; >> 59 class ITuple; >> 60 class IHistogram1D; >> 61 } >> 62 56 class Histo 63 class Histo 57 { 64 { 58 public: << 59 Histo(); << 60 65 61 ~Histo(); << 66 public: >> 67 Histo(G4int ver); >> 68 >> 69 ~Histo(); >> 70 >> 71 void book(); >> 72 // Book predefined histogramms >> 73 >> 74 void save(); >> 75 // Save histogramms to file >> 76 >> 77 void reset(); >> 78 // Reset histogramms >> 79 >> 80 void add1D(const G4String&, const G4String&, G4int nb=100, G4double x1=0., >> 81 G4double x2=1., G4double u=1.); >> 82 // In this method histogramms are predefined >> 83 >> 84 void setHisto1D(G4int, G4int, G4double, G4double, G4double); >> 85 // It change bins and boundaries >> 86 >> 87 void fill(G4int, G4double, G4double); >> 88 // Histogramms are filled >> 89 >> 90 void scale(G4int, G4double); >> 91 >> 92 void addTuple(const G4String&, const G4String&, const G4String&); >> 93 // In this method nTuple is booked >> 94 >> 95 void fillTuple(const G4String&, G4double); >> 96 // Fill nTuple parameter >> 97 >> 98 void addRow(); >> 99 // Save tuple event >> 100 >> 101 void setFileName(const G4String& nam) {histName = nam;}; >> 102 >> 103 void setFileOption(const G4String& nam) {option = nam;}; >> 104 >> 105 void setFileType(const G4String& nam); >> 106 >> 107 void print(G4int i); >> 108 >> 109 void setVerbose(G4int val) {verbose = val;}; >> 110 >> 111 private: 62 112 63 // Book predefined histogramms << 113 G4String histName; 64 void Book(); << 114 G4String histType; >> 115 G4String option; >> 116 G4String tupleName; >> 117 G4String tupleId; >> 118 G4String tupleList; >> 119 G4int nHisto; >> 120 G4int verbose; >> 121 G4int defaultAct; 65 122 66 // Save histogramms to file << 123 std::vector<AIDA::IHistogram1D*> histo; 67 void Save(); << 124 AIDA::ITuple* ntup; 68 << 125 AIDA::ITree* tree; 69 // In this method 1-D histogramms are pred << 126 HistoMessenger* messenger; 70 void Add1D(const G4String&, const G4String << 127 std::vector<G4int> active; 71 G4double u = 1.); << 128 std::vector<G4int> bins; 72 << 129 std::vector<G4double> xmin; 73 // It change bins and boundaries << 130 std::vector<G4double> xmax; 74 void SetHisto1D(G4int, G4int, G4double, G4 << 131 std::vector<G4double> unit; 75 << 132 std::vector<G4String> ids; 76 // Histogram activation/deactivation << 133 std::vector<G4String> tittles; 77 void Activate(G4int, G4bool); << 78 << 79 // Histogramms are filled << 80 void Fill(G4int, G4double, G4double); << 81 << 82 // Histogramms are scaled << 83 void ScaleH1(G4int, G4double); << 84 << 85 // In this method nTuple is booked << 86 void AddTuple(const G4String&); << 87 << 88 // In this method nTuple is booked << 89 void AddTupleI(const G4String&); << 90 void AddTupleF(const G4String&); << 91 void AddTupleD(const G4String&); << 92 << 93 // Fill nTuple parameter << 94 void FillTupleI(G4int, G4int); << 95 void FillTupleF(G4int, G4float); << 96 void FillTupleD(G4int, G4double); << 97 << 98 // Save tuple event << 99 void AddRow(); << 100 << 101 // Set output file << 102 void SetFileName(const G4String&); << 103 << 104 inline void SetVerbose(G4int val) { fVerbo << 105 << 106 inline G4bool IsActive() const { return fH << 107 << 108 private: << 109 G4RootAnalysisManager* fManager; << 110 HistoMessenger* fMessenger; << 111 << 112 G4String fHistName; << 113 G4String fTupleName; << 114 G4String fTupleTitle; << 115 G4int fNHisto; << 116 G4int fVerbose; << 117 G4bool fDefaultAct; << 118 G4bool fHistoActive; << 119 G4bool fNtupleActive; << 120 << 121 std::vector<G4int> fHisto; << 122 std::vector<G4int> fTupleI; << 123 std::vector<G4int> fTupleF; << 124 std::vector<G4int> fTupleD; << 125 std::vector<G4int> fBins; << 126 std::vector<G4bool> fActive; << 127 std::vector<G4double> fXmin; << 128 std::vector<G4double> fXmax; << 129 std::vector<G4double> fUnit; << 130 std::vector<G4String> fIds; << 131 std::vector<G4String> fTitles; << 132 std::vector<G4String> fNtupleI; << 133 std::vector<G4String> fNtupleF; << 134 std::vector<G4String> fNtupleD; << 135 }; 134 }; 136 135 137 #endif 136 #endif 138 137