RhoKalmanVtxFitter Class Reference

#include <RhoKalmanVtxFitter.h>

Inheritance diagram for RhoKalmanVtxFitter:

Public Member Functions
	RhoKalmanVtxFitter (RhoCandidate *b)
virtual	~RhoKalmanVtxFitter ()
void	SetExpansionPoint (TVector3 P)
void	SetNIterations (int i)
double	FitVertexFast (TVector3 &vtx, TMatrixD &cov, bool skipcov=false, int niter=2)
void	SetDebug (bool db=true)
void	SetSilent ()
Bool_t	Fit ()
Bool_t	FitAll ()
double	Chi2Contribution (const RhoCandidate *)
double	GetChi2 () const
int	GetNdf () const
double	GetProb () const
void	SetVerbose (Bool_t v=kTRUE)

Protected Member Functions
RhoCandidate &	HeadOfTree () const
RhoCandidate *	CopyCand (RhoCandidate *)
	uppermost particle composite in tree More...
RhoCandidate *	CopyTree (RhoCandidate *)
void	InsertChi2 (const RhoCandidate *bc, const double chi2)
void	SetDaugthersFromComposite (RhoCandidate *cand)
void	FindAndAddFinalStateDaughters (RhoCandidate *cand)
void	SetFourMomentumByDaughters (RhoCandidate *composite)
void	SetDecayVertex (RhoCandidate *composite, const TVector3 &vtx, const TMatrixD &CovVV)

Protected Attributes
Bool_t	fVerbose
RhoCandidate *	fHeadOfTree
std::vector< RhoCandidate * >	fDaughters
double	fChiSquare
int	fNDegreesOfFreedom

Private Member Functions
Bool_t	FitNode (RhoCandidate *b)
Bool_t	Calculate ()
Bool_t	CalcPrgParams (RhoCandidate *cand, TVector3 expansionpoint)
Bool_t	CalcPrgParams (RhoCandidate *cand)

Private Attributes
RhoCandidate *	fCurrentHead
bool	fDebug
int	fNIterations
TVector3	fExpansionPoint
Float_t	fPrgParams [5]
TMatrixD	fPrgCov
TMatrixD	fJacobian

Detailed Description

Definition at line 21 of file RhoKalmanVtxFitter.h.

Constructor & Destructor Documentation

RhoKalmanVtxFitter::RhoKalmanVtxFitter

(

RhoCandidate *

b

)

Definition at line 16 of file RhoKalmanVtxFitter.cxx.

References fCurrentHead, RhoFitterBase::fHeadOfTree, and RhoFitterBase::fVerbose.

                                                       :
  RhoFitterBase(b),
  fDebug(false),
  fNIterations(2),
  fExpansionPoint(0.,0.,0.),  // Assume Zero origin as default
  fPrgCov(5,5),
  fJacobian(5,7)
{
  fVerbose=false;
  fCurrentHead=fHeadOfTree;
}

RhoKalmanVtxFitter::~RhoKalmanVtxFitter ( )

virtual

Definition at line 28 of file RhoKalmanVtxFitter.cxx.

{
}

Member Function Documentation

Bool_t RhoKalmanVtxFitter::CalcPrgParams ( RhoCandidate *  cand, TVector3  expansionpoint  )

private

Definition at line 33 of file RhoKalmanVtxFitter.cxx.

References Bool_t, RhoCandidate::Cov7(), Double_t, fJacobian, fPrgCov, fPrgParams, RhoCandidate::GetCharge(), RhoCandidate::GetPosition(), mom, RhoCandidate::P4(), RhoCalculationTools::P7toPRG(), and pos.

Referenced by Calculate(), and FitVertexFast().

{
  // calculate helix and its covariance in the preigee representation

  TLorentzVector mom = tcand->P4();
  TVector3 pos= tcand->GetPosition();
  Double_t charge = tcand->GetCharge();
  TMatrixD cov77 = tcand->Cov7();
  Bool_t test = RhoCalculationTools::P7toPRG(pos, mom, charge, cov77, expansionpoint, fPrgParams, fPrgCov, fJacobian, false);
  return test;
}

Bool_t RhoKalmanVtxFitter::CalcPrgParams ( RhoCandidate *  cand )

inlineprivate

Definition at line 40 of file RhoKalmanVtxFitter.h.

{return CalcPrgParams(cand,fExpansionPoint);}; // calculate helix and its covariance

Bool_t RhoKalmanVtxFitter::Calculate ( )

private

Definition at line 192 of file RhoKalmanVtxFitter.cxx.

References Bool_t, CalcPrgParams(), RhoCandidate::Daughter(), Double_t, RhoCandidate::E(), RhoFitterBase::fChiSquare, fCurrentHead, fDebug, fExpansionPoint, fJacobian, RhoFitterBase::fNDegreesOfFreedom, fNIterations, fPrgCov, fPrgParams, RhoFitterBase::fVerbose, RhoCandidate::GetCharge(), RhoCandidate::GetMomentum(), RhoCandidate::GetPosition(), i, RhoFitterBase::InsertChi2(), RhoCandidate::NDaughters(), Print(), printf(), RhoCandidate::SetCov7(), RhoFitterBase::SetDecayVertex(), RhoFitterBase::SetFourMomentumByDaughters(), RhoCandidate::SetP3(), RhoCandidate::SetPos(), T, and W.

Referenced by FitNode().

{
  // We'll do a vertex finding and Fitting in empty space at a constant magnetic field along z-direction
  // We don't dive into the particle tree, it's just the daughters of this candidate

  int nTrk = fCurrentHead->NDaughters();
  fNDegreesOfFreedom=2*nTrk-3;

  TVector3 vtx(0.,0.,0.);
  std::vector<TVector3> momenta;
  TMatrixD tmpcov(3,3);
  Double_t chisquare=500000;//some high value
  Double_t determinant = 0.;
  TMatrixD qip(5,nTrk); // starting parameters of all tracks
  std::vector<TMatrixD> B;
  std::vector<TMatrixD> GI; // only Gi^-1 is used
  std::vector<TMatrixD> BGI; // helper to save on multiplications (more memory needed, though)
  std::vector<TMatrixD> D; // helper
  std::vector<TMatrixD> E; //helper
  TMatrixD A(3,3);
  TMatrixD T(3,1);
  std::vector<TMatrixD> U;
  std::vector<TMatrixD> W;
  std::vector<TMatrixD> dq; // helper
  TMatrixD CovFitFull(3+3*nTrk,3+3*nTrk); //3 vtx coordinates, nTrk*3 momenta
  TMatrixD CovVV(3,3);

  // Find a nice seed value
  bool verbi=fVerbose;
  fVerbose=false;
  bool debi=fDebug;
  fDebug=false; // temporary silence
  //FitVertexFast(vtx,tmpcov,false);//fast vertex seed (with expansion point used)
  fDebug=debi;
  fVerbose=verbi;
  if(fVerbose) {std::cout<<" #$# Fit #$# Vertex after fast prefit: "; vtx.Print();}
  //  fExpansionPoint = vtx; // move to vincinity
  for (int iteration = 0 ; iteration < fNIterations ; iteration++) {
    if(fVerbose) {
      std::cout<<"Iteration "<<iteration<<std::endl;
      std::cout<<" #$# Fit #$# Begin iteration "<<iteration<<" #$#$#$# "<<std::endl;
    }
    //cleanup
    B.clear();
    GI.clear(); // only Gi^-1 is used
    BGI.clear(); // helper to save on multiplications (more memory needed, though)
    D.clear(); // helper
    E.clear(); //helper
    U.clear();
    W.clear();
    dq.clear(); // helper
    A.Zero();
    T.Zero();
    CovFitFull.Zero();
    CovVV.Zero();

    for(int i=0; i<nTrk; i++) {
      if(fVerbose) {std::cout<<" #$# Fit #$# track "<<i<<" #$#$#$# "<<std::endl;}
      if(fVerbose) {std::cout<<" #$# Fit #$# Expansion point        "; fExpansionPoint.Print();}
      if(fVerbose) {std::cout<<" #$# Fit #$# Vertex before (global) "; (vtx).Print();}
      RhoCandidate* tcand = fCurrentHead->Daughter(i);

      const Double_t charge = tcand->GetCharge();

      if(iteration==0) {
        //On the first run fetch the starting values (being updated)
        //Bool_t testprg = CalcPrgParams(tcand,fExpansionPoint);
        Bool_t testprg = CalcPrgParams(tcand,vtx);
        if(!testprg) {printf("#$# Fit #$# CANNOT CALCULATE TRACK PARAMETERS (momenta at origin)"); return kFALSE;}
        qip[0][i]=0.;//fPrgParams[0];// the initial guess is AT the guess vertex. No epsilon_0
        qip[1][i]=0.;//fPrgParams[1];// the initial guess is AT the guess vertex. No z_0
        qip[2][i]=fPrgParams[2];//
        qip[3][i]=fPrgParams[3];//phi0
        qip[4][i]=fPrgParams[4];

        if(fDebug) {printf(" ##$# Fit #$# Measurement p: "); tcand->GetMomentum().Print();}
        if(fDebug) {printf(" ##$# Fit #$# Measurement x: "); tcand->GetPosition().Print();}
        momenta.push_back(tcand->GetMomentum());
        momenta[i].SetPhi(fPrgParams[3]); // set phi while keeping magnitude and theta constant
        if(fDebug) {printf(" ##$# Fit #$# Momentum at Expansionpoint (phi=%g): ",fPrgParams[3]); momenta[i].Print();}
      }

      // get parameters closest to current vertex estimate
      Bool_t testvtx = CalcPrgParams(tcand,vtx);
      if(!testvtx) { return kFALSE; }
      if (fVerbose) {
        printf("#$# Fit #$# Helix params:\n\t epsilon = %.4g cm\n\t Z0\t = %.4g cm\n\t theta\t = %.4g\n\t phi0\t = %.4g\n\t rho\t = %.4g 1/cm\n\t charge\t = %g e\n",fPrgParams[0],fPrgParams[1],fPrgParams[2],fPrgParams[3],fPrgParams[4],charge);
        std::cout<<" #$# Fit #$# Helix cov: ";
        fPrgCov.Print();
      }

      TMatrixD qiv(5,1); // parameters describing relative to expansion point the vertex seed
      qiv[0][0]=fPrgParams[0];//epsilon AT vertex zero
      qiv[1][0]=fPrgParams[1];//z0 AT vertex zero
      qiv[2][0]=fPrgParams[2];//
      qiv[3][0]=fPrgParams[3];//phi0
      qiv[4][0]=fPrgParams[4];//
      //dqi = q_measured - F_0 == div-qip;
      TMatrixD qipi(5,1);
      for (int kk=0; kk<5; kk++) { qipi[kk][0]=qip[kk][i]; }
      TMatrixD dqi(qiv,TMatrixD::kMinus,qipi);
      dq.push_back(dqi);
      if(fDebug) {std::cout<<" #$# Fit #$# dqi  "; dqi.Print();}

      TMatrixD COVi(5,5); // track parameter cov for (epsilon,z_p,Phi_p)
      //Double_t r=pocai.Perp();

      COVi=fPrgCov;
      if(fVerbose) {std::cout<<" #$# Fit #$# COVi  "; COVi.Print();}

      //Take the derivatives at the current vertex estimate (seed) (cf. Avery1 p.4)
      TMatrixD Di(5,3); // Derivative in V
      Di[0][0]=fJacobian[0][0]; // dEpsilon  / dvx
      Di[0][1]=fJacobian[0][1]; // dEpsilon   /dvy
      Di[0][2]=fJacobian[0][2]; // dEpsilon   /dvz
      Di[1][0]=fJacobian[1][0]; // dZ0   /dvx
      Di[1][1]=fJacobian[1][1]; // dZ0   /dvy
      Di[1][2]=fJacobian[1][2]; // dZ0   /dvz
      Di[2][0]=fJacobian[2][0]; //
      Di[2][1]=fJacobian[2][1]; //
      Di[2][2]=fJacobian[2][2]; //
      Di[3][0]=fJacobian[3][0]; // dPhi0 /dvx
      Di[3][1]=fJacobian[3][1]; // dPhi0 /dvy
      Di[3][2]=fJacobian[3][2]; // dPhi0 /dvz
      Di[4][0]=fJacobian[4][0]; //
      Di[4][1]=fJacobian[4][1]; //
      Di[4][2]=fJacobian[4][2]; //
      //Di.T();
      D.push_back(Di);
      if(fDebug) {std::cout<<" #$# Fit #$# Di  "; Di.Print();}

      TMatrixD Ei(5,3); // Derivative in p
      Ei[0][0]=fJacobian[0][3]; // dEpsilon   /dpx
      Ei[0][1]=fJacobian[0][4]; // dEpsilon   /dpy
      Ei[0][2]=fJacobian[0][5]; // dEpsilon   /dpz
      Ei[1][0]=fJacobian[1][3]; // dZ0   /dpx
      Ei[1][1]=fJacobian[1][4]; // dZ0   /dpy
      Ei[1][2]=fJacobian[1][5]; // dZ0   /dpz
      Ei[2][0]=fJacobian[2][3]; //
      Ei[2][1]=fJacobian[2][4]; //
      Ei[2][2]=fJacobian[2][5]; //
      Ei[3][0]=fJacobian[3][3]; // dPhi0 /dpx
      Ei[3][1]=fJacobian[3][4]; // dPhi0 /dpy
      Ei[3][2]=fJacobian[3][5]; // dPhi0 /dpz
      Ei[4][0]=fJacobian[4][3]; //
      Ei[4][1]=fJacobian[4][4]; //
      Ei[4][2]=fJacobian[4][5]; //
      //Ei.T();
      E.push_back(Ei);
      if(fDebug) {std::cout<<" #$# Fit #$# Ei  "; Ei.Print();}

      //det(aA) = a^n det(A)
      TMatrixD Wi(COVi);
      determinant=0.;
      Wi.InvertFast(&determinant);
      if (determinant==0) {
        std::cout<<"RhoKalmanVtxFitter: COVi Inversion failed, abort fit."<<std::endl;
        return kFALSE;
      }
      W.push_back(Wi);
      if(fDebug) {
        std::cout<<" #$# Fit #$# Wi (det(cov) = "<<determinant<<")"<<std::endl;
        Wi.Print();
        for(int sdsd=0; sdsd<3; sdsd++) for(int asas=sdsd; asas<3; asas++) {
            printf("W{%i}[%i][%i] = %6.9f\n",i,sdsd,asas,Wi[sdsd][asas]);
            printf("W{%i}[%i][%i] = %6.9f\n",i,asas,sdsd,Wi[asas][sdsd]);
          }
      }

      // calculate matrices to be stored: A, Bi, Gi, T, Ui
      TMatrixD DitWi(Di,TMatrixD::kTransposeMult,Wi);
      if(fDebug) {std::cout<<" #$# Fit #$# DitWi  "; DitWi.Print();}

      TMatrixD EitWi(Ei,TMatrixD::kTransposeMult,Wi);
      if(fDebug) {std::cout<<" #$# Fit #$# EitWi  "; EitWi.Print();}

      TMatrixD Ai(DitWi,TMatrixD::kMult,Di);
      A += Ai;
      if(fDebug) {std::cout<<" #$# Fit #$# Ai ("<<i<<") = DitWiDi "; Ai.Print();}

      TMatrixD Bi(DitWi,TMatrixD::kMult,Ei);
      B.push_back(Bi);
      if(fDebug) {std::cout<<" #$# Fit #$# Bi = DitWiEi "; Bi.Print();}

      TMatrixD Gi(EitWi,TMatrixD::kMult,Ei);
      if(fDebug) {std::cout<<" #$# Fit #$# Gi = EitWiEi "; Gi.Print();}

      TMatrixD GIi(Gi);
      determinant=0.;
      GIi.InvertFast(&determinant);
      if (determinant==0) {
        std::cout<<"RhoKalmanVtxFitter: GIi Inversion failed, abort fit."<<std::endl;
        return kFALSE;
      }
      GI.push_back(GIi);
      if(fDebug) {std::cout<<" #$# Fit #$# GIi = Gi^-1 (det="<<determinant<<")"; GIi.Print();}

      TMatrixD BiGIi(Bi,TMatrixD::kMult,GIi);
      BGI.push_back(BiGIi);
      if(fDebug) {std::cout<<" #$# Fit #$# BiGIi "; BiGIi.Print();}

      TMatrixD Ti(DitWi,TMatrixD::kMult,dqi);
      T += Ti;
      if(fDebug) {std::cout<<" #$# Fit #$# Ti = DitWi*dqi "; Ti.Print();}

      TMatrixD Ui(EitWi,TMatrixD::kMult,dqi);
      U.push_back(Ui);
      if(fDebug) {std::cout<<" #$# Fit #$# Ui = EitWi*dqi"; Ui.Print();}

    } // loop tracks
    if(fDebug) {std::cout<<" #$# Fit #$# A = "; A.Print();}

    // Calculate Vertex update
    TMatrixD WV(A);
    if(fDebug) {std::cout<<" #$# Fit #$# WV = "; WV.Print();}
    TMatrixD Vpre(T);
    if(fDebug) {std::cout<<" #$# Fit #$# Vpre "; Vpre.Print();}
    for(int i=0; i<nTrk; i++) {
      TMatrixD BiGIi(BGI[i]);
      TMatrixD BiGIiUi(BiGIi,TMatrixD::kMult,U[i]);
      Vpre-=BiGIiUi;
      if(fDebug) {std::cout<<" #$# Fit #$# GIi ("<<i<<") "; GI[i].Print();}
      if(fDebug) {std::cout<<" #$# Fit #$# Bi ("<<i<<") "; B[i].Print();}
      if(fDebug) {std::cout<<" #$# Fit #$# BiGIi ("<<i<<") "; BiGIi.Print();}
      TMatrixD BiGIiBti(BiGIi,TMatrixD::kMultTranspose,B[i]);
      if(fDebug) {std::cout<<" #$# Fit #$# BiGIiBti ("<<i<<") = "; BiGIiBti.Print();}
      WV-=BiGIiBti;
      if(fDebug) {std::cout<<" #$# Fit #$# WV = "; WV.Print();}
      if(fDebug) {std::cout<<" #$# Fit #$# Vpre "; Vpre.Print();}
    }
    CovVV=WV;
    //CovVV*=1e9;//catch numerics
    CovVV.InvertFast(&determinant);
    //CovVV*=1e9;//catch numerics back
    if (determinant==0) {
      std::cout<<"RhoKalmanVtxFitter: WV Inversion failed, abort fit."<<std::endl;
      return kFALSE;
    }
    if(fDebug) {std::cout<<" #$# Fit #$# WV= "; WV.Print();}
    if(fDebug) {std::cout<<" #$# Fit #$# CovVV  (det="<<determinant<<") ="; CovVV.Print();}
    //TMatrixD CovVV(TMatrixD::kInverted,WV); // no determinant returned -> No check possible
    TMatrixD uV(CovVV,TMatrixD::kMult,Vpre);
    vtx.SetXYZ(vtx.X()+uV[0][0],vtx.Y()+uV[1][0],vtx.Z()+uV[2][0]);
    if(fDebug) {std::cout<<" #$# Fit #$# Vertex update = "; uV.Print();}
    if(fVerbose) {std::cout<<" #$# Fit #$# Vertex after          "; vtx.Print();}
    if(fVerbose) {std::cout<<" #$# Fit #$# Vertex after (global) "; vtx.Print();}

//     //skip momentim update until the very last iteration
//     if(iteration < fNIterations-1){
//       fExpansionPoint=vtx; //
//       continue;
//     }

    // Calculate Momentum updates & Chi^2
    std::vector<TMatrixD> uq;
    //dV.T();
    double chiq=0;
    for(int i=0; i<nTrk; i++) {
      TMatrixD BtiV(B[i],TMatrixD::kTransposeMult,uV);
      TMatrixD uPi(GI[i],TMatrixD::kMult,U[i]-BtiV);
      // update momenta
      momenta[i]+=TVector3(uPi[0][0],uPi[1][0],uPi[2][0]);
      if(fVerbose) {std::cout<<" #$# Fit #$# Momenum update on particle %i"<<i; uPi.Print();}
      //qip[3][i]=momenta[i].Phi();
      // update track parameters at current vertex (phi_0)
      // track param residuals for chisquare
      TMatrixD uvi(D[i],TMatrixD::kMult,uV);
      TMatrixD upi(E[i],TMatrixD::kMult,uPi);
      dq[i]-=uvi;
      dq[i]-=upi;
      //chisquare
      TMatrixD dqitWi(dq[i],TMatrixD::kTransposeMult,W[i]);
      TMatrixD chis(dqitWi,TMatrixD::kMult,dq[i]);
      chiq+=chis[0][0];
          if(fDebug){std::cout<<" #$# Fit #$# Insert Chisquare"<<std::endl;}
      InsertChi2(fCurrentHead->Daughter(i),chisquare);
    }
    if(chiq>0 && chiq<10000) { chisquare=chiq; }
    else { chisquare = -20; }
    // TODO: end iteration loop at satisfying chisq?
  }// end of iteratoin loop
  fChiSquare=chisquare;

  // calculating other cov matrices

  // POS COV
  for(int k=0; k<3; k++) for(int l=0; l<3; l++) {
      CovFitFull[k][l]=CovVV[k][l];
    }

  // MOM-POS COV
  for(int i=0; i<nTrk; i++) {
    TMatrixD CovVPi(CovVV,TMatrixD::kMult,BGI[i]);
    CovVPi*=-1.;
    for(int k=0; k<3; k++) for(int l=0; l<3; l++) {
        CovFitFull[k+3*(i+1)][l]=CovVPi[k][l]; ;
        CovFitFull[l][k+3*(i+1)]=CovVPi[k][l]; ;
      }
    // MOM-MOM COV
    for(int j=0; j<nTrk; j++) {
      // Caution: Identity of tracks i and j switch, to use CovVPi from above.
      TMatrixD CovPPijtmp(GI[j],TMatrixD::kMultTranspose,B[j]);
      TMatrixD CovPPij(CovPPijtmp,TMatrixD::kMult,CovVPi);
      CovPPij*=-1.;
      if(i==j) { CovPPij+=GI[j]; }
      for(int k=0; k<3; k++) for(int l=0; l<3; l++) {
          CovFitFull[k+3*(j+1)][l+3*(i+1)]=CovPPij[k][l]; // momentum cov
          CovFitFull[l+3*(i+1)][k+3*(j+1)]=CovPPij[k][l]; // momentum cov
        }
    }
  }
  if(fDebug) {std::cout<<" #$# Fit #$# CovFitFull: "; CovFitFull.Print();}


  //TODO... Continue checking

  // update particle candidates
  TMatrixD CovP7(7,7);
  for(int k=0; k<3; k++) for(int l=0; l<3; l++) {
      CovP7[k][l]=CovVV[k][l];
    }

  for(int i=0; i<nTrk; i++) {
    // The candidate tree is the to-be-fitted copy
    // It's safe to modify!
    RhoCandidate* tcand = fCurrentHead->Daughter(i);
    // set fitted values
    tcand->SetPos(vtx);
    if(fDebug) {std::cout<<" #$# Fit #$# Final Momentum :("<<i<<") "; momenta[i].Print();}
    tcand->SetP3(momenta[i]);

    for(int k=0; k<3; k++) for(int l=0; l<3; l++) {
        CovP7[k+3][l+3]=CovFitFull[3*(i+1)+k][3*(i+1)+l]; // momentum cov
        CovP7[k+3][l]=CovFitFull[3*(i+1)+k][l];   // momentum-position cov
        CovP7[l][k+3]=CovFitFull[l][3*(i+1)+k];   // momentum-position cov
      }
    // Calculate energy covariances:
    double invE = 1./tcand->E();
    CovP7[3][6] = CovP7[6][3] = (momenta[i].X()*CovP7[3][3]+momenta[i].Y()*CovP7[3][4]+momenta[i].Z()*CovP7[3][5])*invE;
    CovP7[4][6] = CovP7[6][4] = (momenta[i].X()*CovP7[4][3]+momenta[i].Y()*CovP7[4][4]+momenta[i].Z()*CovP7[4][5])*invE;
    CovP7[5][6] = CovP7[6][5] = (momenta[i].X()*CovP7[5][3]+momenta[i].Y()*CovP7[5][4]+momenta[i].Z()*CovP7[5][5])*invE;

    CovP7[6][6] = (momenta[i].X()*momenta[i].X()*CovP7[3][3]+momenta[i].Y()*momenta[i].Y()*CovP7[4][4]+momenta[i].Z()*momenta[i].Z()*CovP7[5][5]
                       +2.0*momenta[i].X()*momenta[i].Y()*CovP7[3][4]
                       +2.0*momenta[i].X()*momenta[i].Z()*CovP7[3][5]
                       +2.0*momenta[i].Y()*momenta[i].Z()*CovP7[4][5])*invE*invE;

    CovP7[6][0] = CovP7[0][6] = (momenta[i].X()*CovP7[3][0]+momenta[i].Y()*CovP7[4][0]+momenta[i].Z()*CovP7[5][0])*invE;
    CovP7[6][1] = CovP7[1][6] = (momenta[i].X()*CovP7[3][1]+momenta[i].Y()*CovP7[4][1]+momenta[i].Z()*CovP7[5][1])*invE;
    CovP7[6][2] = CovP7[2][6] = (momenta[i].X()*CovP7[3][2]+momenta[i].Y()*CovP7[4][2]+momenta[i].Z()*CovP7[5][2])*invE;

    tcand->SetCov7(CovP7);
    //tcand->Lock(); //?? autolocking
    //if(fDebug) {std::cout<<" #$# Fit #$# CovP7 with mom for traj "<<i<<": "; CovP7.Print();}
  }
  fCurrentHead->SetPos(vtx); // that's where the P4 is defined
  SetFourMomentumByDaughters(fCurrentHead);//Cov7 updated by daughters, too!
  SetDecayVertex(fCurrentHead,vtx,CovVV);
  return kTRUE; // all good now!
}

Double_t RhoFitterBase::Chi2Contribution ( const RhoCandidate *  b )

inherited

access to the fitted candidates

Definition at line 86 of file RhoFitterBase.cxx.

References Double_t, RhoFitterBase::fChi2Map, uid(), and RhoCandidate::Uid().

{
  if(!b) return -999.;
  Int_t uid = b->Uid();
  Double_t chi2=fChi2Map[uid];
  return chi2 >=0.0 ? chi2 : -1.;
}

RhoCandidate * RhoFitterBase::CopyCand ( RhoCandidate *  b )

protectedinherited

uppermost particle composite in tree

Definition at line 51 of file RhoFitterBase.cxx.

References RhoFactory::Instance(), RhoFactory::NewCandidate(), RhoCandidate::RemoveAssociations(), and RhoCandidate::SetFit().

Referenced by RhoFitterBase::CopyTree().

{
  RhoCandidate* newCand = RhoFactory::Instance()->NewCandidate ( b );
  newCand->RemoveAssociations();
  b->SetFit(newCand);//ready to be modified
  return newCand;
}

RhoCandidate * RhoFitterBase::CopyTree ( RhoCandidate *  head )

protectedinherited

Definition at line 61 of file RhoFitterBase.cxx.

References RhoFitterBase::CopyCand(), RhoCandidate::Daughter(), i, RhoCandidate::IsComposite(), RhoCandidate::NDaughters(), and RhoCandidate::SetMotherLink().

Referenced by RhoFitterBase::RhoFitterBase().

{
  //std::cout<<"\n\tcopy tree "<<head->Uid()<<" "<<&head<<" "<<head->PdgCode()<<" "<<head->NDaughters()<<"...";
  RhoCandidate* headcopy=CopyCand(head);
  RhoCandidate* daucopy=0;
  RhoCandidate* dau=0;
  for(Int_t i=0;i<head->NDaughters();i++)
  {
    dau=head->Daughter(i);
    //std::cout<<"  daugter "<<dau->Uid()<<" "<<i<<" "<<dau->PdgCode()<<" at "<<dau<<"  ";
    if(dau == head) {
      std::cout<<endl<<"*** Candidate is its own mother???  *** \n"<<std::endl;
      std::cout<<"  print:   "<<*head<<std::endl;;
    }
    assert(dau != head);
    if(dau->IsComposite()) daucopy=CopyTree(dau);
        else daucopy=CopyCand(dau);
        //std::cout<<"CopyTree: copied candidate "<<dau->Uid()<<std::endl;
    daucopy->SetMotherLink(headcopy); //daughter link is set automatically, too
  }
  return headcopy;
}

void RhoFitterBase::FindAndAddFinalStateDaughters ( RhoCandidate *  cand )

protectedinherited

Definition at line 149 of file RhoFitterBase.cxx.

References RhoCandidate::Daughter(), RhoFitterBase::fDaughters, RhoCandidate::IsComposite(), RhoCandidate::IsLocked(), and RhoCandidate::NDaughters().

Referenced by Rho4CFitter::Fit(), RhoKinFitter::Fit(), and Rho4CFitter::FitConserveMasses().

{
  RhoCandidate* tc;
  for (int k=0;k<cand->NDaughters();k++) {
    tc=cand->Daughter(k);
    if (!tc->IsComposite() || tc->IsLocked()) { fDaughters.push_back(tc); }
    else { FindAndAddFinalStateDaughters(tc); }
  }
  return;
}

Bool_t RhoFitterBase::Fit ( )

inherited

This function has to be overridden by the real algoritms implementations

Definition at line 95 of file RhoFitterBase.cxx.

References RhoFitterBase::fChi2Map, RhoFitterBase::fHeadOfTree, and RhoFitterBase::FitNode().

Referenced by ana_check(), ana_check_psi(), ana_complete(), ana_day1(), ana_invariantmass_2pi_tpc(), ana_multi(), anaDMesonsCharged(), anaideal_complete(), anatut_psi2s(), PndTripleAnaTask::DsDs2317Analysis(), PndPmtTask::Exec(), PndProdAnaTask::Exec(), PndLLbarAnaTask::Exec(), PndMyAnalysisTask::Exec(), PndSimpleCombinerTask::Exec(), DecayTreeFitter::Fitter::FitAll(), PndTripleAnaTask::JpsiAnalysis(), newana_check_eta(), poormantracks(), run_ana_eta_c_stt(), run_ana_eta_c_stt_v2(), run_ana_eta_c_tpc(), run_ana_invariantmass_2pi_stt(), run_ana_invariantmass_2pi_tpc(), run_ana_invariantmass_4pi_stt(), run_ana_invariantmass_4pi_tpc(), run_ana_mertens_evt7(), tut_ana(), tut_ana_fast(), tut_ana_fit(), tut_ana_ntp(), tut_ana_ntp_noqa(), and tut_ana_ntp_qa().

{
  fChi2Map.clear();
  return FitNode(fHeadOfTree);
}

Bool_t RhoFitterBase::FitAll ( )

inherited

Definition at line 101 of file RhoFitterBase.cxx.

References RhoFitterBase::fHeadOfTree, RhoCandidate::IsLocked(), and RhoFitterBase::IterateAndFit().

{
  if(fHeadOfTree->IsLocked())
  {
    Warning("RhoFitterBase::FitAll","You tried to fit a locked candidate. Retuning kFALSE now.");
    return kFALSE;
  }
  return IterateAndFit(fHeadOfTree);
}

Bool_t RhoKalmanVtxFitter::FitNode ( RhoCandidate *  b )

privatevirtual

Reimplemented from RhoFitterBase.

Definition at line 46 of file RhoKalmanVtxFitter.cxx.

References b, Bool_t, Calculate(), RhoFitterBase::fChiSquare, and fCurrentHead.

{
  fCurrentHead=b;
  Bool_t check = Calculate();
  if(check == kFALSE) {
    fChiSquare=-22222;
    Error("RhoKalmanVtxFitter::FitNode()","Fit failed for composite %p. Set chisquare to %f.",b,fChiSquare);
  }
  return check;
}

Double_t RhoKalmanVtxFitter::FitVertexFast	(	TVector3 &	vtx,
		TMatrixD &	cov,
		bool	skipcov = false,
		int	niter = 2
	)

Definition at line 75 of file RhoKalmanVtxFitter.cxx.

References Bool_t, c, CalcPrgParams(), cos(), RhoCandidate::Daughter(), Double_t, RhoFitterBase::fChiSquare, fCurrentHead, fDebug, fExpansionPoint, fJacobian, RhoFitterBase::fNDegreesOfFreedom, fPrgCov, fPrgParams, RhoFitterBase::fVerbose, i, RhoCandidate::NDaughters(), niter, printf(), s, sin(), t, and CAMath::Tan().

Referenced by PndPmtTask::Exec(), poormantracks(), and PndRhoTupleQA::qaPRG().

{
  // Calculate a vertex of n tracks without considering the changes in momentum vector
  // the variables vtx & cov (3x3) are written and the Chi^2 is returned.
  // the current head candidate is set by FitNode during tree navigation or is the composite
  // which the fitter was initialized with.
  int nTrk = fCurrentHead->NDaughters();
  fNDegreesOfFreedom=2*nTrk-3;
  for (int iterate=0;iterate<niter;iterate++){
    fExpansionPoint=vtx;
    std::vector<TMatrixD> w;
    std::vector<TMatrixD> xp;
    TMatrixD sumw(3,3);
    TMatrixD sumwx(3,1);
    Double_t determinant = 0.;
    for(int i=0; i<nTrk; i++) {
      RhoCandidate* tcand = fCurrentHead->Daughter(i);
      // get Helix in PRG notation from Billoir paper
      // (epsilon,z0,theta,phi0,rho)
      Bool_t testprg = CalcPrgParams(tcand,fExpansionPoint);
      if(!testprg) { return -2; }

      if (fVerbose) {
        printf("#$# Fast #$# Helix params:\n\t epsilon\t = %.4g cm\n\t Z0\t = %.4g cm\n\t theta\t = %.4g\n\t phi0\t = %.4g\n\t rho\t = %.4g\n",fPrgParams[0],fPrgParams[1],fPrgParams[2],fPrgParams[3],fPrgParams[4]);
        std::cout<<" #$# Fast #$# Helix cov: ";
        fPrgCov.Print();
      }

      Double_t s = sin(fPrgParams[3]);
      Double_t c = cos(fPrgParams[3]);
      Double_t t = TMath::Tan(fPrgParams[2]);
      if(t!=0.) { t = 1/t; }
      else { return -3; }   // t = cot(theta)

        TMatrixD COVi(2,2); // track parameter cov for (epsilon,z0)
        COVi[0][0]=fPrgCov[0][0];
        COVi[0][1]=fPrgCov[0][1];
        COVi[1][0]=fPrgCov[1][0];
        COVi[1][1]=fPrgCov[1][1];
        if(fVerbose) {std::cout<<" #$# Fast #$# COVi  "; COVi.Print();}

        TMatrixD Di(2,3); // Derivative in V, 2x3
        Di[0][0]=fJacobian[0][0]; // dEpsilon  / dvx
        Di[0][1]=fJacobian[0][1]; // dEpsilon   /dvy
        Di[0][2]=fJacobian[0][2]; // dEpsilon   /dvz
        Di[1][0]=fJacobian[1][0]; // dZ0   /dvx
        Di[1][1]=fJacobian[1][1]; // dZ0   /dvy
        Di[1][2]=fJacobian[1][2]; // dZ0   /dvz
        if(fDebug) {std::cout<<" #$# Fast #$# Di  "; Di.Print();}

        TMatrixD xpi(1,3);
        xpi[0][0]=s*fPrgParams[0];
        xpi[0][1]=-c*fPrgParams[0];
        xpi[0][2]=fPrgParams[1];
        xp.push_back(xpi); // save for later use
        if(fVerbose) {std::cout<<" #$# Fast #$# xpi "; xpi.Print();}

        //propagate to Vertex
        TMatrixD Wi(COVi);
        Wi.InvertFast(&determinant);
        if (determinant == 0.) {
          std::cout<<"RhoKalmanVtxFitter: COVi Inversion failed, abort fit -888"<<std::endl;
          return -888;
        }
        //TMatrixD Wi(TMatrixD::kInverted,COVi); // no determinant returned -> No check possible
        if(fDebug) {std::cout<<" #$# Fast #$# Wi"<<std::endl; Wi.Print();}
        TMatrixD wi(TMatrixD(Di,TMatrixD::kTransposeMult,Wi),TMatrixD::kMult,Di);
        w.push_back(wi); // save for later use
        if(fDebug) {std::cout<<" #$# Fast #$# wi  "; wi.Print();}
        TMatrixD wixpi(wi,TMatrixD::kMultTranspose,xpi);
        if(fDebug) {std::cout<<" #$# Fast #$# wixpi  "; wixpi.Print();}

        sumw+=wi; // sum up weights
        if(fDebug) {std::cout<<" #$# Fast #$# sumw  "; sumw.Print();}
        sumwx+=wixpi; // sum up weighted positions
        if(fDebug) {std::cout<<" #$# Fast #$# sumwx  "; sumwx.Print();}
    }

    TMatrixD cV(sumw);
    cV.InvertFast(&determinant);
    if (determinant==0) {
      std::cout<<"RhoKalmanVtxFitter: sumw Inversion failed, retunring -777"<<std::endl;
      return -777;
    }
    //TMatrixD cV(TMatrixD::kInverted,sumw); // no determinant returned -> No check possible
    if(fVerbose) {std::cout<<" #$# Fast #$# cV  "; cV.Print();}
    TMatrixD V(cV,TMatrixD::kMult,sumwx);// result vertex
    V.T(); // make it a row vector
    if(fVerbose) {std::cout<<" #$# Fast #$# V  "; V.Print();}
    vtx.SetXYZ(V[0][0],V[0][1],V[0][2]);
    vtx += fExpansionPoint; // move back to origin
    if(fVerbose) {std::cout<<" #$# Fast #$# Vtx: "; vtx.Print();}

    if (skipcov) {
      fChiSquare = -111.;// skip chisqare calculation
    } else {
      Double_t chisq=0.; // calculate chisquare
      for(int i=0; i<nTrk; i++) {
        if(fDebug) {std::cout<<" #$# Fast #$# V  "; V.Print();}
        TMatrixD resid(xp[i],TMatrixD::kMinus,V);
        if(fDebug) {std::cout<<" #$# Fast #$# resid  "; resid.Print();}
        if(fDebug) {std::cout<<" #$# Fast #$# w["<<i<<"]  "; w[i].Print();}
        TMatrixD chisqi(TMatrixD(resid,TMatrixD::kMult,w[i]),TMatrixD::kMultTranspose,resid);
        if(fDebug) {std::cout<<" #$# Fast #$# chisqi  "; chisqi.Print();}
        chisq+=chisqi[0][0];
      }
      if(fVerbose) {std::cout<<" #$# Fast #$# chisq = "<<chisq<<std::endl;}
      cov=cV;
      fChiSquare=chisq;
    }
  }
  return fChiSquare;

}

double RhoFitterBase::GetChi2 ( ) const

inlineinherited

Definition at line 48 of file RhoFitterBase.h.

References RhoFitterBase::fChiSquare.

Referenced by ana_check(), ana_check_psi(), ana_complete(), ana_day1(), ana_multi(), anaideal_complete(), anatut_psi2s(), PndTripleAnaTask::DsDs2317Analysis(), PndPmtTask::Exec(), PndProdAnaTask::Exec(), PndLLbarAnaTask::Exec(), PndAnaWithTrigger::Exec(), PndMyAnalysisTask::Exec(), PndScrutAnaTask::Exec(), PndSimpleCombinerTask::Exec(), PndTripleAnaTask::JpsiAnalysis(), newana_check_eta(), poormantracks(), PndRhoTupleQA::qaFitter(), run_ana_eta_c_stt(), run_ana_eta_c_stt_v2(), run_ana_eta_c_tpc(), PndTripleAnaTask::ThreePiAnalysis(), tut_ana(), tut_ana_fast(), tut_ana_fit(), tut_ana_ntp(), tut_ana_ntp_noqa(), and tut_ana_ntp_qa().

{return fChiSquare;};

int RhoFitterBase::GetNdf ( ) const

inlineinherited

Definition at line 49 of file RhoFitterBase.h.

References RhoFitterBase::fNDegreesOfFreedom.

Referenced by poormantracks(), and PndRhoTupleQA::qaFitter().

{return fNDegreesOfFreedom;};

double RhoFitterBase::GetProb ( ) const

inlineinherited

Definition at line 50 of file RhoFitterBase.h.

References RhoFitterBase::fChiSquare, and RhoFitterBase::fNDegreesOfFreedom.

Referenced by ana_complete(), ana_day1(), ana_multi(), anaideal_complete(), PndPmtTask::Exec(), PndLLbarAnaTask::Exec(), poormantracks(), PndRhoTupleQA::qaFitter(), tut_ana(), tut_ana_fast(), tut_ana_fit(), tut_ana_ntp(), tut_ana_ntp_noqa(), and tut_ana_ntp_qa().

{return TMath::Prob(fChiSquare,fNDegreesOfFreedom);};

RhoCandidate& RhoFitterBase::HeadOfTree ( ) const

inlineprotectedinherited

Definition at line 57 of file RhoFitterBase.h.

References RhoFitterBase::fHeadOfTree.

{ return *fHeadOfTree; }

void RhoFitterBase::InsertChi2 ( const RhoCandidate *  bc, const double  chi2  )

inlineprotectedinherited

Definition at line 66 of file RhoFitterBase.h.

References RhoFitterBase::fChi2Map, and RhoCandidate::Uid().

Referenced by Calculate().

{fChi2Map[ bc->Uid()] = chi2;}

void RhoFitterBase::SetDaugthersFromComposite ( RhoCandidate *  cand )

protectedinherited

Definition at line 136 of file RhoFitterBase.cxx.

References RhoCandidate::Daughter(), RhoFitterBase::fDaughters, RhoCandidate::IsComposite(), and RhoCandidate::NDaughters().

Referenced by RhoKinHyperonFitter::Fit(), RhoKinVtxFitter::FitNode(), RhoKinHyperonVtxFitter::FitNode(), and RhoKinHyperonFitter::SetMatrices().

{
  fDaughters.clear();
  if(cand->IsComposite()){
    RhoCandidate* tc;
    for (int k=0;k<cand->NDaughters();k++) {
      tc=cand->Daughter(k);
      fDaughters.push_back(tc);
    }
  }
  return;
}

void RhoKalmanVtxFitter::SetDebug ( bool  db = true )

inline

Definition at line 32 of file RhoKalmanVtxFitter.h.

References fDebug, and RhoFitterBase::fVerbose.

{if(db) { fVerbose=true; } fDebug=db;};

void RhoFitterBase::SetDecayVertex ( RhoCandidate *  composite, const TVector3 &  vtx, const TMatrixD &  CovVV  )

protectedinherited

Definition at line 178 of file RhoFitterBase.cxx.

References RhoCandidate::SetDecayVtx().

Referenced by Calculate(), RhoKinVtxFitter::SetOutput(), and RhoKinHyperonVtxFitter::SetOutput().

{
  RhoError decaypointcov(CovVV);
  RhoVector3Err decayvertex(vtx,decaypointcov);
  composite->SetDecayVtx(decayvertex);
}

void RhoKalmanVtxFitter::SetExpansionPoint ( TVector3  P )

inline

Definition at line 26 of file RhoKalmanVtxFitter.h.

References fExpansionPoint, and P.

Referenced by PndPmtTask::Exec().

{fExpansionPoint=P;}

void RhoFitterBase::SetFourMomentumByDaughters ( RhoCandidate *  composite )

protectedinherited

Definition at line 160 of file RhoFitterBase.cxx.

References RhoCandidate::Cov7(), RhoCandidate::Daughter(), RhoCandidate::IsComposite(), RhoCandidate::IsLocked(), RhoCandidate::NDaughters(), RhoCandidate::P4(), RhoCandidate::SetCov7(), and RhoCandidate::SetP4().

Referenced by Calculate(), Rho4CFitter::Fit(), Rho4CFitter::FitConserveMasses(), RhoKinVtxFitter::SetOutput(), and RhoKinFitter::SetOutput().

{
  RhoCandidate* tc;
  TLorentzVector tmpLV;
  TMatrixD tmpCov(7,7);
  // Sum daughter fourmomenta, dive into nodes if necessary
  for (int k=0;k<composite->NDaughters();k++) {
    tc=composite->Daughter(k);
    if (tc->IsComposite() && !tc->IsLocked()) { SetFourMomentumByDaughters(tc); }
    tmpLV += tc->P4();
    tmpCov = tmpCov + tc->Cov7();
  }
  composite->SetP4(tmpLV);
  composite->SetCov7(tmpCov);
  //std::cout<<" Base fitter cov7 from tc "<<tc->Uid()<<"/"<<tc->Charge()<<"/"<<tc->PdgCode()<<": ";tmpCov.Print();
  return;
}

void RhoKalmanVtxFitter::SetNIterations ( int  i )

inline

Definition at line 27 of file RhoKalmanVtxFitter.h.

References fNIterations, and i.

Referenced by PndPmtTask::Exec().

{fNIterations=i;};

void RhoKalmanVtxFitter::SetSilent ( )

inline

Definition at line 33 of file RhoKalmanVtxFitter.h.

References fDebug, and RhoFitterBase::fVerbose.

Referenced by PndPmtTask::Exec(), poormantracks(), and PndRhoTupleQA::qaPRG().

{fVerbose=false; fDebug=false;};

void RhoFitterBase::SetVerbose ( Bool_t  v = kTRUE )

inlineinherited

Definition at line 52 of file RhoFitterBase.h.

References RhoFitterBase::fVerbose, and v.

{fVerbose=v;}

Member Data Documentation

double RhoFitterBase::fChiSquare

protectedinherited

Definition at line 74 of file RhoFitterBase.h.

Referenced by Calculate(), DecayTreeFitter::Fitter::chiSquare(), RhoKinVtxFitter::Compute(), RhoKinHyperonVtxFitter::Compute(), Rho4CFitter::Do4CFit(), Rho4CFitter::Do4CFitWithMassConservation(), FitNode(), FitVertexFast(), RhoFitterBase::GetChi2(), RhoFitterBase::GetProb(), Rho4CFitter::Rho4CFitter(), RhoKinFitter::Solve(), and RhoKinHyperonFitter::Solve().

RhoCandidate* RhoKalmanVtxFitter::fCurrentHead

private

Definition at line 36 of file RhoKalmanVtxFitter.h.

Referenced by Calculate(), FitNode(), FitVertexFast(), and RhoKalmanVtxFitter().

std::vector<RhoCandidate*> RhoFitterBase::fDaughters

protectedinherited

Definition at line 69 of file RhoFitterBase.h.

Referenced by RhoKinVtxFitter::Compute(), RhoKinHyperonVtxFitter::Compute(), Rho4CFitter::Do4CFit(), Rho4CFitter::Do4CFitWithMassConservation(), RhoFitterBase::FindAndAddFinalStateDaughters(), Rho4CFitter::Fit(), RhoKinFitter::Fit(), RhoKinHyperonFitter::Fit(), Rho4CFitter::FitConserveMasses(), RhoKinVtxFitter::GetCovariance(), RhoKinHyperonVtxFitter::GetCovariance(), RhoKinFitter::Read4MomKinMatrix(), RhoKinHyperonFitter::Read4MomKinMatrix(), RhoKinVtxFitter::ReadKinMatrix(), RhoKinHyperonVtxFitter::ReadKinMatrix(), RhoKinVtxFitter::ReadMassKinMatrix(), RhoKinFitter::ReadMassKinMatrix(), RhoKinHyperonVtxFitter::ReadMassKinMatrix(), RhoKinVtxFitter::ReadMatrix(), RhoKinFitter::ReadMatrix(), RhoKinHyperonVtxFitter::ReadMatrix(), RhoKinFitter::ReadMomKinMatrix(), RhoKinHyperonFitter::ReadMomKinMatrix(), RhoKinFitter::ReadTotEKinMatrix(), RhoKinHyperonFitter::ReadTotEKinMatrix(), RhoKinFitter::ReadTotMomKinMatrix(), RhoKinHyperonFitter::ReadTotMomKinMatrix(), RhoFitterBase::SetDaugthersFromComposite(), RhoKinVtxFitter::SetMatrices(), RhoKinFitter::SetMatrices(), RhoKinHyperonFitter::SetMatrices(), RhoKinHyperonVtxFitter::SetMatrices(), RhoKinVtxFitter::SetOutput(), RhoKinFitter::SetOutput(), RhoKinHyperonFitter::SetOutput(), RhoKinHyperonVtxFitter::SetOutput(), RhoKinHyperonVtxFitter::TransportToPoca(), RhoKinVtxFitter::TransportToVertex(), and RhoKinHyperonVtxFitter::TransportToVertex().

bool RhoKalmanVtxFitter::fDebug

private

Definition at line 40 of file RhoKalmanVtxFitter.h.

Referenced by Calculate(), FitVertexFast(), SetDebug(), and SetSilent().

TVector3 RhoKalmanVtxFitter::fExpansionPoint

private

Definition at line 44 of file RhoKalmanVtxFitter.h.

Referenced by Calculate(), FitVertexFast(), and SetExpansionPoint().

RhoCandidate* RhoFitterBase::fHeadOfTree

protectedinherited

Definition at line 62 of file RhoFitterBase.h.

Referenced by Rho4CFitter::Fit(), RhoKinFitter::Fit(), RhoKinHyperonFitter::Fit(), RhoFitterBase::Fit(), RhoFitterBase::FitAll(), Rho4CFitter::FitConserveMasses(), RhoFitterBase::HeadOfTree(), RhoKinHyperonFitter::ReadMatrix(), RhoFitterBase::RhoFitterBase(), RhoKalmanVtxFitter(), RhoKinHyperonFitter::SetMatrices(), RhoKinFitter::SetOutput(), and RhoKinHyperonFitter::SetOutput().

TMatrixD RhoKalmanVtxFitter::fJacobian

private

Definition at line 47 of file RhoKalmanVtxFitter.h.

Referenced by CalcPrgParams(), Calculate(), and FitVertexFast().

int RhoFitterBase::fNDegreesOfFreedom

protectedinherited

Definition at line 75 of file RhoFitterBase.h.

Referenced by Calculate(), RhoKinVtxFitter::Compute(), RhoKinHyperonVtxFitter::Compute(), Rho4CFitter::Do4CFit(), Rho4CFitter::Do4CFitWithMassConservation(), FitVertexFast(), RhoFitterBase::GetNdf(), RhoFitterBase::GetProb(), RhoKinFitter::SetOutput(), and RhoKinHyperonFitter::SetOutput().

int RhoKalmanVtxFitter::fNIterations

private

Definition at line 43 of file RhoKalmanVtxFitter.h.

Referenced by Calculate(), and SetNIterations().

TMatrixD RhoKalmanVtxFitter::fPrgCov

private

Definition at line 46 of file RhoKalmanVtxFitter.h.

Referenced by CalcPrgParams(), Calculate(), and FitVertexFast().

Float_t RhoKalmanVtxFitter::fPrgParams[5]

private

Definition at line 45 of file RhoKalmanVtxFitter.h.

Referenced by CalcPrgParams(), Calculate(), and FitVertexFast().

Bool_t RhoFitterBase::fVerbose

protectedinherited

Definition at line 58 of file RhoFitterBase.h.

Referenced by Calculate(), RhoKinVtxFitter::Compute(), RhoKinHyperonVtxFitter::Compute(), FitVertexFast(), RhoKalmanVtxFitter(), SetDebug(), RhoKinVtxFitter::SetOutput(), RhoKinHyperonVtxFitter::SetOutput(), SetSilent(), and RhoFitterBase::SetVerbose().

---

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

• RhoKalmanVtxFitter.h
• RhoKalmanVtxFitter.cxx