CbmRichPmt.h

Go to the documentation of this file.

/* Copyright (C) 2015-2020 GSI/JINR-LIT, Darmstadt/Dubna
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Semen Lebedev, Andrey Lebedev [committer], Florian Uhlig */
 
/*
 * CbmRichPmt.h
 *
 *  Created on: Dec 16, 2015
 *      Author: slebedev
 */
 
#ifndef RICH_DETECTOR_CBMRICHPMT_H_
#define RICH_DETECTOR_CBMRICHPMT_H_
 
#include "CbmRichPmtType.h"  // for CbmRichPmtTypeEnum
 
#include <Rtypes.h>      // for THashConsistencyHolder, ClassDef
#include <RtypesCore.h>  // for Double_t, Bool_t
#include <TObject.h>     // for TObject
 
#include <map>     // for map
#include <vector>  // for vector
 
class TGraph;
 
class CbmRichPmtQEData {
public:
  CbmRichPmtTypeEnum fDetectorType;   // pmt type
  Double_t fLambdaMin;                // minimum wavelength in QE table
  Double_t fLambdaMax;                // maximum wavelength in QE table
  Double_t fLambdaStep;               // wavelaength in QE table
  std::vector<Double_t> fEfficiency;  // Array of QE
 
  CbmRichPmtQEData() : fDetectorType(), fLambdaMin(0.), fLambdaMax(0.), fLambdaStep(0.), fEfficiency() {}
};
 
class CbmRichPmt : public TObject {
public:
  CbmRichPmt();
 
  virtual ~CbmRichPmt();
 
  static Double_t getLambda(Double_t momentum);
 
  Bool_t isPhotonDetected(CbmRichPmtTypeEnum detType, Double_t momentum);
 
  void SetCollectionEfficiency(Double_t collEff) { fCollectionEfficiency = collEff; }
 
  /*
         * Get QE curve as TGraph for specified detector type
         */
  TGraph* getQEGraph(CbmRichPmtTypeEnum detType);
 
private:
  Double_t fCollectionEfficiency;  // collection efficiency. Final QE = QE * fCollectionEfficiency
  std::map<CbmRichPmtTypeEnum, CbmRichPmtQEData*> fPmtDataMap;  // store QE for specific PMT type
 
  static const Double_t c;    // speed of light
  static const Double_t h;    // Planck constant
  static const Double_t e;    // elementary charge
  Double_t fRefractiveIndex;  // refractive index of the gas. Set in the constructor from CbmRichGeoManager
 
  void InitQE();
 
  void ClearMap();
 
 
  CbmRichPmt(const CbmRichPmt&);
 
  CbmRichPmt& operator=(const CbmRichPmt&);
 
  ClassDef(CbmRichPmt, 1)
};
 
#endif /* RICH_DETECTOR_CBMRICHPMT_H_ */