TRestAxionFieldPropagationProcess.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_TRestAxionFieldPropagationProcess
#define RestCore_TRestAxionFieldPropagationProcess
 
#include "TRestAxionEvent.h"
#include "TRestAxionEventProcess.h"
#include "TRestAxionField.h"
#include "TRestAxionMagneticField.h"
#include "TRestPhysics.h"
#include "TVector3.h"
#include "TVectorD.h"
 
class TRestAxionFieldPropagationProcess : public TRestAxionEventProcess {
   private:
    Double_t fBufferGasAdditionalLength = 0;  //<
 
    Double_t fAccuracy = 0.1;
 
    Int_t fNumIntervals = 100;
 
    Int_t fQawoLevels = 20;
 
    TVector3 fReMap = TVector3(30, 30, 100);
 
    TRestAxionMagneticField* fMagneticField = nullptr;  
 
    TRestAxionField* fAxionField = nullptr;  
 
    TRestAxionBufferGas* fBufferGas = nullptr;  
 
    void Initialize() override;
 
   public:
    void InitProcess() override;
 
    RESTValue GetInputEvent() const override { return fAxionEvent; }
    RESTValue GetOutputEvent() const override { return fAxionEvent; }
 
    TRestEvent* ProcessEvent(TRestEvent* eventInput) override;
 
    void PrintMetadata() override {
        BeginPrintProcess();
 
        RESTMetadata << "Integration parameters" << RESTendl;
        RESTMetadata << "- Integration accuracy/tolerance : " << fAccuracy << RESTendl;
        RESTMetadata << "- Max number of integration intervals : " << fNumIntervals << RESTendl;
        RESTMetadata << "- Number of QAWO levels : " << fQawoLevels << RESTendl;
        RESTMetadata << " " << RESTendl;
        RESTMetadata << "Field re-mapping size : (" << fReMap.X() << ", " << fReMap.Y() << ", " << fReMap.Z()
                     << ")" << RESTendl;
        RESTMetadata << " " << RESTendl;
        RESTMetadata << "Buffer gas additional length : " << fBufferGasAdditionalLength * units("m") << " m"
                     << RESTendl;
 
        EndPrintProcess();
    }
 
    const char* GetProcessName() const override { return "axionFieldPropagation"; }
 
    // Constructor
    TRestAxionFieldPropagationProcess();
    TRestAxionFieldPropagationProcess(char* cfgFileName);
 
    // Destructor
    ~TRestAxionFieldPropagationProcess();
 
    ClassDefOverride(TRestAxionFieldPropagationProcess, 2);
};
#endif