PixelTripletLowPtGenerator Class Reference

#include <PixelTripletLowPtGenerator.h>

Inheritance diagram for PixelTripletLowPtGenerator:

Public Member Functions
virtual void	hitTriplets (const TrackingRegion &region, OrderedHitTriplets &trs, const edm::Event &ev, const edm::EventSetup &es)
void	init (const HitPairGenerator &pairs, LayerCacheType *layerCache) override
const HitPairGenerator &	pairGenerator () const
	PixelTripletLowPtGenerator (const edm::ParameterSet &cfg, edm::ConsumesCollector &iC)
void	setSeedingLayers (SeedingLayerSetsHits::SeedingLayerSet pairLayers, std::vector< SeedingLayerSetsHits::SeedingLayer > thirdLayers) override
virtual	~PixelTripletLowPtGenerator ()
Public Member Functions inherited from HitTripletGeneratorFromPairAndLayers
virtual	~HitTripletGeneratorFromPairAndLayers ()
Public Member Functions inherited from HitTripletGenerator
virtual void	clear ()
	HitTripletGenerator (unsigned int size=500)
virtual void	hitTriplets (const TrackingRegion &reg, OrderedHitTriplets &prs, const edm::EventSetup &es)
virtual const OrderedHitTriplets &	run (const TrackingRegion &region, const edm::Event &ev, const edm::EventSetup &es)
virtual	~HitTripletGenerator ()
Public Member Functions inherited from OrderedHitsGenerator
	OrderedHitsGenerator ()
virtual	~OrderedHitsGenerator ()

Private Member Functions
GlobalPoint	getGlobalPosition (const TrackingRecHit *recHit)
void	getTracker (const edm::EventSetup &es)

Private Attributes
std::string	builderName
bool	checkClusterShape
bool	checkMultipleScattering
double	maxAngleRatio
double	nSigMultipleScattering
edm::ParameterSet	ps
double	rzTolerance
edm::EDGetTokenT < SiPixelClusterShapeCache >	theClusterShapeCacheToken
TripletFilter *	theFilter
LayerCacheType *	theLayerCache
std::vector < SeedingLayerSetsHits::SeedingLayer >	theLayers
HitPairGenerator *	thePairGenerator
const TrackerGeometry *	theTracker

Additional Inherited Members
Public Types inherited from HitTripletGeneratorFromPairAndLayers
typedef LayerHitMapCache	LayerCacheType
Public Attributes inherited from OrderedHitsGenerator
unsigned int	theMaxElement

Detailed Description

Definition at line 24 of file PixelTripletLowPtGenerator.h.

Constructor & Destructor Documentation

PixelTripletLowPtGenerator::PixelTripletLowPtGenerator	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &	iC
	)

Definition at line 25 of file PixelTripletLowPtGenerator.cc.

                                                                                                           :
  theTracker(nullptr), theFilter(nullptr), ps(cfg), thePairGenerator(nullptr), theLayerCache(nullptr),
  theClusterShapeCacheToken(iC.consumes<SiPixelClusterShapeCache>(cfg.getParameter<edm::InputTag>("clusterShapeCacheSrc")))
{}

virtual PixelTripletLowPtGenerator::~PixelTripletLowPtGenerator ( )

inlinevirtual

Definition at line 31 of file PixelTripletLowPtGenerator.h.

References theFilter, and thePairGenerator.

{ delete thePairGenerator; delete theFilter; }

Member Function Documentation

GlobalPoint PixelTripletLowPtGenerator::getGlobalPosition ( const TrackingRecHit *  recHit )

private

Definition at line 73 of file PixelTripletLowPtGenerator.cc.

References TrackingRecHit::geographicalId(), and TrackingRecHit::localPosition().

Referenced by hitTriplets().

{
  DetId detId = recHit->geographicalId();

  return
    theTracker->idToDet(detId)->toGlobal(recHit->localPosition());
}

void PixelTripletLowPtGenerator::getTracker ( const edm::EventSetup &  es )

private

Definition at line 54 of file PixelTripletLowPtGenerator.cc.

References edm::EventSetup::get(), edm::ESHandle< class >::product(), and patCandidatesForDimuonsSequences_cff::tracker.

Referenced by hitTriplets().

{
  if(theTracker == 0)
  {
    // Get tracker geometry
    edm::ESHandle<TrackerGeometry> tracker;
    es.get<TrackerDigiGeometryRecord>().get(tracker);

    theTracker = tracker.product();
  }

  if(theFilter == 0)
  {
    theFilter = new TripletFilter(es); 
  }
}

void PixelTripletLowPtGenerator::hitTriplets ( const TrackingRegion &  region, OrderedHitTriplets &  trs, const edm::Event &  ev, const edm::EventSetup &  es  )

virtual

Implements HitTripletGenerator.

Definition at line 82 of file PixelTripletLowPtGenerator.cc.

References builderName, dtNoiseDBValidation_cfg::cerr, checkClusterShape, checkMultipleScattering, TripletFilter::checkTrack(), edm::EventSetup::get(), edm::Event::getByToken(), getGlobalPosition(), HitInfo::getInfo(), ThirdHitPrediction::getRanges(), getTracker(), HitPairGenerator::hitPairs(), RecHitsSortedInPhi::hits(), i, ThirdHitPrediction::isCompatibleWithMultipleScattering(), DetLayer::location(), PixelRecoRange< T >::max(), maxAngleRatio, PixelRecoRange< T >::min(), nSigMultipleScattering, phi, edm::ESHandle< class >::product(), rzTolerance, OrderedHitPairs::size(), findQualityFiles::size, DetLayer::subDetector(), theClusterShapeCacheToken, theFilter, theLayerCache, theLayers, and thePairGenerator.

{

  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHand;
  es.get<IdealGeometryRecord>().get(tTopoHand);
  const TrackerTopology *tTopo=tTopoHand.product();

  edm::Handle<SiPixelClusterShapeCache> clusterShapeCache;
  ev.getByToken(theClusterShapeCacheToken, clusterShapeCache);

  // Generate pairs
  OrderedHitPairs pairs; pairs.reserve(30000);
  thePairGenerator->hitPairs(region,pairs,ev,es);

  if (pairs.size() == 0) return;

  int size = theLayers.size(); 

  // Set aliases
  const RecHitsSortedInPhi **thirdHitMap = new const RecHitsSortedInPhi*[size]; 
  for(int il=0; il<size; il++)
    thirdHitMap[il] = &(*theLayerCache)(theLayers[il], region, ev, es);

  // Get tracker
  getTracker(es);

  // Look at all generated pairs
  for(OrderedHitPairs::const_iterator ip = pairs.begin();
                                      ip!= pairs.end(); ip++)
  {
    // Fill rechits and points
    vector<const TrackingRecHit*> recHits(3);
    vector<GlobalPoint> points(3);

    recHits[0] = (*ip).inner()->hit();
    recHits[1] = (*ip).outer()->hit();

#ifdef Debug
    cerr << " RecHits " + HitInfo::getInfo(*recHits[0]) +
                          HitInfo::getInfo(*recHits[1]) << endl;
#endif

    for(int i=0; i<2; i++)
      points[i] = getGlobalPosition(recHits[i]);

    // Initialize helix prediction
    ThirdHitPrediction
      thePrediction(region,
                    points[0],points[1], es,
                    nSigMultipleScattering,maxAngleRatio,builderName);

    // Look at all layers
    for(int il=0; il<size; il++)
    {
      const DetLayer * layer = theLayers[il].detLayer();

#ifdef Debug
      cerr << "  check layer " << layer->subDetector()
                        << " " << layer->location() << endl;
#endif

      // Get ranges for the third hit
      float phi[2],rz[2];
      thePrediction.getRanges(layer, phi,rz);

      PixelRecoRange<float> phiRange(phi[0]              , phi[1]             );
      PixelRecoRange<float>  rzRange( rz[0] - rzTolerance, rz[1] + rzTolerance);

      // Get third hit candidates from cache
      typedef RecHitsSortedInPhi::Hit Hit;
      vector<Hit> thirdHits = thirdHitMap[il]->hits(phiRange.min(),phiRange.max());
      typedef vector<Hit>::const_iterator IH;

      for (IH th=thirdHits.begin(), eh=thirdHits.end(); th < eh; ++th) 
      {
        // Fill rechit and point
        recHits[2] = (*th)->hit();
        points[2]  = getGlobalPosition(recHits[2]);

#ifdef Debug
        cerr << "  third hit " + HitInfo::getInfo(*recHits[2]) << endl;
#endif

        // Check if third hit is compatible with multiple scattering
        vector<GlobalVector> globalDirs;
        if(thePrediction.isCompatibleWithMultipleScattering
             (points[2], recHits, globalDirs, es) == false)
        {
#ifdef Debug
          cerr << "  not compatible: multiple scattering" << endl;
#endif
          if(checkMultipleScattering) continue;
        }

        // Convert to localDirs
/*
        vector<LocalVector> localDirs;
        vector<GlobalVector>::const_iterator globalDir = globalDirs.begin();
        for(vector<const TrackingRecHit *>::const_iterator
                                            recHit  = recHits.begin();
                                            recHit != recHits.end(); recHit++)
        {
          localDirs.push_back(theTracker->idToDet(
                             (*recHit)->geographicalId())->toLocal(*globalDir));
          globalDir++;
        }
*/

        // Check if the cluster shapes are compatible with thrusts
        if(checkClusterShape)
        {
          if(! theFilter->checkTrack(recHits,globalDirs,tTopo, *clusterShapeCache))
          {
#ifdef Debug
            cerr << "  not compatible: cluster shape" << endl;
#endif
            continue;
          }
        }

        // All checks passed, put triplet back
        result.push_back(OrderedHitTriplet((*ip).inner(),(*ip).outer(),*th));
      }
    }
  } 
  delete [] thirdHitMap;

  return;
}

void PixelTripletLowPtGenerator::init ( const HitPairGenerator &  pairs, LayerCacheType *  layerCache  )

overridevirtual

Implements HitTripletGeneratorFromPairAndLayers.

Definition at line 31 of file PixelTripletLowPtGenerator.cc.

References builderName, checkClusterShape, checkMultipleScattering, HitPairGenerator::clone(), edm::ParameterSet::getParameter(), maxAngleRatio, nSigMultipleScattering, ps, rzTolerance, theLayerCache, and thePairGenerator.

{
  thePairGenerator = pairs.clone();
  theLayerCache    = layerCache;

  checkMultipleScattering = ps.getParameter<bool>("checkMultipleScattering");
  nSigMultipleScattering  = ps.getParameter<double>("nSigMultipleScattering");
  checkClusterShape       = ps.getParameter<bool>("checkClusterShape"); 
  rzTolerance             = ps.getParameter<double>("rzTolerance");
  maxAngleRatio           = ps.getParameter<double>("maxAngleRatio");
  builderName             = ps.getParameter<string>("TTRHBuilder");
}

const HitPairGenerator& PixelTripletLowPtGenerator::pairGenerator ( ) const

inline

Definition at line 40 of file PixelTripletLowPtGenerator.h.

References thePairGenerator.

{ return *thePairGenerator; }

void PixelTripletLowPtGenerator::setSeedingLayers ( SeedingLayerSetsHits::SeedingLayerSet  pairLayers, std::vector< SeedingLayerSetsHits::SeedingLayer >  thirdLayers  )

overridevirtual

Implements HitTripletGeneratorFromPairAndLayers.

Definition at line 46 of file PixelTripletLowPtGenerator.cc.

References HitPairGenerator::setSeedingLayers(), theLayers, and thePairGenerator.

                                                                                                           {
  thePairGenerator->setSeedingLayers(pairLayers);
  theLayers = thirdLayers;
}

Member Data Documentation

std::string PixelTripletLowPtGenerator::builderName

private

Definition at line 60 of file PixelTripletLowPtGenerator.h.

Referenced by hitTriplets(), and init().

bool PixelTripletLowPtGenerator::checkClusterShape

private

Definition at line 62 of file PixelTripletLowPtGenerator.h.

Referenced by hitTriplets(), and init().

bool PixelTripletLowPtGenerator::checkMultipleScattering

private

Definition at line 61 of file PixelTripletLowPtGenerator.h.

Referenced by hitTriplets(), and init().

double PixelTripletLowPtGenerator::maxAngleRatio

private

Definition at line 58 of file PixelTripletLowPtGenerator.h.

Referenced by hitTriplets(), and init().

double PixelTripletLowPtGenerator::nSigMultipleScattering

private

Definition at line 56 of file PixelTripletLowPtGenerator.h.

Referenced by hitTriplets(), and init().

edm::ParameterSet PixelTripletLowPtGenerator::ps

private

Definition at line 49 of file PixelTripletLowPtGenerator.h.

Referenced by init().

double PixelTripletLowPtGenerator::rzTolerance

private

Definition at line 57 of file PixelTripletLowPtGenerator.h.

Referenced by hitTriplets(), and init().

edm::EDGetTokenT<SiPixelClusterShapeCache> PixelTripletLowPtGenerator::theClusterShapeCacheToken

private

Definition at line 54 of file PixelTripletLowPtGenerator.h.

Referenced by hitTriplets().

TripletFilter* PixelTripletLowPtGenerator::theFilter

private

Definition at line 47 of file PixelTripletLowPtGenerator.h.

Referenced by hitTriplets(), and ~PixelTripletLowPtGenerator().

LayerCacheType* PixelTripletLowPtGenerator::theLayerCache

private

Definition at line 52 of file PixelTripletLowPtGenerator.h.

Referenced by hitTriplets(), and init().

std::vector<SeedingLayerSetsHits::SeedingLayer> PixelTripletLowPtGenerator::theLayers

private

Definition at line 51 of file PixelTripletLowPtGenerator.h.

Referenced by hitTriplets(), and setSeedingLayers().

HitPairGenerator* PixelTripletLowPtGenerator::thePairGenerator

private

Definition at line 50 of file PixelTripletLowPtGenerator.h.

Referenced by hitTriplets(), init(), pairGenerator(), setSeedingLayers(), and ~PixelTripletLowPtGenerator().

const TrackerGeometry* PixelTripletLowPtGenerator::theTracker

private

Definition at line 46 of file PixelTripletLowPtGenerator.h.