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

// -------------------------------------------------------------------------
// -----                      PndLinTrack header file                  -----
// -----                       Created by M. Michel                    -----
// -----                       Modifided by A.Karavdina                -----
// -------------------------------------------------------------------------

#ifndef PNDLINTRACK_H
#define PNDLINTRACK_H

#include "TMatrixDSym.h"
#include "TMatrixTSym.h"
#include "TString.h"
#include "TVector3.h"

#include <stdio.h>
#include <iostream>

class PndLinTrack : public TObject {
 public:
  PndLinTrack();

  /* PndLinTrack(TString detName, Double_t z0, Double_t p0, Double_t p1,
   * Double_t p2, Double_t p3,  */
  /*          Double_t p0err, Double_t p1err, Double_t p2err, Double_t p3err, */
  /*          Double_t chi, Int_t first, Int_t last, Int_t cand); */
  PndLinTrack(TString detName, Double_t p0, Double_t p1, Double_t p2, Double_t p3, Double_t p4, Double_t p5, Double_t chi, Int_t first, Int_t last, Int_t cand);

  virtual ~PndLinTrack();

  void SetDetName(TString name) { fDetName = name; }
  void SetCovarianceMatrix(TMatrixDSym covmatrix)
  {
    for (size_t i = 0; i < 6; i++) {
      for (size_t j = 0; j < 6; j++) {
        fCovMatrix[i][j] = covmatrix[i][j];
      }
    }
  }

  TString GetDetName() const { return fDetName; }
  void GetPar(Double_t *par) const
  {
    par[0] = fP0;
    par[1] = fP1;
    par[2] = fP2;
    par[3] = fP3;
    par[4] = fP4;
    par[5] = fP5;
  }
  /* TVector3   GetStartVec()       const { return TVector3(fP0,fP2,0.);} */
  TVector3 GetStartVec() const { return TVector3(fP0, fP2, fP4); }
  // TVector3   GetStartVec();
  TVector3 GetDirectionVec() const { return TVector3(fP1, fP3, fP5); }
  // TVector3   GetDirectionVec();
  TVector3 GetStartErrVec();
  TVector3 GetDirectionErrVec();
  /* TVector3      GetParErr(Double_t* errpar)  */
  /*   const {
   * errpar[0]=fP0err;errpar[1]=fP1err;errpar[2]=fP2err;errpar[3]=fP3err;errpar[4]=fP4err;}
   */
  Double_t GetChiSquare() const { return fChiSq; }
  Int_t GetFirstHit() const { return fFirst; }
  Int_t GetLastHit() const { return fLast; }
  Int_t GetTCandID() const { return fCandId; }

  void GetParErr(Double_t *errpar);
  TMatrixDSym GetCovarianceMatrix() { return fCovMatrix; }

 private:
  TString fDetName; // Detector name
  //  Double_t fz0;// z-position of first LMD plane
  Double_t fP0, fP1, fP2, fP3, fP4, fP5; // fit-parameter
  // start (P0, P2, P4), direction (P1, P3, P5) straight line
  TMatrixDSym fCovMatrix;
  // Double_t fP0err, fP1err, fP2err, fP3err, fP4err; // errors of parameters
  Double_t fChiSq; // Chi-Square of Fit
  Int_t fFirst, fLast;
  Int_t fCandId; // TrackCand id

  ClassDef(PndLinTrack, 10);
};

#endif