GENFIT
Rev:NoNumberAvailable
|
Go to the documentation of this file.
8 #ifndef GBLTRAJECTORY_H_
9 #define GBLTRAJECTORY_H_
16 #include "TMatrixDSymEigen.h"
28 GblTrajectory(
const std::vector<GblPoint> &aPointList,
bool flagCurv =
true,
29 bool flagU1dir =
true,
bool flagU2dir =
true);
30 GblTrajectory(
const std::vector<GblPoint> &aPointList,
unsigned int aLabel,
31 const TMatrixDSym &aSeed,
bool flagCurv =
true,
bool flagU1dir =
32 true,
bool flagU2dir =
true);
34 const std::vector<std::pair<std::vector<GblPoint>, TMatrixD> > &aPointaAndTransList);
36 const std::vector<std::pair<std::vector<GblPoint>, TMatrixD> > &aPointaAndTransList,
37 const TMatrixD &extDerivatives,
const TVectorD &extMeasurements,
38 const TVectorD &extPrecisions);
42 unsigned int getResults(
int aSignedLabel, TVectorD &localPar,
43 TMatrixDSym &localCov)
const;
44 unsigned int getMeasResults(
unsigned int aLabel,
unsigned int &numRes,
45 TVectorD &aResiduals, TVectorD &aMeasErrors, TVectorD &aResErrors,
46 TVectorD &aDownWeights);
47 unsigned int getScatResults(
unsigned int aLabel,
unsigned int &numRes,
48 TVectorD &aResiduals, TVectorD &aMeasErrors, TVectorD &aResErrors,
49 TVectorD &aDownWeights);
50 void getLabels(std::vector<unsigned int> &aLabelList);
51 void getLabels(std::vector<std::vector< unsigned int> > &aLabelList);
52 unsigned int fit(
double &Chi2,
int &Ndf,
double &lostWeight,
53 std::string optionList =
"");
86 std::pair<std::vector<unsigned int>, TMatrixD>
getJacobian(
87 int aSignedLabel)
const;
90 unsigned int nJacobian = 1)
const;
100 void getResAndErr(
unsigned int aData,
double &aResidual,
101 double &aMeadsError,
double &aResError,
double &aDownWeight);
unsigned int externalPoint
Label of external point (or 0)
GblTrajectory(const std::vector< GblPoint > &aPointList, bool flagCurv=true, bool flagU1dir=true, bool flagU2dir=true)
Create new (simple) trajectory from list of points.
void getLabels(std::vector< unsigned int > &aLabelList)
Get (list of) labels of points on (simple) trajectory.
Millepede-II (binary) record.
void calcJacobians()
Calculate Jacobians to previous/next scatterer from point to point ones.
std::vector< unsigned int > scatDataIndex
mapping points to data blocks from scatterers
unsigned int numMeasurements
Total number of measurements.
void printTrajectory(unsigned int level=0)
Print GblTrajectory.
unsigned int numParameters
Number of fit parameters.
std::pair< std::vector< unsigned int >, TMatrixD > getJacobian(int aSignedLabel) const
Get jacobian for transformation from fit to track parameters at point.
void getResAndErr(unsigned int aData, double &aResidual, double &aMeadsError, double &aResError, double &aDownWeight)
Get residual and errors from data block.
unsigned int numOffsets
Number of (points with) offsets on trajectory.
VVector theVector
Vector of linear equation system.
std::vector< TMatrixD > innerTransformations
Transformations at innermost points of.
void milleOut(MilleBinary &aMille)
Write valid trajectory to Millepede-II binary file.
std::vector< GblData > theData
List of data blocks.
void buildLinearEquationSystem()
Build linear equation system from data (blocks).
bool isValid() const
Retrieve validity of trajectory.
ROOT::Math::SMatrix< double, 2, 7 > SMatrix27
Namespace for the general broken lines package.
TVectorD externalMeasurements
unsigned int numInnerTrans
Number of inner transformations to external parameters.
unsigned int getNumPoints() const
Retrieve number of point from trajectory.
unsigned int numCurvature
Number of curvature parameters (0 or 1) or external parameters.
unsigned int getMeasResults(unsigned int aLabel, unsigned int &numRes, TVectorD &aResiduals, TVectorD &aMeasErrors, TVectorD &aResErrors, TVectorD &aDownWeights)
Get residuals at point from measurement.
std::vector< std::vector< GblPoint > > thePoints
(list of) List of points on trajectory
TVectorD externalPrecisions
(Symmetric) Bordered Band Matrix.
TMatrixD externalDerivatives
void construct()
Construct trajectory from list of points.
std::vector< unsigned int > theDimension
List of active dimensions (0=u1, 1=u2) in fit.
unsigned int getResults(int aSignedLabel, TVectorD &localPar, TMatrixDSym &localCov) const
Get fit results at point.
void prepare()
Prepare fit for simple or composed trajectory.
ROOT::Math::SMatrix< double, 5, 5 > SMatrix55
void predict()
Calculate predictions for all points.
std::vector< unsigned int > numPoints
Number of points on (sub)trajectory.
double downWeight(unsigned int aMethod)
Down-weight all points.
TMatrixDSym externalSeed
Precision (inverse covariance matrix) of external seed.
void defineOffsets()
Define offsets from list of points.
unsigned int numLocals
Total number of (additional) local parameters.
unsigned int numAllPoints
Number of all points on trajectory.
std::vector< unsigned int > measDataIndex
mapping points to data blocks from measurements
void getFitToLocalJacobian(std::vector< unsigned int > &anIndex, SMatrix55 &aJacobian, const GblPoint &aPoint, unsigned int measDim, unsigned int nJacobian=1) const
Get (part of) jacobian for transformation from (trajectory) fit to track parameters at point.
unsigned int fit(double &Chi2, int &Ndf, double &lostWeight, std::string optionList="")
Perform fit of trajectory.
void printData()
Print GblData blocks for trajectory.
void getFitToKinkJacobian(std::vector< unsigned int > &anIndex, SMatrix27 &aJacobian, const GblPoint &aPoint) const
Get jacobian for transformation from (trajectory) fit to kink parameters at point.
unsigned int numTrajectories
Number of trajectories (in composed trajectory)
BorderedBandMatrix theMatrix
(Bordered band) matrix of linear equation system
Simple Vector based on std::vector<double>
unsigned int getScatResults(unsigned int aLabel, unsigned int &numRes, TVectorD &aResiduals, TVectorD &aMeasErrors, TVectorD &aResErrors, TVectorD &aDownWeights)
Get (kink) residuals at point from scatterer.
bool fitOK
Trajectory has been successfully fitted (results are valid)
void printPoints(unsigned int level=0)
Print GblPoints on trajectory.
bool constructOK
Trajectory has been successfully constructed (ready for fit/output)