Trk::ForCFT Struct Reference

#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	nmcnst
double	wm [vkalNTrkM]
double	wmfit [vkalMaxNMassCnst]
double	localbmag
double	Ap
double	Bp
double	Cp
double	Dp
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 122 of file PrCFit.cxx.

                       {
  nmcnst=0;
  useMassCnst=0; usePhiCnst=0; useThetaCnst=0; usePointingCnst=0; usePlaneCnst=0;
  useAprioriVrt=0; usePassNear=0;
  Ap=Bp=Dp=Cp=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] = 139.57018;
     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 43 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 178 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 140 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 147 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 161 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 154 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 60 of file ForCFT.h.

{ useAprioriVrt = 1;}

vksetUseMassCnst()

void Trk::ForCFT::vksetUseMassCnst ( )

noexcept

Definition at line 168 of file PrCFit.cxx.

{ useMassCnst = 1;}

vksetUsePassNear()

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

inlinenoexcept

Definition at line 62 of file ForCFT.h.

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

vksetUsePhiCnst()

void Trk::ForCFT::vksetUsePhiCnst ( )

noexcept

Definition at line 169 of file PrCFit.cxx.

{ usePhiCnst = 1;}

vksetUsePlaneCnst()

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

noexcept

Definition at line 170 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 61 of file ForCFT.h.

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

vksetUseThetaCnst()

void Trk::ForCFT::vksetUseThetaCnst ( )

inlinenoexcept

Definition at line 59 of file ForCFT.h.

{ useThetaCnst = 1;}

Member Data Documentation

Ap

double Trk::ForCFT::Ap

Definition at line 30 of file ForCFT.h.

Bp

double Trk::ForCFT::Bp

Definition at line 30 of file ForCFT.h.

covvrt

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

Definition at line 39 of file ForCFT.h.

Cp

double Trk::ForCFT::Cp

Definition at line 30 of file ForCFT.h.

Dp

double Trk::ForCFT::Dp

Definition at line 30 of file ForCFT.h.

indtrkmc

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

Definition at line 46 of file ForCFT.h.

irob

int Trk::ForCFT::irob

Definition at line 43 of file ForCFT.h.

IterationNumber

int Trk::ForCFT::IterationNumber

Definition at line 47 of file ForCFT.h.

IterationPrecision

double Trk::ForCFT::IterationPrecision

Definition at line 48 of file ForCFT.h.

localbmag

double Trk::ForCFT::localbmag

Definition at line 28 of file ForCFT.h.

nmcnst

int Trk::ForCFT::nmcnst

Definition at line 25 of file ForCFT.h.

robres

double Trk::ForCFT::robres[vkalNTrkM]

Definition at line 45 of file ForCFT.h.

RobustScale

double Trk::ForCFT::RobustScale

Definition at line 44 of file ForCFT.h.

useAprioriVrt

int Trk::ForCFT::useAprioriVrt

Definition at line 19 of file ForCFT.h.

useMassCnst

int Trk::ForCFT::useMassCnst

Definition at line 15 of file ForCFT.h.

usePassNear

int Trk::ForCFT::usePassNear

Definition at line 20 of file ForCFT.h.

usePhiCnst

int Trk::ForCFT::usePhiCnst

Definition at line 16 of file ForCFT.h.

usePlaneCnst

int Trk::ForCFT::usePlaneCnst

Definition at line 21 of file ForCFT.h.

usePointingCnst

int Trk::ForCFT::usePointingCnst

Definition at line 18 of file ForCFT.h.

useThetaCnst

int Trk::ForCFT::useThetaCnst

Definition at line 17 of file ForCFT.h.

vrt

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

Definition at line 39 of file ForCFT.h.

vrtstp

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

Definition at line 42 of file ForCFT.h.

wgtvrt

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

Definition at line 39 of file ForCFT.h.

wm

double Trk::ForCFT::wm[vkalNTrkM]

Definition at line 26 of file ForCFT.h.

wmfit

double Trk::ForCFT::wmfit[vkalMaxNMassCnst]

Definition at line 27 of file ForCFT.h.

---

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

• ForCFT.h
• PrCFit.cxx