TRestGeant4AnalysisProcess.h

/*************************************************************************
 * This file is part of the REST software framework.                     *
 *                                                                       *
 * Copyright (C) 2016 GIFNA/TREX (University of Zaragoza)                *
 * For more information see http://gifna.unizar.es/trex                  *
 *                                                                       *
 * REST is free software: you can redistribute it and/or modify          *
 * it under the terms of the GNU General Public License as published by  *
 * the Free Software Foundation, either version 3 of the License, or     *
 * (at your option) any later version.                                   *
 *                                                                       *
 * REST is distributed in the hope that it will be useful,               *
 * but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the          *
 * GNU General Public License for more details.                          *
 *                                                                       *
 * You should have a copy of the GNU General Public License along with   *
 * REST in $REST_PATH/LICENSE.                                           *
 * If not, see http://www.gnu.org/licenses/.                             *
 * For the list of contributors see $REST_PATH/CREDITS.                  *
 *************************************************************************/
 
#ifndef RestCore_TRestGeant4AnalysisProcess
#define RestCore_TRestGeant4AnalysisProcess
 
#include <TRestGeant4Event.h>
#include <TRestGeant4Metadata.h>
 
#include "TRestEventProcess.h"
 
class TRestGeant4AnalysisProcess : public TRestEventProcess {
   private:
    TRestGeant4Event* fInputG4Event;  
 
    TRestGeant4Event* fOutputG4Event;  
 
    TRestGeant4Metadata* fG4Metadata;  
 
    std::vector<std::string> fEnergyInObservables;  
 
    std::vector<Int_t> fVolumeID;  
 
    std::vector<std::string> fVolumeName;  
 
    std::vector<std::string> fMeanPosObservables;  
 
    std::vector<Int_t> fVolumeID2;  
 
    std::vector<std::string> fDirID;  
 
    std::vector<std::string> fProcessObservables;  
 
    std::vector<Int_t> fVolumeID3;  
 
    std::vector<std::string> fProcessVolumeName;  
 
    std::vector<std::string> fProcessName;  
 
    std::vector<std::string> fContainsProcessObservables;  
 
    std::vector<std::string> fContainsProcessG4Name;  
 
    std::vector<std::string> fPerProcessObservables;  
 
    std::vector<std::vector<std::string>> fPerProcessG4Names;  
 
    std::vector<std::string> fPerProcessParticle;  
 
    std::vector<std::string> fTrackCounterObservables;  
 
    std::vector<std::string> fParticleTrackCounter;  
 
    std::vector<std::string> fTracksEDepObservables;  
 
    std::vector<std::string> fParticleTrackEdep;  
 
    Bool_t fPerProcessSensitiveEnergy = false;
    Bool_t fPerProcessSensitiveEnergyNorm = false;
 
    void Initialize() override;
 
    void LoadDefaultConfig();
 
    void PrintNotActiveVolumeWarning(const TString& volumeName);
 
   protected:
    std::map<std::string, std::string> GetAliasObservableNameToVolume();
 
   public:
    static const std::map<std::string, std::string>& GetRestToGeant4ProcessNameMap();
    std::string GetGeant4ProcessName(const std::string& name) const;
 
    RESTValue GetInputEvent() const override { return fInputG4Event; }
    RESTValue GetOutputEvent() const override { return fOutputG4Event; }
 
    void InitProcess() override;
    TRestEvent* ProcessEvent(TRestEvent* inputEvent) override;
    void EndProcess() override;
 
    void LoadConfig(const std::string& configFilename, const std::string& name = "");
 
    void PrintMetadata() override {
        BeginPrintProcess();
 
        EndPrintProcess();
    }
 
    TRestEventProcess* Maker() { return new TRestGeant4AnalysisProcess; }
 
    const char* GetProcessName() const override { return "Geant4Analysis"; }
 
    TRestGeant4AnalysisProcess();
    TRestGeant4AnalysisProcess(const char* configFilename);
    ~TRestGeant4AnalysisProcess();
 
    ClassDefOverride(TRestGeant4AnalysisProcess, 3);
};
#endif