TRestComponentDataSet.h

/*************************************************************************
 * This file is part of the REST software framework.                     *
 *                                                                       *
 * Copyright (C) 2016 GIFNA/TREX (University of Zaragoza)                *
 * For more information see https://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 https://www.gnu.org/licenses/.                            *
 * For the list of contributors see $REST_PATH/CREDITS.                  *
 *************************************************************************/

#ifndef REST_TRestComponentDataSet
#define REST_TRestComponentDataSet

#include <THn.h>

#include "TRestComponent.h"
#include "TRestDataSet.h"

class TRestComponentDataSet : public TRestComponent {
   private:
    std::vector<std::string> fWeights;  //<

    std::vector<Int_t> fNSimPerNode;  //<

    std::vector<std::string> fDataSetFileNames;  //<

    std::vector<THnD*> fNodeDensity;  //<



    TRestDataSet fDataSet;  

    Bool_t fDataSetLoaded = false;  

    Bool_t ValidDataSet();

   protected:
    std::vector<Double_t> ExtractParameterizationNodes();
    std::vector<Int_t> ExtractNodeStatistics();
    void FillHistograms();

    Bool_t VariablesOk();
    Bool_t WeightsOk();

   public:
    Bool_t LoadDataSets();

    Bool_t IsDataSetLoaded() { return fDataSetLoaded; }

    Double_t GetRate(std::vector<Double_t> point) override;
    Double_t GetTotalRate() override;

    TCanvas* DrawComponent(std::vector<std::string> drawVariables, std::vector<std::string> scanVariables,
                           Int_t binScanSize = 1, TString drawOption = "");

    THnD* GetDensityForNode(Double_t value);
    THnD* GetDensityForActiveNode();
    THnD* GetDensity() { return GetDensityForActiveNode(); }

    TH1D* GetHistogram(Double_t node, std::string varName);
    TH2D* GetHistogram(Double_t node, std::string varName1, std::string varName2);
    TH3D* GetHistogram(Double_t node, std::string varName1, std::string varName2, std::string varName3);

    TH1D* GetHistogram(std::string varName);
    TH2D* GetHistogram(std::string varName1, std::string varName2);
    TH3D* GetHistogram(std::string varName1, std::string varName2, std::string varName3);

    void PrintStatistics();

    void PrintMetadata() override;
    void Initialize() override;
    void InitFromConfigFile() override;

    TRestComponentDataSet();
    TRestComponentDataSet(const char* cfgFileName, const std::string& name);
    ~TRestComponentDataSet();

    ClassDefOverride(TRestComponentDataSet, 1);
};
#endif