da/d38/ForCFT_8h_source.html

/*

  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration

*/


#ifndef _TrkVKalVrtCore_ForCFT_H

#define _TrkVKalVrtCore_ForCFT_H


#include "TrkVKalVrtCore/CommonPars.h"


namespace Trk {


  struct ForCFT

  {

//

    int useMassCnst{};

    int usePhiCnst{};

    int useThetaCnst{};

    int usePointingCnst{};

    int useAprioriVrt{};

    int usePassNear{};

    int usePlaneCnst{};

    int useRadiusCnst{};


//

// For several(up to 8) mass constraints

    int nmcnst;

    double wm[vkalNTrkM];

    double wmfit[vkalMaxNMassCnst];

    double localbmag;

// Since 20/09/2009 - Vertex in plane constraint

    double Ap{};

    double Bp{};

    double Cp{};

    double Dp{};

    double RC{};

    double radiusRefP[2]{};


//

// Vertex, covariance and it's weight (inverse covariance)

// Used for pointing and a priori vertex constraints

// Since 15/03/2009:

//    For "passing near" vrt+covvrt+charge from ForVrtClose structure is used.

//    For pointing cnst vertex from constraint object itself is used

//---------------------------------------------------------------------------

    double vrt[3]{}, covvrt[6]{}, wgtvrt[6]{};

//

// temporary vertex in global ref.frame

    double vrtstp[3]{};

    int irob;

    double RobustScale;

    double robres[vkalNTrkM];

    int indtrkmc[vkalMaxNMassCnst][vkalNTrkM];

    int IterationNumber;

    double IterationPrecision;


    void prcfit( long int ntrk, double  *wm, const double  *wmfit, double  bmag, const double  *vrt, const double  *vrte) noexcept;


    void vksetIterationNum(long int Iter) noexcept;

    void vksetIterationPrec(double Prec) noexcept;

    void vksetRobustScale(double Scale) noexcept;

    void vksetRobustness(long int Rob) noexcept;

    void vksetUseMassCnst() noexcept;

    void vksetUsePhiCnst() noexcept;

    void vksetUsePlaneCnst(double a, double b, double c, double d) noexcept;

    void vksetUseThetaCnst() noexcept { useThetaCnst = 1;}

    void vksetUseAprioriVrt() noexcept { useAprioriVrt = 1;}

    void vksetUsePointingCnst(int iType = 1 ) noexcept { usePointingCnst = iType<2 ? 1 : 2 ;}

    void vksetUsePassNear(int iType = 1 )  noexcept { usePassNear = iType<2 ? 1 : 2 ;}

    void setmasscnst_(long int ncnsttrk, long int *indextrk, double  wmcnst) noexcept;


    ForCFT() noexcept;//Body is in PrCFit.cxx

    ~ForCFT() = default;


  };


}

#endif