CbmTrdParModGas.h

Go to the documentation of this file.

/* Copyright (C) 2018-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Florian Uhlig [committer], Alexandru Bercuci */
 
#ifndef CBMTRDPARMODGAS_H
#define CBMTRDPARMODGAS_H
 
#include "CbmTrdParMod.h"  // for CbmTrdParMod
 
#include <Rtypes.h>      // for CLRBIT, SETBIT, TESTBIT, THashConsistencyH...
#include <RtypesCore.h>  // for Double_t, Int_t, Float_t, Char_t, UShort_t
#include <TString.h>     // for TString
 
class TDirectory;
class TH2F;
 
#define NSHELLS 3  // no of atomic shells to be considered in the PE effect
 
class CbmTrdParModGas : public CbmTrdParMod {
public:
  enum CbmTrdParModGasDef
  {
    kDetType = 0 
      ,
    kNobleGasType 
      ,
    kPID 
  };
  CbmTrdParModGas(const char* title = "TRD gas properties definition");
  virtual ~CbmTrdParModGas();
 
  Double_t GetCO2() const { return fPercentCO2; }
  const Char_t* GetDetName() const { return (GetDetType() ? "MB" : "GSI"); }
  Int_t GetDetType() const { return TESTBIT(fConfig, kDetType); }
  Double_t GetDriftTime(Double_t y0, Double_t z0) const;
  void GetElectricPotential(Int_t& ua, Int_t& ud) const
  {
    ua = fUa;
    ud = fUd;
  }
  Float_t GetBindingEnergy(const Char_t shell = 'K', Bool_t main = kTRUE) const;
  Float_t EkevFC(Float_t ekev) const;
  Float_t EfCkeV(Float_t efc) const;
  TString GetFileName() const { return fFileNamePID; }
  Double_t GetNobleGas() const { return 1. - fPercentCO2; }
  const Char_t* GetNobleGasName() const { return (GetNobleGasType() - 1 ? "Ar" : "Xe"); }
  Int_t GetNobleGasType() const { return TESTBIT(fConfig, kNobleGasType) ? 2 : 1; }
  Float_t GetNonIonizingBR(const Char_t shell = 'K') const;
  Char_t GetPEshell(Float_t Ex) const;
  Int_t GetPidType() const { return TESTBIT(fConfig, kPID); }
  Int_t GetShellId(const Char_t shell) const;
  UShort_t GetUanode() const { return fUa; }
  UShort_t GetUdrift() const { return fUd; }
  virtual void Print(Option_t* opt = "") const;
  Double_t ScanDriftTime(Double_t y0, Double_t z0, Double_t dzdy, Double_t dy) const;
 
  void SetDriftMap(TH2F* hm, TDirectory* d);
  void SetDw(Double_t dw) { fDw = dw; }
  void SetGasThick(Double_t sz) { fGasThick = sz; }
  void SetDetType(Int_t gsi = 0) { gsi ? SETBIT(fConfig, kDetType) : CLRBIT(fConfig, kDetType); }
  void SetCO2(Double_t p) { fPercentCO2 = p; }
  void SetFileName(const Char_t* fn) { fFileNamePID = fn; }
  void SetNobleGas(Double_t p) { fPercentCO2 = 1 - p; }
  void SetNobleGasType(Int_t ar = 1) { ar ? SETBIT(fConfig, kNobleGasType) : CLRBIT(fConfig, kNobleGasType); }
  void SetPidType(Int_t like = 1) { like ? SETBIT(fConfig, kPID) : CLRBIT(fConfig, kPID); }
  void SetElectricPotential(Int_t ua, Int_t ud)
  {
    fUa = ua;
    fUd = ud;
  }
 
private:
  CbmTrdParModGas(const CbmTrdParModGas& ref);
  const CbmTrdParModGas& operator=(const CbmTrdParModGas& ref);
 
  Int_t fConfig;         
  UShort_t fUa;          
  UShort_t fUd;          
  Double_t fDw;          
  Double_t fGasThick;    
  Double_t fPercentCO2;  
  TH2F* fDriftMap;       
  TString fFileNamePID;  
 
  static Float_t fgkBindingEnergy[2][NSHELLS];  
  static Float_t fgkBR[2][NSHELLS - 1];         
  static Float_t fgkWi[3];                      
  static Float_t fgkGGainUaPar[2];              
  static Float_t fgkE0;                         
  ClassDef(CbmTrdParModGas,
           1)  // Definition of gas parameters for one TRD module
};
 
#endif