LayerCrossingSide.h

Go to the documentation of this file.

#ifndef TkDetLayers_LayerCrossingSide_h
#define TkDetLayers_LayerCrossingSide_h

// temporary solution
#include "TrackingTools/GeomPropagators/interface/Propagator.h"

class LayerCrossingSide {
public:

  int barrelSide(const TrajectoryStateOnSurface& startingState, const Propagator& prop) const {
    GlobalPoint pos = startingState.globalPosition();
    GlobalVector radial(pos.x(), pos.y(), 0);
    if (startingState.globalMomentum().dot( radial) > 0) {  // momentum points outwards
      return (prop.propagationDirection() == alongMomentum ? 0 : 1);
    }
    else {  // momentum points inwards
      return (prop.propagationDirection() == oppositeToMomentum ? 0 : 1);
    }
  }

  int endcapSide(const TrajectoryStateOnSurface& startingState, const Propagator& prop) const {
    float zpos = startingState.globalPosition().z();
    if (startingState.globalMomentum().z() * zpos > 0) {  // momentum points outwards
      return (prop.propagationDirection() == alongMomentum ? 0 : 1);
    }
    else {  // momentum points inwards
      return (prop.propagationDirection() == oppositeToMomentum ? 0 : 1);
    }
  }

};

#endif