#include <ForCFT.h>

Collaboration diagram for Trk::ForCFT:

Public Member Functions
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
void	vksetUseAprioriVrt () noexcept
void	vksetUsePointingCnst (int iType=1) noexcept
void	vksetUsePassNear (int iType=1) noexcept
void	setmasscnst_ (long int ncnsttrk, long int *indextrk, double wmcnst) noexcept
	ForCFT () noexcept
	~ForCFT ()=default

Public Attributes
int	useMassCnst {}
int	usePhiCnst {}
int	useThetaCnst {}
int	usePointingCnst {}
int	useAprioriVrt {}
int	usePassNear {}
int	usePlaneCnst {}
int	useRadiusCnst {}
int	nmcnst
double	wm [vkalNTrkM]
double	wmfit [vkalMaxNMassCnst]
double	localbmag
double	Ap {}
double	Bp {}
double	Cp {}
double	Dp {}
double	RC {}
double	radiusRefP [2] {}
double	vrt [3] {}
double	covvrt [6] {}
double	wgtvrt [6] {}
double	vrtstp [3] {}
int	irob
double	RobustScale
double	robres [vkalNTrkM]
int	indtrkmc [vkalMaxNMassCnst][vkalNTrkM]
int	IterationNumber
double	IterationPrecision

Detailed Description

Definition at line 12 of file ForCFT.h.

Constructor & Destructor Documentation

◆ ForCFT()

Trk::ForCFT::ForCFT ( )

noexcept

Definition at line 123 of file PrCFit.cxx.

                       {
  nmcnst=0;
  useMassCnst=0; usePhiCnst=0; useThetaCnst=0; usePointingCnst=0; usePlaneCnst=0;
  useAprioriVrt=0; usePassNear=0;
  Ap=Bp=Dp=Cp=RC=0.;
  IterationNumber = 50;
  IterationPrecision=1.e-3;
  RobustScale = 1.; irob=0;
  for (int ic=0; ic<vkalMaxNMassCnst; ++ic) wmfit[ic] = -10000.;
  for (int it=0; it<vkalNTrkM; ++it) {
     wm[it] = ParticleConstants::chargedPionMassInMeV;
     robres[it] = 1.;
     for(int ic=0; ic<vkalMaxNMassCnst; ic++) indtrkmc[ic][it]=0;
  }
  localbmag=1.997;   // Safety: standard magnetic field in ID
}

◆ ~ForCFT()

Trk::ForCFT::~ForCFT ( )

default

Member Function Documentation

◆ prcfit()

void Trk::ForCFT::prcfit	(	long int	ntrk,
		double *	wm,
		const double *	wmfit,
		double	bmag,
		const double *	vrt,
		const double *	vrte )

noexcept

Definition at line 44 of file PrCFit.cxx.

{
    long int i__1;
    double   summ;
 
 
/*------------------------------------------------------------------------------*/
/*        SETTING OF INITIAL PARAMETERS FOR C-FIT               */
/*        (DEPENDING ON TYPE OF FIT)                        */
/*          INPUT:                              */
/*         NTRK        : number of tracks               */
/*                 WM(JTK)     : track masses ( if negative - track doesn't */
/*                 : participate in mass constraint )       */
/*                 WMFIT       : masses for mass constraints            */
/*                 VRT(3)      : point for directional constraint       */
/*                                        (usually primary vertex)      */
/*                 VRTE(6)     : error matrix for VRT ( for IFLAG=6)        */
/*                 BMAG        : magnetic field (in Tesla)          */
/*    Magnetic field is obligatory , other parameters MUST be set       */
/*    when used in constraints                          */
/*     General mass constraint for all tracks is set in PRCFIT through WMFIT    */
/*    if WMFIT > Sum_of_track_masses                        */
/*       Any number of additional mass constraints can be set up with       */
/*     SetMassCnst(NCnstTrk,IndexTrk,WMCNST)     where              */
/*               INPUT:                             */
/*                    NCnstTrk           :track number of track subset for cnst.*/
/*                     IndexTrk(NCnstTrk) :index of each track from subset in   */
/*                                         :general track set for vertex fit.   */
/*                       WmCnst             :Mass for given constraint      */
/* Author: V.Kostioukhine (1995 ,2004,2005)                     */
/*                                          */
/*      All functions here must be called AFTER(!!!) prcfit.        */
/*                                      */
/*------------------------------------------------------------------------------*/
 
    this->localbmag = bmag;
    this->nmcnst = 0;
    for (int i=0; i<8; ++i) this->wmfit[i] = -10000.;
    summ = 0.;
    i__1 = ntrk<vkalNTrkM ? ntrk: vkalNTrkM;
    for (int i=0; i<i__1; ++i) {
      this->wm[i] =  std::abs(wm[i]);
      summ += wm[i];
    }
    if ((*wmfit) > summ) {
/*  Set general mass constraint based on ALL tracks */
    this->nmcnst = 1;
    for (int i = 0; i < vkalNTrkM; ++i) {
        indtrkmc[0][i] = 0;
        if (i < ntrk) {indtrkmc[0][i] = 1;}
    }
    this->wmfit[0] = (*wmfit);
    }
 
    this->vrt[0] = vrt[0];
    this->vrt[1] = vrt[1];
    this->vrt[2] = vrt[2];
    this->covvrt[0] = vrte[0];
    this->covvrt[1] = vrte[1];
    this->covvrt[2] = vrte[2];
    this->covvrt[3] = vrte[3];
    this->covvrt[4] = vrte[4];
    this->covvrt[5] = vrte[5];
    this->irob = 0;
    this->IterationNumber    = 50;
    this->IterationPrecision = 1.e-3;
    for (int i = 0; i < vkalNTrkM; ++i) this->robres[i]=1.; //Safety
//
// Reset all constraints
//
this->useMassCnst = 0;
this->usePhiCnst = 0;
this->useThetaCnst = 0;
this->useAprioriVrt = 0;
this->usePointingCnst = 0;
this->usePassNear = 0;
//forcft_1.usePlaneCnst = 0;   //Used only on demand=> must NOT be reset here!!!
}

◆ setmasscnst_()

void Trk::ForCFT::setmasscnst_	(	long int	ncnsttrk,
		long int *	indextrk,
		double	wmcnst )

noexcept

Definition at line 179 of file PrCFit.cxx.

{
    if (indextrk==nullptr) return;  //Protection!  Track indices start from 1 (not 0)!
    --indextrk;
 
    ++nmcnst;
    if (nmcnst > 8) return ;
    for (int i = 0; i < vkalNTrkM; ++i) {
    indtrkmc[nmcnst-1][i] = 0;
    }
    for (int i = 0; i < ncnsttrk;  ++i) {
    if (indextrk[i] > 0 && indextrk[i] < vkalNTrkM) {
        indtrkmc[nmcnst-1][indextrk[i]] = 1;
    }
    }
    wmfit[nmcnst - 1] = wmcnst;
 }

◆ vksetIterationNum()

void Trk::ForCFT::vksetIterationNum ( long int Iter )

noexcept

Definition at line 141 of file PrCFit.cxx.

{
    if (Iter<3)   Iter=3;
    if (Iter>100) Iter=100;
    IterationNumber    = Iter;
}

◆ vksetIterationPrec()

void Trk::ForCFT::vksetIterationPrec ( double Prec )

noexcept

Definition at line 148 of file PrCFit.cxx.

{
    if (Prec<1.e-5)   Prec=1.e-5;
    if (Prec>1.e-1)   Prec=1.e-1;
    IterationPrecision = Prec;
}

◆ vksetRobustness()

void Trk::ForCFT::vksetRobustness ( long int Rob )

noexcept

Definition at line 162 of file PrCFit.cxx.

{
    if (Rob<0)   Rob=0;
    if (Rob>7)   Rob=7;
    irob = Rob;
}

◆ vksetRobustScale()

void Trk::ForCFT::vksetRobustScale ( double Scale )

noexcept

Definition at line 155 of file PrCFit.cxx.

{
    if (Scale<0.01)   Scale=0.01;
    if (Scale>100.)   Scale=100.;
    RobustScale = Scale;
}

◆ vksetUseAprioriVrt()

void Trk::ForCFT::vksetUseAprioriVrt ( )

inlinenoexcept

Definition at line 66 of file ForCFT.h.

66{ useAprioriVrt = 1;}

◆ vksetUseMassCnst()

void Trk::ForCFT::vksetUseMassCnst ( )

noexcept

Definition at line 169 of file PrCFit.cxx.

169{ useMassCnst = 1;}

◆ vksetUsePassNear()

void Trk::ForCFT::vksetUsePassNear ( int iType = 1 )

inlinenoexcept

Definition at line 68 of file ForCFT.h.

68{ usePassNear = iType<2 ? 1 : 2 ;}

◆ vksetUsePhiCnst()

void Trk::ForCFT::vksetUsePhiCnst ( )

noexcept

Definition at line 170 of file PrCFit.cxx.

170{ usePhiCnst = 1;}

◆ vksetUsePlaneCnst()

void Trk::ForCFT::vksetUsePlaneCnst	(	double	a,
		double	b,
		double	c,
		double	d )

noexcept

Definition at line 171 of file PrCFit.cxx.

                                                                               {
  if(a+b+c+d == 0.){  usePlaneCnst = 0;
  }else{              usePlaneCnst = 1; }
  Ap = a; Bp = b; Cp = c; Dp = d;
}

◆ vksetUsePointingCnst()

void Trk::ForCFT::vksetUsePointingCnst ( int iType = 1 )

inlinenoexcept

Definition at line 67 of file ForCFT.h.

67{ usePointingCnst = iType<2 ? 1 : 2 ;}

◆ vksetUseThetaCnst()

void Trk::ForCFT::vksetUseThetaCnst ( )

inlinenoexcept

Definition at line 65 of file ForCFT.h.

65{ useThetaCnst = 1;}

Member Data Documentation

◆ Ap

double Trk::ForCFT::Ap {}

Definition at line 31 of file ForCFT.h.

31{};

◆ Bp

double Trk::ForCFT::Bp {}

Definition at line 32 of file ForCFT.h.

32{};

◆ covvrt

double Trk::ForCFT::covvrt[6] {}

Definition at line 45 of file ForCFT.h.

45{}, covvrt[6]{}, wgtvrt[6]{};

Trk::ForCFT::wgtvrt

double wgtvrt[6]

Definition ForCFT.h:45

◆ Cp

double Trk::ForCFT::Cp {}

Definition at line 33 of file ForCFT.h.

33{};

◆ Dp

double Trk::ForCFT::Dp {}

Definition at line 34 of file ForCFT.h.

34{};

◆ indtrkmc

int Trk::ForCFT::indtrkmc[vkalMaxNMassCnst][vkalNTrkM]

Definition at line 52 of file ForCFT.h.

◆ irob

int Trk::ForCFT::irob

Definition at line 49 of file ForCFT.h.

◆ IterationNumber

int Trk::ForCFT::IterationNumber

Definition at line 53 of file ForCFT.h.

◆ IterationPrecision

double Trk::ForCFT::IterationPrecision

Definition at line 54 of file ForCFT.h.

◆ localbmag

double Trk::ForCFT::localbmag

Definition at line 29 of file ForCFT.h.

◆ nmcnst

int Trk::ForCFT::nmcnst

Definition at line 26 of file ForCFT.h.

◆ radiusRefP

double Trk::ForCFT::radiusRefP[2] {}

Definition at line 36 of file ForCFT.h.

36{};

◆ RC

double Trk::ForCFT::RC {}

Definition at line 35 of file ForCFT.h.

35{};

◆ robres

double Trk::ForCFT::robres[vkalNTrkM]

Definition at line 51 of file ForCFT.h.

◆ RobustScale

double Trk::ForCFT::RobustScale

Definition at line 50 of file ForCFT.h.

◆ useAprioriVrt

int Trk::ForCFT::useAprioriVrt {}

Definition at line 19 of file ForCFT.h.

19{};

◆ useMassCnst

int Trk::ForCFT::useMassCnst {}

Definition at line 15 of file ForCFT.h.

15{};

◆ usePassNear

int Trk::ForCFT::usePassNear {}

Definition at line 20 of file ForCFT.h.

20{};

◆ usePhiCnst

int Trk::ForCFT::usePhiCnst {}

Definition at line 16 of file ForCFT.h.

16{};

◆ usePlaneCnst

int Trk::ForCFT::usePlaneCnst {}

Definition at line 21 of file ForCFT.h.

21{};

◆ usePointingCnst

int Trk::ForCFT::usePointingCnst {}

Definition at line 18 of file ForCFT.h.

18{};

◆ useRadiusCnst

int Trk::ForCFT::useRadiusCnst {}

Definition at line 22 of file ForCFT.h.

22{};

◆ useThetaCnst

int Trk::ForCFT::useThetaCnst {}

Definition at line 17 of file ForCFT.h.

17{};

◆ vrt

double Trk::ForCFT::vrt[3] {}

Definition at line 45 of file ForCFT.h.

45{}, covvrt[6]{}, wgtvrt[6]{};

◆ vrtstp

double Trk::ForCFT::vrtstp[3] {}

Definition at line 48 of file ForCFT.h.

48{};

◆ wgtvrt

double Trk::ForCFT::wgtvrt[6] {}

Definition at line 45 of file ForCFT.h.

45{}, covvrt[6]{}, wgtvrt[6]{};

◆ wm

double Trk::ForCFT::wm[vkalNTrkM]

Definition at line 27 of file ForCFT.h.

◆ wmfit

double Trk::ForCFT::wmfit[vkalMaxNMassCnst]

Definition at line 28 of file ForCFT.h.

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

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ForCFT()

◆ ~ForCFT()

Member Function Documentation

◆ prcfit()

◆ setmasscnst_()

◆ vksetIterationNum()

◆ vksetIterationPrec()

◆ vksetRobustness()

◆ vksetRobustScale()

◆ vksetUseAprioriVrt()

◆ vksetUseMassCnst()

◆ vksetUsePassNear()

◆ vksetUsePhiCnst()

◆ vksetUsePlaneCnst()

◆ vksetUsePointingCnst()

◆ vksetUseThetaCnst()

Member Data Documentation

◆ Ap

◆ Bp

◆ covvrt

◆ Cp

◆ Dp

◆ indtrkmc

◆ irob

◆ IterationNumber

◆ IterationPrecision

◆ localbmag

◆ nmcnst

◆ radiusRefP

◆ RC

◆ robres

◆ RobustScale

◆ useAprioriVrt

◆ useMassCnst

◆ usePassNear

◆ usePhiCnst

◆ usePlaneCnst

◆ usePointingCnst

◆ useRadiusCnst

◆ useThetaCnst

◆ vrt

◆ vrtstp

◆ wgtvrt

◆ wm

◆ wmfit