AnalyticalCurvilinearJacobian.h

Go to the documentation of this file.

#ifndef AnalyticalCurvilinearJacobian_H
#define AnalyticalCurvilinearJacobian_H
 
#include "DataFormats/Math/interface/AlgebraicROOTObjects.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
 
class GlobalTrajectoryParameters;
 
class AnalyticalCurvilinearJacobian {
public:
  AnalyticalCurvilinearJacobian() : theJacobian(ROOT::Math::SMatrixNoInit()) {}
 
  AnalyticalCurvilinearJacobian(const GlobalTrajectoryParameters& globalParameters,
                                const GlobalPoint& x,
                                const GlobalVector& p,
                                const double& s);
  AnalyticalCurvilinearJacobian(const GlobalTrajectoryParameters& globalParameters,
                                const GlobalPoint& x,
                                const GlobalVector& p,
                                const GlobalVector& theFieldInInverseGeV,
                                const double& s);
 
public:
  void computeFullJacobian(
      const GlobalTrajectoryParameters&, const GlobalPoint&, const GlobalVector&, const GlobalVector&, const double& s);
  void computeInfinitesimalJacobian(
      const GlobalTrajectoryParameters&, const GlobalPoint&, const GlobalVector&, const GlobalVector&, const double& s);
  void computeStraightLineJacobian(const GlobalTrajectoryParameters&,
                                   const GlobalPoint&,
                                   const GlobalVector&,
                                   const double& s);
 
  const AlgebraicMatrix55& jacobian() const { return theJacobian; }
 
private:
  AlgebraicMatrix55 theJacobian;
};
 
#endif