/data/refman/pasoursint/CMSSW_5_3_1/src/TrackingTools/AnalyticalJacobians/interface/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(AlgebraicMatrixID()){}

  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