CaCloneMerger.h

Go to the documentation of this file.

/* Copyright (C) 2022 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Sergei Zharko [committer], Maksym Zyzak */
 
 
#pragma once  // include this header only once per compilation unit
 
#include "CaHit.h"  // For ca::HitIndex_t
#include "CaInputData.h"
#include "CaParameters.h"
#include "CaSimd.h"  // TEMPORARY FOR fvec, fscal
#include "CaTrack.h"
#include "CaVector.h"
#include "CaWindowData.h"
 
namespace cbm::algo::ca
{
  class Track;
  class Framework;
 
  class CloneMerger {
   public:
    CloneMerger(const ca::Parameters<fvec>& pars, const fscal mass);
 
    ~CloneMerger();
 
    CloneMerger(const CloneMerger&) = default;
 
    CloneMerger(CloneMerger&&) = default;
 
    CloneMerger& operator=(const CloneMerger&) = delete;
 
    CloneMerger& operator=(CloneMerger&&) = delete;
 
 
    //void Exec(Vector<Track>& extTracks, Vector<ca::HitIndex_t>&, const ca::InputData& input);
    void Exec(const ca::InputData& input, WindowData& wData);
 
   private:
    // ***************
    // ** Functions **
    // ***************
 
    static void InvertCholesky(fvec a[15]);
 
    static void MultiplySS(fvec const C[15], fvec const V[15], fvec K[5][5]);
 
    static void MultiplyMS(fvec const C[5][5], fvec const V[15], fvec K[15]);
 
    static void MultiplySR(fvec const C[15], fvec const r_in[5], fvec r_out[5]);
 
    static void FilterTracks(fvec const r[5], fvec const C[15], fvec const m[5], fvec const V[15], fvec R[5],
                             fvec W[15], fvec* chi2);
 
 
    // ***************
    // ** Variables **
    // ***************
 
    Vector<unsigned short> fTrackFirstStation{"CloneMerger::fTrackFirstStation"};
 
    Vector<unsigned short> fTrackLastStation{"CloneMerger::fTrackLastStation"};
 
    Vector<ca::HitIndex_t> fTrackFirstHit{"CloneMerger::fTrackFirstHit"};
 
    Vector<ca::HitIndex_t> fTrackLastHit{"CloneMerger::fTrackLastHit"};
 
    Vector<unsigned short> fTrackNeighbour{"CloneMerger::fTrackNeighbour"};
 
    Vector<fscal> fTrackChi2{"CloneMerger::fTrackChi2"};
 
    Vector<char> fTrackIsStored{"CloneMerger::fTrackIsStored"};
 
    Vector<char> fTrackIsDownstreamNeighbour{"CloneMerger::fTrackIsDownstreamNeighbour"};
 
    Vector<Track> fTracksNew{"CaCloneMerger::fTracksNew"};  
 
    Vector<ca::HitIndex_t> fRecoHitsNew{"CaCloneMerger::fRecoHitsNew"};  
 
    const Parameters<fvec>& fParameters;            
    fscal fDefaultMass{constants::phys::MuonMass};  
  };
 
}  // namespace cbm::algo::ca