#include <HITrackingRegionProducer.h>

Inheritance diagram for HITrackingRegionProducer:

Public Member Functions
int	estimateMultiplicity (const edm::Event &ev, const edm::EventSetup &es) const

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

std::vector< std::unique_ptr< TrackingRegion > >	regions (const edm::Event &ev, const edm::EventSetup &es) const override

	~HITrackingRegionProducer () override

Public Member Functions inherited from TrackingRegionProducer
virtual	~TrackingRegionProducer ()

Private Attributes
GlobalVector	theDirection

GlobalPoint	theOrigin

double	theOriginHalfLength

double	theOriginRadius

bool	thePrecise

double	thePtMin

edm::InputTag	theSiPixelRecHits

edm::EDGetTokenT< SiPixelRecHitCollection >	theSiPixelRecHitsToken

Detailed Description

Definition at line 20 of file HITrackingRegionProducer.h.

Constructor & Destructor Documentation

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

inline

Definition at line 22 of file HITrackingRegionProducer.h.

References edm::ParameterSet::getParameter(), theDirection, theOrigin, theOriginHalfLength, theOriginRadius, thePrecise, thePtMin, theSiPixelRecHits, and theSiPixelRecHitsToken.

                                                                                   {
     edm::ParameterSet regionPSet = cfg.getParameter<edm::ParameterSet>("RegionPSet");
 
     thePtMin = regionPSet.getParameter<double>("ptMin");
     theOriginRadius = regionPSet.getParameter<double>("originRadius");
     theOriginHalfLength = regionPSet.getParameter<double>("originHalfLength");
     double xPos = regionPSet.getParameter<double>("originXPos");
     double yPos = regionPSet.getParameter<double>("originYPos");
     double zPos = regionPSet.getParameter<double>("originZPos");
     double xDir = regionPSet.getParameter<double>("directionXCoord");
     double yDir = regionPSet.getParameter<double>("directionYCoord");
     double zDir = regionPSet.getParameter<double>("directionZCoord");
     thePrecise = regionPSet.getParameter<bool>("precise");
     theSiPixelRecHits = regionPSet.getParameter<edm::InputTag>("siPixelRecHits");
     theSiPixelRecHitsToken = iC.consumes<SiPixelRecHitCollection>(theSiPixelRecHits);
     theOrigin = GlobalPoint(xPos, yPos, zPos);
     theDirection = GlobalVector(xDir, yDir, zDir);
   }

HITrackingRegionProducer::~HITrackingRegionProducer ( )

inlineoverride

Definition at line 41 of file HITrackingRegionProducer.h.

41 {}

Member Function Documentation

int HITrackingRegionProducer::estimateMultiplicity	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)		const

inline

Definition at line 43 of file HITrackingRegionProducer.h.

References edmNew::DetSetVector< T >::begin(), DetId::det(), edmNew::DetSet< T >::detId(), edmNew::DetSetVector< T >::end(), edm::EventSetup::get(), edm::Event::getByToken(), hfClusterShapes_cfi::hits, TrackerTopology::pxbLayer(), edmNew::DetSet< T >::size(), DetId::subdetId(), and theSiPixelRecHitsToken.

Referenced by regions().

                                                                               {
     //rechits
     edm::Handle<SiPixelRecHitCollection> recHitColl;
     ev.getByToken(theSiPixelRecHitsToken, recHitColl);
 
     //Retrieve tracker topology from geometry
     edm::ESHandle<TrackerTopology> tTopo;
     es.get<IdealGeometryRecord>().get(tTopo);
 
     int numRecHits = 0;
     //FIXME: this can be optimized quite a bit by looping only on the per-det 'items' of DetSetVector
     for (SiPixelRecHitCollection::const_iterator recHitIdIterator = recHitColl->begin(),
                                                  recHitIdIteratorEnd = recHitColl->end();
          recHitIdIterator != recHitIdIteratorEnd;
          recHitIdIterator++) {
       SiPixelRecHitCollection::DetSet hits = *recHitIdIterator;
       DetId detId = DetId(hits.detId());      // Get the Detid object
       unsigned int detType = detId.det();     // det type, tracker=1
       unsigned int subid = detId.subdetId();  //subdetector type, barrel=1, fpix=2
 
       unsigned int layer = 0;
       layer = tTopo->pxbLayer(detId);
       if (detType == 1 && subid == 1 && layer == 1) {
         numRecHits += hits.size();
       }
     }
     return numRecHits;
   }

std::vector<std::unique_ptr<TrackingRegion> > HITrackingRegionProducer::regions	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)		const

inlineoverridevirtual

Implements TrackingRegionProducer.

Definition at line 72 of file HITrackingRegionProducer.h.

References estimateMultiplicity(), RectangularEtaPhiTrackingRegion::kNever, LogTrace, Pi, mps_fire::result, theDirection, theOrigin, theOriginHalfLength, theOriginRadius, thePrecise, and thePtMin.

                                                                                                 {
     int estMult = estimateMultiplicity(ev, es);
 
     // fit from MC information
     float aa = 1.90935e-04;
     float bb = -2.90167e-01;
     float cc = 3.86125e+02;
 
     float estTracks = aa * estMult * estMult + bb * estMult + cc;
 
     LogTrace("heavyIonHLTVertexing") << "[HIVertexing]";
     LogTrace("heavyIonHLTVertexing") << " [HIVertexing: hits in the 1. layer:" << estMult << "]";
     LogTrace("heavyIonHLTVertexing") << " [HIVertexing: estimated number of tracks:" << estTracks << "]";
 
     float regTracking = 400.;  //if we have more tracks -> regional tracking
     float etaB = 10.;
     float phiB = TMath::Pi() / 2.;
 
     float decEta = estTracks / 600.;
     etaB = 2.5 / decEta;
 
     if (estTracks > regTracking) {
       LogTrace("heavyIonHLTVertexing") << " [HIVertexing: Regional Tracking]";
       LogTrace("heavyIonHLTVertexing") << "  [Regional Tracking: eta range: -" << etaB << ", " << etaB << "]";
       LogTrace("heavyIonHLTVertexing") << "  [Regional Tracking: phi range: -" << phiB << ", " << phiB << "]";
       LogTrace("heavyIonHLTVertexing") << "  [Regional Tracking: factor of decrease: " << decEta * 2.
                                        << "]";  // 2:from phi
     }
 
     // tracking region selection
     std::vector<std::unique_ptr<TrackingRegion> > result;
     if (estTracks > regTracking) {  // regional tracking
       result.push_back(std::make_unique<RectangularEtaPhiTrackingRegion>(
           theDirection,
           theOrigin,
           thePtMin,
           theOriginRadius,
           theOriginHalfLength,
           etaB,
           phiB,
           RectangularEtaPhiTrackingRegion::UseMeasurementTracker::kNever,
           thePrecise));
     } else {  // global tracking
       LogTrace("heavyIonHLTVertexing") << " [HIVertexing: Global Tracking]";
       result.push_back(std::make_unique<GlobalTrackingRegion>(
           thePtMin, theOrigin, theOriginRadius, theOriginHalfLength, thePrecise));
     }
     return result;
   }