CbmVMCSettings.h

Go to the documentation of this file.

/* Copyright (C) 2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Florian Uhlig [committer] */
 
#ifndef CBMVMCSETTINGS_H
#define CBMVMCSETTINGS_H 1
 
#include <Logger.h>
 
#include <TObject.h>
 
template<typename T>
void CheckValueInRange(const T& value, const T& low, const T& high, std::string functionName)
{
  if ((value < low) || (value > high)) {
    LOG(fatal) << "You try to set the value " << value << "which is out of the bounds(" << low << "," << high
               << ") for the function " << functionName;
  }
}
 
class TVirtualMC;
 
class CbmVMCSettings : public TObject {
 
public:
  CbmVMCSettings()                      = default;
  ~CbmVMCSettings()                     = default;
  CbmVMCSettings(const CbmVMCSettings&) = delete;
  CbmVMCSettings& operator=(const CbmVMCSettings&) = delete;
 
 
  void Init(TVirtualMC*);
 
  void SetProcessPairProduction(Int_t val)
  {
    CheckValueInRange(val, 0, 2, "SetProcessPairProduction");
    fProcessPairProduction = val;
  }
 
  void SetProcessComptonScattering(Int_t val)
  {
    CheckValueInRange(val, 0, 2, "SetProcessComptonScattering");
    fProcessComptonScattering = val;
  }
 
  void SetProcessPhotoEffect(Int_t val)
  {
    CheckValueInRange(val, 0, 2, "SetProcessPhotoEffect");
    fProcessPhotoEffect = val;
  }
 
  void SetProcessPhotoFission(Int_t val)
  {
    CheckValueInRange(val, 0, 2, "SetProcessPhotoFission");
    fProcessPhotoFission = val;
  }
 
  void SetProcessDeltaRay(Int_t val)
  {
    CheckValueInRange(val, 0, 2, "SetProcessDeltaRay");
    fProcessDeltaRay = val;
  }
 
  void SetProcessAnnihilation(Int_t val)
  {
    CheckValueInRange(val, 0, 2, "SetProcessAnnihilation");
    fProcessAnnihilation = val;
  }
 
  void SetProcessBremsstrahlung(Int_t val)
  {
    CheckValueInRange(val, 0, 2, "SetProcessBremsstrahlung");
    fProcessBremsstrahlung = val;
  }
 
  void SetProcessHadronicInteraction(Int_t val)
  {
    CheckValueInRange(val, 0, 5, "SetProcessHadronicInteraction");
    fProcessHadronicInteraction = val;
  }
 
  void SetProcessMuonNuclearInteraction(Int_t val)
  {
    CheckValueInRange(val, 0, 2, "SetProcessMuonNuclearInteraction");
    fProcessMuonNuclearInteraction = val;
  }
 
  void SetProcessDecay(Int_t val)
  {
    CheckValueInRange(val, 0, 2, "SetProcessDecay");
    fProcessDecay = val;
  }
 
  void SetProcessEnergyLossModel(Int_t val)
  {
    CheckValueInRange(val, 0, 4, "SetProcessEnergyLossModel");
    fProcessEnergyLossModel = val;
  }
 
  void SetProcessMultipleScattering(Int_t val)
  {
    CheckValueInRange(val, 0, 3, "SetProcessMultipleScattering");
    fProcessMultipleScattering = val;
  }
 
  void SetEnergyCutGammas(Double_t val)
  {
    CheckValueInRange(val, 0., 100., "SetEnergyCutGammas");
    fEnergyCutGammas = val;
  }
 
  void SetEnergyCutElectrons(Double_t val)
  {
    CheckValueInRange(val, 0., 100., "SetEnergyCutElectrons");
    fEnergyCutElectrons = val;
  }
 
  void SetEnergyCutNeutralHadrons(Double_t val)
  {
    CheckValueInRange(val, 0., 100., "SetEnergyCutNeutralHadrons");
    fEnergyCutNeutralHadrons = val;
  }
 
  void SetEnergyCutChargedHadrons(Double_t val)
  {
    CheckValueInRange(val, 0., 100., "SetEnergyCutChargedHadrons");
    fEnergyCutChargedHadrons = val;
  }
 
  void SetEnergyCutMuons(Double_t val)
  {
    CheckValueInRange(val, 0., 100., "SetEnergyCutMuons");
    fEnergyCutMuons = val;
  }
 
  void SetEnergyCutElectronBremsstrahlung(Double_t val)
  {
    CheckValueInRange(val, 0., 100., "SetEnergyCutElectronBremsstrahlung");
    fEnergyCutElectronBremsstrahlung = val;
  }
 
  void SetEnergyCutMuonHadronBremsstrahlung(Double_t val)
  {
    CheckValueInRange(val, 0., 100., "SetEnergyCutMuonHadronBremsstrahlung");
    fEnergyCutMuonHadronBremsstrahlung = val;
  }
 
  void SetEnergyCutElectronDeltaRay(Double_t val)
  {
    CheckValueInRange(val, 0., 100., "SetEnergyCutElectronDeltaRay");
    fEnergyCutElectronDeltaRay = val;
  }
 
  void SetEnergyCutMuonDeltaRay(Double_t val)
  {
    CheckValueInRange(val, 0., 100., "SetEnergyCutMuonDeltaRay");
    fEnergyCutMuonDeltaRay = val;
  }
 
  void SetEnergyCutMuonPairProduction(Double_t val)
  {
    CheckValueInRange(val, 0., 100., "SetEnergyCutMuonPairProduction");
    fEnergyCutMuonPairProduction = val;
  }
 
  void SetTimeCutTof(Double_t val)
  {
    CheckValueInRange(val, 0., 2., "SetTimeCutTof");
    fTimeCutTof = val;
  }
 
private:
  Int_t fProcessPairProduction {1};
  Int_t fProcessComptonScattering {1};
  Int_t fProcessPhotoEffect {1};
  Int_t fProcessPhotoFission {0};
  Int_t fProcessDeltaRay {1};
  Int_t fProcessAnnihilation {1};
  Int_t fProcessBremsstrahlung {1};
  Int_t fProcessHadronicInteraction {1};
  Int_t fProcessMuonNuclearInteraction {1};
  Int_t fProcessDecay {1};
  Int_t fProcessEnergyLossModel {1};
  Int_t fProcessMultipleScattering {1};
 
  Double_t fEnergyCutGammas {1.e-3};                    // GeV
  Double_t fEnergyCutElectrons {1.e-3};                 // GeV
  Double_t fEnergyCutNeutralHadrons {1.e-3};            // GeV
  Double_t fEnergyCutChargedHadrons {1.e-3};            // GeV
  Double_t fEnergyCutMuons {1.e-3};                     // GeV
  Double_t fEnergyCutElectronBremsstrahlung {1.e-3};    // GeV
  Double_t fEnergyCutMuonHadronBremsstrahlung {1.e-3};  // GeV
  Double_t fEnergyCutElectronDeltaRay {1.e-3};          // GeV
  Double_t fEnergyCutMuonDeltaRay {1.e-3};              // GeV
  Double_t fEnergyCutMuonPairProduction {1.e-3};        // GeV
 
  Double_t fTimeCutTof {1.0};  // s
 
  ClassDef(CbmVMCSettings, 1);
};
 
#endif /* CBMVMCSETTINGS_H */