Files
2025-12-05 08:54:02 +01:00

159 lines
5.0 KiB
C++

//
// ********************************************************************
// * License and Disclaimer *
// * *
// * The Geant4 software is copyright of the Copyright Holders of *
// * the Geant4 Collaboration. It is provided under the terms and *
// * conditions of the Geant4 Software License, included in the file *
// * LICENSE and available at http://cern.ch/geant4/license . These *
// * include a list of copyright holders. *
// * *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work make any representation or warranty, express or implied, *
// * regarding this software system or assume any liability for its *
// * use. Please see the license in the file LICENSE and URL above *
// * for the full disclaimer and the limitation of liability. *
// * *
// * This code implementation is the result of the scientific and *
// * technical work of the GEANT4 collaboration. *
// * By using, copying, modifying or distributing the software (or *
// * any work based on the software) you agree to acknowledge its *
// * use in resulting scientific publications, and indicate your *
// * acceptance of all terms of the Geant4 Software license. *
// ********************************************************************
//
/// \file HistoManager.hh
/// \brief Definition of the HistoManager class
#ifndef HistoManager_h
#define HistoManager_h 1
//---------------------------------------------------------------------------
//
// ClassName: HistoManager
//
// Description: Singleton class to hold Emc parameters and build histograms.
// User cannot access to the constructor.
// The pointer of the only existing object can be got via
// HistoManager::GetPointer() static method.
// The first invokation of this static method makes
// the singleton object.
//
// Author: V.Ivanchenko 27/09/00
//
//----------------------------------------------------------------------------
//
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
#include "G4DataVector.hh"
#include "G4DynamicParticle.hh"
#include "G4Track.hh"
#include "G4VPhysicalVolume.hh"
#include "globals.hh"
#include <vector>
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
class Histo;
class HistoManager
{
public:
// With description
static HistoManager* GetPointer();
private:
HistoManager();
public: // Without description
~HistoManager();
void BookHisto();
void BeginOfRun();
void EndOfRun(G4int runID);
void BeginOfEvent();
void EndOfEvent();
void ScoreNewTrack(const G4Track* aTrack);
void AddEnergy(G4double edep, G4int idx, G4int copyNo);
void AddDeltaElectron(const G4DynamicParticle*);
void AddPhoton(const G4DynamicParticle*);
inline void ResetTrackLength() { fTrackLength = 0.0, fTrackAbs = true; };
inline void AddPositron(const G4DynamicParticle*) { ++fPosit; };
inline void SetVerbose(G4int val) { fVerbose = val; };
inline G4int GetVerbose() const { return fVerbose; };
inline void SetHistoNumber(G4int val) { fNHisto = val; };
inline void SetFirstEventToDebug(G4int val) { fEvt1 = val; };
inline G4int FirstEventToDebug() const { return fEvt1; };
inline void SetLastEventToDebug(G4int val) { fEvt2 = val; };
inline G4int LastEventToDebug() const { return fEvt2; };
inline void SetMaxEnergy(G4double val) { fMaxEnergy = val; };
inline G4double GetMaxEnergy() const { return fMaxEnergy; };
inline void AddStep() { fStep += 1.0; };
// Acceptance parameters
inline void SetBeamEnergy(G4double val) { fBeamEnergy = val; };
void SetEdepAndRMS(G4int, const G4ThreeVector&);
private:
// MEMBERS
static HistoManager* fManager;
const G4ParticleDefinition* fGamma;
const G4ParticleDefinition* fElectron;
const G4ParticleDefinition* fPositron;
Histo* fHisto;
G4int fNHisto;
G4int fVerbose;
G4int fEvt1;
G4int fEvt2;
G4double fBeamEnergy;
G4double fMaxEnergy;
G4double fMaxEnergyAbs;
G4double fTrackLength;
G4double fStep;
G4bool fTrackAbs; // Track is in absorber
G4int fEvt;
G4int fElec;
G4int fPosit;
G4int fGam;
G4int fLowe;
G4int fBinsE, fBinsEA, fBinsED;
G4double fEabs1, fEabs2, fEabs3, fEabs4;
G4double fE[25];
G4DataVector fEvertex;
G4DataVector fNvertex;
G4DataVector fBrem;
G4DataVector fPhot;
G4DataVector fComp;
G4DataVector fConv;
G4double fEdeptrue[3];
G4double fRmstrue[3];
G4double fLimittrue[3];
G4double fEdep[6];
G4double fErms[6];
G4double fEdeptr[6];
G4double fErmstr[6];
G4int fStat[6];
G4int fNmax;
};
#endif