CAHitTripletGenerator Class Reference

#include <CAHitTripletGenerator.h>

Classes
class	QuantityDependsPt
class	QuantityDependsPtEval

Public Types
typedef LayerHitMapCache	LayerCacheType
typedef OrderedHitSeeds	ResultType

Public Member Functions
	CAHitTripletGenerator (const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)
	CAHitTripletGenerator (const edm::ParameterSet &cfg, edm::ConsumesCollector &iC)
void	hitNtuplets (const IntermediateHitDoublets &regionDoublets, std::vector< OrderedHitSeeds > &result, const edm::EventSetup &es, const SeedingLayerSetsHits &layers)
void	initEvent (const edm::Event &ev, const edm::EventSetup &es)
	~CAHitTripletGenerator ()=default

Static Public Member Functions
static void	fillDescriptions (edm::ParameterSetDescription &desc)
static const char *	fillDescriptionsLabel ()

Static Public Attributes
static constexpr unsigned int	minLayers = 3

Private Attributes
const float	caHardPtCut = 0.f
CACut	caPhiCut
CACut	caThetaCut
const float	extraHitRPhitolerance
const QuantityDependsPt	maxChi2
std::unique_ptr< SeedComparitor >	theComparitor
LayerCacheType	theLayerCache
const bool	useBendingCorrection

Detailed Description

Definition at line 30 of file CAHitTripletGenerator.h.

Member Typedef Documentation

LayerCacheType

typedef LayerHitMapCache CAHitTripletGenerator::LayerCacheType

Definition at line 32 of file CAHitTripletGenerator.h.

ResultType

typedef OrderedHitSeeds CAHitTripletGenerator::ResultType

Definition at line 35 of file CAHitTripletGenerator.h.

Constructor & Destructor Documentation

CAHitTripletGenerator() [1/2]

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

inline

Definition at line 38 of file CAHitTripletGenerator.h.

: CAHitTripletGenerator(cfg, iC) {}

CAHitTripletGenerator() [2/2]

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

Definition at line 28 of file CAHitTripletGenerator.cc.

    : extraHitRPhitolerance(cfg.getParameter<double>(
          "extraHitRPhitolerance")),  //extra window in ThirdHitPredictionFromCircle range (divide by R to get phi)
      maxChi2(cfg.getParameter<edm::ParameterSet>("maxChi2")),
      useBendingCorrection(cfg.getParameter<bool>("useBendingCorrection")),
      caThetaCut(cfg.getParameter<double>("CAThetaCut"),
                 cfg.getParameter<std::vector<edm::ParameterSet>>("CAThetaCut_byTriplets")),
      caPhiCut(cfg.getParameter<double>("CAPhiCut"),
               cfg.getParameter<std::vector<edm::ParameterSet>>("CAPhiCut_byTriplets")),
      caHardPtCut(cfg.getParameter<double>("CAHardPtCut")) {
  edm::ParameterSet comparitorPSet = cfg.getParameter<edm::ParameterSet>("SeedComparitorPSet");
  std::string comparitorName = comparitorPSet.getParameter<std::string>("ComponentName");
  if (comparitorName != "none") {
    theComparitor = SeedComparitorFactory::get()->create(comparitorName, comparitorPSet, iC);
  }
}

References looper::cfg, get, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, and theComparitor.

~CAHitTripletGenerator()

CAHitTripletGenerator::~CAHitTripletGenerator ( )

default

Member Function Documentation

fillDescriptions()

void CAHitTripletGenerator::fillDescriptions ( edm::ParameterSetDescription &  desc )

static

Definition at line 46 of file CAHitTripletGenerator.cc.

                                                                             {
  desc.add<double>("extraHitRPhitolerance", 0.06);
  desc.add<bool>("useBendingCorrection", false);
  desc.add<double>("CAThetaCut", 0.00125);
  desc.add<double>("CAPhiCut", 0.1);
 
  edm::ParameterSetDescription validatorCACut;
  validatorCACut.add<string>("seedingLayers", "BPix1+BPix2+BPix3");
  validatorCACut.add<double>("cut", 0.00125);
  std::vector<edm::ParameterSet> defaultCACutVector;
  edm::ParameterSet defaultCACut;
  defaultCACut.addParameter<string>("seedingLayers", "");
  defaultCACut.addParameter<double>("cut", -1.);
  defaultCACutVector.push_back(defaultCACut);
  desc.addVPSet("CAThetaCut_byTriplets", validatorCACut, defaultCACutVector);
  desc.addVPSet("CAPhiCut_byTriplets", validatorCACut, defaultCACutVector);
 
  desc.add<double>("CAHardPtCut", 0);
 
  edm::ParameterSetDescription descMaxChi2;
  descMaxChi2.add<double>("pt1", 0.8);
  descMaxChi2.add<double>("pt2", 2);
  descMaxChi2.add<double>("value1", 50);
  descMaxChi2.add<double>("value2", 8);
  descMaxChi2.add<bool>("enabled", true);
  desc.add<edm::ParameterSetDescription>("maxChi2", descMaxChi2);
 
  edm::ParameterSetDescription descComparitor;
  descComparitor.add<std::string>("ComponentName", "none");
  descComparitor.setAllowAnything();  // until we have moved SeedComparitor too to EDProducers
  desc.add<edm::ParameterSetDescription>("SeedComparitorPSet", descComparitor);
}

References edm::ParameterSetDescription::add(), edm::ParameterSet::addParameter(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

fillDescriptionsLabel()

static const char* CAHitTripletGenerator::fillDescriptionsLabel ( )

inlinestatic

Definition at line 44 of file CAHitTripletGenerator.h.

{ return "caHitTriplet"; }

hitNtuplets()

void CAHitTripletGenerator::hitNtuplets	(	const IntermediateHitDoublets &	regionDoublets,
		std::vector< OrderedHitSeeds > &	result,
		const edm::EventSetup &	es,
		const SeedingLayerSetsHits &	layers
	)

Definition at line 159 of file CAHitTripletGenerator.cc.

                                                                            {
  CAGraph g;
 
  std::vector<const HitDoublets*> hitDoublets;
 
  std::vector<CACell::CAntuplet> foundTriplets;
 
  int index = 0;
  for (const auto& regionLayerPairs : regionDoublets) {
    const TrackingRegion& region = regionLayerPairs.region();
    hitDoublets.clear();
    foundTriplets.clear();
 
    if (index == 0) {
      createGraphStructure(layers, g);
      caThetaCut.setCutValuesByLayerIds(g);
      caPhiCut.setCutValuesByLayerIds(g);
    } else {
      clearGraphStructure(layers, g);
    }
    fillGraph(layers, regionLayerPairs, g, hitDoublets);
    CellularAutomaton ca(g);
    ca.findTriplets(hitDoublets, foundTriplets, region, caThetaCut, caPhiCut, caHardPtCut);
 
    auto& allCells = ca.getAllCells();
 
    const QuantityDependsPtEval maxChi2Eval = maxChi2.evaluator(es);
 
    // re-used thoughout, need to be vectors because of RZLine interface
    std::array<float, 3> bc_r;
    std::array<float, 3> bc_z;
    std::array<float, 3> bc_errZ2;
    std::array<GlobalPoint, 3> gps;
    std::array<GlobalError, 3> ges;
    std::array<bool, 3> barrels;
 
    unsigned int numberOfFoundTriplets = foundTriplets.size();
    for (unsigned int tripletId = 0; tripletId < numberOfFoundTriplets; ++tripletId) {
      OrderedHitTriplet tmpTriplet(allCells[foundTriplets[tripletId][0]].getInnerHit(),
                                   allCells[foundTriplets[tripletId][0]].getOuterHit(),
                                   allCells[foundTriplets[tripletId][1]].getOuterHit());
 
      auto isBarrel = [](const unsigned id) -> bool { return id == PixelSubdetector::PixelBarrel; };
      for (unsigned int i = 0; i < 2; ++i) {
        auto const& ahit = allCells[foundTriplets[tripletId][i]].getInnerHit();
        gps[i] = ahit->globalPosition();
        ges[i] = ahit->globalPositionError();
        barrels[i] = isBarrel(ahit->geographicalId().subdetId());
      }
 
      auto const& ahit = allCells[foundTriplets[tripletId][1]].getOuterHit();
      gps[2] = ahit->globalPosition();
      ges[2] = ahit->globalPositionError();
      barrels[2] = isBarrel(ahit->geographicalId().subdetId());
 
      PixelRecoLineRZ line(gps[0], gps[2]);
      ThirdHitPredictionFromCircle predictionRPhi(gps[0], gps[2], extraHitRPhitolerance);
      const float curvature = predictionRPhi.curvature(ThirdHitPredictionFromCircle::Vector2D(gps[1].x(), gps[1].y()));
      const float abscurv = std::abs(curvature);
      const float thisMaxChi2 = maxChi2Eval.value(abscurv);
      float chi2 = std::numeric_limits<float>::quiet_NaN();
      // TODO: Do we have any use case to not use bending correction?
 
      if (useBendingCorrection) {
        // Following PixelFitterByConformalMappingAndLine
        const float simpleCot = (gps.back().z() - gps.front().z()) / (gps.back().perp() - gps.front().perp());
        const float pt = 1.f / PixelRecoUtilities::inversePt(abscurv, es);
        for (int i = 0; i < 3; ++i) {
          const GlobalPoint& point = gps[i];
          const GlobalError& error = ges[i];
          bc_r[i] = sqrt(sqr(point.x() - region.origin().x()) + sqr(point.y() - region.origin().y()));
          bc_r[i] += pixelrecoutilities::LongitudinalBendingCorrection(pt, es)(bc_r[i]);
          bc_z[i] = point.z() - region.origin().z();
          bc_errZ2[i] = (barrels[i]) ? error.czz() : error.rerr(point) * sqr(simpleCot);
        }
        RZLine rzLine(bc_r, bc_z, bc_errZ2, RZLine::ErrZ2_tag());
        chi2 = rzLine.chi2();
      } else {
        RZLine rzLine(gps, ges, barrels);
        chi2 = rzLine.chi2();
      }
 
      if (edm::isNotFinite(chi2) || chi2 > thisMaxChi2) {
        continue;
      }
 
      if (theComparitor) {
        if (!theComparitor->compatible(tmpTriplet)) {
          continue;
        }
      }
      result[index].emplace_back(tmpTriplet);
    }
    index++;
  }
}

References funct::abs(), caHardPtCut, caPhiCut, caThetaCut, hltPixelTracks_cff::chi2, RZLine::chi2(), ThirdHitPredictionFromCircle::curvature(), PixelRecoUtilities::curvature(), relativeConstraints::error, CAHitTripletGenerator::QuantityDependsPt::evaluator(), extraHitRPhitolerance, CellularAutomaton::findTriplets(), g, CellularAutomaton::getAllCells(), mps_fire::i, triggerObjects_cff::id, PixelRecoUtilities::inversePt(), PixelPluginsPhase0_cfi::isBarrel, edm::isNotFinite(), hgcalTopologyTester_cfi::layers, mps_splice::line, maxChi2, PixelSubdetector::PixelBarrel, point, DiDispStaMuonMonitor_cfi::pt, ihd::RegionLayerSets< T >::region(), HLT_FULL_cff::region, mps_fire::result, CACut::setCutValuesByLayerIds(), sqr(), mathSSE::sqrt(), theComparitor, useBendingCorrection, CAHitTripletGenerator::QuantityDependsPtEval::value(), x, and y.

initEvent()

void CAHitTripletGenerator::initEvent	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)

Definition at line 79 of file CAHitTripletGenerator.cc.

                                                                                 {
  if (theComparitor)
    theComparitor->init(ev, es);
}

References ev, and theComparitor.

Member Data Documentation

caHardPtCut

const float CAHitTripletGenerator::caHardPtCut = 0.f

private

Definition at line 130 of file CAHitTripletGenerator.h.

Referenced by hitNtuplets().

caPhiCut

CACut CAHitTripletGenerator::caPhiCut

private

Definition at line 129 of file CAHitTripletGenerator.h.

Referenced by hitNtuplets().

caThetaCut

CACut CAHitTripletGenerator::caThetaCut

private

Definition at line 128 of file CAHitTripletGenerator.h.

Referenced by hitNtuplets().

extraHitRPhitolerance

const float CAHitTripletGenerator::extraHitRPhitolerance

private

Definition at line 123 of file CAHitTripletGenerator.h.

Referenced by hitNtuplets().

maxChi2

const QuantityDependsPt CAHitTripletGenerator::maxChi2

private

Definition at line 125 of file CAHitTripletGenerator.h.

Referenced by hitNtuplets().

minLayers

constexpr unsigned int CAHitTripletGenerator::minLayers = 3

staticconstexpr

Definition at line 34 of file CAHitTripletGenerator.h.

theComparitor

std::unique_ptr<SeedComparitor> CAHitTripletGenerator::theComparitor

private

Definition at line 56 of file CAHitTripletGenerator.h.

Referenced by CAHitTripletGenerator(), hitNtuplets(), and initEvent().

theLayerCache

LayerCacheType CAHitTripletGenerator::theLayerCache

private

Definition at line 54 of file CAHitTripletGenerator.h.

useBendingCorrection

const bool CAHitTripletGenerator::useBendingCorrection

private

Definition at line 126 of file CAHitTripletGenerator.h.

Referenced by hitNtuplets().