#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=default

Public Member Functions inherited from TrackingRegionProducer
virtual	~TrackingRegionProducer ()

Detailed Description

Definition at line 23 of file HITrackingRegionProducer.h.

Constructor & Destructor Documentation

◆ HITrackingRegionProducer()

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

inline

Definition at line 25 of file HITrackingRegionProducer.h.

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

       : theTtopoToken(iC.esConsumes()), theFieldToken(iC.esConsumes()) {
     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);
     if (thePrecise) {
       theMSMakerToken = iC.esConsumes();
     }
   }

◆ ~HITrackingRegionProducer()

HITrackingRegionProducer::~HITrackingRegionProducer ( )

overridedefault

Member Function Documentation

◆ estimateMultiplicity()

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

inline

Definition at line 50 of file HITrackingRegionProducer.h.

References edmNew::DetSetVector< T >::begin(), DetId::det(), edmNew::DetSetVector< T >::end(), makeMEIFBenchmarkPlots::ev, edm::EventSetup::getData(), hfClusterShapes_cfi::hits, pixelTopology::layer, DetId::subdetId(), theSiPixelRecHitsToken, and theTtopoToken.

Referenced by regions().

                                                                               {
     //rechits
     edm::Handle<SiPixelRecHitCollection> recHitColl;
     ev.getByToken(theSiPixelRecHitsToken, recHitColl);
 
     //Retrieve tracker topology from geometry
     const auto& tTopo = es.getData(theTtopoToken);
 
     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;
   }

◆ regions()

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

inlineoverridevirtual

Implements TrackingRegionProducer.

Definition at line 78 of file HITrackingRegionProducer.h.

References gpuPixelDoublets::cc, estimateMultiplicity(), makeMEIFBenchmarkPlots::ev, edm::EventSetup::getData(), LogTrace, Pi, mps_fire::result, theDirection, theFieldToken, theMSMakerToken, 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
       const auto& field = es.getData(theFieldToken);
       const MultipleScatteringParametrisationMaker* msmaker = nullptr;
       if (thePrecise) {
         msmaker = &es.getData(theMSMakerToken);
       }
 
       result.push_back(std::make_unique<RectangularEtaPhiTrackingRegion>(theDirection,
                                                                          theOrigin,
                                                                          thePtMin,
                                                                          theOriginRadius,
                                                                          theOriginHalfLength,
                                                                          etaB,
                                                                          phiB,
                                                                          field,
                                                                          msmaker,
                                                                          thePrecise));
     } else {  // global tracking
       LogTrace("heavyIonHLTVertexing") << " [HIVertexing: Global Tracking]";
       result.push_back(std::make_unique<GlobalTrackingRegion>(
           thePtMin, theOrigin, theOriginRadius, theOriginHalfLength, thePrecise));
     }
     return result;
   }

Member Data Documentation

◆ theDirection

GlobalVector HITrackingRegionProducer::theDirection

private

Definition at line 146 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

◆ theFieldToken

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> HITrackingRegionProducer::theFieldToken

private

Definition at line 139 of file HITrackingRegionProducer.h.

Referenced by regions().

◆ theMSMakerToken

edm::ESGetToken<MultipleScatteringParametrisationMaker, TrackerMultipleScatteringRecord> HITrackingRegionProducer::theMSMakerToken

private

Definition at line 140 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

◆ theOrigin

GlobalPoint HITrackingRegionProducer::theOrigin

private

Definition at line 142 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

◆ theOriginHalfLength

double HITrackingRegionProducer::theOriginHalfLength

private

Definition at line 144 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

◆ theOriginRadius

double HITrackingRegionProducer::theOriginRadius

private

Definition at line 143 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

◆ thePrecise

bool HITrackingRegionProducer::thePrecise

private

Definition at line 145 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

◆ thePtMin

double HITrackingRegionProducer::thePtMin

private

Definition at line 141 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer(), and regions().

◆ theSiPixelRecHits

edm::InputTag HITrackingRegionProducer::theSiPixelRecHits

private

Definition at line 136 of file HITrackingRegionProducer.h.

Referenced by HITrackingRegionProducer().

◆ theSiPixelRecHitsToken

edm::EDGetTokenT<SiPixelRecHitCollection> HITrackingRegionProducer::theSiPixelRecHitsToken

private

Definition at line 137 of file HITrackingRegionProducer.h.

Referenced by estimateMultiplicity(), and HITrackingRegionProducer().

◆ theTtopoToken

edm::ESGetToken<TrackerTopology, IdealGeometryRecord> HITrackingRegionProducer::theTtopoToken

private

Definition at line 138 of file HITrackingRegionProducer.h.

Referenced by estimateMultiplicity().

Private Attributes
GlobalVector	theDirection

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	theFieldToken

edm::ESGetToken< MultipleScatteringParametrisationMaker, TrackerMultipleScatteringRecord >	theMSMakerToken

GlobalPoint	theOrigin

double	theOriginHalfLength

double	theOriginRadius

bool	thePrecise

double	thePtMin

edm::InputTag	theSiPixelRecHits

edm::EDGetTokenT< SiPixelRecHitCollection >	theSiPixelRecHitsToken

edm::ESGetToken< TrackerTopology, IdealGeometryRecord >	theTtopoToken

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HITrackingRegionProducer()

◆ ~HITrackingRegionProducer()

Member Function Documentation

◆ estimateMultiplicity()

◆ regions()

Member Data Documentation

◆ theDirection

◆ theFieldToken

◆ theMSMakerToken

◆ theOrigin

◆ theOriginHalfLength

◆ theOriginRadius

◆ thePrecise

◆ thePtMin

◆ theSiPixelRecHits

◆ theSiPixelRecHitsToken

◆ theTtopoToken