TRestDetectorElectronDiffusionProcess.h

 
#ifndef RestCore_TRestDetectorElectronDiffusionProcess
#define RestCore_TRestDetectorElectronDiffusionProcess
 
#include <TRandom3.h>
#include <TRestDetectorGas.h>
#include <TRestDetectorHitsEvent.h>
#include <TRestDetectorReadout.h>
#include <TRestEventProcess.h>
 
class TRestDetectorElectronDiffusionProcess : public TRestEventProcess {
   private:
    TRestDetectorHitsEvent* fInputHitsEvent;   
    TRestDetectorHitsEvent* fOutputHitsEvent;  
 
    TRestDetectorGas* fGas;          
    TRestDetectorReadout* fReadout;  
 
    TRandom3* fRandom;  
 
    void InitFromConfigFile() override;
 
    void Initialize() override;
 
    void LoadDefaultConfig();
 
   protected:
    Double_t fElectricField;
    Double_t fAttachment;
    Double_t fGasPressure;
    Double_t fWValue;
    Double_t fLongitudinalDiffusionCoefficient;
    Double_t fTransversalDiffusionCoefficient;
    Bool_t fPoissonElectronExcitation;
    Bool_t fUnitElectronEnergy;
    UInt_t fMaxHits;
    Double_t fSeed = 0;
    Bool_t fCheckIsInside = true;
 
   public:
    RESTValue GetInputEvent() const override { return fInputHitsEvent; }
    RESTValue GetOutputEvent() const override { return fOutputHitsEvent; }
 
    void InitProcess() override;
    TRestEvent* ProcessEvent(TRestEvent* inputEvent) override;
    void EndProcess() override;
 
    void LoadConfig(const std::string& configFilename, const std::string& name = "");
 
    void PrintMetadata() override {
        BeginPrintProcess();
 
        RESTMetadata << " eField : " << fElectricField * units("V/cm") << " V/cm" << RESTendl;
        RESTMetadata << " attachment coeficient : " << fAttachment << " V/cm" << RESTendl;
        RESTMetadata << " gas pressure : " << fGasPressure << " atm" << RESTendl;
        RESTMetadata << " longitudinal diffusion coefficient : " << fLongitudinalDiffusionCoefficient
                     << " cm^1/2" << RESTendl;
        RESTMetadata << " transversal diffusion coefficient : " << fTransversalDiffusionCoefficient
                     << " cm^1/2" << RESTendl;
        RESTMetadata << " W value : " << fWValue << " eV" << RESTendl;
 
        RESTMetadata << " Maximum number of hits : " << fMaxHits << RESTendl;
 
        RESTMetadata << " seed : " << fSeed << RESTendl;
 
        EndPrintProcess();
    }
 
    TRestMetadata* GetProcessMetadata() const { return fGas; }
 
    const char* GetProcessName() const override { return "electronDiffusion"; }
 
    inline Double_t GetElectricField() const { return fElectricField; }
    inline Double_t GetAttachmentCoefficient() const { return fAttachment; }
    inline Double_t GetGasPressure() const { return fGasPressure; }
 
    // Constructor
    TRestDetectorElectronDiffusionProcess();
    TRestDetectorElectronDiffusionProcess(const char* configFilename);
    // Destructor
    ~TRestDetectorElectronDiffusionProcess();
 
    ClassDefOverride(TRestDetectorElectronDiffusionProcess, 4);  // Template for a REST "event process" class
                                                                 // inherited from TRestEventProcess
};
#endif