TRestAxionHelioscopeSignal.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_TRestAxionHelioscopeSignal
#define REST_TRestAxionHelioscopeSignal
 
#include <TRestAxionBufferGas.h>
#include <TRestAxionField.h>
#include <TRestAxionSolarFlux.h>
 
#include "TRestComponent.h"
 
class TRestAxionHelioscopeSignal : public TRestComponent {
   private:
    std::string fConversionType = "IAXO";
 
    Int_t fBores = 2;
 
    Double_t fMagnetRadius = 350;
 
    Double_t fMagnetLength = 10000;
 
    Double_t fMagnetStrength = 2;
 
    Double_t fOpticsEfficiency = 1;
 
    Double_t fWindowEfficiency = 1;
 
    TRestAxionBufferGas* fGas = nullptr;
 
    TRestAxionSolarFlux* fFlux = nullptr;
 
    TRestAxionField* fField = nullptr;
 
   protected:
    void InitFromConfigFile() override;
 
    void FillHistograms() override;
 
   public:
    Double_t GetSignalRate(std::vector<Double_t> point, Double_t mass = 0);
    Double_t GetSignalRate(Double_t mass, Double_t Eo, Double_t Ef);
 
    TRestAxionBufferGas* GetGas() { return fGas; }
    TRestAxionSolarFlux* GetFlux() { return fFlux; }
 
    void PrintMetadata() override;
 
    void SetNumberOfBores(const Int_t& bores) { fBores = bores; }
    void SetMagnetRadius(const Double_t& radius) { fMagnetRadius = radius; }
    void SetMagnetLength(const Double_t& length) { fMagnetLength = length; }
    void SetMagnetStrength(const Double_t& strength) { fMagnetStrength = strength; }
    void SetType(const std::string& type) { fConversionType = type; }
 
    Int_t GetNumberOfBores() const { return fBores; }
    Double_t GetMagnetRadius() const { return fMagnetRadius; }
    Double_t GetMagnetLength() const { return fMagnetLength; }
    Double_t GetMagnetStrength() const { return fMagnetStrength; }
    std::string GetType() const { return fConversionType; }
 
    void Initialize() override;
    TRestAxionHelioscopeSignal(const char* cfgFileName, const std::string& name);
    TRestAxionHelioscopeSignal();
    ~TRestAxionHelioscopeSignal();
 
    ClassDefOverride(TRestAxionHelioscopeSignal, 1);
};
#endif