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