LowPtClusterShapeSeedComparitor Class Reference

#include <LowPtClusterShapeSeedComparitor.h>

Inheritance diagram for LowPtClusterShapeSeedComparitor:

Public Member Functions
virtual bool	compatible (const SeedingHitSet &hits, const edm::EventSetup &es)
	LowPtClusterShapeSeedComparitor (const edm::ParameterSet &ps)
virtual	~LowPtClusterShapeSeedComparitor ()
Public Member Functions inherited from SeedComparitor
virtual	~SeedComparitor ()

Private Member Functions
float	areaParallelogram (const Global2DVector &a, const Global2DVector &b)
std::vector< GlobalVector >	getGlobalDirs (const std::vector< GlobalPoint > &globalPoss)
std::vector< GlobalPoint >	getGlobalPoss (const TransientTrackingRecHit::ConstRecHitContainer &recHits)

Private Attributes
const ClusterShapeHitFilter *	theFilter

Detailed Description

Definition at line 23 of file LowPtClusterShapeSeedComparitor.h.

Constructor & Destructor Documentation

LowPtClusterShapeSeedComparitor::LowPtClusterShapeSeedComparitor

(

const edm::ParameterSet &

ps

)

Definition at line 23 of file LowPtClusterShapeSeedComparitor.cc.

{

}

LowPtClusterShapeSeedComparitor::~LowPtClusterShapeSeedComparitor ( )

virtual

Definition at line 29 of file LowPtClusterShapeSeedComparitor.cc.

{
}

Member Function Documentation

float LowPtClusterShapeSeedComparitor::areaParallelogram ( const Global2DVector &  a, const Global2DVector &  b  )

private

Definition at line 35 of file LowPtClusterShapeSeedComparitor.cc.

References PV2DBase< T, PVType, FrameType >::x(), and PV2DBase< T, PVType, FrameType >::y().

{  
  return a.x() * b.y() - a.y() * b.x();
}

bool LowPtClusterShapeSeedComparitor::compatible ( const SeedingHitSet &  hits, const edm::EventSetup &  es  )

virtual

Implements SeedComparitor.

Definition at line 120 of file LowPtClusterShapeSeedComparitor.cc.

References SeedingHitSet::container(), edm::EventSetup::get(), HitInfo::getInfo(), i, TrackingRecHit::isValid(), LogDebug, LogTrace, convertSQLiteXML::ok, and edm::ESHandle< class >::product().

{

  // Get cluster shape hit filter
  edm::ESHandle<ClusterShapeHitFilter> shape;
  es.get<CkfComponentsRecord>().get("ClusterShapeHitFilter",shape);
  theFilter = shape.product();

  const TransientTrackingRecHit::ConstRecHitContainer & thits = hits.container();

  // Get global positions
  vector<GlobalPoint>  globalPoss = getGlobalPoss(thits);

  // Get global directions
  vector<GlobalVector> globalDirs = getGlobalDirs(globalPoss);

  bool ok = true;

  // Check whether shape of pixel cluster is compatible
  // with local track direction

  if (globalDirs.size()!=globalPoss.size() || globalDirs.size()!=thits.size())
    {
      LogDebug("LowPtClusterShapeSeedComparitor")<<"not enough global dir calculated:"
                         <<"\nnHits: "<<thits.size()
                         <<"\nnPos: "<<globalPoss.size()
                         <<"\nnDir: "<<globalDirs.size()
                         <<" will say the seed is good anyway.";
      return true;
    }

  for(int i = 0; i < 3; i++)
  {
    const SiPixelRecHit* pixelRecHit =
      dynamic_cast<const SiPixelRecHit *>(thits[i]->hit());

    if (!pixelRecHit){
      edm::LogError("LowPtClusterShapeSeedComparitor")<<"this is not a pixel cluster";
      ok = false; break;
    }

    if(!pixelRecHit->isValid())
    { 
      ok = false; break; 
    }
    
    LogDebug("LowPtClusterShapeSeedComparitor")<<"about to compute compatibility."
                           <<"hit ptr: "<<pixelRecHit
                           <<"global direction:"<< globalDirs[i];


    if(! theFilter->isCompatible(*pixelRecHit, globalDirs[i]) )
    {
      LogTrace("LowPtClusterShapeSeedComparitor")
         << " clusShape is not compatible"
         << HitInfo::getInfo(*thits[i]->hit());

      ok = false; break;
    }
  }

  return ok;
}

vector< GlobalVector > LowPtClusterShapeSeedComparitor::getGlobalDirs ( const std::vector< GlobalPoint > &  globalPoss )

private

Definition at line 42 of file LowPtClusterShapeSeedComparitor.cc.

References trackerHits::c, CircleFromThreePoints::center(), PixelRecoUtilities::curvature(), CircleFromThreePoints::curvature(), dir, LogDebug, M_PI, mag2(), L1TEmulatorMonitor_cff::p, slope, funct::sqr(), mathSSE::sqrt(), v, PV2DBase< T, PVType, FrameType >::x(), Basic2DVector< T >::x(), PV2DBase< T, PVType, FrameType >::y(), and Basic2DVector< T >::y().

{
  // Get 2d points
  vector<Global2DVector> p;
  for(vector<GlobalPoint>::const_iterator ig = g.begin();
                                          ig!= g.end(); ig++)
     p.push_back( Global2DVector(ig->x(), ig->y()) );

  //
  vector<GlobalVector> globalDirs;

  // Determine circle
  CircleFromThreePoints circle(g[0],g[1],g[2]);

  if(circle.curvature() != 0.)
  {
    Global2DVector c (circle.center().x(), circle.center().y());

    float rad2 = (p[0] - c).mag2();
    float area = fabsf(areaParallelogram(p[1] - p[0], p[1] - c));

    float a12;
    if(area >= rad2) a12 = M_PI/2;
                else a12 = asin(area / rad2);

    float slope = (g[1].z() - g[0].z()) / a12;

    float cotTheta = slope * circle.curvature(); // == sinhEta
    float coshEta  = sqrt(1 + sqr(cotTheta));    // == 1/sinTheta

    // Calculate globalDirs
    float sinTheta =       1. / coshEta;
    float cosTheta = cotTheta * sinTheta;

    int dir;
    if(areaParallelogram(p[0] - c, p[1] - c) > 0) dir = 1; else dir = -1;

    float curvature = circle.curvature();

    for(vector<Global2DVector>::const_iterator ip = p.begin();
                                               ip!= p.end(); ip++)
    {
      Global2DVector v = (*ip - c)*curvature*dir;
      globalDirs.push_back(GlobalVector(-v.y()*sinTheta,
                                         v.x()*sinTheta,
                                               cosTheta));
    }
  }
  else{
    LogDebug("LowPtClusterShapeSeedComparitor")<<"the curvature is null:"
                           <<"\n point1: "<<g[0]
                           <<"\n point2: "<<g[1]
                           <<"\n point3: "<<g[2];
  }
  return globalDirs;
}

vector< GlobalPoint > LowPtClusterShapeSeedComparitor::getGlobalPoss ( const TransientTrackingRecHit::ConstRecHitContainer &  recHits )

private

Definition at line 101 of file LowPtClusterShapeSeedComparitor.cc.

{
  vector<GlobalPoint> globalPoss;

  for(TransientTrackingRecHit::ConstRecHitContainer::const_iterator recHit = thits.begin();
      recHit!= thits.end();
      recHit++)
  {
    DetId detId = (*recHit)->hit()->geographicalId();

    GlobalPoint gpos = (*recHit)->globalPosition();

    globalPoss.push_back(gpos);
  }

  return globalPoss;
}

Member Data Documentation

const ClusterShapeHitFilter* LowPtClusterShapeSeedComparitor::theFilter

private

Definition at line 38 of file LowPtClusterShapeSeedComparitor.h.