HITrackingRegionForPrimaryVtxProducer Class Reference

#include <HITrackingRegionForPrimaryVtxProducer.h>

Inheritance diagram for HITrackingRegionForPrimaryVtxProducer:

Public Member Functions
int	estimateMultiplicity (const edm::Event &ev, const edm::EventSetup &es) const
	HITrackingRegionForPrimaryVtxProducer (const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)
std::vector< std::unique_ptr< TrackingRegion > >	regions (const edm::Event &ev, const edm::EventSetup &es) const override
	~HITrackingRegionForPrimaryVtxProducer () override=default
Public Member Functions inherited from TrackingRegionProducer
virtual	~TrackingRegionProducer ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
bool	doVariablePtMin
edm::InputTag	theBeamSpotTag
edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotToken
GlobalVector	theDirection
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	theFieldToken
double	theFixedError
edm::ESGetToken< MultipleScatteringParametrisationMaker, TrackerMultipleScatteringRecord >	theMSMakerToken
double	theNSigmaZ
double	theOriginRadius
bool	thePrecise
double	thePtMin
double	theSigmaZVertex
edm::InputTag	theSiPixelRecHits
edm::EDGetTokenT< SiPixelRecHitCollection >	theSiPixelRecHitsToken
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	theTtopoToken
bool	theUseFixedError
bool	theUseFoundVertices
edm::InputTag	vertexCollName
edm::EDGetTokenT< reco::VertexCollection >	vertexCollToken

Detailed Description

Definition at line 34 of file HITrackingRegionForPrimaryVtxProducer.h.

Constructor & Destructor Documentation

HITrackingRegionForPrimaryVtxProducer()

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

inline

Definition at line 36 of file HITrackingRegionForPrimaryVtxProducer.h.

References looper::cfg, doVariablePtMin, edm::ParameterSet::getParameter(), theBeamSpotTag, theBeamSpotToken, theDirection, theFixedError, theMSMakerToken, theNSigmaZ, theOriginRadius, thePrecise, thePtMin, theSigmaZVertex, theSiPixelRecHits, theSiPixelRecHitsToken, theUseFixedError, theUseFoundVertices, vertexCollName, and vertexCollToken.

      : theTtopoToken(iC.esConsumes()), theFieldToken(iC.esConsumes()) {
    edm::ParameterSet regionPSet = cfg.getParameter<edm::ParameterSet>("RegionPSet");
    thePtMin = regionPSet.getParameter<double>("ptMin");
    theOriginRadius = regionPSet.getParameter<double>("originRadius");
    theNSigmaZ = regionPSet.getParameter<double>("nSigmaZ");
    theBeamSpotTag = regionPSet.getParameter<edm::InputTag>("beamSpot");
    theBeamSpotToken = iC.consumes<reco::BeamSpot>(theBeamSpotTag);
    thePrecise = regionPSet.getParameter<bool>("precise");
    theSiPixelRecHits = regionPSet.getParameter<edm::InputTag>("siPixelRecHits");
    theSiPixelRecHitsToken = iC.consumes<SiPixelRecHitCollection>(theSiPixelRecHits);
    doVariablePtMin = regionPSet.getParameter<bool>("doVariablePtMin");
    double xDir = regionPSet.getParameter<double>("directionXCoord");
    double yDir = regionPSet.getParameter<double>("directionYCoord");
    double zDir = regionPSet.getParameter<double>("directionZCoord");
    theDirection = GlobalVector(xDir, yDir, zDir);

    // for using vertex instead of beamspot
    theSigmaZVertex = regionPSet.getParameter<double>("sigmaZVertex");
    theFixedError = regionPSet.getParameter<double>("fixedError");
    theUseFoundVertices = regionPSet.getParameter<bool>("useFoundVertices");
    theUseFixedError = regionPSet.getParameter<bool>("useFixedError");
    vertexCollName = regionPSet.getParameter<edm::InputTag>("VertexCollection");
    vertexCollToken = iC.consumes<reco::VertexCollection>(vertexCollName);

    if (thePrecise) {
      theMSMakerToken = iC.esConsumes();
    }
  }

~HITrackingRegionForPrimaryVtxProducer()

HITrackingRegionForPrimaryVtxProducer::~HITrackingRegionForPrimaryVtxProducer ( )

overridedefault

Member Function Documentation

estimateMultiplicity()

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

inline

Definition at line 96 of file HITrackingRegionForPrimaryVtxProducer.h.

References edmNew::copyDetSetRange(), makeMEIFBenchmarkPlots::ev, edm::EventSetup::getData(), theSiPixelRecHitsToken, and theTtopoToken.

Referenced by regions().

                                                                              {
    //rechits
    edm::Handle<SiPixelRecHitCollection> recHitColl;
    ev.getByToken(theSiPixelRecHitsToken, recHitColl);

    auto const& ttopo = es.getData(theTtopoToken);

    std::vector<const TrackingRecHit*> theChosenHits;
    edmNew::copyDetSetRange(*recHitColl, theChosenHits, ttopo.pxbDetIdLayerComparator(1));
    return theChosenHits.size();
  }

fillDescriptions()

static void HITrackingRegionForPrimaryVtxProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

inlinestatic

Definition at line 68 of file HITrackingRegionForPrimaryVtxProducer.h.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                           {
    edm::ParameterSetDescription desc;

    desc.add<double>("ptMin", 0.7);
    desc.add<bool>("doVariablePtMin", true);
    desc.add<double>("originRadius", 0.2);
    desc.add<double>("nSigmaZ", 3.0);
    desc.add<edm::InputTag>("beamSpot", edm::InputTag("offlineBeamSpot"));
    desc.add<bool>("precise", true);
    desc.add<bool>("useMultipleScattering", false);
    desc.add<bool>("useFakeVertices", false);
    desc.add<edm::InputTag>("siPixelRecHits", edm::InputTag("siPixelRecHits"));
    desc.add<double>("directionXCoord", 1.0);
    desc.add<double>("directionYCoord", 1.0);
    desc.add<double>("directionZCoord", 0.0);
    desc.add<bool>("useFoundVertices", true);
    desc.add<edm::InputTag>("VertexCollection", edm::InputTag("hiPixelClusterVertex"));
    desc.add<bool>("useFixedError", true);
    desc.add<double>("fixedError", 3.0);
    desc.add<double>("sigmaZVertex", 3.0);

    // Only for backwards-compatibility
    edm::ParameterSetDescription descRegion;
    descRegion.add<edm::ParameterSetDescription>("RegionPSet", desc);

    descriptions.add("hiTrackingRegionFromClusterVtx", descRegion);
  }

regions()

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

inlineoverridevirtual

Implements TrackingRegionProducer.

Definition at line 108 of file HITrackingRegionForPrimaryVtxProducer.h.

References cms::cuda::bs, gpuPixelDoublets::cc, doVariablePtMin, estimateMultiplicity(), makeMEIFBenchmarkPlots::ev, edm::EventSetup::getData(), edm::HandleBase::isValid(), LogTrace, HiEvtPlane_cfi::minpt, Pi, mps_fire::result, theBeamSpotToken, theDirection, theFieldToken, theFixedError, theMSMakerToken, theNSigmaZ, theOriginRadius, thePrecise, thePtMin, theSigmaZVertex, theUseFixedError, theUseFoundVertices, spclusmultinvestigator_cfi::vertexCollection, and vertexCollToken.

                                                                                                {
    int estMult = estimateMultiplicity(ev, es);

    // from MC relating first layer pixel hits to "findable" sim tracks with pt>1 GeV
    float cc = -38.6447;
    float bb = 0.0581765;
    float aa = 1.34306e-06;

    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 = 60.;  //if we have more tracks -> regional tracking
    float etaB = 10.;
    float phiB = TMath::Pi() / 2.;

    float decEta = estTracks / 90.;
    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
    }

    float minpt = thePtMin;
    float varPtCutoff = 1500;  //cutoff
    if (doVariablePtMin && estMult < varPtCutoff) {
      minpt = 0.075;
      if (estMult > 0)
        minpt += estMult * (thePtMin - 0.075) / varPtCutoff;  // lower ptMin linearly with pixel hit multiplicity
    }

    // tracking region selection
    std::vector<std::unique_ptr<TrackingRegion> > result;
    double halflength;
    GlobalPoint origin;
    edm::Handle<reco::BeamSpot> bsHandle;
    ev.getByToken(theBeamSpotToken, bsHandle);
    if (bsHandle.isValid()) {
      const reco::BeamSpot& bs = *bsHandle;
      origin = GlobalPoint(bs.x0(), bs.y0(), bs.z0());
      halflength = theNSigmaZ * bs.sigmaZ();

      if (theUseFoundVertices) {
        edm::Handle<reco::VertexCollection> vertexCollection;
        ev.getByToken(vertexCollToken, vertexCollection);

        for (reco::VertexCollection::const_iterator iV = vertexCollection->begin(); iV != vertexCollection->end();
             iV++) {
          if (iV->isFake() || !iV->isValid())
            continue;
          origin = GlobalPoint(bs.x0(), bs.y0(), iV->z());
          halflength = (theUseFixedError ? theFixedError : (iV->zError()) * theSigmaZVertex);
        }
      }

      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, origin, thePtMin, theOriginRadius, halflength, etaB, phiB, field, msmaker, thePrecise));
      } else {  // global tracking
        LogTrace("heavyIonHLTVertexing") << " [HIVertexing: Global Tracking]";
        result.push_back(
            std::make_unique<GlobalTrackingRegion>(minpt, origin, theOriginRadius, halflength, thePrecise));
      }
    }
    return result;
  }

Member Data Documentation

doVariablePtMin

bool HITrackingRegionForPrimaryVtxProducer::doVariablePtMin

private

Definition at line 198 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

theBeamSpotTag

edm::InputTag HITrackingRegionForPrimaryVtxProducer::theBeamSpotTag

private

Definition at line 192 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer().

theBeamSpotToken

edm::EDGetTokenT<reco::BeamSpot> HITrackingRegionForPrimaryVtxProducer::theBeamSpotToken

private

Definition at line 193 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

theDirection

GlobalVector HITrackingRegionForPrimaryVtxProducer::theDirection

private

Definition at line 195 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

theFieldToken

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

private

Definition at line 207 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by regions().

theFixedError

double HITrackingRegionForPrimaryVtxProducer::theFixedError

private

Definition at line 201 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

theMSMakerToken

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

private

Definition at line 208 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

theNSigmaZ

double HITrackingRegionForPrimaryVtxProducer::theNSigmaZ

private

Definition at line 191 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

theOriginRadius

double HITrackingRegionForPrimaryVtxProducer::theOriginRadius

private

Definition at line 190 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

thePrecise

bool HITrackingRegionForPrimaryVtxProducer::thePrecise

private

Definition at line 194 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

thePtMin

double HITrackingRegionForPrimaryVtxProducer::thePtMin

private

Definition at line 189 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

theSigmaZVertex

double HITrackingRegionForPrimaryVtxProducer::theSigmaZVertex

private

Definition at line 200 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

theSiPixelRecHits

edm::InputTag HITrackingRegionForPrimaryVtxProducer::theSiPixelRecHits

private

Definition at line 196 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer().

theSiPixelRecHitsToken

edm::EDGetTokenT<SiPixelRecHitCollection> HITrackingRegionForPrimaryVtxProducer::theSiPixelRecHitsToken

private

Definition at line 197 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by estimateMultiplicity(), and HITrackingRegionForPrimaryVtxProducer().

theTtopoToken

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> HITrackingRegionForPrimaryVtxProducer::theTtopoToken

private

Definition at line 206 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by estimateMultiplicity().

theUseFixedError

bool HITrackingRegionForPrimaryVtxProducer::theUseFixedError

private

Definition at line 203 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

theUseFoundVertices

bool HITrackingRegionForPrimaryVtxProducer::theUseFoundVertices

private

Definition at line 202 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().

vertexCollName

edm::InputTag HITrackingRegionForPrimaryVtxProducer::vertexCollName

private

Definition at line 204 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer().

vertexCollToken

edm::EDGetTokenT<reco::VertexCollection> HITrackingRegionForPrimaryVtxProducer::vertexCollToken

private

Definition at line 205 of file HITrackingRegionForPrimaryVtxProducer.h.

Referenced by HITrackingRegionForPrimaryVtxProducer(), and regions().