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