PndDrcHitProducerReal.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                   PndDrcHitProducer header file               -----
// -----               Created 30/10/09  by Dipanwita Dutta
// -----                                                               -----
// -----                                                               -----
// -------------------------------------------------------------------------

#ifndef PNDDRCHITPRODUCERREAL_H
#define PNDDRCHITPRODUCERREAL_H

#include <PndPersistencyTask.h>
#include "TClonesArray.h"
#include "PndMCTrack.h"
#include "PndDrcHit.h"
#include "PndDrcPDHit.h"
#include "PndStack.h"
#include "FairBaseParSet.h"
#include "PndGeoDrcPar.h"
#include "PndGeoDrc.h"
#include "TString.h"
#include <string>
#include "TH1.h"
#include "TH2.h"

#ifndef ROOT_TParticlePDG
#include "TParticlePDG.h"
#endif
#ifndef ROOT_TDatabasePDG
#include "TDatabasePDG.h"
#endif

class PndDrcHitProducerReal : public PndPersistencyTask {

 public:
  PndDrcHitProducerReal();

  PndDrcHitProducerReal(Int_t verbose, Int_t det_type);

  virtual ~PndDrcHitProducerReal();

  //  virtual void SetParContainers();
  virtual InitStatus Init();

  virtual void Exec(Option_t *option);

  virtual void Finish();

  PndDrcHit *AddHit(Int_t detID, TVector3 posHit, TVector3 dPosHit, Double_t thetaC, Double_t errThetaC, Int_t index);

  PndDrcPDHit *AddPDHit(Int_t detID, TVector3 posPDHit, TVector3 dPosPDHit, Double_t time, Double_t timeThreshold, Int_t index1);

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

  void SetPhotonDetEff(Double_t &lambda_min, Double_t &lambda_max, Double_t &lambda_step, Double_t efficiency[]);
  void SetFakeDetEff(Double_t &lambda_min, Double_t &lambda_max, Double_t &lambda_step, Double_t efficiency[]);
  void SetPhotonDetEffOld(Double_t &lambda_min, Double_t &lambda_max, Double_t &lambda_step, Double_t efficiency[]);
  void SetPhotonDetEffNew(Double_t &lambda_min, Double_t &lambda_max, Double_t &lambda_step, Double_t efficiency[]);

  void SetPhotonTransportEff(Double_t &, Double_t &, Double_t &, Double_t &, Int_t &, Double_t fEfficiency[]);
  void NumberOfBounces(TVector3, TVector3, Int_t *, Int_t *, Double_t *, Double_t *);
  Double_t FindPhiRot(Double_t, Double_t);
  Double_t FindOutPoint(Double_t, Double_t, Double_t, Double_t *, Bool_t);

  void SetIsDetEfficiency(Bool_t isDetEff) { fisDetEff = isDetEff; }
  void SetIsPixelization(Bool_t isPixel) { fisPixel = isPixel; }
  void SetTransportEfficiency(Bool_t isTran = 0) { fisTransportEff = isTran; }

  //########################################
  void SetTilt(Double_t tilt = 0.) { ftilt = tilt; } // degrees
  //########################################

  void DrawDetEfficiency(Bool_t dr = kFALSE) { fDrawing = dr; }

 protected:
 private:
  void ProcessPhotonPoint();
  void ProcessBarPoint();
  void Smear(Double_t &time, Double_t sigt);

  // basic parameters of DIRC
  Double_t fpi;
  Double_t fzup;
  Double_t fzdown;
  Double_t fradius;
  Double_t fhthick;
  Double_t fpipehAngle;
  Double_t fbbGap;
  Double_t fbbnum;
  Double_t fbarnum;
  Double_t fphi0;
  Double_t fdphi;
  Double_t flside;
  Double_t fbarwidth;

  Bool_t fisDetEff;
  Bool_t fisPixel;
  Bool_t fisTransportEff;
  Bool_t fDrawing;
  Double_t ftilt;
  Int_t fDetectorID;
  TVector3 fPosHit;
  TVector3 fDPosHit;
  TVector3 fPosPDHit;
  TVector3 fDPosPDHit;
  Double_t fThetaC, fErrThetaC;
  Double_t fTime, fTimeThreshold;
  Int_t fRefIndex;
  Int_t fPDRefIndex;

  TClonesArray *fBarPointArray; // DRC MC points in the bars
  TClonesArray *fPDPointArray;  // DRC MC points in the photon plane
  TClonesArray *fHitArray;      // DRC hits
  TClonesArray *fPDHitArray;    // DRC Photon Detectorhits
  TClonesArray *fMCArray;       // DRC Hits in the photon detector

  // TObjArray *fVolumeArray;
  // TClonesArray* fListStack;     // Tracks

  PndGeoDrcPar *fPar;

  PndGeoDrc *fGeo; 

  Int_t fVerbose;

  Int_t fDetType; // detector type

  Double_t nRefrac;          // Refractive index of photon detector
  Double_t fPixelDim;        // Pixel Diemsion of photocathode
  Double_t fSigmat;          // Time Resolution in ps
  Double_t fCollectionEff;   // Collection Efficiency
  Double_t fPackingFraction; // Packing Fraction or Active Area Ratio
  Double_t fRoughness;       // Surface roughness (bars)
  Int_t detection;
  Int_t nevents;

  void SetParameters();

  ClassDef(PndDrcHitProducerReal, 1)
};

#endif