CbmRich.h

Go to the documentation of this file.

 
#ifndef CBM_RICH
#define CBM_RICH
 
 
#include "FairDetector.h"
 
#include "Rtypes.h"
#include "TGeoMatrix.h"
#include "TString.h"
#include "TVector3.h"
#include <TClonesArray.h>  // for ROOTCLING
#include <map>
 
class CbmRichRefPlanePoint;
class CbmRichPoint;
class CbmRichMirrorPoint;
class FairVolume;
class TGeoMatrix;
class TGeoNode;
class TGeoMedium;
class CbmRich : public FairDetector {
 
public:
  CbmRich();
 
  CbmRich(
    const char* name,
    Bool_t active,
    Double_t px = 0.,
    Double_t py = 0.,
    Double_t pz =
      258.75,  // Z coordinate for v16a = 270, for v17a = 258.75, for v18a = 0
    Double_t rx = 0.,
    Double_t ry = 0.,
    Double_t rz = 0.);
 
 
  virtual ~CbmRich();
 
 
  virtual void Initialize();
 
 
  virtual Bool_t ProcessHits(FairVolume* vol = 0);
 
  virtual void EndOfEvent();
 
 
  virtual void Register();
 
 
  virtual TClonesArray* GetCollection(Int_t iColl) const;
 
 
  virtual void Print(Option_t*) const;
 
 
  virtual void Reset();
 
 
  virtual void CopyClones(TClonesArray* cl1, TClonesArray* cl2, Int_t offset);
 
 
  virtual void ConstructGeometry();
 
 
  void ConstructGdmlGeometry(TGeoMatrix* geoMatrix);
 
  void ExpandNodeForGdml(TGeoNode* node);
 
 
  void ConstructOpGeometry();
 
  void SetRichGlassPropertiesForGeant4();
 
 
  virtual Bool_t CheckIfSensitive(std::string name);
 
  /*
    * \brief set fRegisterPhotonsOnSensitivePlane parameter
    */
  void SetRegisterPhotonsOnSensitivePlane(Bool_t b) {
    fRegisterPhotonsOnSensitivePlane = b;
  }
 
private:
  Int_t fPosIndex;
  // set to true if you want to register photons onto the sensitive gas plane,
  // if false then only charged particles are registered
  Bool_t fRegisterPhotonsOnSensitivePlane;
 
  TClonesArray* fRichPoints;          // MC points onto the photodetector plane
  TClonesArray* fRichRefPlanePoints;  // points on the reference plane
  TClonesArray* fRichMirrorPoints;    // mirror points
 
 
  // GDML geometry
  static std::map<TString, TGeoMedium*>
    fFixedMedia;            // List of media "repaired" after importing GMDL
  TGeoRotation* fRotation;  // Rotation matrix of the RICH detector
  TGeoCombiTrans*
    fPositionRotation;  // Full combined matrix for position and rotation of the RICH detector
 
  CbmRichPoint* AddHit(Int_t trackID,
                       Int_t detID,
                       TVector3 pos,
                       TVector3 mom,
                       Double_t time,
                       Double_t length,
                       Double_t eLoss);
 
  CbmRichPoint* AddRefPlaneHit(Int_t trackID,
                               Int_t detID,
                               TVector3 pos,
                               TVector3 mom,
                               Double_t time,
                               Double_t length,
                               Double_t eLoss);
 
  CbmRichPoint* AddMirrorHit(Int_t trackID,
                             Int_t detID,
                             TVector3 pos,
                             TVector3 mom,
                             Double_t time,
                             Double_t length,
                             Double_t eLoss);
 
  CbmRich(const CbmRich&);
 
  CbmRich& operator=(const CbmRich&);
 
  ClassDef(CbmRich, 2)
};
 
#endif