PndApollonius::TripletSolution Struct Reference

#include <PndApolloniusTriplet.h>

Public Types
enum	detID {   MVDpixel, MVDstrip, STT, GEM,   NOTDEFINED }

Public Member Functions
	TripletSolution ()
	TripletSolution (TVector3 track)
void	AddHits (detID detector, std::vector< FairHit *> data)
void	AddHits (std::vector< PndSttHit *> sttHits)
void	AddHit (detID detector, FairHit *data)
int	GetNHits ()
TVector2	HitOnTrack (FairHit *hit)
void	IsClockwise ()
void	SortHits (TVector2 &firstHit, detID detector)
void	SortStt (PndSttGeometryMap *fGeometryMap)
FairHit *	GetFirstHit (PndSttGeometryMap *fGeometryMap)
void	SortAllHits (PndSttGeometryMap *fGeometryMap)

Public Attributes
std::vector< std::string >	detID_Names {"MVD_Pixel", "MVD_Strip", "STT", "GEM", "NOTDEFINED"}
Triplet	fTriplet {nullptr, nullptr, nullptr}
TVector3	fTrack
std::map< detID, std::vector< FairHit * > >	fHits
std::vector< FairHit * >	fAllHits
double	fMeanSquare = 1000000.
bool	fClockwise = true

Friends
std::ostream &	operator<< (std::ostream &output, const TripletSolution &t)

Detailed Description

Definition at line 53 of file PndApolloniusTriplet.h.

Member Enumeration Documentation

detID

enum PndApollonius::TripletSolution::detID

Enumerator
MVDpixel
MVDstrip
STT
GEM
NOTDEFINED

Definition at line 55 of file PndApolloniusTriplet.h.

{ MVDpixel, MVDstrip, STT, GEM, NOTDEFINED };

Constructor & Destructor Documentation

TripletSolution() [1/2]

PndApollonius::TripletSolution::TripletSolution ( )

inline

Definition at line 58 of file PndApolloniusTriplet.h.

{};

TripletSolution() [2/2]

PndApollonius::TripletSolution::TripletSolution ( TVector3  track )

inline

Definition at line 59 of file PndApolloniusTriplet.h.

: fTrack(track){};

Member Function Documentation

AddHit()

void PndApollonius::TripletSolution::AddHit ( detID  detector, FairHit *  data  )

inline

Definition at line 65 of file PndApolloniusTriplet.h.

  {
    if (fHits.count(detector) == 0) {
      fHits[detector];
    }
    if (std::find(fHits[detector].begin(), fHits[detector].end(), data) == fHits[detector].end()) { // add only if hit is not in dataset
      fHits[detector].push_back(data);
    }
  }

AddHits() [1/2]

void PndApollonius::TripletSolution::AddHits ( detID  detector, std::vector< FairHit *>  data  )

inline

Definition at line 61 of file PndApolloniusTriplet.h.

Referenced by PndApollonius::ApolloniusTripletFunctions::FindHitsCloseToCircle().

{ fHits[detector] = data; }

AddHits() [2/2]

void PndApollonius::TripletSolution::AddHits ( std::vector< PndSttHit *>  sttHits )

inline

Definition at line 63 of file PndApolloniusTriplet.h.

{ fHits[STT].insert(fHits[STT].end(), sttHits.begin(), sttHits.end()); }

GetFirstHit()

FairHit* PndApollonius::TripletSolution::GetFirstHit ( PndSttGeometryMap *  fGeometryMap )

inline

Definition at line 147 of file PndApolloniusTriplet.h.

References GEM.

  {
    std::vector<FairHit *> mvdPixelHits = fHits[MVDpixel];
    std::vector<FairHit *> mvdStripHits = fHits[MVDstrip];
    std::vector<FairHit *> sttHits = fHits[STT];
    std::vector<FairHit *> gemHits = fHits[GEM];

    if (mvdPixelHits.size() > 0) {
      double dSquare = 100000.;
      FairHit *first = nullptr;
      for (auto hit : mvdPixelHits) {
        if (((hit->GetX() * hit->GetX()) + (hit->GetY() * hit->GetY() + (hit->GetZ() * hit->GetZ()))) < dSquare) {
          dSquare = ((hit->GetX() * hit->GetX()) + (hit->GetY() * hit->GetY()) + (hit->GetZ() * hit->GetZ()));
          first = hit;
        }
      }
      return first;
    } else if (mvdStripHits.size() > 0) {
      double dSquare = 100000.;
      FairHit *first = nullptr;
      for (auto hit : mvdStripHits) {
        if (((hit->GetX() * hit->GetX()) + (hit->GetY() * hit->GetY() + (hit->GetZ() * hit->GetZ()))) < dSquare) {
          dSquare = ((hit->GetX() * hit->GetX()) + (hit->GetY() * hit->GetY()) + (hit->GetZ() * hit->GetZ()));
          first = hit;
        }
      }
      return first;

    } else if (sttHits.size() > 0) {
      return sttHits.front();
    } else if (gemHits.size() > 0) {
      double dSquare = 100000.;
      FairHit *first = nullptr;
      for (auto hit : gemHits) {
        if (((hit->GetX() * hit->GetX()) + (hit->GetY() * hit->GetY())) < dSquare) {
          dSquare = ((hit->GetX() * hit->GetX()) + (hit->GetY() * hit->GetY()));
          first = hit;
        }
      }
      return first;
    } else
      return nullptr;
  }

GetNHits()

int PndApollonius::TripletSolution::GetNHits ( )

inline

Definition at line 75 of file PndApolloniusTriplet.h.

References GEM.

Referenced by PndApollonius::ApolloniusTripletFunctions::CombineIdenticalSolutions(), and PndApollonius::ApolloniusTripletFunctions::MeanSquareDistance().

{ return fHits[MVDpixel].size() + fHits[MVDstrip].size() + fHits[STT].size() + fHits[GEM].size(); }

HitOnTrack()

TVector2 PndApollonius::TripletSolution::HitOnTrack ( FairHit *  hit )

inline

Definition at line 77 of file PndApolloniusTriplet.h.

{ return TVector2(hit->GetX() - fTrack.X(), hit->GetY() - fTrack.Y()); }

IsClockwise()

void PndApollonius::TripletSolution::IsClockwise ( )

inline

Definition at line 79 of file PndApolloniusTriplet.h.

  {
    std::vector<FairHit *> sttHits = fHits[STT];
    if (sttHits.size() > 2) {

      fClockwise = true;
      std::sort(sttHits.begin(), sttHits.end(), [](FairHit *first, FairHit *second) { return (*(static_cast<PndSttHit *>(first)) < *(static_cast<PndSttHit *>(second))); });
      TVector2 firstHit(HitOnTrack(sttHits.front()));
      TVector2 secondHit(HitOnTrack(sttHits.back()));
      if (firstHit.DeltaPhi(secondHit) < 0) {
        fClockwise = false;
      }
    }
  }

SortAllHits()

void PndApollonius::TripletSolution::SortAllHits ( PndSttGeometryMap *  fGeometryMap )

inline

Definition at line 191 of file PndApolloniusTriplet.h.

References GEM.

  {
    std::vector<FairHit *> sttHits = fHits[STT];
    fAllHits.clear();
    if (sttHits.size() > 2) {
      FairHit *first = GetFirstHit(fGeometryMap);

      TVector2 firstHit(HitOnTrack(first));
      TVector2 secondHit(HitOnTrack(sttHits.back()));
      bool hasChanged = false;
      if (firstHit.DeltaPhi(secondHit) < 0 && fClockwise) {
        hasChanged = true;
        fClockwise = false;
      } else if (firstHit.DeltaPhi(secondHit) >= 0 && !fClockwise) {
        hasChanged = true;
        fClockwise = true;
      }

      if (hasChanged) {
        SortHits(firstHit, MVDpixel);
        SortHits(firstHit, MVDstrip);
        SortHits(firstHit, STT);
        SortHits(firstHit, GEM);
        fAllHits.insert(fAllHits.end(), fHits[MVDpixel].begin(), fHits[MVDpixel].end());
        fAllHits.insert(fAllHits.end(), fHits[MVDstrip].begin(), fHits[MVDstrip].end());
        fAllHits.insert(fAllHits.end(), fHits[STT].begin(), fHits[STT].end());
        fAllHits.insert(fAllHits.end(), fHits[GEM].begin(), fHits[GEM].end());
      } else {
        SortHits(firstHit, MVDpixel);
        SortHits(firstHit, MVDstrip);
        fAllHits.insert(fAllHits.end(), fHits[MVDpixel].begin(), fHits[MVDpixel].end());
        fAllHits.insert(fAllHits.end(), fHits[MVDstrip].begin(), fHits[MVDstrip].end());
        fAllHits.insert(fAllHits.end(), fHits[STT].begin(), fHits[STT].end());
        SortHits(firstHit, GEM);
        fAllHits.insert(fAllHits.end(), fHits[GEM].begin(), fHits[GEM].end());
      }
    }
  }

SortHits()

void PndApollonius::TripletSolution::SortHits ( TVector2 &  firstHit, detID  detector  )

inline

Definition at line 94 of file PndApolloniusTriplet.h.

References CAMath::Pi().

  {
    if (fClockwise) {
      std::sort(fHits[detector].begin(), fHits[detector].end(), [&](FairHit *a, FairHit *b) {
        double phi_a = HitOnTrack(a).DeltaPhi(firstHit);
        double phi_b = HitOnTrack(b).DeltaPhi(firstHit);
        if (phi_a <= 0)
          phi_a += 2 * TMath::Pi();
        if (phi_b <= 0)
          phi_b += 2 * TMath::Pi();

        return (phi_a > phi_b);
      });
    } else {
      std::sort(fHits[detector].begin(), fHits[detector].end(), [&](FairHit *a, FairHit *b) {
        double phi_a = HitOnTrack(a).DeltaPhi(firstHit);
        double phi_b = HitOnTrack(b).DeltaPhi(firstHit);
        if (phi_a < 0)
          phi_a += 2 * TMath::Pi();
        if (phi_b < 0)
          phi_b += 2 * TMath::Pi();
        return (phi_a < phi_b);
      });
    }
  }

SortStt()

void PndApollonius::TripletSolution::SortStt ( PndSttGeometryMap *  fGeometryMap )

inline

Definition at line 120 of file PndApolloniusTriplet.h.

References PndSttGeometryMap::GetRow().

Referenced by PndApollonius::ApolloniusTripletFunctions::FindHitsCloseToCircle(), and PndApollonius::ApolloniusTripletFunctions::IsContinuous().

  {
    std::vector<FairHit *> sttHits = fHits[STT];
    if (sttHits.size() > 2) {

      IsClockwise();

      std::sort(sttHits.begin(), sttHits.end(), [](FairHit *first, FairHit *second) { return (*(static_cast<PndSttHit *>(first)) < *(static_cast<PndSttHit *>(second))); });
      std::vector<FairHit *> possibleFirstHits;
      int minRow = 100;
      for( auto hit : sttHits){
        if (fGeometryMap->GetRow(static_cast<PndSttHit *>(hit)->GetTubeID()) <= minRow) {
          minRow = fGeometryMap->GetRow(static_cast<PndSttHit *>(hit)->GetTubeID());
          possibleFirstHits.push_back(hit);
        } else
          break;
      }
      TVector2 firstHit;
      if (fClockwise) {
        firstHit = HitOnTrack(possibleFirstHits.back());
      } else {
        firstHit = HitOnTrack(possibleFirstHits.front());
      }
      SortHits(firstHit, STT);
    }
  }

Friends And Related Function Documentation

operator<<

std::ostream& operator<< ( std::ostream &  output, const TripletSolution &  t  )

friend

Definition at line 230 of file PndApolloniusTriplet.h.

  {
    output << "Track x/y/r " << t.fTrack.X() << "/" << t.fTrack.Y() << "/" << t.fTrack.Z() << " ";
    if (t.fTriplet.fTripletHits[0] != nullptr) {
      output << " trip: " << t.fTriplet.fTripletHits[0]->GetTubeID() << "/" << t.fTriplet.fTripletHits[1]->GetTubeID() << "/" << t.fTriplet.fTripletHits[2]->GetTubeID() << " ";
    }
    if (t.fHits.count(TripletSolution::detID::MVDpixel) > 0) {
      output << t.detID_Names.at(TripletSolution::detID::MVDpixel) << " : ";
      for (auto hit : t.fHits.at(TripletSolution::detID::MVDpixel)) {
        output << hit->GetEntryNr() << "/";
      }
      output << " ";
    }
    if (t.fHits.count(TripletSolution::detID::MVDstrip) > 0) {
      output << t.detID_Names.at(TripletSolution::detID::MVDstrip) << " : ";
      for (auto hit : t.fHits.at(TripletSolution::detID::MVDstrip)) {
        output << hit->GetEntryNr() << "/";
      }
      output << " ";
    }
    if (t.fHits.count(TripletSolution::detID::STT) > 0) {
      output << t.detID_Names.at(TripletSolution::detID::STT) << " : ";
      for (auto hit : t.fHits.at(TripletSolution::detID::STT)) {
        PndSttHit *myHit = static_cast<PndSttHit *>(hit);
        output << myHit->GetTubeID() << "/";
      }
      output << " ";
    }
    if (t.fHits.count(TripletSolution::detID::GEM) > 0) {
      output << t.detID_Names.at(TripletSolution::detID::GEM) << " : ";
      for (auto hit : t.fHits.at(TripletSolution::detID::GEM)) {
        output << hit->GetEntryNr() << "/";
      }
      output << " ";
    }

    output << "MSD: " << t.fMeanSquare;
    return output;
  }

Member Data Documentation

detID_Names

std::vector<std::string> PndApollonius::TripletSolution::detID_Names {"MVD_Pixel", "MVD_Strip", "STT", "GEM", "NOTDEFINED"}

Definition at line 56 of file PndApolloniusTriplet.h.

fAllHits

std::vector<FairHit *> PndApollonius::TripletSolution::fAllHits

Definition at line 273 of file PndApolloniusTriplet.h.

fClockwise

bool PndApollonius::TripletSolution::fClockwise = true

Definition at line 275 of file PndApolloniusTriplet.h.

fHits

std::map<detID, std::vector<FairHit *> > PndApollonius::TripletSolution::fHits

Definition at line 272 of file PndApolloniusTriplet.h.

Referenced by PndApollonius::ApolloniusTripletFunctions::ContainsTriplet(), PndApollonius::ApolloniusTripletFunctions::FindBestSolutions(), PndApollonius::ApolloniusTripletFunctions::IsContinuous(), and PndApollonius::ApolloniusTripletFunctions::MeanSquareDistance().

fMeanSquare

double PndApollonius::TripletSolution::fMeanSquare = 1000000.

Definition at line 274 of file PndApolloniusTriplet.h.

Referenced by PndApollonius::ApolloniusTripletFunctions::FindBestSolutions(), and PndApollonius::ApolloniusTripletFunctions::FindHitsCloseToCircle().

fTrack

TVector3 PndApollonius::TripletSolution::fTrack

Definition at line 271 of file PndApolloniusTriplet.h.

Referenced by PndApollonius::ApolloniusTripletFunctions::MeanSquareDistance().

fTriplet

Triplet PndApollonius::TripletSolution::fTriplet {nullptr, nullptr, nullptr}

Definition at line 270 of file PndApolloniusTriplet.h.

Referenced by PndApollonius::ApolloniusTripletFunctions::GenerateTripletTracks().

---

The documentation for this struct was generated from the following file:

• PndApolloniusTriplet.h