CaTriplet.h

Go to the documentation of this file.

/* Copyright (C) 2019-2023 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Valentina Akishina, Sergey Gorbunov[committer] */
 
 
#pragma once  // include this header only once per compilation unit
 
#include "CaHit.h"
#include "CaSimd.h"
 
#include <string>
 
namespace cbm::algo::ca
{
 
  class Triplet {
   public:
    Triplet() = default;
 
    Triplet(ca::HitIndex_t iHitL, ca::HitIndex_t iHitM, ca::HitIndex_t iHitR, unsigned int iStaL, unsigned int iStaM,
            unsigned int iStaR, unsigned char Level, unsigned int firstNeighbour, char nNeighbours, fscal Chi2,
            fscal Qp, fscal Cqp, fscal tx, fscal Ctx, fscal ty, fscal Cty, bool isMomentumFitted)
      : fChi2(Chi2)
      , fQp(Qp)
      , fCqp(Cqp)
      , fTx(tx)
      , fCtx(Ctx)
      , fTy(ty)
      , fCty(Cty)
      , fFirstNeighbour(firstNeighbour)
      , fHitL(iHitL)
      , fHitM(iHitM)
      , fHitR(iHitR)
      , fNneighbours(nNeighbours)
      , fLevel(Level)
      , fSta((iStaL << 4) + ((iStaM - iStaL - 1) << 2) + (iStaR - iStaL - 2))
      , fIsMomentumFitted(isMomentumFitted)
    {
    }
 
 
    void SetLevel(unsigned char Level) { fLevel = Level; }
    unsigned char GetLevel() const { return fLevel; }
 
    ca::HitIndex_t GetLHit() const { return fHitL; }
    ca::HitIndex_t GetMHit() const { return fHitM; }
    ca::HitIndex_t GetRHit() const { return fHitR; }
 
    void SetNNeighbours(int n) { fNneighbours = n; }
    int GetNNeighbours() const { return fNneighbours; }
 
    void SetFNeighbour(unsigned int n) { fFirstNeighbour = n; }
    unsigned int GetFNeighbour() const { return fFirstNeighbour; }
 
    fscal GetQp() const { return fQp; }
    fscal GetChi2() const { return fChi2; }
    fscal GetTime() const { return -111.; }
 
    int GetLSta() const { return fSta >> 4; }
    int GetMSta() const { return ((fSta % 16) >> 2) + GetLSta() + 1; }
    int GetRSta() const { return (fSta % 4) + GetLSta() + 2; }
 
    fscal GetCqp() const { return fCqp; }
    fscal GetTx() const { return fTx; }
    fscal GetCtx() const { return fCtx; }
    fscal GetTy() const { return fTy; }
    fscal GetCty() const { return fCty; }
 
    bool IsMomentumFitted() const { return fIsMomentumFitted; }
    void SetIsMomentumFitted(bool val) { fIsMomentumFitted = val; }
 
    std::string ToString(int indentLevel = 0) const;
 
   private:
 
    fscal fChi2{0.};  
    fscal fQp{0.};    
    fscal fCqp{0.};   
    fscal fTx{0.};    
    fscal fCtx{0.};   
    fscal fTy{0.};    
    fscal fCty{0.};   
 
    unsigned int fFirstNeighbour{0};  
    ca::HitIndex_t fHitL{0};          
    ca::HitIndex_t fHitM{0};          
    ca::HitIndex_t fHitR{0};          
    int fNneighbours{0};              
 
    unsigned char fLevel{0};
 
    unsigned short fSta{0};     
    bool fIsMomentumFitted{0};  
  };
 
}  // namespace cbm::algo::ca