#include <PixelTripletLowPtGenerator.h>

Inheritance diagram for PixelTripletLowPtGenerator:

Public Member Functions
void	hitTriplets (const TrackingRegion &region, OrderedHitTriplets &trs, const edm::Event &ev, const edm::EventSetup &es, const SeedingLayerSetsHits::SeedingLayerSet &pairLayers, const std::vector< SeedingLayerSetsHits::SeedingLayer > &thirdLayers) override

void	hitTriplets (const TrackingRegion &region, OrderedHitTriplets &result, const edm::EventSetup &es, const HitDoublets &doublets, const RecHitsSortedInPhi *thirdHitMap, const std::vector< const DetLayer > &thirdLayerDetLayer, const int nThirdLayers) override

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

	~PixelTripletLowPtGenerator () override

Public Member Functions inherited from HitTripletGeneratorFromPairAndLayers
	HitTripletGeneratorFromPairAndLayers (unsigned int maxElement=0)

	HitTripletGeneratorFromPairAndLayers (const edm::ParameterSet &pset)

void	init (std::unique_ptr< HitPairGeneratorFromLayerPair > &&pairs, LayerCacheType *layerCache)

const HitPairGeneratorFromLayerPair &	pairGenerator () const

virtual	~HitTripletGeneratorFromPairAndLayers ()

Private Member Functions
GlobalPoint	getGlobalPosition (const TrackingRecHit *recHit)

void	getTracker (const edm::EventSetup &es)

Private Attributes
bool	checkClusterShape

bool	checkMultipleScattering

edm::ESGetToken< ClusterShapeHitFilter, CkfComponentsRecord >	m_clusterFilterToken

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	m_geomToken

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	m_magfieldToken

edm::ESGetToken< MultipleScatteringParametrisationMaker, TrackerMultipleScatteringRecord >	m_msmakerToken

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	m_topoToken

edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	m_ttrhBuilderToken

double	maxAngleRatio

double	nSigMultipleScattering

double	rzTolerance

edm::EDGetTokenT< SiPixelClusterShapeCache >	theClusterShapeCacheToken

std::unique_ptr< TripletFilter >	theFilter

const TrackerGeometry *	theTracker

Additional Inherited Members
Public Types inherited from HitTripletGeneratorFromPairAndLayers
typedef LayerHitMapCache	LayerCacheType

Static Public Member Functions inherited from HitTripletGeneratorFromPairAndLayers
static void	fillDescriptions (edm::ParameterSetDescription &desc)

Protected Attributes inherited from HitTripletGeneratorFromPairAndLayers
LayerCacheType *	theLayerCache

const unsigned int	theMaxElement

std::unique_ptr< HitPairGeneratorFromLayerPair >	thePairGenerator

Detailed Description

Definition at line 35 of file PixelTripletLowPtGenerator.h.

Constructor & Destructor Documentation

◆ PixelTripletLowPtGenerator()

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

Definition at line 25 of file PixelTripletLowPtGenerator.cc.

References looper::cfg, checkClusterShape, checkMultipleScattering, maxAngleRatio, nSigMultipleScattering, and rzTolerance.

     : HitTripletGeneratorFromPairAndLayers(),  // no theMaxElement used in this class
       m_geomToken(iC.esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>()),
       m_topoToken(iC.esConsumes<TrackerTopology, TrackerTopologyRcd>()),
       m_magfieldToken(iC.esConsumes()),
       m_ttrhBuilderToken(iC.esConsumes(edm::ESInputTag("", cfg.getParameter<string>("TTRHBuilder")))),
       m_msmakerToken(iC.esConsumes()),
       m_clusterFilterToken(iC.esConsumes(edm::ESInputTag("", "ClusterShapeHitFilter"))),
       theTracker(nullptr),
       theClusterShapeCacheToken(
           iC.consumes<SiPixelClusterShapeCache>(cfg.getParameter<edm::InputTag>("clusterShapeCacheSrc"))) {
   checkMultipleScattering = cfg.getParameter<bool>("checkMultipleScattering");
   nSigMultipleScattering = cfg.getParameter<double>("nSigMultipleScattering");
   checkClusterShape = cfg.getParameter<bool>("checkClusterShape");
   rzTolerance = cfg.getParameter<double>("rzTolerance");
   maxAngleRatio = cfg.getParameter<double>("maxAngleRatio");
 }

◆ ~PixelTripletLowPtGenerator()

PixelTripletLowPtGenerator::~PixelTripletLowPtGenerator ( )

override

Definition at line 44 of file PixelTripletLowPtGenerator.cc.

44 {}

Member Function Documentation

◆ getGlobalPosition()

GlobalPoint PixelTripletLowPtGenerator::getGlobalPosition ( const TrackingRecHit * recHit )

private

Definition at line 61 of file PixelTripletLowPtGenerator.cc.

References hcalRecHitTable_cff::detId, TrackerGeometry::idToDet(), rpcPointValidation_cfi::recHit, theTracker, and GeomDet::toGlobal().

Referenced by hitTriplets().

                                                                                       {
   DetId detId = recHit->geographicalId();
 
   return theTracker->idToDet(detId)->toGlobal(recHit->localPosition());
 }

◆ getTracker()

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

private

Definition at line 49 of file PixelTripletLowPtGenerator.cc.

References edm::EventSetup::getData(), m_clusterFilterToken, m_geomToken, theFilter, and theTracker.

Referenced by hitTriplets().

                                                                    {
   if (theTracker == nullptr) {
     // Get tracker geometry
     theTracker = &es.getData(m_geomToken);
   }
 
   if (!theFilter) {
     theFilter = std::make_unique<TripletFilter>(&es.getData(m_clusterFilterToken));
   }
 }

◆ hitTriplets() [1/2]

void PixelTripletLowPtGenerator::hitTriplets	(	const TrackingRegion &	region,
		OrderedHitTriplets &	trs,
		const edm::Event &	ev,
		const edm::EventSetup &	es,
		const SeedingLayerSetsHits::SeedingLayerSet &	pairLayers,
		const std::vector< SeedingLayerSetsHits::SeedingLayer > &	thirdLayers
	)

overridevirtual

Implements HitTripletGeneratorFromPairAndLayers.

Definition at line 68 of file PixelTripletLowPtGenerator.cc.

References DMR_cfg::cerr, checkClusterShape, checkMultipleScattering, SiPixelPhase1OnlineDQM_cff::clusterShapeCache, makeMEIFBenchmarkPlots::ev, edm::EventSetup::getData(), getGlobalPosition(), HitInfo::getInfo(), ThirdHitPrediction::getRanges(), getTracker(), RecHitsSortedInPhi::hits(), mps_fire::i, ThirdHitPrediction::isCompatibleWithMultipleScattering(), nano_mu_digi_cff::layer, m_magfieldToken, m_msmakerToken, m_topoToken, m_ttrhBuilderToken, volumeBasedMagneticField_160812_cfi::magfield, PixelRecoRange< T >::max(), maxAngleRatio, PixelRecoRange< T >::min(), nSigMultipleScattering, phi, hiPixelPairStep_cff::points, FastTrackerRecHitMaskProducer_cfi::recHits, nano_mu_digi_cff::region, mps_fire::result, rzTolerance, findQualityFiles::size, theClusterShapeCacheToken, theFilter, HitTripletGeneratorFromPairAndLayers::theLayerCache, and HitTripletGeneratorFromPairAndLayers::thePairGenerator.

                                                                                                              {
   //Retrieve tracker topology from geometry
   const TrackerTopology* tTopo = &es.getData(m_topoToken);
   const auto& magfield = es.getData(m_magfieldToken);
   const auto& ttrhBuilder = es.getData(m_ttrhBuilderToken);
   const auto& msmaker = es.getData(m_msmakerToken);
 
   edm::Handle<SiPixelClusterShapeCache> clusterShapeCache;
   ev.getByToken(theClusterShapeCacheToken, clusterShapeCache);
 
   // Generate pairs
   OrderedHitPairs pairs;
   pairs.reserve(30000);
   thePairGenerator->hitPairs(region, pairs, ev, es, pairLayers);
 
   if (pairs.empty())
     return;
 
   int size = thirdLayers.size();
 
   // Set aliases
   const RecHitsSortedInPhi** thirdHitMap = new const RecHitsSortedInPhi*[size];
   for (int il = 0; il < size; il++)
     thirdHitMap[il] = &(*theLayerCache)(thirdLayers[il], region);
 
   // 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], magfield, ttrhBuilder, nSigMultipleScattering, maxAngleRatio);
 
     // Look at all layers
     for (int il = 0; il < size; il++) {
       const DetLayer* layer = thirdLayers[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, msmaker) == 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;
 }

◆ hitTriplets() [2/2]

void PixelTripletLowPtGenerator::hitTriplets	(	const TrackingRegion &	region,
		OrderedHitTriplets &	result,
		const edm::EventSetup &	es,
		const HitDoublets &	doublets,
		const RecHitsSortedInPhi **	thirdHitMap,
		const std::vector< const DetLayer *> &	thirdLayerDetLayer,
		const int	nThirdLayers
	)

overridevirtual

Implements HitTripletGeneratorFromPairAndLayers.

Definition at line 193 of file PixelTripletLowPtGenerator.cc.

References Exception.

                                                                      {
   throw cms::Exception("Error") << "PixelTripletLowPtGenerator::hitTriplets is not implemented \n";
 }