PndPhoGunShortP.h

Go to the documentation of this file.

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

#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:
  PndPhoGunShortP();

  PndPhoGunShortP(Int_t verbose);

  virtual ~PndPhoGunShortP();

  virtual InitStatus Init();

  virtual void Exec(Option_t *option);

  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; 
  PndGeoHandling *fGeoH;

  void SetDefaultParameters();

  Double_t InBarCoordSyst(TVector3);
  Double_t FindReflectionType(Double_t, Double_t, Double_t, TString);

  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