CAHitTripletGenerator Class Reference

#include <CAHitTripletGenerator.h>

Inheritance diagram for CAHitTripletGenerator:

Classes
class	QuantityDependsPt
class	QuantityDependsPtEval

Public Types
typedef LayerHitMapCache	LayerCacheType
typedef OrderedHitTriplets	ResultType

Public Member Functions
	CAHitTripletGenerator (const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC, bool needSeedingLayerSetsHits=true)
	CAHitTripletGenerator (const edm::ParameterSet &cfg, edm::ConsumesCollector &iC, bool needSeedingLayerSetsHits=true)
void	hitNtuplets (const IntermediateHitDoublets::RegionLayerSets &regionLayerPairs, OrderedHitTriplets &result, const edm::EventSetup &es, const SeedingLayerSetsHits &layers)
virtual void	hitTriplets (const TrackingRegion &reg, OrderedHitTriplets &triplets, const edm::Event &ev, const edm::EventSetup &es)
	from base class More...
void	initEvent (const edm::Event &ev, const edm::EventSetup &es)
virtual	~CAHitTripletGenerator ()
Public Member Functions inherited from HitTripletGenerator
virtual void	clear () final
	HitTripletGenerator (unsigned int size=500)
	HitTripletGenerator (HitTripletGenerator const &other)
virtual void	hitTriplets (const TrackingRegion &reg, OrderedHitTriplets &prs, const edm::EventSetup &es)
virtual const OrderedHitTriplets &	run (const TrackingRegion &region, const edm::Event &ev, const edm::EventSetup &es) final
virtual	~HitTripletGenerator ()
Public Member Functions inherited from OrderedHitsGenerator
	OrderedHitsGenerator ()
virtual	~OrderedHitsGenerator ()

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

Static Public Attributes
static unsigned int	minLayers = 3

Private Member Functions
void	hitTriplets (const TrackingRegion &reg, OrderedHitTriplets &result, std::vector< const HitDoublets * > &hitDoublets, const CAGraph &g, const edm::EventSetup &es)

Private Attributes
const float	caHardPtCut = 0.f
const float	caPhiCut = 1.f
const float	caThetaCut = 0.00125f
const float	extraHitRPhitolerance
const QuantityDependsPt	maxChi2
std::unique_ptr< SeedComparitor >	theComparitor
LayerCacheType	theLayerCache
edm::EDGetTokenT < SeedingLayerSetsHits >	theSeedingLayerToken
const bool	useBendingCorrection

Additional Inherited Members
Public Attributes inherited from OrderedHitsGenerator
unsigned int	theMaxElement

Detailed Description

Definition at line 31 of file CAHitTripletGenerator.h.

Member Typedef Documentation

typedef LayerHitMapCache CAHitTripletGenerator::LayerCacheType

Definition at line 33 of file CAHitTripletGenerator.h.

typedef OrderedHitTriplets CAHitTripletGenerator::ResultType

Definition at line 36 of file CAHitTripletGenerator.h.

Constructor & Destructor Documentation

CAHitTripletGenerator::CAHitTripletGenerator ( const edm::ParameterSet &  cfg, edm::ConsumesCollector &&  iC, bool  needSeedingLayerSetsHits = true  )

inline

Definition at line 40 of file CAHitTripletGenerator.h.

: CAHitTripletGenerator(cfg, iC, needSeedingLayerSetsHits) {}

CAHitTripletGenerator::CAHitTripletGenerator	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &	iC,
		bool	needSeedingLayerSetsHits = true
	)

Definition at line 36 of file CAHitTripletGenerator.cc.

References edm::ConsumesCollector::consumes(), beamerCreator::create(), edm::ParameterSet::exists(), reco::get(), edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, theComparitor, and theSeedingLayerToken.

                                                                            :
        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")),
        caPhiCut(cfg.getParameter<double>("CAPhiCut")),
        caHardPtCut(cfg.getParameter<double>("CAHardPtCut"))
{

    if(needSeedingLayerSetsHits)
      theSeedingLayerToken = iC.consumes<SeedingLayerSetsHits>(cfg.getParameter<edm::InputTag>("SeedingLayers"));

    if (cfg.exists("SeedComparitorPSet"))
    {
        edm::ParameterSet comparitorPSet = cfg.getParameter < edm::ParameterSet
                > ("SeedComparitorPSet");
        std::string comparitorName = comparitorPSet.getParameter < std::string
                > ("ComponentName");
        if (comparitorName != "none")
        {
            theComparitor.reset(
                    SeedComparitorFactory::get()->create(comparitorName,
                            comparitorPSet, iC));
        }
    }
}

CAHitTripletGenerator::~CAHitTripletGenerator ( )

virtual

Definition at line 65 of file CAHitTripletGenerator.cc.

{
}

Member Function Documentation

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

static

Definition at line 69 of file CAHitTripletGenerator.cc.

References edm::ParameterSetDescription::add(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                             {
  desc.add<double>("extraHitRPhitolerance", 0.06);
  desc.add<bool>("useBendingCorrection", false);
  desc.add<double>("CAThetaCut", 0.00125);
  desc.add<double>("CAPhiCut", 0.1);
  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);
}

static const char* CAHitTripletGenerator::fillDescriptionsLabel ( )

inlinestatic

Definition at line 46 of file CAHitTripletGenerator.h.

{ return "caHitTriplet"; }

void CAHitTripletGenerator::hitNtuplets	(	const IntermediateHitDoublets::RegionLayerSets &	regionLayerPairs,
		OrderedHitTriplets &	result,
		const edm::EventSetup &	es,
		const SeedingLayerSetsHits &	layers
	)

Definition at line 196 of file CAHitTripletGenerator.cc.

References ihd::RegionLayerSets< T >::begin(), ihd::RegionLayerSets< T >::end(), newFWLiteAna::found, g, hitTriplets(), SeedingLayerSetsHits::SeedingLayer::index(), SurfaceOrientation::inner, IntermediateHitDoublets::LayerPairHitDoublets::innerLayerIndex(), SurfaceOrientation::outer, IntermediateHitDoublets::LayerPairHitDoublets::outerLayerIndex(), ihd::RegionLayerSets< T >::region(), and mps_fire::result.

                                                                            {
  CAGraph g;

  std::vector<const HitDoublets *> hitDoublets;

  auto layerPairEqual = [](const IntermediateHitDoublets::LayerPairHitDoublets& pair,
                           SeedingLayerSetsHits::LayerIndex inner,
                           SeedingLayerSetsHits::LayerIndex outer) {
    return pair.innerLayerIndex() == inner && pair.outerLayerIndex() == outer;
  };
  fillGraph(layers, g, hitDoublets, [&](const SeedingLayerSetsHits::SeedingLayer& inner, const SeedingLayerSetsHits::SeedingLayer& outer, std::vector<const HitDoublets *>& hitDoublets) {
      using namespace std::placeholders;
      auto found = std::find_if(regionLayerPairs.begin(), regionLayerPairs.end(), std::bind(layerPairEqual, _1, inner.index(), outer.index()));
      if(found != regionLayerPairs.end()) {
        hitDoublets.emplace_back(&(found->doublets()));
        return true;
      }
      return false;
    });

  hitTriplets(regionLayerPairs.region(), result, hitDoublets, g, es);
}

void CAHitTripletGenerator::hitTriplets ( const TrackingRegion &  reg, OrderedHitTriplets &  triplets, const edm::Event &  ev, const edm::EventSetup &  es  )

virtual

from base class

Implements HitTripletGenerator.

Definition at line 159 of file CAHitTripletGenerator.cc.

References LayerHitMapCache::clear(), HitPairGeneratorFromLayerPair::doublets(), alignCSCRings::e, g, edm::Event::getByToken(), SurfaceOrientation::inner, SeedingLayerSetsHits::numberOfLayersInSet(), SurfaceOrientation::outer, theComparitor, theLayerCache, and theSeedingLayerToken.

Referenced by hitNtuplets().

{
    edm::Handle<SeedingLayerSetsHits> hlayers;
    ev.getByToken(theSeedingLayerToken, hlayers);
    const SeedingLayerSetsHits& layers = *hlayers;
    if (layers.numberOfLayersInSet() != 3)
        throw cms::Exception("Configuration")
                << "CAHitTripletGenerator expects SeedingLayerSetsHits::numberOfLayersInSet() to be 3, got "
                << layers.numberOfLayersInSet();

    CAGraph g;

    std::vector<HitDoublets> hitDoublets;



    HitPairGeneratorFromLayerPair thePairGenerator(0, 1, &theLayerCache);
        fillGraph(layers, g, hitDoublets, [&](const SeedingLayerSetsHits::SeedingLayer& inner, const SeedingLayerSetsHits::SeedingLayer& outer, std::vector<HitDoublets>& hitDoublets) {
            hitDoublets.emplace_back(thePairGenerator.doublets(region, ev, es, inner, outer));
            return true;
          });


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

        std::vector<const HitDoublets *> hitDoubletsPtr;
        hitDoubletsPtr.reserve(hitDoublets.size());
        for(const auto& e: hitDoublets)
          hitDoubletsPtr.emplace_back(&e);

        hitTriplets(region, result, hitDoubletsPtr, g, es);
        theLayerCache.clear();
}

void CAHitTripletGenerator::hitTriplets ( const TrackingRegion &  reg, OrderedHitTriplets &  result, std::vector< const HitDoublets * > &  hitDoublets, const CAGraph &  g, const edm::EventSetup &  es  )

private

Definition at line 223 of file CAHitTripletGenerator.cc.

References funct::abs(), caHardPtCut, caPhiCut, caThetaCut, RZLine::chi2(), beam_dqm_sourceclient-live_cfg::chi2, LayerHitMapCache::clear(), ThirdHitPredictionFromCircle::curvature(), PixelRecoUtilities::curvature(), GlobalErrorBase< T, ErrorWeightType >::czz(), relativeConstraints::error, CAHitTripletGenerator::QuantityDependsPt::evaluator(), extraHitRPhitolerance, CellularAutomaton::findTriplets(), i, PixelRecoUtilities::inversePt(), GeomDetEnumerators::isBarrel(), edm::isNotFinite(), geometryCSVtoXML::line, maxChi2, TrackingRegion::origin(), PixelSubdetector::PixelBarrel, point, EnergyCorrector::pt, GlobalErrorBase< T, ErrorWeightType >::rerr(), funct::sqr(), mathSSE::sqrt(), theComparitor, theLayerCache, useBendingCorrection, CAHitTripletGenerator::QuantityDependsPtEval::value(), x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                 {
    std::vector<CACell::CAntuplet> foundTriplets;
    CellularAutomaton ca(g);

    ca.findTriplets(hitDoublets, foundTriplets, region, caThetaCut, caPhiCut,
            caHardPtCut);

    unsigned int numberOfFoundTriplets = foundTriplets.size();

    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;

    for (unsigned int tripletId = 0; tripletId < numberOfFoundTriplets;
            ++tripletId)
    {

        OrderedHitTriplet tmpTriplet(foundTriplets[tripletId][0]->getInnerHit(),
                foundTriplets[tripletId][0]->getOuterHit(),
                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 = foundTriplets[tripletId][i]->getInnerHit();
            gps[i] = ahit->globalPosition();
            ges[i] = ahit->globalPositionError();
            barrels[i] = isBarrel(ahit->geographicalId().subdetId());
        }

        auto const& ahit = 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.push_back(tmpTriplet);

    }
    theLayerCache.clear();
}

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

Definition at line 90 of file CAHitTripletGenerator.cc.

References theComparitor.

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

Member Data Documentation

const float CAHitTripletGenerator::caHardPtCut = 0.f

private

Definition at line 144 of file CAHitTripletGenerator.h.

Referenced by hitTriplets().

const float CAHitTripletGenerator::caPhiCut = 1.f

private

Definition at line 143 of file CAHitTripletGenerator.h.

Referenced by hitTriplets().

const float CAHitTripletGenerator::caThetaCut = 0.00125f

private

Definition at line 142 of file CAHitTripletGenerator.h.

Referenced by hitTriplets().

const float CAHitTripletGenerator::extraHitRPhitolerance

private

Definition at line 137 of file CAHitTripletGenerator.h.

Referenced by hitTriplets().

const QuantityDependsPt CAHitTripletGenerator::maxChi2

private

Definition at line 139 of file CAHitTripletGenerator.h.

Referenced by hitTriplets().

unsigned int CAHitTripletGenerator::minLayers = 3

static

Definition at line 35 of file CAHitTripletGenerator.h.

std::unique_ptr<SeedComparitor> CAHitTripletGenerator::theComparitor

private

Definition at line 71 of file CAHitTripletGenerator.h.

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

LayerCacheType CAHitTripletGenerator::theLayerCache

private

Definition at line 69 of file CAHitTripletGenerator.h.

Referenced by hitTriplets().

edm::EDGetTokenT<SeedingLayerSetsHits> CAHitTripletGenerator::theSeedingLayerToken

private

Definition at line 67 of file CAHitTripletGenerator.h.

Referenced by CAHitTripletGenerator(), and hitTriplets().

const bool CAHitTripletGenerator::useBendingCorrection

private

Definition at line 140 of file CAHitTripletGenerator.h.

Referenced by hitTriplets().