HITrackingRegionProducer Class Reference

#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 ()

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

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, phase1PixelTopology::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 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().