Hadr01/include/HistoManager.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 /// \file hadronic/Hadr01/include/HistoManager.hh
0027 /// \brief Definition of the HistoManager class
0028 //
0029 //---------------------------------------------------------------------------
0030 //
0031 // ClassName:   HistoManager
0032 //
0033 // Description: Singleton class to hold parameters and build histograms.
0034 //              User cannot access to the constructor.
0035 //              The pointer of the only existing object can be got via
0036 //              HistoManager::GetPointer() static method.
0037 //              The first invokation of this static method makes
0038 //              the singleton object.
0039 //
0040 // Author:      V.Ivanchenko 27/09/00
0041 //
0042 // Modified:
0043 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchenko)
0044 // 03.10.2006 Add csFlag (V.Ivanchenko)
0045 // 16.11.2006 Add beamFlag (V.Ivanchenko)
0046 //
0047 //----------------------------------------------------------------------------
0048 //
0049
0050 #ifndef HistoManager_h
0051 #define HistoManager_h 1
0052
0053 #include "G4Element.hh"
0054 #include "G4Material.hh"
0055 #include "globals.hh"
0056
0057 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
0058
0059 class Histo;
0060 class G4Track;
0061 class G4Step;
0062 class G4ParticleDefinition;
0063
0064 class HistoManager
0065 {
0066   public:
0067     static HistoManager* GetPointer();
0068
0069   private:
0070     HistoManager();
0071
0072   public:
0073     ~HistoManager();
0074
0075     void BookHisto();
0076
0077     void BeginOfRun();
0078     void EndOfRun();
0079     void BeginOfEvent();
0080     void EndOfEvent();
0081     void Fill(G4int id, G4double x, G4double w);
0082
0083     void ScoreNewTrack(const G4Track*);
0084     void AddTargetStep(const G4Step*);
0085     void AddLeakingParticle(const G4Track*);
0086
0087     void SetVerbose(G4int val);
0088
0089     inline void SetTargetRadius(G4double val)
0090     {
0091       fRadius = val;
0092       fR2 = val * val;
0093     };
0094     inline void SetTargetLength(G4double val) { fLength = val; };
0095     inline void SetNumberOfSlices(G4int val) { fNSlices = val; };
0096     inline void SetNumberOfBinsE(G4int val) { fNBinsE = val; };
0097     inline void SetDefaultBeamPositionFlag(G4bool f) { fBeamFlag = f; };
0098     inline void SetMaxEnergyDeposit(G4double val) { fEdepMax = val; };
0099
0100     inline G4double Radius() const { return fRadius; };
0101     inline G4double Length() const { return fLength; };
0102     inline G4bool DefaultBeamPosition() const { return fBeamFlag; };
0103     inline G4int NumberOfSlices() const { return fNSlices; };
0104     inline G4int GetVerbose() const { return fVerbose; };
0105
0106     inline G4int PrintBertiniXS() const { return fPrintBertiniXS; };
0107     inline void SetPrintBertiniXS(G4int key) { fPrintBertiniXS = key; };
0108
0109   private:
0110     static HistoManager* fManager;
0111
0112     const G4ParticleDefinition* fPrimaryDef;
0113     const G4ParticleDefinition* fNeutron;
0114
0115     G4double fR2 = 0.0;
0116     G4double fRadius;
0117     G4double fLength;
0118     G4double fEdepMax;
0119     G4double fEdepEvt = 0.0;
0120     G4double fEdepEM = 0.0;
0121     G4double fEdepPI = 0.0;
0122     G4double fEdepP = 0.0;
0123     G4double fEdepSum = 0.0;
0124     G4double fEdepSum2 = 0.0;
0125     G4double fAbsZ0 = 0.0;
0126     G4double fPrimaryKineticEnergy = 0.0;
0127
0128     G4int fVerbose = 0;
0129     G4int fNBinsE = 100;
0130     G4int fNSlices = 300;
0131
0132     G4int fNevt = 0;
0133     G4int fNelec = 0;
0134     G4int fNposit = 0;
0135     G4int fNgam = 0;
0136     G4int fNprot_leak = 0;
0137     G4int fNpiofNleak = 0;
0138     G4int fNcpions = 0;
0139     G4int fNpi0 = 0;
0140     G4int fNkaons = 0;
0141     G4int fNmuons = 0;
0142     G4int fNions = 0;
0143     G4int fNdeut = 0;
0144     G4int fNalpha = 0;
0145     G4int fNneutron = 0;
0146     G4int fNproton = 0;
0147     G4int fNaproton = 0;
0148     G4int fNneu_forw = 0;
0149     G4int fNneu_leak = 0;
0150     G4int fNneu_back = 0;
0151     G4int fNstep = 0;
0152     G4int fNHisto = 28;
0153     G4int fPrintBertiniXS = -1;  // 0 - all
0154
0155     G4bool fBeamFlag = true;
0156     G4bool fHistoBooked = false;
0157
0158     Histo* fHisto;
0159 };
0160
0161 #endif