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

// -------------------------------------------------------------------------
// -----               PndLmdTrackQ header file                  -----
// -----       Created 18/06/2013 by A.Karavdina          -----
// -------------------------------------------------------------------------

#ifndef PNDLMDTRACKQ_H
#define PNDLMDTRACKQ_H

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

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

//#include "PndTrack.h"

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

  PndLmdTrackQ(double Plab);

  ~PndLmdTrackQ();

  void SetTrkRecStatus(int st) { fTrkRecStatus = st; }

  int GetTrkRecStatus() const { return fTrkRecStatus; }

  void SetPDGcode(int pdg) { fPDGcode = pdg; }

  int GetPDGcode() const { return fPDGcode; }

  void SetLMDpoint(double x, double y, double z)
  {
    fXrecLMD = x;
    fYrecLMD = y;
    fZrecLMD = z;
  }

  void GetLMDpoint(TVector3 &lmdp) { lmdp.SetXYZ(fXrecLMD, fYrecLMD, fZrecLMD); }

  void SetLMDdir(double theta, double phi)
  {
    fThetarecLMD = theta;
    fPhirecLMD = phi;
  }

  double GetLMDtheta() const { return fThetarecLMD; }

  double GetLMDphi() const { return fPhirecLMD; }

  void SetLMDchi2(double chi2) { fchi2LMD = chi2; }

  double GetLMDchi2() const { return fchi2LMD; }

  void SetIPpoint(double x, double y, double z)
  {
    fXrec = x;
    fYrec = y;
    fZrec = z;
  }

  void GetIPpoint(TVector3 &lmdp) { lmdp.SetXYZ(fXrec, fYrec, fZrec); }

  void SetIPmom(double theta, double phi, double mom)
  {
    fThetarec = theta;
    fPhirec = phi;
    fMomrec = mom;
  }

  double GetIPtheta() const { return fThetarec; }

  double GetIPphi() const { return fPhirec; }

  double GetIPmom() const { return fMomrec; }

  void SetIPerrpoint(double errx, double erry, double errz)
  {
    ferrXrec = errx;
    ferrYrec = erry;
    ferrZrec = errz;
  }

  void GetIPerrpoint(TVector3 &lmdp) { lmdp.SetXYZ(ferrXrec, ferrYrec, ferrZrec); }

  void SetIPerrmom(double errtheta, double errphi, double errmom)
  {
    ferrThetarec = errtheta;
    ferrPhirec = errphi;
    ferrMomrec = errmom;
  }

  void GetIPerrmom(TVector3 &lmerrdp) { lmerrdp.SetXYZ(ferrThetarec, ferrPhirec, ferrMomrec); }

  void SetMCpoint(double x, double y, double z)
  {
    fXmc = x;
    fYmc = y;
    fZmc = z;
  }

  void GetMCpoint(TVector3 &lmdp) { lmdp.SetXYZ(fXmc, fYmc, fZmc); }

  void SetMCmom(double theta, double phi, double mom)
  {
    fThetamc = theta;
    fPhimc = phi;
    fMommc = mom;
  }

  double GetMCtheta() const { return fThetamc; }

  double GetMCphi() const { return fPhimc; }

  double GetMCmom() const { return fMommc; }

  void GetMCpointLMD(TVector3 &lmdp) { lmdp.SetXYZ(fXmcLMD, fYmcLMD, fZmcLMD); }

  void SetMCpointLMD(double x, double y, double z)
  {
    fXmcLMD = x;
    fYmcLMD = y;
    fZmcLMD = z;
  }

  void SetMCmomLMD(double theta, double phi, double mom)
  {
    fThetamcLMD = theta;
    fPhimcLMD = phi;
    fMommcLMD = mom;
  }

  double GetMCthetaLMD() const { return fThetamcLMD; }

  double GetMCphiLMD() const { return fPhimcLMD; }

  double GetMCmomLMD() const { return fMommcLMD; }

  void SetSecondary(int sec) { fSecondary = sec; }

  int GetSecondary() const { return fSecondary; }

  void SetNumMChits(int num) { fnumMChits = num; }

  int GetNumMChits() const { return fnumMChits; }

  void SetNumDoubleMChits(int num) { fnumMCdoublehits = num; }

  int GetNumDoubleMChits() const { return fnumMCdoublehits; }

  void SetSumEvPDG(int sumid) { fsumID = sumid; }

  int GetSumEvPDG() const { return fsumID; }

  /* void SetTotEvCharge(int totch){ ftotCharge = totch;} */
  /* int GetTotEvCharge() const{return ftotCharge;} */

  void SetEvMCMulti(int tot) { fnumEvparticles = tot; }

  int GetEvMCMulti() const { return fnumEvparticles; }

  void SetEvRECMulti(int tot) { fnumEvRECparticles = tot; }

  int GetEvRECMulti() const { return fnumEvRECparticles; }

  void SetEvTime(double evtm) { fEvTime = evtm; }
  void SetTrkTime(double trktm) { fTrkTime = trktm; }
  double GetEvTime() const { return fEvTime; }
  double GetTrkTime() const { return fTrkTime; }

  void SetModule(int mod) { fModule = mod; }
  void SetHalf(int hf) { fHalf = hf; }
  int GetModule() const { return fModule; }
  int GetHalf() const { return fHalf; }

 private:
  //  Int_t ftotCharge;
  Int_t fsumID, fnumEvparticles, fnumEvRECparticles;
  Double_t fXrecLMD, fYrecLMD, fZrecLMD, fThetarecLMD, fPhirecLMD;
  Double_t fXrec, fYrec, fZrec, fThetarec, fPhirec, fMomrec;
  Double_t ferrXrec, ferrYrec, ferrZrec, ferrThetarec, ferrPhirec, ferrMomrec;
  Double_t fXmc, fYmc, fZmc, fThetamc, fPhimc, fMommc;
  Double_t fXmcLMD, fYmcLMD, fZmcLMD, fThetamcLMD, fPhimcLMD, fMommcLMD;
  Double_t fchi2LMD;
  int fTrkRecStatus;    // 0 = good,
  int fSecondary;       // 0=primary, 1=secondary
  int fPDGcode;         // PDG code particle made the main contribution (65 % of hits) to the rec.track, -9999=n/a
  int fnumMChits;       // number of MC hits for MC trk
  int fnumMCdoublehits; // number of doubled [diff side of one plane] MC hits for MC trk
                        // >0 = ghost: 1 = hit mixture, 2 = doubled trk
                        // <0 = missed: -1=trk-search losses,-2=little amount of hits,-3=no MC hits, -10=trk wasn't back-propag
  double fEvTime;       // global time for event
  double fTrkTime;      // time for 1st hit in the Trk = ToF (flight time from IP to LMD)

  int fModule; // module number of 1st hit
  int fHalf;   // half of the plane number of 1st hit
  ClassDef(PndLmdTrackQ, 11);
};

#endif