TrajectoryStateOnSurface.h

Go to the documentation of this file.

#ifndef TrajectoryStateOnSurface_H
#define TrajectoryStateOnSurface_H

#include "TrackingTools/TrajectoryState/interface/BasicTrajectoryState.h"
#include "TrackingTools/TrajectoryState/interface/SurfaceSideDefinition.h"

#include <iosfwd>

class TrajectoryStateOnSurface : private  BasicTrajectoryState::Proxy 
{
  typedef BasicTrajectoryState::SurfaceSide SurfaceSide;
  typedef BasicTrajectoryState::Proxy             Base;

public:
  // construct
  TrajectoryStateOnSurface() {}
  TrajectoryStateOnSurface( BasicTrajectoryState* p) : Base(p) {}
  TrajectoryStateOnSurface( const FreeTrajectoryState& fts,
                const Surface& aSurface,
                const SurfaceSide side = SurfaceSideDefinition::atCenterOfSurface);
  TrajectoryStateOnSurface( const GlobalTrajectoryParameters& gp,
                const Surface& aSurface,
                const SurfaceSide side = SurfaceSideDefinition::atCenterOfSurface);
  TrajectoryStateOnSurface( const GlobalTrajectoryParameters& gp,
                const CartesianTrajectoryError& err,
                const Surface& aSurface,
                const SurfaceSide side = SurfaceSideDefinition::atCenterOfSurface);
  TrajectoryStateOnSurface( const GlobalTrajectoryParameters& gp,
                const CurvilinearTrajectoryError& err,
                const Surface& aSurface,
                const SurfaceSide side = SurfaceSideDefinition::atCenterOfSurface, 
                double weight = 1.);
  TrajectoryStateOnSurface( const GlobalTrajectoryParameters& gp,
                const CurvilinearTrajectoryError& err,
                const Surface& aSurface,
                double weight);
  TrajectoryStateOnSurface( const LocalTrajectoryParameters& p,
                const Surface& aSurface,
                const MagneticField* field,
                const SurfaceSide side = SurfaceSideDefinition::atCenterOfSurface);
  TrajectoryStateOnSurface( const LocalTrajectoryParameters& p,
                const LocalTrajectoryError& err,
                const Surface& aSurface,
                const MagneticField* field,
                const SurfaceSide side = SurfaceSideDefinition::atCenterOfSurface, 
                double weight = 1.);
  TrajectoryStateOnSurface( const LocalTrajectoryParameters& p,
                const LocalTrajectoryError& err,
                const Surface& aSurface,
                const MagneticField* field,
                double weight);

  virtual ~TrajectoryStateOnSurface() {}

#if defined( __GXX_EXPERIMENTAL_CXX0X__)

  TrajectoryStateOnSurface(TrajectoryStateOnSurface && rh) :
    Base(std::forward<Base>(rh)){}
    
  TrajectoryStateOnSurface & operator=(TrajectoryStateOnSurface && rh) {
    Base::swap(rh);
    return *this;
  }

#endif

  void swap(TrajectoryStateOnSurface & rh) {
    Base::swap(rh);
  }


  bool isValid() const {
    return Base::isValid() && data().isValid();
  }

  bool hasError() const {
    return data().hasError();
  }

  FreeTrajectoryState* freeState(bool withErrors = true) const {
    return data().freeTrajectoryState();
  }

  FreeTrajectoryState* freeTrajectoryState(bool withErrors = true) const {
    return freeState();
  }

  const MagneticField *magneticField() const { return data().magneticField(); }

  const GlobalTrajectoryParameters& globalParameters() const {
    return data().globalParameters();
  }
  GlobalPoint globalPosition() const {
    return data().globalPosition();
  }
  GlobalVector globalMomentum() const {
    return data().globalMomentum();
  }
  GlobalVector globalDirection() const {
    return data().globalDirection();
  }
  TrackCharge charge() const {
    return data().charge();
  }
  double signedInverseMomentum() const {
    return data().signedInverseMomentum();
  }
  double transverseCurvature() const {
    return data().transverseCurvature();
  }
  const CartesianTrajectoryError& cartesianError() const {
    return data().cartesianError();
  }
  const CurvilinearTrajectoryError& curvilinearError() const {
    return data().curvilinearError();
  }
  const LocalTrajectoryParameters& localParameters() const {
    return data().localParameters();
  }
  LocalPoint localPosition() const {
    return data().localPosition();
  }
  LocalVector localMomentum() const {
    return data().localMomentum();
  }
  LocalVector localDirection() const {
    return data().localDirection();
  }
  const LocalTrajectoryError& localError() const {
    return data().localError();
  }
  const Surface& surface() const {
    return data().surface();
  }

  double weight() const {return data().weight();} 

  void rescaleError(double factor) {
    unsharedData().rescaleError(factor);
  }

  std::vector<TrajectoryStateOnSurface> components() const {
    return data().components();
  }
  /*
  std::vector<TrajectoryStateOnSurface> components() const {
    std::vector<BasicTrajectoryState::RCPtr> c( data().components());
    std::vector<TrajectoryStateOnSurface> result; 
    result.reserve(c.size());
    for (std::vector<BasicTrajectoryState::RCPtr>::iterator i=c.begin();
     i != c.end(); i++) result.push_back(&(**i));
    return result;
  }
  */

  SurfaceSide surfaceSide() const {
    return data().surfaceSide();
  }

  void update( const LocalTrajectoryParameters& p,
           const Surface& aSurface,
           const MagneticField* field,
           const SurfaceSide side = SurfaceSideDefinition::atCenterOfSurface);
  void update( const LocalTrajectoryParameters& p,
           const LocalTrajectoryError& err,
               const Surface& aSurface,
               const MagneticField* field,
               const SurfaceSide side = SurfaceSideDefinition::atCenterOfSurface, 
               double weight = 1.);

};

inline void swap( TrajectoryStateOnSurface & rh,  TrajectoryStateOnSurface & lh) {
  // use base swap
  rh.swap(lh);
}

std::ostream& operator<<(std::ostream& os, const TrajectoryStateOnSurface& tsos);
#endif