/data/refman/pasoursint/CMSSW_4_1_8_patch9/src/TrackingTools/AnalyticalJacobians/interface/AnalyticalCurvilinearJacobian.h

Go to the documentation of this file.

#ifndef AnalyticalCurvilinearJacobian_H
#define AnalyticalCurvilinearJacobian_H

#include "TrackingTools/AnalyticalJacobians/interface/CurvilinearJacobian.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h" 
#include "DataFormats/GeometryVector/interface/GlobalVector.h" 

class GlobalTrajectoryParameters;

class AnalyticalCurvilinearJacobian : public CurvilinearJacobian {
 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);
  
  virtual ~AnalyticalCurvilinearJacobian() {}
  
  virtual const AlgebraicMatrix55& jacobian() const {return theJacobian;}
  virtual const AlgebraicMatrix jacobian_old() const {return asHepMatrix(theJacobian);}
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);

 private:
  
  AlgebraicMatrix55 theJacobian;

};  

#endif