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.6 2007-05-24 13:52:31 vnivanch Exp $ 27 /// \brief Definition of the Histo class << 27 // GEANT4 tag $Name: not supported by cvs2svn $ 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 class IAnalysisFactory; >> 62 } >> 63 56 class Histo 64 class Histo 57 { 65 { 58 public: << 59 Histo(); << 60 66 61 ~Histo(); << 67 public: >> 68 Histo(G4int ver); >> 69 >> 70 ~Histo(); >> 71 >> 72 void book(); >> 73 // Book predefined histogramms >> 74 >> 75 void save(); >> 76 // Save histogramms to file >> 77 >> 78 void reset(); >> 79 // Reset histogramms >> 80 >> 81 void add1D(const G4String&, const G4String&, G4int nb=100, G4double x1=0., >> 82 G4double x2=1., G4double u=1.); >> 83 // In this method histogramms are predefined >> 84 >> 85 void setHisto1D(G4int, G4int, G4double, G4double, G4double); >> 86 // It change bins and boundaries >> 87 >> 88 void fill(G4int, G4double, G4double); >> 89 // Histogramms are filled >> 90 >> 91 void scale(G4int, G4double); >> 92 >> 93 void addTuple(const G4String&, const G4String&, const G4String&); >> 94 // In this method nTuple is booked >> 95 >> 96 void fillTuple(const G4String&, G4double); >> 97 // Fill nTuple parameter >> 98 >> 99 void addRow(); >> 100 // Save tuple event >> 101 >> 102 void setFileName(const G4String& nam) {histName = nam;}; >> 103 >> 104 void setFileOption(const G4String& nam) {option = nam;}; >> 105 >> 106 void setFileType(const G4String& nam); >> 107 >> 108 void print(G4int i); >> 109 >> 110 void setVerbose(G4int val) {verbose = val;}; >> 111 >> 112 private: 62 113 63 // Book predefined histogramms << 114 G4String histName; 64 void Book(); << 115 G4String histType; >> 116 G4String option; >> 117 G4String tupleName; >> 118 G4String tupleId; >> 119 G4String tupleList; >> 120 G4int nHisto; >> 121 G4int verbose; >> 122 G4int defaultAct; 65 123 66 // Save histogramms to file << 124 std::vector<AIDA::IHistogram1D*> histo; 67 void Save(); << 125 AIDA::IAnalysisFactory* af; 68 << 126 AIDA::ITuple* ntup; 69 // In this method 1-D histogramms are pred << 127 AIDA::ITree* tree; 70 void Add1D(const G4String&, const G4String << 128 HistoMessenger* messenger; 71 G4double u = 1.); << 129 std::vector<G4int> active; 72 << 130 std::vector<G4int> bins; 73 // It change bins and boundaries << 131 std::vector<G4double> xmin; 74 void SetHisto1D(G4int, G4int, G4double, G4 << 132 std::vector<G4double> xmax; 75 << 133 std::vector<G4double> unit; 76 // Histogram activation/deactivation << 134 std::vector<G4String> ids; 77 void Activate(G4int, G4bool); << 135 std::vector<G4String> titles; 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 }; 136 }; 136 137 137 #endif 138 #endif 138 139