MpdCellTrack.cxx 3.4 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
  1. /// \class MpdCellTrack
  2. ///
  3. /// Cellular automaton track object for the MPD inner tracking system
  4. /// \author Alexander Zinchenko, Maxim Strelchenko (LHEP, JINR, Dubna)
  5. #include "MpdCellTrack.h"
  6. #include <map>
  7. //__________________________________________________________________________
  8. MpdCellTrack::MpdCellTrack()
  9. : TObject(),
  10. fDetectorID(0),
  11. fFlag(1),
  12. fLength(0.),
  13. fIndex(0),
  14. fHitType(kFixedP),
  15. fNofDim(2),
  16. fSignal(0.),
  17. fDist(0.),
  18. fCode("")
  19. {
  20. /// Default constructor
  21. //for (Int_t i = 0; i < 2; ++i) fMeas[i] = fErr[i] = fCosSin[i] = 0;
  22. // fCosSin[0] = 1.;
  23. }
  24. //__________________________________________________________________________
  25. MpdCellTrack::MpdCellTrack(Int_t detID, Int_t nDim, HitType hitType, TVector3 &meas, Double_t *err,
  26. Double_t *cossin, Double_t signal, Double_t dist,
  27. Int_t index, Int_t index1, Int_t trackNo)
  28. : TObject(),
  29. fDetectorID(detID),
  30. fFlag(1),
  31. fLength(0.),
  32. fHitType(hitType),
  33. fNofDim(nDim),
  34. fSignal(signal),
  35. fDist(dist),
  36. fMeas(meas),
  37. fTrackNo(trackNo),
  38. fCode("")
  39. {
  40. /// Constructor
  41. for (Int_t i = 0; i < 2; ++i) {
  42. //fMeas[i] = meas[i];
  43. fErr[i] = err[i];
  44. fCosSin[i] = cossin[i];
  45. }
  46. SetIndex(index);
  47. SetIndex(index1);
  48. }
  49. //__________________________________________________________________________
  50. MpdCellTrack::MpdCellTrack (const MpdCellTrack& track)
  51. : TObject(track),
  52. fDetectorID(track.fDetectorID),
  53. fFlag(track.fFlag),
  54. fLength(track.fLength),
  55. fIndex(track.fIndex),
  56. fHitType(track.fHitType),
  57. fNofDim(track.fNofDim),
  58. fSignal(track.fSignal),
  59. fDist(track.fDist),
  60. fMeas(track.fMeas),
  61. fTrackNo(track.fTrackNo),
  62. fCode(track.fCode)
  63. {
  64. ///copy constructor
  65. for (Int_t i = 0; i < 2; ++i) {
  66. //fMeas[i] = hit.fMeas[i];
  67. fErr[i] = track.fErr[i];
  68. fCosSin[i] = track.fCosSin[i];
  69. }
  70. }
  71. //__________________________________________________________________________
  72. MpdCellTrack::~MpdCellTrack()
  73. {
  74. /// Destructor
  75. }
  76. //__________________________________________________________________________
  77. Int_t MpdCellTrack::Compare(const TObject* hit) const
  78. {
  79. /// "Compare" function to sort in descending order in fDist
  80. MpdCellTrack *kHit = (MpdCellTrack*) hit;
  81. if (kHit->GetType() == GetType()) {
  82. // Check layers
  83. if (GetLayer() < kHit->GetLayer()) return 1;
  84. else if (GetLayer() > kHit->GetLayer()) return -1;
  85. if (GetType() == kFixedP) {
  86. // Sort according to sector number
  87. if (GetDetectorID() % 1000000 < kHit->GetDetectorID() % 1000000) return -1;
  88. else if (GetDetectorID() % 1000000 > kHit->GetDetectorID() % 1000000) return 1;
  89. }
  90. }
  91. if (TMath::Abs(GetDist()) < TMath::Abs(kHit->GetDist())) return 1;
  92. else if (TMath::Abs(GetDist()) > TMath::Abs(kHit->GetDist())) return -1;
  93. return 0;
  94. }
  95. //__________________________________________________________________________
  96. void MpdCellTrack::Print(Option_t *opt)
  97. {
  98. /// Print hit info
  99. }
  100. //__________________________________________________________________________
  101. Double_t MpdCellTrack::GetPos() const
  102. {
  103. /// Distance to (0,0,0)
  104. /*
  105. if (fHitType != kFixedP) return fDist;
  106. printf(" !!! Not implemented for kFixedP hits. Exit. \n");
  107. exit(0);
  108. */
  109. return fDist;
  110. }
  111. //__________________________________________________________________________
  112. void MpdCellTrack::SetIndex(Int_t indx)
  113. {
  114. /// Add point index
  115. Int_t size = fIndex.GetSize();
  116. fIndex.Set (size + 1);
  117. fIndex[size] = indx;
  118. }
  119. ClassImp(MpdCellTrack);