/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/TrackingTools/TrajectoryState/interface/TrajectoryStateClosestToBeamLine.h

Go to the documentation of this file.

#ifndef TrajectoryStateClosestToBeamLine_H
#define TrajectoryStateClosestToBeamLine_H

#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/Measurement1D.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"

class TrajectoryStateClosestToBeamLine
{
public:

  typedef FreeTrajectoryState           FTS;

  TrajectoryStateClosestToBeamLine() : valid(false) {}

  TrajectoryStateClosestToBeamLine
    (const FTS& stateAtPCA, const GlobalPoint & pointOnBeamLine,
     const reco::BeamSpot & beamSpot);

  ~TrajectoryStateClosestToBeamLine(){}

  FTS const & trackStateAtPCA() const {
    if (!isValid()) throw TrajectoryStateException(
      "TrajectoryStateClosestToBeamLine is invalid.");
    return theFTS;
  }

  GlobalPoint const & beamLinePCA() const {
    if (!isValid()) throw TrajectoryStateException(
      "TrajectoryStateClosestToBeamLine is invalid.");
    return thePointOnBeamLine;
  }

  Measurement1D transverseImpactParameter() const;

  reco::BeamSpot const & beamSpot() {
    if (!isValid()) throw TrajectoryStateException(
      "TrajectoryStateClosestToBeamLine is invalid.");
    return theBeamSpot;
  }

  inline bool isValid() const {return valid;}

private:

  bool valid;
  FTS theFTS;
  GlobalPoint thePointOnBeamLine;
  reco::BeamSpot theBeamSpot;

};
#endif