CbmRichHitProducer.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file CbmRichHitProducer.h

#ifndef CBM_RICH_HIT_PRODUCER
#define CBM_RICH_HIT_PRODUCER

#include "FairTask.h"

class CbmGeoRichPar;
class TClonesArray;
class TVector3;

class CbmRichHitProducer : public FairTask
{
public:
   CbmRichHitProducer();

   virtual ~CbmRichHitProducer();

   virtual void SetParContainers();

   virtual InitStatus Init();

   virtual void Exec(
        Option_t* option);

   virtual void Finish();

   void SetDetectorType(Int_t detType){ fDetType = detType;}

   void SetNofNoiseHits(Int_t noise) {fNofNoiseHits = noise;}

   void SetCollectionEfficiency(Double_t collEff){fCollectionEfficiency = collEff;}

   void SetSigmaMirror(Double_t sigMirror) {fSigmaMirror = sigMirror;}

   void SetCrossTalkHitProb(Double_t crosstalk) {fCrossTalkHitProb = crosstalk;}

   void AddCrossTalkHits(
         Double_t x,
         Double_t y,
         Int_t pointInd,
         Int_t RichDetID);

   void AddHit(
         TVector3 &posHit,
         TVector3 &posHitErr,
         Int_t address,
         Int_t pmtID,
         Double_t ampl,
         Int_t index);

   void FindRichHitPositionSinglePMT(
         Double_t xPoint,
         Double_t yPoint,
         Double_t& xHit,
         Double_t& yHit,
         Int_t& pmtID);

   void FindRichHitPositionMAPMT(
         Double_t sigma,
         Double_t xPoint,
         Double_t yPoint,
         Double_t& xHit,
         Double_t& yHit,
         Int_t & pmtID);

   void FindRichHitPositionCsI(
         Double_t xPoint,
         Double_t yPoint,
         Double_t& xHit,
         Double_t& yHit,
         Int_t& pmtID);

   Double_t OnePhotonAmplitude(
        Double_t x);

   Double_t GetAmplitude();

   void SetPhotoDetPar(
         Int_t det_type,
         Double_t& lambda_min,
         Double_t& lambda_max,
         Double_t& lambda_step,
         Double_t efficiency[]);

   static void TiltPoint(
         TVector3 *inPos,
         TVector3 *outPos,
         Double_t phi,
         Double_t theta,
         Double_t detZOrig,
         Bool_t noTilting = false);

private:

   TClonesArray* fRichPoints; // RICH MC points
   TClonesArray* fRichHits; // RICH hits
   //TClonesArray* fMcTracks; // Monte-Carlo tracks

   Int_t fNHits; // Number of hits
   Int_t fNDoubleHits; // Number of double hits

   static constexpr Double_t c = 2.998E8; // speed of light
   static constexpr Double_t h = 6.626E-34; // Planck constant
   static constexpr Double_t e = 1.6022E-19; // elementary charge

   Double_t fNRefrac; // refraction index
   Int_t fDetection; // flag for detection
   Int_t fNEvents; // event number

   Double_t fDetX; // X-coordinate of photodetector
   Double_t fDetY; // Y-coordinate of photodetector
   Double_t fDetZ; // Z-coordinate of photodetector
   Double_t fDetZ_org; // X-coordinate of photodetector (original from parameter file)
   Double_t fDetWidthX; // width of photodetector in x
   Double_t fDetWidthY; // width of photodetector in y

   TObjArray *fSensNodes;
   TObjArray *fPassNodes;
   CbmGeoRichPar *fPar;

   // Parameters of photodetector
   Double_t fPhotomulRadius; // radius of photomultiplier [cm], set according to detector type
   Double_t fPhotomulDist; // distance between PMT tube [cm], set according to detector type
   Int_t fDetType; // detector type
   Int_t fNofNoiseHits; // number of noise hits
   Double_t fCollectionEfficiency; // collection efficiency for photoelectrons in PMT optics
   Double_t fSigmaMirror; // additinal smearing of MC Points due to light scattering in mirror

   Double_t fTheta; // angle by which photodetector was tilted around x-axis
   Double_t fPhi; // angle by which photodetector was tilted around y-axis

   Double_t fCrossTalkHitProb;// probability of cross talk hit
   Int_t fNofCrossTalkHits; // Number of cross talk hits

   CbmRichHitProducer(const CbmRichHitProducer&);

   CbmRichHitProducer& operator=(const CbmRichHitProducer&);

   ClassDef(CbmRichHitProducer,2)
};

#endif