FairRoot/PandaRoot
Classes | Functions
PndSttHelixTrackFitter.h File Reference
#include "PndSttTrack.h"
#include "PndSttTrackFitter.h"
#include "TList.h"
#include <map>

Go to the source code of this file.

Classes

class  PndSttHelixTrackFitter
 

Functions

void fcnHelix (Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag)
 
void fcnHelix2 (Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag)
 

Function Documentation

void fcnHelix ( Int_t &  npar,
Double_t gin,
Double_t f,
Double_t par,
Int_t  iflag 
)

Definition at line 1593 of file PndSttHelixTrackFitter.cxx.

References Double_t, i, and sqrt().

Referenced by PndSttHelixTrackFitter::MinuitFit().

1594 {
1595 
1596  TMatrixT<Double_t> *mama = (TMatrixT<Double_t> *)gMinuit->GetObjectFit();
1597 
1598  Double_t chisq = 0;
1599  Double_t delta = 0;
1600  Int_t hitcounter = mama->GetNrows();
1601  for (Int_t i = 0; i < hitcounter; i++)
1602  {
1603  delta =sqrt((mama[0][i][0]-par[0])*(mama[0][i][0]-par[0])+(mama[0][i][1]-par[1])*(mama[0][i][1]-par[1])) -par[2] ;
1604  if(mama[0][i][2] == 0) chisq += (delta * delta * 12.);
1605  else chisq += (delta*delta)/(mama[0][i][2] * mama[0][i][2] / 12.);
1606  }
1607  f = chisq;
1608 }
Int_t i
Definition: run_full.C:25
friend F32vec4 sqrt(const F32vec4 &a)
Definition: P4_F32vec4.h:29
Double_t par[3]
Double_t
TFile * f
Definition: bump_analys.C:12
void fcnHelix2 ( Int_t &  npar,
Double_t gin,
Double_t f,
Double_t par,
Int_t  iflag 
)

Definition at line 1610 of file PndSttHelixTrackFitter.cxx.

References Double_t, i, and sqrt().

Referenced by PndSttHelixTrackFitter::MinuitFit().

1611 {
1612  TMatrixT<Double_t> *mama = (TMatrixT<Double_t> *)gMinuit->GetObjectFit();
1613 
1614  Double_t chisq = 0;
1615  Double_t delta = 0;
1616  Int_t hitcounter = mama->GetNrows();
1617  for (Int_t i = 0; i < hitcounter; i++)
1618  {
1619  delta =sqrt((mama[0][i][0]-par[0])*(mama[0][i][0]-par[0])+(mama[0][i][1]-par[1])*(mama[0][i][1]-par[1])) -par[2] ;
1620  if(mama[0][i][2] == 0) chisq += (delta * delta * 12.);
1621  else chisq += (delta*delta)/(pow(mama[0][i][3],2));
1622 
1623  }
1624 
1625  f = chisq;
1626 }
Int_t i
Definition: run_full.C:25
friend F32vec4 sqrt(const F32vec4 &a)
Definition: P4_F32vec4.h:29
Double_t par[3]
Double_t
TFile * f
Definition: bump_analys.C:12