PndTrkTrackFinder.cxx File Reference

#include "PndTrkTrackFinder.h"
#include "PndSttHit.h"
#include "PndSttTube.h"
#include "PndSttMapCreator.h"
#include "PndSdsHit.h"
#include "PndTrackCand.h"
#include "PndTrackCandHit.h"
#include "PndTrack.h"
#include "PndTrkConformalHit.h"
#include "PndTrkConformalHitList.h"
#include "PndTrkConformalTransform.h"
#include "PndTrkTools.h"
#include "PndTrkSkewHit.h"
#include "PndTrkFitter.h"
#include "FairRootManager.h"
#include "FairRunAna.h"
#include "FairRuntimeDb.h"
#include "TClonesArray.h"
#include "TVector3.h"
#include "TArc.h"
#include "TLine.h"
#include "TMarker.h"
#include "TSpectrum2.h"
#include "TSpectrum.h"
#include "TStopwatch.h"
#include "TH3D.h"
#include "TMinuit.h"
#include "PndTrkClusterList.h"
#include "PndTrkTrackList.h"
#include "PndTrkClean.h"
#include "PndTrkNeighboringMap.h"
#include "PndTrkIndivisibleHit.h"
#include <iostream>

Go to the source code of this file.

Functions
Double_t	fit_distance (float x, float y, Double_t *par)
void	Chi2Calculation (Int_t &, Double_t *, Double_t &f, Double_t *par, Int_t)
void	Chi2Calculation2 (Int_t &, Double_t *, Double_t &f, Double_t *par, Int_t)

Function Documentation

void Chi2Calculation	(	Int_t &	,
		Double_t *	,
		Double_t &	f,
		Double_t *	par,
		Int_t
	)

Definition at line 65 of file PndTrkTrackFinder.cxx.

References Double_t, and fit_distance().

Referenced by PndTrkTrackFinder::MinuitFit().

{


  TMatrixT<Double_t> *objtofit = (TMatrixT<Double_t> *) gMinuit->GetObjectFit();
  
  Double_t chi2 = 0;
  Int_t hitcounter = objtofit->GetNrows();

  for (Int_t ihit = 0; ihit < hitcounter; ihit++)
    { 
      double r_reco = fit_distance(objtofit[0][ihit][0], objtofit[0][ihit][1], par);
      double delta  = (objtofit[0][ihit][2] - r_reco)/objtofit[0][ihit][3];

      cout << "reco iso " << r_reco << endl;
      cout << "drift " << objtofit[0][ihit][2] << endl;
      cout << "delta " << delta << endl;
      chi2 += delta * delta;
  }
  f = chi2;
  return;
}

void Chi2Calculation2	(	Int_t &	,
		Double_t *	,
		Double_t &	f,
		Double_t *	par,
		Int_t
	)

Definition at line 88 of file PndTrkTrackFinder.cxx.

References Double_t, and CAMath::Sqrt().

Referenced by PndTrkTrackFinder::MinuitFit2().

{

  TMatrixT<Double_t> *objtofit = (TMatrixT<Double_t> *) gMinuit->GetObjectFit();
  
  Double_t chi2 = 0;
  Int_t hitcounter = objtofit->GetNrows();

  for (Int_t ihit = 0; ihit < hitcounter; ihit++)
    { 

      double delta = TMath::Sqrt((objtofit[0][ihit][0] - par[0])*(objtofit[0][ihit][0] - par[0])+(objtofit[0][ihit][1] - par[1])*(objtofit[0][ihit][1] - par[1])) - par[2]+ par[3] *  objtofit[0][ihit][2];
      if(objtofit[0][ihit][2] == 0) chi2 += (delta * delta * 12.);
      else chi2 += (delta*delta)/(pow(objtofit[0][ihit][3],2));
    }
 f = chi2;
  return;
}

Double_t fit_distance	(	float	x,
		float	y,
		Double_t *	par
	)

Definition at line 58 of file PndTrkTrackFinder.cxx.

References CAMath::Sqrt().

Referenced by Chi2Calculation().

{
  double value=(par[0] * x - y + par[1]) / TMath::Sqrt(par[0] * par[0] + 1);
  return value;
}