TRestAxionBufferGas.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 _TRestAxionBufferGas
#define _TRestAxionBufferGas
 
#include <TRestMetadata.h>
 
class TRestAxionBufferGas : public TRestMetadata {
   private:
    void Initialize();
    void Clear();
 
    void InitFromConfigFile();
 
    std::vector<TString> fBufferGasName;  //->
 
    std::vector<Double_t> fBufferGasDensity;  //->
 
    std::vector<std::vector<Double_t> > fAbsEnergy;  //->
 
    std::vector<std::vector<Double_t> > fGasAbsCoefficient;  //->
 
    std::vector<std::vector<Double_t> > fFactorEnergy;  //->
 
    std::vector<std::vector<Double_t> > fGasFormFactor;  //->
 
    void ReadGasData(TString gasName);
 
    Int_t FindGasIndex(TString gName);
    Int_t GetEnergyIndex(std::vector<Double_t> enVector, Double_t energy);
 
   public:
    void SetGasDensity(TString gasName, Double_t density);
    Double_t GetGasDensity(TString gasName);
 
    void SetGasMixture(TString gasMixture, TString gasDensities = "0");
 
    Int_t GetNumberOfGases() { return (Int_t)fBufferGasName.size(); }
 
    Double_t GetAbsorptionCoefficient(TString gasName, Double_t energy);
    Double_t GetFormFactor(TString gasName, Double_t energy);
 
    Double_t GetPhotonAbsorptionLength(Double_t energy);
 
    Double_t GetPhotonAbsorptionLengthIneV(Double_t energy);
 
    Double_t cmToeV(double l_Inv);
 
    Double_t GetPhotonMass(double en);
 
    Double_t GetDensityForMass(double m_gamma, double en = 4.2);
 
    void PrintAbsorptionGasData(TString gasName);
    void PrintFormFactorGasData(TString gasName);
 
    void PrintMetadata();
 
    TRestAxionBufferGas();
    TRestAxionBufferGas(const char* cfgFileName, std::string name = "");
 
    ~TRestAxionBufferGas();
 
    ClassDef(TRestAxionBufferGas, 1);
};
#endif