/data/refman/pasoursint/CMSSW_5_2_7_hltpatch2/src/TrackingTools/TrajectoryState/interface/TrajectoryStateClosestToPoint.h

Go to the documentation of this file.

#ifndef TrajectoryStateClosestToPoint_H
#define TrajectoryStateClosestToPoint_H

#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
#include "TrackingTools/TrajectoryParametrization/interface/PerigeeTrajectoryParameters.h"
#include "TrackingTools/TrajectoryParametrization/interface/PerigeeTrajectoryError.h"
#include "TrackingTools/TrajectoryState/interface/PerigeeConversions.h"
#include "TrackingTools/TrajectoryParametrization/interface/TrajectoryStateExceptions.h"
#include <vector>

class TrajectoryStateClosestToPoint
{
  typedef TrajectoryStateOnSurface      TSOS;
  typedef FreeTrajectoryState           FTS;

public:

  TrajectoryStateClosestToPoint():
    valid(false), theFTSavailable(false), errorIsAvailable(false) {}

  TrajectoryStateClosestToPoint(const PerigeeTrajectoryParameters& perigeeParameters, double pt,
                                const GlobalPoint& referencePoint, const MagneticField* field);

  TrajectoryStateClosestToPoint(const PerigeeTrajectoryParameters& perigeeParameters, double pt,
    const PerigeeTrajectoryError& perigeeError, const GlobalPoint& referencePoint,
    const MagneticField* field);


  const GlobalPoint & referencePoint() const {
    if(!isValid()) throw TrajectoryStateException(
        "TrajectoryStateClosestToPoint is invalid and cannot return any parameters");
    return theRefPoint;
  }


  const PerigeeTrajectoryParameters & perigeeParameters() const {
    if(!isValid()) throw TrajectoryStateException(
        "TrajectoryStateClosestToPoint is invalid and cannot return any parameters");
    return theParameters;
  }

  double pt() const {
    if(!isValid()) throw TrajectoryStateException(
        "TrajectoryStateClosestToPoint is invalid and cannot return any parameters");
    return thePt;
  }

  const PerigeeTrajectoryError & perigeeError() const {
    if(!isValid()) throw TrajectoryStateException(
        "TrajectoryStateClosestToPoint is invalid and cannot return any parameters");
    if (!errorIsAvailable) throw TrajectoryStateException(
      "TrajectoryStateClosestToPoint: attempt to access errors when none available");
    return thePerigeeError;
  }

  GlobalPoint position() const {
    if(!isValid()) throw TrajectoryStateException(
        "TrajectoryStateClosestToPoint is invalid and cannot return any parameters");
    return perigeeConversions.positionFromPerigee(theParameters, theRefPoint);
  }


  GlobalVector momentum() const {
    if(!isValid()) throw TrajectoryStateException(
        "TrajectoryStateClosestToPoint is invalid and cannot return any parameters");
    return perigeeConversions.momentumFromPerigee(theParameters, thePt, theRefPoint);
  }


  TrackCharge charge() const {
    if(!isValid()) throw TrajectoryStateException(
        "TrajectoryStateClosestToPoint is invalid and cannot return any parameters");
    return theParameters.charge();
  }


  const FreeTrajectoryState & theState() const {
    if(!isValid()) throw TrajectoryStateException(
        "TrajectoryStateClosestToPoint is invalid and cannot return any parameters");
    if (!theFTSavailable) calculateFTS();
    return theFTS;
  }


  bool hasError() const {
    if(!isValid()) throw TrajectoryStateException(
        "TrajectoryStateClosestToPoint is invalid and cannot return any parameters");
    return errorIsAvailable;
  }

  bool isValid() const {
    return valid;
  }


private:

  friend class TrajectoryStateClosestToPointBuilder;
  friend class PerigeeConversions;

  TrajectoryStateClosestToPoint(const FTS& originalFTS, const GlobalPoint& referencePoint);

  void calculateFTS() const;

  bool valid;
  const MagneticField* theField;

  mutable FTS theFTS;
  mutable bool theFTSavailable;
  
  GlobalPoint theRefPoint;
  PerigeeTrajectoryParameters theParameters;
  double thePt;
  PerigeeTrajectoryError thePerigeeError;
  bool errorIsAvailable;
  PerigeeConversions perigeeConversions;
  
};
#endif