PndDrcLutReco.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

// -----------------------------------------
// PndDrcLutReco.h
//
// Created on: 13.07.2013
// Author: R.Dzhygadlo at gsi.de
// -----------------------------------------
// Class for reconstruction in DIRC using look-up table method

#ifndef PNDDRCLUTRECO_H
#define PNDDRCLUTRECO_H

#include "FairTask.h"

#include "PndMCTrack.h"
#include "PndDrcBarPoint.h"
#include "PndDrcEVPoint.h"
#include "PndDrcPDPoint.h"
#include "PndDrcHit.h"
#include "PndDrcPDHit.h"
#include "PndDrcDigi.h"
#include "PndGeoDrc.h"
#include "PndDrcTrackInfo.h"
#include "PndDrcPhotonInfo.h"
#include "PndDrcAmbiguityInfo.h"
#include "PndPersistencyTask.h"

#include "TClonesArray.h"
#include "TString.h"
#include "TFile.h"
#include "TH1.h"
#include "TF1.h"
#include "TSpectrum.h"
#include "TTree.h"

class PndDrcLutReco : public PndPersistencyTask {

 public:
  // Default constructor
  PndDrcLutReco();

  // Standard constructors
  PndDrcLutReco(Int_t verbose);
  PndDrcLutReco(Int_t verbose, TString infilename);

  // Destructor
  virtual ~PndDrcLutReco();

  virtual InitStatus Init();

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

  // Finish task
  virtual void Finish();

  void SetOutputFile(TString infilename = "luttab.root") { fInputFile = infilename; }

 private:
  void LoopOverMcTracks();
  void FillAmbiguities(PndDrcPhotonInfo *photoninfo, Int_t barId, Int_t recalculatedSensorId, Double_t directz, Double_t barHitTime);
  void DetermineCherenkov(PndDrcTrackInfo *trackinfo, Int_t boxId);
  void DetermineBarId(Double_t phi, Double_t &boxPhi, Int_t &boxId, Int_t &barId);
  Double_t FindPeak();
  Int_t FindPdg(Double_t mom, Double_t cangle);
  PndGeoDrc *fGeo;
  Int_t fDetectorID;
  Double_t fBboxNum, fPipehAngle, fDphi, fBarPhi;

  TClonesArray *fMCArray; // DRC MCPoints in the photon detector
  TClonesArray *fBarPointArray;
  TClonesArray *fEVPointArray;
  TClonesArray *fPDPointArray; // DRC points in the photon detector
  TClonesArray *fDigiArray;
  TClonesArray *fPDHitArray; // DRC Hits in the photon detector
  TClonesArray *fLut[5];
  TClonesArray *fDrcTrackInfoArray;

  TFile *fFile;
  TTree *fTree;

  PndMCTrack *fMCTrack;
  PndDrcBarPoint *fBarPoint;
  PndDrcEVPoint *fEVPoint;
  PndDrcPDPoint *fPDPoint;
  PndDrcDigi *fDigi;
  PndDrcPDHit *fPDHit;

  // Set the parameters to the default values.
  void SetDefaultParameters();

  // Verbosity level
  Int_t fVerbose;
  Int_t nevents;
  TString fInputFile;
  Int_t fEvType, fRadType, fLensType;
  TH1F *fHist;
  TH1F *fHist2;
  TF1 *fFit;
  TSpectrum *fSpect;

  ClassDef(PndDrcLutReco, 1)
};

#endif