// -------------------------------------------------------------------------
// -----                   PndPhoGunShortP header file                  -----
// -----               Created 11/09/13 by Harphool Kumawat            -----
// -----                                                               -----
// -----                                                               -----
// -------------------------------------------------------------------------

/**  PhoGun.h

 **
 ** Class for Analysising DRC Cherenkov Photons fired in the DIRC barrel directly


 **/
 
#ifndef PNDPHOGUNSHORTP_H
#define PNDPHOGUNSHORTP_H
 
#include "FairTask.h"
#include "TClonesArray.h"
#include "PndMCTrack.h"
#include "PndDrcHit.h"
#include "PndDrcPDHit.h"
#include "PndStack.h"
#include "FairBaseParSet.h"
#include "PndGeoDrcPar.h"
#include "TString.h"
#include <string>

#include "PndDrcBarPoint.h"
#include "PndDrcPDPoint.h"
#include "PndDrcEVPoint.h"
#include "PndDrcHit.h"
#include "PndDrcPDHit.h"
#include "PndDrcDigi.h"
#include "PndGeoHandling.h"

#include "PndDrcLutNode.h"

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

class PndPhoGunShortP : public FairTask {

public:

  /** Default constructor **/
  PndPhoGunShortP();

  /** Constructor with verbosity **/
  PndPhoGunShortP(Int_t verbose);

  /** Destructor **/
  virtual ~PndPhoGunShortP();

  virtual InitStatus Init();

  /** Executed task **/ 
  virtual void Exec(Option_t* option);

  /** Finish task **/ 
  virtual void Finish();
 
  void SetOutputFile(TString outName = "lutnode.root"){fOutputName = outName;}
 
private:

  void ProcessPhotonHit();
  void InitLut();
 
  
  Int_t fDetectorID;  

  TClonesArray* fBarPointArray; // DRC MC points in the bars
  TClonesArray* fPDPointArray; // DRC points in the photon detector
  TClonesArray* fPDHitArray; // DRC Hits in the photon detector
  TClonesArray* fMCArray; // DRC MCPoints in the photon detector
  TClonesArray* fEVPointArray; // DRC MCPoints in the photon detector
  TClonesArray *fDigiArray;

  //TClonesArray* fChPhoArray; // output array of PndChPho

  TList* fHistoList;
  TFile*    froot;
  
  PndDrcDigi* fDigi;
  PndGeoDrc* fGeo;                 //!< Basic geometry data of barrel DRC.
  PndGeoHandling* fGeoH;
  
  /** Set the parameters to the default values. **/
  void SetDefaultParameters();
  
  Double_t InBarCoordSyst(TVector3);
  Double_t FindReflectionType(Double_t, Double_t, Double_t, TString);
  
  /** Verbosity level **/
  Int_t fVerbose;
  
  Int_t nevents;
  Int_t ambiguity;
  
  TString fOutputName;
  TString fTableName;
  Double_t ftilt;
  Int_t	   fBarId; 
  
  // parameters:
  Double_t fR;
  Double_t fRBottom;
  Double_t fHThick;
  Double_t fBboxNum;
  Double_t fPipehAngle;
  Double_t fBarBoxGap;
  Double_t fLength;
  Double_t fDphi;
  Double_t fEVlen;
  Double_t fpi;
  Double_t fEVdrop;
 
  Double_t fNoDD, fNoU, fNoB,  fNoR, fNoL;
  Double_t fNoBL, fNoUR, fNoUL, fNoRL, fNoLR, fNoRB, fNoBU, fNoLB, fNoBR,fNoRU,fNoLU;
  Double_t fNoBLU, fNoBRL, fNoBRU, fNoBLR, fNoBUR, fNoBUL, fNoRBL, fNoRLU,fNoRLB, fNoRUL, fNoRBU, fNoURL, fNoULR, fNoTotal, fNoLRB, fNoLRU, fNoLBR,fNoLBU, fNoLUR, fNoLUB, fNoURB, fNoULB, fNoUBR,fNoUBL;
  Double_t fNoBLRU,fNoBRLU,fNoBULB,fNoBURB,fNoUBLR,fNoULBR,fNoURBL,fNoULRB,fNoURLB,fNoUBRL,fNoRULB,fNoRBUL,fNoRBLU,fNoLBUR,fNoLBRU,fNoLURB,fNoLUBR,fNoLRUL,fNoLRLU,fNoLURL,fNoRULR,fNoBULR,fNoBURL,fNoBLUR,fNoBRUL;
  Double_t fNoRBULR,fNoRBLUR,fNoRBULB,fNoRBRUL,fNoBRUBL,fNoBURBL,fNoBRULB,fNoBLRUL,fNoBLUBR,fNoBULBR,fNoBURLB,fNoBLURB,fNoBULRB,fNoLBURL,fNoLBLUR;

  Double_t fNweirdPhotons;
  Double_t flambdah; 
  Double_t fPixIndex; 
  Double_t ftime;
  Double_t fPhiRot;   
  TVector3 fPphoInit; 
  TVector3 fStartVertex;
  TVector3 fPDSec,fEVSec,fEvSec;
  Double_t fPDPhi,fEVPhi,fPhiRotEV;
  Double_t Ang_pipe, Rout1, Rin1, Rin2, PlanB[9], PlanU[9], PlanR[6];  
  Double_t determint1,determint2,determint3,determint4;
  TString ReflectionType,ReflName;
  TArrayD   fmatrixdata;
  Double_t fZin;
  Double_t fLowZ;
  
  Double_t fkxBar, fkyBar, fkzBar;  
  TVector3 fPphoB; // kBar
  
  
  PndDrcEVPoint* EVpt;
  PndDrcEVPoint* EVt;
  PndDrcPDPoint* Ppt;
  PndDrcPDHit* pdhit;
  PndMCTrack* tr;
  TClonesArray *fLut[5];
  
  TFile *fFile; 
  TTree *fTree;
  
  ClassDef(PndPhoGunShortP,1)

};

#endif