/home/mweber/Dokumente/CMS/Alignment/kalmanalignment/trunk/simulation/Track.h

Go to the documentation of this file.

#ifndef _Track_h
#define _Track_h

#include <vector>

#include <CLHEP/Matrix/Vector.h>
#include <CLHEP/Matrix/Matrix.h>
#include <CLHEP/Matrix/SymMatrix.h>

#include "Det.h"
#include "Hit.h"

class Track {
protected:
  // track parameters on a given surface
  CLHEP::HepVector    fOrigin;
  CLHEP::HepVector    fDestination;
  CLHEP::HepSymMatrix fCovariance;
  double fMomentum;
  double fChi2;
  int fNdof;
  std::vector<Hit *> fHits;

public:
  Track();
  Track(CLHEP::HepVector origin, CLHEP::HepVector destination, double momentum);
  virtual ~Track();
  CLHEP::HepVector GetOrigin() const;
  CLHEP::HepVector GetDestination() const;
  double GetMomentum() const { return fMomentum; }
  double GetMass() const { return 0.105685; }
  int GetCharge() const { return 1; }
  double GetChi2() const { return fChi2; }
  int GetNdof() const { return fNdof; }
  void SetChi2Ndof(double chi2, int ndof) { fChi2 = chi2; fNdof = ndof; }
  Hit * MisalignedIntersection(const Det & adet) const;
  Hit * NominalIntersection(const Det & adet) const;
  CLHEP::HepVector * IntersectionWithPlane(const ReferenceFrame & frame) const;
  // Parameters
  // format is (tx, ty, x, y) 
  // where tx = dx/dz = tan(theta)*cos(phi)
  // where ty = dy/dz = tan(theta)*sin(phi)
  // x, y global x and y coordinates at z=0
  CLHEP::HepVector GetGlobalParameters() const;
  void SetGlobalParameters(CLHEP::HepVector p);
  void SetGlobalCovariance(CLHEP::HepSymMatrix covariance) { fCovariance = covariance; };
  CLHEP::HepSymMatrix GetGlobalCovariance() const { return fCovariance; }
  CLHEP::HepVector GetLocalParameters(const ReferenceFrame & frame) const;
  CLHEP::HepSymMatrix GetLocalCovariance(const ReferenceFrame & frame);
  void SetLocalCovariance(const ReferenceFrame & frame, const CLHEP::HepSymMatrix & covariance);
  // derivatives from local to global parameters
  CLHEP::HepMatrix Jacobian(const ReferenceFrame & frame) const;
  void SetLocalParameters(const ReferenceFrame & frame, CLHEP::HepVector p);
  const std::vector<Hit *> & GetHits() const { return fHits; };
  std::vector<Hit *> & GetHits() { return fHits; };
};

#endif