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