OMSimTools.hh

Go to the documentation of this file.

 
#pragma once
 
#include "OMSimLogger.hh"
#include "OMSimDetectorComponent.hh"
#include "OMSimCommandArgsTable.hh"
#include <G4LogicalSkinSurface.hh>
 
#include "globals.hh"
#include <TGraph.h>
#include <TH2D.h>
#include <variant>
#include <optional>
 
#include <G4ThreeVector.hh>
#include <G4RotationMatrix.hh>
#include <G4LogicalVolume.hh>
#include <CADMesh.hh>
#include <G4MultiUnion.hh>
 
namespace Tools
{
    std::vector<double> arange(double start, double stop, double step);
    TGraph* create1dInterpolator(const std::vector<double>& pX, const std::vector<double>& pY, const std::string& pName);
    TGraph* create1dInterpolator(const std::string& pFileName);
    TH2D *create2DHistogramFromDataFile(const std::string &pFilePath);
 
    void ensureDirectoryExists(const std::string &filepath);
    std::pair<std::vector<double>, std::vector<double>> histogram(
        const std::vector<double> &data,
        const std::variant<int, std::vector<double>> &bins = 10,
        const std::optional<std::pair<double, double>> &range = std::nullopt,
        const std::vector<double> &weights = std::vector<double>());
 
    G4String getThreadIDStr();
    std::vector<std::vector<double>> loadtxt(const std::string &pFilePath,
                                             bool pUnpack = true,
                                             size_t pSkipRows = 0,
                                             char pDelimiter = ' ',
                                             char pComments = '#');
    std::vector<double> linspace(double start, double end, int num_points, bool endpoint = true);
    std::vector<double> logspace(double start, double end, int num_points, double base = 10.0, bool endpoint = true);
    void sortVectorByReference(std::vector<G4double> &referenceVector, std::vector<G4double> &sortVector);
    double median(std::vector<double> p_vec);
    double mean(const std::vector<double> &p_vec, const std::vector<double> &p_weights = {});
    double std(const std::vector<double>& vec,  const std::vector<double> &p_weights = {});
    
    void throwError(const G4String& message);
    std::vector<G4String> splitStringByDelimiter(G4String const &p_string, char p_delim);
    std::vector<G4String> splitStringByDelimiter(char *p_char, char p_delim);
    extern std::string visualisationURL;
 
    void AppendCADComponent(OMSimDetectorComponent* component,
                            const G4double scaleCAD,
                            const G4ThreeVector& offsetCAD,
                            const G4RotationMatrix& rotation,
                            const std::string& filename,
                            const std::string& componentname,
                            G4Material* material,
                            G4VisAttributes visAttributes,
                            G4OpticalSurface* opticalSurface = nullptr); 
};