CompareHitPairsY Class Reference

#include <CosmicHitPairGeneratorFromLayerPair.h>

Public Member Functions
	CompareHitPairsY (const edm::EventSetup &iSetup)
bool	operator() (const OrderedHitPair &h1, const OrderedHitPair &h2)

Private Attributes
edm::ESHandle< TrackerGeometry >	tracker

Detailed Description

Definition at line 14 of file CosmicHitPairGeneratorFromLayerPair.h.

Constructor & Destructor Documentation

CompareHitPairsY::CompareHitPairsY ( const edm::EventSetup &  iSetup )

inline

Definition at line 16 of file CosmicHitPairGeneratorFromLayerPair.h.

References edm::EventSetup::get(), and tracker.

{ iSetup.get<TrackerDigiGeometryRecord>().get(tracker); };

Member Function Documentation

bool CompareHitPairsY::operator() ( const OrderedHitPair &  h1, const OrderedHitPair &  h2  )

inline

Definition at line 17 of file CosmicHitPairGeneratorFromLayerPair.h.

References funct::abs(), TrackingRecHit::geographicalId(), TrackerGeometry::idToDet(), OrderedHitPair::inner(), TrackingRecHit::localPosition(), OrderedHitPair::outer(), GeomDet::toGlobal(), tracker, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                      {
    const TrackingRecHit* trh1i = h1.inner()->hit();
    const TrackingRecHit* trh2i = h2.inner()->hit();
    const TrackingRecHit* trh1o = h1.outer()->hit();
    const TrackingRecHit* trh2o = h2.outer()->hit();
    GlobalPoint in1p = tracker->idToDet(trh1i->geographicalId())->surface().toGlobal(trh1i->localPosition());
    GlobalPoint in2p = tracker->idToDet(trh2i->geographicalId())->surface().toGlobal(trh2i->localPosition());
    GlobalPoint ou1p = tracker->idToDet(trh1o->geographicalId())->surface().toGlobal(trh1o->localPosition());
    GlobalPoint ou2p = tracker->idToDet(trh2o->geographicalId())->surface().toGlobal(trh2o->localPosition());
    if (ou1p.y() * ou2p.y() < 0)
      return ou1p.y() > ou2p.y();
    else {
      float dist1 = 100 * std::abs(ou1p.z() - in1p.z()) - std::abs(ou1p.y()) - 0.1 * std::abs(in1p.y());
      float dist2 = 100 * std::abs(ou2p.z() - in2p.z()) - std::abs(ou2p.y()) - 0.1 * std::abs(in2p.y());
      return dist1 < dist2;
    }
  }

Member Data Documentation

edm::ESHandle<TrackerGeometry> CompareHitPairsY::tracker

private

Definition at line 36 of file CosmicHitPairGeneratorFromLayerPair.h.

Referenced by CompareHitPairsY(), and operator()().