PndHoughUtilities.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".              *
//****************************************************************************

//  PndHoughUtilities
//  This class contains different methods to handle PndTracks vs Circles

#ifndef PndHoughUtilities_H_
#define PndHoughUtilities_H_
#include "FairTask.h"
#include "FairRootManager.h"
//#include "PndGeoSttPar.h"
//#include "PndSttGeometryMap.h"
//#include "PndSttStrawMap.h"
//#include "PndSttCellTrackFinder.h"
#include "PndSttHit.h"
#include "PndMCTrack.h"

//#include "TEllipse.h"
#include "TH2F.h"
#include "TGraph.h"
//#include "TCanvas.h" // DRAW

#include <vector>
#include <tuple>
//#include "TDatabasePDG.h"
#include "PndHoughData.h"

class PndHoughUtilities {

 public:
  PndHoughUtilities(PndHoughData *fData) : fMapFairLinktoFairHit(fData->GetMapFairLinktoFairHit()), fMapFairLinktoIsochrone(fData->GetMapFairLinktoIsochrone())
  {
    ioman = FairRootManager::Instance();
  };

  PndHoughUtilities() { ioman = FairRootManager::Instance(); };

  virtual ~PndHoughUtilities(){};

  std::vector<double> getCircleFromPndTrack(PndTrack &tr);
  double Pt(double B, double r);
  std::tuple<TVector3, Int_t> getPandqforHit(int i, double B, TVector3 &circle, PndTrackCand &myCand);
  FairTrackParP getTrackParPForHit(Int_t i, Double_t B, PndTrackCand &myCand, TVector3 &circle);
  PndTrack getPndTrack(Double_t B, PndTrackCand &myCand, TVector3 &circle);

  TVector2 calcPointOnCircle(Int_t i, PndTrackCand &myCand, std::vector<double> &circle);
  TVector2 findCorrectIntersectionPoint(std::vector<double> &intersectionPoints, TVector2 &point);
  std::vector<double> calcIntersectionPointCircleLine(std::vector<double> &circle, double m, double b, double Ax);

  double calcDistanceTrackCandToTrack(PndTrackCand &trackCand, std::vector<double> &circle);
  double calcDistanceFromIsochroneToCircle(double &xi, double &yi, double &ri, double &xT, double &yT, double &rT);

 private:
  FairRootManager *ioman = nullptr;

  std::map<FairLink, FairHit *> fMapFairLinktoFairHit;
  std::map<FairLink, double> fMapFairLinktoIsochrone;

  ClassDef(PndHoughUtilities, 1);
};

#endif /*PndHoughUtilities_H_*/