PndRich.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

#ifndef PNDRICH_H
#define PNDRICH_H

#include "FairDetector.h"
#include "PndGeoHandling.h"
#include "PndRichGeo.h"

#include "TVector3.h"
#include "TLorentzVector.h"

class PndRichPDPoint;
class PndRichBarPoint;
class FairVolume;
class TClonesArray;

class PndRich : public FairDetector {

 public:
  PndRich(const char *Name, Bool_t Active);

  PndRich();

  virtual ~PndRich();

  virtual void Initialize();

  virtual Bool_t ProcessHits(FairVolume *v = 0);

  virtual void Register();

  virtual TClonesArray *GetCollection(Int_t iColl) const;

  virtual void Reset();

  void ConstructGeometry();
  void ConstructOpGeometry();

  void SetRunCherenkov(Bool_t const &ch) { fRunCherenkov = ch; };

  void SetGeoVersion(UInt_t const &version) { fGeoVersion = version; };

  void UseProtection(Bool_t const &up) { fUseProtection = up; };

  std::vector<std::string> fListOfSensitives;
  bool CheckIfSensitive(std::string name);

  Double_t lhcbaerindex(Double_t n400, Double_t wl);

  PndRichPDPoint *AddPDPoint(Int_t trackID, Int_t detID, TVector3 pos, TVector3 mom, Double_t time, Double_t length, Double_t eLoss, UInt_t EventId = 0);
  PndRichBarPoint *AddBarPoint(Int_t trackID, Int_t detID, TVector3 pos, TVector3 mom, Double_t time, Double_t length, Int_t pdgCode, Double_t thetaC, Int_t eventID, Double_t mass,
                               TVector3 pos0, TVector3 mom0);

  virtual void CopyClones(TClonesArray *, TClonesArray *, Int_t) { ; } // cl1 cl2  offset //[R.K.03/2017] unused variable(s)
  virtual void SetSpecialPhysicsCuts() { ; }
  virtual void EndOfEvent();
  virtual void FinishPrimary() { ; }
  virtual void FinishRun();
  virtual void BeginPrimary() { ; }
  virtual void PostTrack() { ; }
  virtual void PreTrack() { ; }
  virtual void BeginEvent();

 private:
  Int_t fTrackID;      
  Int_t fVolumeID;     
  TLorentzVector fPos; 
  TLorentzVector fMom; 
  Double_t fTime;      
  Double_t fLength;    
  Double_t fELoss;     

  PndRichGeo *fGeo;
  PndGeoHandling *fGeoH; 

  Bool_t fUseProtection;
  Bool_t fRunCherenkov;
  UInt_t fGeoVersion;

  TClonesArray *fPndRichPDPointCollection;
  TClonesArray *fPndRichBarPointCollection;

  PndRich(const PndRich &);
  PndRich &operator=(const PndRich &);

  void DefGeoVersion();

  std::vector<Double_t> fWlPhoton;
  std::vector<Double_t> fPDE;

  std::map<Int_t, Int_t> trackid;

  Double_t fnOpt;
  Double_t fZabar;

  // TGeoElement *C;
  // TGeoElement *H;
  // TGeoElement *B;
  // TGeoIsotope *B10;
  // TGeoIsotope *B11;
  // TGeoMaterial *matRcihProt;
  // TGeoMedium *med;

  ClassDef(PndRich, 1)
};

#endif // PNDRICH_H