/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_7/src/Geometry/TrackerGeometryBuilder/src/GeomDetLess.cc

Go to the documentation of this file.

#include "Geometry/TrackerGeometryBuilder/interface/GeomDetLess.h"
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
#include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
#include "FWCore/Utilities/interface/Exception.h"

bool GeomDetLess::insideOutLess( const GeomDet* a, const GeomDet* b) const
{
  if (a == b) return false;

  DetId ida(a->geographicalId());
  DetId idb(b->geographicalId());

  if      ( (ida.subdetId() == StripSubdetector::TIB || ida.subdetId() == StripSubdetector::TOB || ida.subdetId() == PixelSubdetector::PixelBarrel) &&
            (idb.subdetId() == StripSubdetector::TIB || idb.subdetId() == StripSubdetector::TOB || idb.subdetId() == PixelSubdetector::PixelBarrel)) {  // barrel with barrel
    return a->surface().position().perp() < b->surface().position().perp();
  }

  if      ( (ida.subdetId() == StripSubdetector::TEC || ida.subdetId() == StripSubdetector::TID || ida.subdetId() == PixelSubdetector::PixelEndcap) &&
            (idb.subdetId() == StripSubdetector::TEC || idb.subdetId() == StripSubdetector::TID || idb.subdetId() == PixelSubdetector::PixelEndcap)) {  // fwd with fwd
    return std::abs(a->surface().position().z()) < std::abs(b->surface().position().z());
  }
  
  //
  //  here I have 1 barrel against one forward
  //

  if      ( (ida.subdetId() == StripSubdetector::TIB || ida.subdetId() == StripSubdetector::TOB || ida.subdetId() == PixelSubdetector::PixelBarrel) &&
            (idb.subdetId() == StripSubdetector::TEC || idb.subdetId() == StripSubdetector::TID || idb.subdetId() == PixelSubdetector::PixelEndcap)) {  // barrel with barrel
    return barrelForwardLess( a, b);
  }else{
    return !barrelForwardLess( b, a);
  }
  
  throw cms::Exception("GeomDetLess") << "GeomDetLess: arguments are not Ok";
  
}

bool GeomDetLess::barrelForwardLess( const GeomDet* bla, 
                                     const GeomDet* flb) const
{
  return std::abs(bla->surface().position().z()) < std::abs( flb->position().z());
}