#include <SeedingRegionByTracks.h>

Inheritance diagram for ticl::SeedingRegionByTracks:

Public Member Functions
void	initialize (const edm::EventSetup &es) override

void	makeRegions (const edm::Event &ev, const edm::EventSetup &es, std::vector< TICLSeedingRegion > &result) override

	SeedingRegionByTracks (const edm::ParameterSet &conf, edm::ConsumesCollector &sumes)

	~SeedingRegionByTracks () override

Public Member Functions inherited from ticl::SeedingRegionAlgoBase
	SeedingRegionAlgoBase (const edm::ParameterSet &conf, edm::ConsumesCollector &sumes)

virtual	~SeedingRegionAlgoBase ()

Static Public Member Functions
static void	fillPSetDescription (edm::ParameterSetDescription &desc)

Static Public Member Functions inherited from ticl::SeedingRegionAlgoBase
static void	fillPSetDescription (edm::ParameterSetDescription &desc)

Private Member Functions
void	buildFirstLayers ()

Private Attributes
edm::ESHandle< MagneticField >	bfield_

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	bfield_token_

const StringCutObjectSelector< reco::Track >	cutTk_

const std::string	detector_

std::unique_ptr< GeomDet >	firstDisk_ [2]

edm::ESGetToken< HGCalDDDConstants, IdealGeometryRecord >	hdc_token_

const HGCalDDDConstants *	hgcons_

std::once_flag	initializeGeometry_

edm::ESHandle< Propagator >	propagator_

edm::ESGetToken< Propagator, TrackingComponentsRecord >	propagator_token_

const std::string	propName_

edm::EDGetTokenT< reco::TrackCollection >	tracks_token_

Additional Inherited Members
Protected Attributes inherited from ticl::SeedingRegionAlgoBase
int	algo_verbosity_

int	algoId_

Detailed Description

Definition at line 31 of file SeedingRegionByTracks.h.

Constructor & Destructor Documentation

◆ SeedingRegionByTracks()

SeedingRegionByTracks::SeedingRegionByTracks	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &	sumes
	)

Definition at line 16 of file SeedingRegionByTracks.cc.

References edm::BeginRun, detector_, edm::ConsumesCollector::esConsumes(), hdc_token_, and AlCaHLTBitMon_QueryRunRegistry::string.

     : SeedingRegionAlgoBase(conf, sumes),
       tracks_token_(sumes.consumes<reco::TrackCollection>(conf.getParameter<edm::InputTag>("tracks"))),
       cutTk_(conf.getParameter<std::string>("cutTk")),
       detector_(conf.getParameter<std::string>("detector")),
       propName_(conf.getParameter<std::string>("propagator")),
       bfield_token_(sumes.esConsumes<MagneticField, IdealMagneticFieldRecord, edm::Transition::BeginRun>()),
       propagator_token_(sumes.esConsumes<Propagator, TrackingComponentsRecord, edm::Transition::BeginRun>(
           edm::ESInputTag("", propName_))) {
   std::string detectorName_ = (detector_ == "HFNose") ? "HGCalHFNoseSensitive" : "HGCalEESensitive";
   hdc_token_ = sumes.esConsumes<HGCalDDDConstants, IdealGeometryRecord, edm::Transition::BeginRun>(
       edm::ESInputTag("", detectorName_));
 }

◆ ~SeedingRegionByTracks()

SeedingRegionByTracks::~SeedingRegionByTracks ( )

override

Definition at line 30 of file SeedingRegionByTracks.cc.

30 {}

Member Function Documentation

◆ buildFirstLayers()

void SeedingRegionByTracks::buildFirstLayers ( )

private

Definition at line 81 of file SeedingRegionByTracks.cc.

References Disk::build(), firstDisk_, hgcons_, HGCalDDDConstants::rangeR(), and HGCalDDDConstants::waferZ().

Referenced by initialize().

                                              {
   float zVal = hgcons_->waferZ(1, true);
   std::pair<double, double> rMinMax = hgcons_->rangeR(zVal, true);
 
   for (int iSide = 0; iSide < 2; ++iSide) {
     float zSide = (iSide == 0) ? (-1. * zVal) : zVal;
     firstDisk_[iSide] =
         std::make_unique<GeomDet>(Disk::build(Disk::PositionType(0, 0, zSide),
                                               Disk::RotationType(),
                                               SimpleDiskBounds(rMinMax.first, rMinMax.second, zSide - 0.5, zSide + 0.5))
                                       .get());
   }
 }

◆ fillPSetDescription()

void SeedingRegionByTracks::fillPSetDescription ( edm::ParameterSetDescription & desc )

static

Definition at line 71 of file SeedingRegionByTracks.cc.

References submitPVResolutionJobs::desc, ticl::SeedingRegionAlgoBase::fillPSetDescription(), ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                 {
   desc.add<edm::InputTag>("tracks", edm::InputTag("generalTracks"));
   desc.add<std::string>("cutTk",
                         "1.48 < abs(eta) < 3.0 && pt > 1. && quality(\"highPurity\") && "
                         "hitPattern().numberOfLostHits(\"MISSING_OUTER_HITS\") < 5");
   desc.add<std::string>("propagator", "PropagatorWithMaterial");
   desc.add<std::string>("detector", "HGCAL");
   SeedingRegionAlgoBase::fillPSetDescription(desc);
 }

◆ initialize()

void SeedingRegionByTracks::initialize ( const edm::EventSetup & es )

overridevirtual

Implements ticl::SeedingRegionAlgoBase.

Definition at line 32 of file SeedingRegionByTracks.cc.

References bfield_, bfield_token_, buildFirstLayers(), edm::EventSetup::getHandle(), hdc_token_, hgcons_, edm::ESHandle< T >::product(), propagator_, and propagator_token_.

                                                               {
   edm::ESHandle<HGCalDDDConstants> hdc = es.getHandle(hdc_token_);
   hgcons_ = hdc.product();
 
   buildFirstLayers();
 
   bfield_ = es.getHandle(bfield_token_);
   propagator_ = es.getHandle(propagator_token_);
 }

◆ makeRegions()

void SeedingRegionByTracks::makeRegions	(	const edm::Event &	ev,
		const edm::EventSetup &	es,
		std::vector< TICLSeedingRegion > &	result
	)

overridevirtual

Implements ticl::SeedingRegionAlgoBase.

Definition at line 42 of file SeedingRegionByTracks.cc.

References a, b, bfield_, cutTk_, reco::TrackBase::eta(), makeMEIFBenchmarkPlots::ev, firstDisk_, mps_fire::i, edm::HandleBase::id(), createfilelist::int, beamSpotPI::nTracks, trajectoryStateTransform::outerFreeState(), edm::ESHandle< T >::product(), Propagator::propagate(), mps_fire::result, jetUpdater_cfi::sort, and tracks_token_.

                                                                               {
   edm::Handle<reco::TrackCollection> tracks_h;
   ev.getByToken(tracks_token_, tracks_h);
   edm::ProductID trkId = tracks_h.id();
   auto bFieldProd = bfield_.product();
   const Propagator &prop = (*propagator_);
 
   int nTracks = tracks_h->size();
   for (int i = 0; i < nTracks; ++i) {
     const reco::Track &tk = (*tracks_h)[i];
     if (!cutTk_((tk))) {
       continue;
     }
 
     FreeTrajectoryState fts = trajectoryStateTransform::outerFreeState((tk), bFieldProd);
     int iSide = int(tk.eta() > 0);
     TrajectoryStateOnSurface tsos = prop.propagate(fts, firstDisk_[iSide]->surface());
     if (tsos.isValid()) {
       result.emplace_back(tsos.globalPosition(), tsos.globalMomentum(), iSide, i, trkId);
     }
   }
   // sorting seeding region by descending momentum
   std::sort(result.begin(), result.end(), [](const TICLSeedingRegion &a, const TICLSeedingRegion &b) {
     return a.directionAtOrigin.perp2() > b.directionAtOrigin.perp2();
   });
 }

Member Data Documentation

◆ bfield_

edm::ESHandle<MagneticField> ticl::SeedingRegionByTracks::bfield_

private

Definition at line 52 of file SeedingRegionByTracks.h.

Referenced by initialize(), and makeRegions().

◆ bfield_token_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> ticl::SeedingRegionByTracks::bfield_token_

private

Definition at line 55 of file SeedingRegionByTracks.h.

Referenced by initialize().

◆ cutTk_

const StringCutObjectSelector<reco::Track> ticl::SeedingRegionByTracks::cutTk_

private

Definition at line 48 of file SeedingRegionByTracks.h.

Referenced by makeRegions().

◆ detector_

const std::string ticl::SeedingRegionByTracks::detector_

private

Definition at line 49 of file SeedingRegionByTracks.h.

Referenced by SeedingRegionByTracks().

◆ firstDisk_

std::unique_ptr<GeomDet> ticl::SeedingRegionByTracks::firstDisk_[2]

private

Definition at line 53 of file SeedingRegionByTracks.h.

Referenced by buildFirstLayers(), and makeRegions().

◆ hdc_token_

edm::ESGetToken<HGCalDDDConstants, IdealGeometryRecord> ticl::SeedingRegionByTracks::hdc_token_

private

Definition at line 54 of file SeedingRegionByTracks.h.

Referenced by initialize(), and SeedingRegionByTracks().

◆ hgcons_

const HGCalDDDConstants* ticl::SeedingRegionByTracks::hgcons_

private

Definition at line 47 of file SeedingRegionByTracks.h.

Referenced by buildFirstLayers(), and initialize().

◆ initializeGeometry_

std::once_flag ticl::SeedingRegionByTracks::initializeGeometry_

private

Definition at line 46 of file SeedingRegionByTracks.h.

◆ propagator_

edm::ESHandle<Propagator> ticl::SeedingRegionByTracks::propagator_

private

Definition at line 50 of file SeedingRegionByTracks.h.

Referenced by initialize().

◆ propagator_token_

edm::ESGetToken<Propagator, TrackingComponentsRecord> ticl::SeedingRegionByTracks::propagator_token_

private

Definition at line 56 of file SeedingRegionByTracks.h.

Referenced by initialize().

◆ propName_

const std::string ticl::SeedingRegionByTracks::propName_

private

Definition at line 51 of file SeedingRegionByTracks.h.

◆ tracks_token_

edm::EDGetTokenT<reco::TrackCollection> ticl::SeedingRegionByTracks::tracks_token_

private

Definition at line 45 of file SeedingRegionByTracks.h.

Referenced by makeRegions().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SeedingRegionByTracks()

◆ ~SeedingRegionByTracks()

Member Function Documentation

◆ buildFirstLayers()

◆ fillPSetDescription()

◆ initialize()

◆ makeRegions()

Member Data Documentation

◆ bfield_

◆ bfield_token_

◆ cutTk_

◆ detector_

◆ firstDisk_

◆ hdc_token_

◆ hgcons_

◆ initializeGeometry_

◆ propagator_

◆ propagator_token_

◆ propName_

◆ tracks_token_