PndDrcLutFill Class Reference

#include <PndDrcLutFill.h>

Inheritance diagram for PndDrcLutFill:

Public Member Functions
	PndDrcLutFill ()
	PndDrcLutFill (Int_t verbose)
	PndDrcLutFill (Int_t verbose, TString outfilename)
virtual	~PndDrcLutFill ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *option)
virtual void	Finish ()
void	InitLut ()
void	SetOutputFile (TString outfilename="luttab.root")

Private Member Functions
void	ProcessPhotonHit ()
void	SetDefaultParameters ()

Private Attributes
PndGeoDrc *	fGeo
Int_t	fDetectorID
Double_t	fBboxNum
Double_t	fPipehAngle
Double_t	fDphi
TClonesArray *	fMCArray
TClonesArray *	fBarPointArray
TClonesArray *	fPDPointArray
TClonesArray *	fEVPointArray
TClonesArray *	fDigiArray
TClonesArray *	fPDHitArray
TClonesArray *	fLut [5]
TFile *	fFile
TTree *	fTree
PndMCTrack *	fMCTrack
PndDrcBarPoint *	fBarPoint
PndDrcPDPoint *	fPDPoint
PndDrcDigi *	fDigi
PndDrcPDHit *	fPDHit
PndDrcEVPoint *	fEVPoint
Int_t	fVerbose
Int_t	nevents
TString	fOutputFile

Detailed Description

Definition at line 26 of file PndDrcLutFill.h.

Constructor & Destructor Documentation

PndDrcLutFill::PndDrcLutFill

(

)

Definition at line 22 of file PndDrcLutFill.cxx.

References fOutputFile.

:FairTask("PndDrcLutFill")
{
  fOutputFile = "luttab.root";
}

PndDrcLutFill::PndDrcLutFill

(

Int_t

verbose

)

Definition at line 29 of file PndDrcLutFill.cxx.

References fOutputFile, fVerbose, and verbose.

  :FairTask("PndDrcLutFill",verbose)
{
  fVerbose = verbose;  
  fOutputFile = "luttab.root";
}

PndDrcLutFill::PndDrcLutFill	(	Int_t	verbose,
		TString	outfilename
	)

Definition at line 36 of file PndDrcLutFill.cxx.

References fOutputFile, fVerbose, and verbose.

  :FairTask("PndDrcLutFill",verbose)
{
  fVerbose = verbose;  
  fOutputFile = outfilename;
}

PndDrcLutFill::~PndDrcLutFill ( )

virtual

Definition at line 44 of file PndDrcLutFill.cxx.

{

}

Member Function Documentation

void PndDrcLutFill::Exec ( Option_t *  option )

virtual

Definition at line 133 of file PndDrcLutFill.cxx.

References fDetectorID, nevents, and ProcessPhotonHit().

{
  nevents++;
  // if(fVerbose>0 && nevents%1000==0)
    std::cout<<"Event # "<< nevents<<std::endl;  
  fDetectorID = 0;
  ProcessPhotonHit();
}

void PndDrcLutFill::Finish ( )

virtual

Definition at line 235 of file PndDrcLutFill.cxx.

References fFile, fLut, and fTree.

{
    fTree->Fill();
    fTree->Write();
    fFile->Write();
 
    for(Int_t l=0; l<5; l++){
      fLut[l]->Clear(); 
    }
    cout << "-I- PndDrcLutFill: Finish" << endl; 
}

InitStatus PndDrcLutFill::Init ( )

virtual

Definition at line 50 of file PndDrcLutFill.cxx.

References PndGeoDrc::BBoxNum(), Double_t, fBarPointArray, fBboxNum, fDigiArray, fDphi, fEVPointArray, fFile, fGeo, fLut, fMCArray, fOutputFile, fPDHitArray, fPDPointArray, fPipehAngle, fTree, InitLut(), nevents, and PndGeoDrc::PipehAngle().

{
  cout << " ---------- INITIALIZATION ------------" << endl;
  nevents = 0;
  // Get RootManager
  FairRootManager* ioman = FairRootManager::Instance();
  if ( ! ioman ) {
    cout << "-E- PndDrcLutFill::Init: " << "RootManager not instantiated!" << endl;
    return kFATAL;
  }

  // Get input array
  fMCArray = (TClonesArray*) ioman->GetObject("MCTrack");
  if ( ! fMCArray ) {
    cout << "-W- PndDrcLutFill::Init: " << "No MCTrack array!" << endl;
    return kERROR;
  } 

  // Get bar points array
  fBarPointArray = (TClonesArray*) ioman->GetObject("DrcBarPoint");
  if ( ! fBarPointArray ) {
    cout << "-W- PndDrcLutReco::Init: " << "No DrcBarPoint array!" << endl;
    return kERROR;
  }

  // Get Photon point array
  fPDPointArray = (TClonesArray*) ioman->GetObject("DrcPDPoint");
  if ( ! fPDPointArray ) {
    cout << "-W- PndDrcLutFill::Init: " << "No DrcPDPoint array!" << endl;
    return kERROR;
  }  

  // Get ev points array
  fEVPointArray = (TClonesArray*) ioman->GetObject("DrcEVPoint");
  if ( ! fEVPointArray ) {
    cout << "-W- PndDrcLutFill::Init: " << "No DrcEVPoint array!" << endl;
    return kERROR;
  }

  // Get digi array
  fDigiArray = (TClonesArray*) ioman->GetObject("DrcDigi");
  if ( ! fDigiArray ) {
    cout << "-W- PndDrcLutFill::Init: " << "No DrcDigi array!" << endl;
    return kERROR;
  }
  // Get input array
  fPDHitArray = (TClonesArray*) ioman->GetObject("DrcPDHit");
  if ( ! fPDHitArray ) {
    cout << "-W- PndDrcLutFill::Init: " << "No DrcPDHit array!" << endl;
    return kERROR;
  }
 
  fFile = TFile::Open(fOutputFile,"RECREATE");
  fTree = new TTree("dircsim","Look-up table for DIRC"); 
  for(Int_t l=0; l<5; l++){
    fLut[l] = new TClonesArray("PndDrcLutNode");
    fTree->Branch(Form("LUT%d",l),&fLut[l],256000,0); 
  }
 
  InitLut();

  fGeo = new PndGeoDrc();
  fBboxNum    = fGeo->BBoxNum();
  fPipehAngle = fGeo->PipehAngle();
  fDphi       = 2.*(180. - 2*fPipehAngle)/(Double_t)fGeo->BBoxNum();

  cout << "-I- PndDrcLutFill: Intialization successfull" << endl;
  return kSUCCESS;
   
}

void PndDrcLutFill::InitLut

(

)

Definition at line 121 of file PndDrcLutFill.cxx.

References fLut, and n.

Referenced by Init().

{
  Int_t Nnodes = 30000;
  for(Int_t l=0; l<5; l++){
    TClonesArray &fLuta = *fLut[l]; 
    for (Long64_t n=0; n<Nnodes; n++) {
      new((fLuta)[n]) PndDrcLutNode(-1);
    }
  }
}

void PndDrcLutFill::ProcessPhotonHit ( )

private

Definition at line 143 of file PndDrcLutFill.cxx.

References At, Double_t, fBarPoint, fBarPointArray, fEVPoint, fEVPointArray, fLut, fMCArray, fMCTrack, fPDHit, fPDHitArray, fPDPoint, fPDPointArray, PndDrcBarPoint::GetBarId(), PndDrcPDPoint::GetBarPointID(), PndMCTrack::GetMomentum(), PndDrcEVPoint::GetNormal(), PndDrcPDHit::GetPosition(), PndDrcPDHit::GetSensorId(), i, Pi, and vec.

Referenced by Exec().

{
  // LUT was generated for:
  //Int_t  lutboxId=3; //[R.K. 01/2017] unused variable?
  Double_t  lutboxPhi=10.825;

  //Int_t nofChPho = 0; //[R.K. 01/2017] unused variable?
  Double_t path;//, barPhi; //[R.K. 01/2017] unused variable?
  TVector3 dir, dirm, vec, posInBar;

  // Loop over PndDrcPDHits
  for(Int_t k=0; k<fPDHitArray->GetEntriesFast(); k++) {
    fPDHit = (PndDrcPDHit*)fPDHitArray->At(k);
    Int_t pointID = fPDHit->GetLink(1).GetIndex();
    Int_t sensorId = fPDHit->GetSensorId();
    if(pointID==-1) continue;
    
    // Int_t digiID= fPDHit->GetRefIndex();
    // fDigi = (PndDrcDigi*) fDigiArray->At(digiID);
    // Int_t sensorId = fDigi->GetSensorId();
    // Int_t pointID= fDigi->GetIndex(0);
     
    fPDPoint = (PndDrcPDPoint*)fPDPointArray->At(pointID);
    fBarPoint= (PndDrcBarPoint*)fBarPointArray->At(fPDPoint->GetBarPointID());
    Int_t barId = fBarPoint->GetBarId();
 
    Int_t trackID = fPDPoint->GetTrackID();
    Double_t time = fPDPoint->GetTime();
    Int_t nev=0; 
    path=0;
    for(int i=0; i<fEVPointArray->GetEntriesFast(); i++){
      fEVPoint = (PndDrcEVPoint*)fEVPointArray->At(i);
      if(trackID == fEVPoint->GetTrackID()){
        // if(fEVPoint->GetDetectorID()==9375) {
        //   // get the direction of the phothon at the ent of the bar
        //   fEVPoint->Momentum(dirm);
        // }
        nev++;
        vec = fEVPoint->GetNormal();
        path += (vec.X()+vec.Y()*10 + vec.Z()*100)*1000*nev;
      }
    }
    
    fMCTrack = (PndMCTrack*)fMCArray->At(trackID);
    dir =  fMCTrack->GetMomentum().Unit();
    dir.RotateZ(-lutboxPhi/180.*TMath::Pi());

    // sensorId =  (sensorId/100 + lutboxId*17)*100 + sensorId%100; 
 
    if(sensorId>30000 || sensorId<0) {
      std::cout<<"WTQ  fPDHit->GetDetectorID()   "<<sensorId <<std::endl;
      continue;
    }
    
    
    // //======================
    // posInBar = fMCTrack->GetStartVertex();
    // Double_t phi = posInBar.Phi()/TMath::Pi()*180;
    // if(phi < 0) phi = 360 + phi;
    // if(phi >= 0 && phi < 90) barPhi = TMath::Floor(phi/fDphi) *fDphi + fDphi/2.;
    // if(phi >= 90 && phi < 270) barPhi = 90  + fPipehAngle + TMath::Floor((phi-90-fPipehAngle)/fDphi) *fDphi + fDphi/2.;
    // if(phi >= 270 && phi < 360) barPhi = 270 + fPipehAngle + TMath::Floor((phi-270-fPipehAngle)/fDphi) *fDphi + fDphi/2.;

    // dirm.RotateZ(-barPhi/180.*TMath::Pi());
    // TVector3 fnX1 = TVector3 (1,0,0);   
    // TVector3 fnY1 = TVector3( 0,1,0); 
    // double criticalAngle = asin(1.00028/fGeo->nQuartz()); 
    // for(int u=0; u<8; u++){
    //   if(u == 0) dir = dirm;
    //   if(u == 1) dir.SetXYZ( dirm.X(), dirm.Y(),-dirm.Z());
    //   if(u == 2) dir.SetXYZ( dirm.X(),-dirm.Y(), dirm.Z());
    //   if(u == 3) dir.SetXYZ(-dirm.X(), dirm.Y(), dirm.Z());
    //   if(u == 4) dir.SetXYZ(-dirm.X(),-dirm.Y(), dirm.Z());
    //   if(u == 5) dir.SetXYZ(-dirm.X(), dirm.Y(),-dirm.Z());
    //   if(u == 6) dir.SetXYZ( dirm.X(),-dirm.Y(),-dirm.Z());
    //   if(u == 7) dir = -dirm;

    //   if(dir.Angle(fnX1) < criticalAngle || dir.Angle(fnY1) < criticalAngle){
    //  //std::cout<<"Wrong combination "<<std::endl;
    //  continue;
    //   }
    //   ((PndDrcLutNode*)(fLut->At(sensorId)))->AddEntry(dir);
    //   ((PndDrcLutNode*)(fLut->At(sensorId)))->AddPathId(pathid);
    // }
    // //======================
    
    ((PndDrcLutNode*)(fLut[barId]->At(sensorId)))->AddEntry(sensorId, dir,path,0,time,fPDHit->GetPosition()); //fDigi->GetDetectorId()
  }
}

void PndDrcLutFill::SetDefaultParameters ( )

private

void PndDrcLutFill::SetOutputFile ( TString  outfilename = "luttab.root" )

inline

Definition at line 51 of file PndDrcLutFill.h.

References fOutputFile.

{fOutputFile = outfilename;}

Member Data Documentation

PndDrcBarPoint* PndDrcLutFill::fBarPoint

private

Definition at line 73 of file PndDrcLutFill.h.

Referenced by ProcessPhotonHit().

TClonesArray* PndDrcLutFill::fBarPointArray

private

Definition at line 62 of file PndDrcLutFill.h.

Referenced by Init(), and ProcessPhotonHit().

Double_t PndDrcLutFill::fBboxNum

private

Definition at line 59 of file PndDrcLutFill.h.

Referenced by Init().

Int_t PndDrcLutFill::fDetectorID

private

Definition at line 58 of file PndDrcLutFill.h.

Referenced by Exec().

PndDrcDigi* PndDrcLutFill::fDigi

private

Definition at line 75 of file PndDrcLutFill.h.

TClonesArray* PndDrcLutFill::fDigiArray

private

Definition at line 65 of file PndDrcLutFill.h.

Referenced by Init().

Double_t PndDrcLutFill::fDphi

private

Definition at line 59 of file PndDrcLutFill.h.

Referenced by Init().

PndDrcEVPoint* PndDrcLutFill::fEVPoint

private

Definition at line 77 of file PndDrcLutFill.h.

Referenced by ProcessPhotonHit().

TClonesArray* PndDrcLutFill::fEVPointArray

private

Definition at line 64 of file PndDrcLutFill.h.

Referenced by Init(), and ProcessPhotonHit().

TFile* PndDrcLutFill::fFile

private

Definition at line 69 of file PndDrcLutFill.h.

Referenced by Finish(), and Init().

PndGeoDrc* PndDrcLutFill::fGeo

private

Definition at line 57 of file PndDrcLutFill.h.

Referenced by Init().

TClonesArray* PndDrcLutFill::fLut[5]

private

Definition at line 67 of file PndDrcLutFill.h.

Referenced by Finish(), Init(), InitLut(), and ProcessPhotonHit().

TClonesArray* PndDrcLutFill::fMCArray

private

Definition at line 61 of file PndDrcLutFill.h.

Referenced by Init(), and ProcessPhotonHit().

PndMCTrack* PndDrcLutFill::fMCTrack

private

Definition at line 72 of file PndDrcLutFill.h.

Referenced by ProcessPhotonHit().

TString PndDrcLutFill::fOutputFile

private

Definition at line 85 of file PndDrcLutFill.h.

Referenced by Init(), PndDrcLutFill(), and SetOutputFile().

PndDrcPDHit* PndDrcLutFill::fPDHit

private

Definition at line 76 of file PndDrcLutFill.h.

Referenced by ProcessPhotonHit().

TClonesArray* PndDrcLutFill::fPDHitArray

private

Definition at line 66 of file PndDrcLutFill.h.

Referenced by Init(), and ProcessPhotonHit().

PndDrcPDPoint* PndDrcLutFill::fPDPoint

private

Definition at line 74 of file PndDrcLutFill.h.

Referenced by ProcessPhotonHit().

TClonesArray* PndDrcLutFill::fPDPointArray

private

Definition at line 63 of file PndDrcLutFill.h.

Referenced by Init(), and ProcessPhotonHit().

Double_t PndDrcLutFill::fPipehAngle

private

Definition at line 59 of file PndDrcLutFill.h.

Referenced by Init().

TTree* PndDrcLutFill::fTree

private

Definition at line 70 of file PndDrcLutFill.h.

Referenced by Finish(), and Init().

Int_t PndDrcLutFill::fVerbose

private

Definition at line 83 of file PndDrcLutFill.h.

Referenced by PndDrcLutFill().

Int_t PndDrcLutFill::nevents

private

Definition at line 84 of file PndDrcLutFill.h.

Referenced by Exec(), and Init().

---

The documentation for this class was generated from the following files:

• PndDrcLutFill.h
• PndDrcLutFill.cxx