TRestAxionSolarQCDFlux.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 _TRestAxionSolarQCDFlux
#define _TRestAxionSolarQCDFlux
 
#include <TRestAxionSolarFlux.h>
#include <TRestAxionSolarModel.h>
 
class TRestAxionSolarQCDFlux : public TRestAxionSolarFlux {
   private:
    std::string fFluxDataFile = "";  //<
 
    std::string fFluxSptFile = "";  //<
 
    Double_t fBinSize = 0;  //<
 
    Double_t fPeakSigma = 0;  //<
 
    std::vector<TH1F*> fFluxTable;  
 
    std::map<Double_t, TH1F*> fFluxLines;  
 
    std::vector<Double_t> fFluxTableIntegrals;  
 
    std::vector<Double_t> fFluxLineIntegrals;  
 
    Double_t fTotalMonochromaticFlux = 0;  
 
    Double_t fTotalContinuumFlux = 0;  
 
    Double_t fFluxRatio = 0;  
 
    TH1F* fContinuumHist = nullptr;  
 
    TH1F* fMonoHist = nullptr;  
 
    TH1F* fTotalHist = nullptr;  
 
    void ReadFluxFile();
    void LoadContinuumFluxTable();
    void LoadMonoChromaticFluxTable();
    void IntegrateSolarFluxes();
 
   public:
    Bool_t isSolarTableLoaded() { return fFluxTable.size() > 0; }
 
    Bool_t isSolarSpectrumLoaded() { return fFluxLines.size() > 0; }
 
    Double_t IntegrateFluxInRange(TVector2 eRange = TVector2(-1, -1), Double_t mass = 0) override;
 
    std::pair<Double_t, Double_t> GetRandomEnergyAndRadius(TVector2 eRange = TVector2(-1, -1),
                                                           Double_t mass = 0) override;
 
    Bool_t LoadTables() override;
 
    Double_t GetTotalFlux(Double_t mass = 0) override {
        return fTotalContinuumFlux + fTotalMonochromaticFlux;
    }
 
    TH1F* GetEnergySpectrum(Double_t m = 0) override { return GetTotalSpectrum(); }
 
    TH1F* GetContinuumSpectrum();
    TH1F* GetMonochromaticSpectrum();
    TH1F* GetTotalSpectrum();
 
    virtual TCanvas* DrawSolarFlux() override;
 
    void InitializeSolarTable(TRestAxionSolarModel* model) {
        // TOBE implemented
        // This method should initialize the tables fFluxTable and fFluxLines
    }
 
    void ExportTables(Bool_t ascii = false) override;
 
    void PrintMetadata() override;
 
    void PrintContinuumSolarTable();
    void PrintIntegratedRingFlux();
    void PrintMonoChromaticFlux();
 
    TRestAxionSolarQCDFlux();
    TRestAxionSolarQCDFlux(const char* cfgFileName, std::string name = "");
    ~TRestAxionSolarQCDFlux();
 
    ClassDefOverride(TRestAxionSolarQCDFlux, 1);
};
#endif