SimTrack.h

Go to the documentation of this file.

#ifndef SimTrack_H
#define SimTrack_H
 
#include "SimDataFormats/Track/interface/CoreSimTrack.h"
#include "DataFormats/Math/interface/Vector3D.h"
#include "DataFormats/Math/interface/LorentzVector.h"
#include "FWCore/Utilities/interface/Exception.h"
 
class SimTrack : public CoreSimTrack {
public:
  typedef CoreSimTrack Core;
 
  SimTrack();
  SimTrack(int ipart, const math::XYZTLorentzVectorD& p);
 
  SimTrack(int ipart, const math::XYZTLorentzVectorD& p, int iv, int ig);
 
  SimTrack(int ipart,
           const math::XYZTLorentzVectorD& p,
           int iv,
           int ig,
           const math::XYZVectorD& tkp,
           const math::XYZTLorentzVectorD& tkm);
 
  SimTrack(const CoreSimTrack& t, int iv, int ig);
 
  int vertIndex() const { return ivert; }
  bool noVertex() const { return ivert == -1; }
 
  int genpartIndex() const { return igenpart; }
  bool noGenpart() const { return igenpart == -1; }
 
  const math::XYZVectorD& trackerSurfacePosition() const { return tkposition; }
 
  const math::XYZTLorentzVectorD& trackerSurfaceMomentum() const { return tkmomentum; }
 
  inline void setTkPosition(const math::XYZVectorD& pos) { tkposition = pos; }
 
  inline void setTkMomentum(const math::XYZTLorentzVectorD& mom) { tkmomentum = mom; }
 
  inline void setVertexIndex(const int v) { ivert = v; }
 
  void setCrossedBoundaryVars(bool crossedBoundary,
                              int idAtBoundary,
                              math::XYZTLorentzVectorF positionAtBoundary,
                              math::XYZTLorentzVectorF momentumAtBoundary) {
    crossedBoundary_ = crossedBoundary;
    idAtBoundary_ = idAtBoundary;
    positionAtBoundary_ = positionAtBoundary;
    momentumAtBoundary_ = momentumAtBoundary;
  }
  bool crossedBoundary() const { return crossedBoundary_; }
  const math::XYZTLorentzVectorF& getPositionAtBoundary() const { return positionAtBoundary_; }
  const math::XYZTLorentzVectorF& getMomentumAtBoundary() const { return momentumAtBoundary_; }
  int getIDAtBoundary() const { return idAtBoundary_; }
 
private:
  int ivert;
  int igenpart;
 
  math::XYZVectorD tkposition;
  math::XYZTLorentzVectorD tkmomentum;
 
  bool crossedBoundary_;
  int idAtBoundary_;
  math::XYZTLorentzVectorF positionAtBoundary_;
  math::XYZTLorentzVectorF momentumAtBoundary_;
};
 
#include <iosfwd>
std::ostream& operator<<(std::ostream& o, const SimTrack& t);
 
#endif