LitTrackFinderNNVecMuon.h

Go to the documentation of this file.

/* Copyright (C) 2009-2012 GSI/JINR-LIT, Darmstadt/Dubna
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Andrey Lebedev [committer] */
 
#ifndef LITTRACKFINDERNNVECMUON_H_
#define LITTRACKFINDERNNVECMUON_H_
 
#include "LitDetectorLayoutMuon.h"
#include "LitHitDataMuon.h"
 
#include <map>
#include <vector>
 
namespace lit
{
  namespace parallel
  {
 
    class LitTrackFinderNNVecMuon {
     public:
      LitTrackFinderNNVecMuon();
 
      virtual ~LitTrackFinderNNVecMuon();
 
      virtual void DoFind(const PixelHitArray& hits, const TrackArray& trackSeeds, TrackArray& tracks);
 
      void SetDetectorLayout(const LitDetectorLayoutMuon<fvec>& layout)
      {
        fLayout = layout;
        fHitData.SetDetectorLayout(layout);
      }
 
     private:
      void ArrangeHits(const PixelHitArray& hits);
 
      /* Initializes the track seeds
     * @param itBegin Iterator to the first track seed.
     * @param itEnd iterator to the last track seed. */
      void InitTrackSeeds(const TrackArray& trackSeeds);
 
      /* Follows tracks through the detector
     * @param itBegin Iterator to the first track.
     * @param itEnd Iterator to the last track. */
      void FollowTracks();
 
      /* TODO: Add comments
     *
     */
      void PropagateThroughAbsorber(const std::vector<unsigned int>& tracksId1, const LitAbsorber<fvec>& absorber);
 
      /* TODO: Add comments
     *
     */
      inline void PropagateThroughAbsorber(const TrackArray& tracks, const LitAbsorber<fvec>& absorber);
 
      /*
     * TODO Add comments
     */
      void CollectHits(std::vector<LitTrackParamScal>& par, LitScalTrack* track, unsigned char stationGroup,
                       unsigned char station, unsigned char nofSubstations);
 
      /*
     * TODO Add comments
     */
      void ProcessStation(const std::vector<unsigned int>& tracksId1, unsigned char stationGroup,
                          unsigned char station);
 
      /* TODO: Add comment
     *
     */
      inline void ProcessStation(const TrackArray& tracks, unsigned char stationGroup, unsigned char station);
 
      /* Adds the nearest hit to the track.
     * @param track Pointer to the track.
     * @param hits Reference to hit vector from which the nearest hit will be attached.
     * @return True if the hit was attached to track.
     */
      bool AddNearestHit(LitScalTrack* track, const PixelHitArray& hits, const std::vector<LitTrackParamScal*>& pars,
                         unsigned int nofHits);
 
      /*
     *
     */
      void MinMaxIndex(const LitTrackParamScal* par, const PixelHitArray& hits, fscal maxErr,
                       PixelHitConstIterator& first, PixelHitConstIterator& last);
 
     private:
      TrackArray fTracks;                   // Local copy of tracks
      LitDetectorLayoutMuon<fvec> fLayout;  // Detector geometry
      LitHitDataMuon<fvec> fHitData;        // Arranged hits
      unsigned char fMaxNofMissingHits;     // Maximum number of missing hits
      bool
        fIsProcessSubstationsTogether;  //  If true than hits from all substations are gathered together and hit selection is done over all this hits at a time
      fscal fSigmaCoef;                 // Sigma coefficient for fast hit search
      fscal fMaxCovSq;                  // Maximum covariance value
      fvec fChiSqPixelHitCut;           // Chi square cut for pixel hits
    };
 
  }  // namespace parallel
}  // namespace lit
#endif /* LITTRACKFINDERNNVECMUON_H_ */