Inheritance diagram for HGCalTrackCollectionProducer:

Public Member Functions
	HGCalTrackCollectionProducer (const edm::ParameterSet &)

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
void	beginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &) override

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::ESHandle< MagneticField >	bField_

float	diskInnerRadius_

float	diskOuterRadius_

const std::vector< double >	DPtovPtCut_

std::array< std::string, 1 >	hgc_names_

std::array< edm::ESHandle< HGCalGeometry >, 1 >	hgcGeometries_

std::unique_ptr< PropagatorWithMaterial >	mat_prop_

std::array< std::vector< ReferenceCountingPointer< BoundDisk > >, 1 >	minusSurface_

const std::vector< unsigned >	NHitCut_

std::array< std::vector< ReferenceCountingPointer< BoundDisk > >, 1 >	plusSurface_

edm::EDGetTokenT< edm::View< reco::PFRecTrack > >	src_

edm::ESHandle< TrackerGeometry >	tkGeom_

const reco::TrackBase::TrackQuality	trackQuality_

const bool	useIterTracking_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 43 of file HGCalTrackCollectionProducer.cc.

Constructor & Destructor Documentation

HGCalTrackCollectionProducer::HGCalTrackCollectionProducer ( const edm::ParameterSet & iConfig )

Definition at line 75 of file HGCalTrackCollectionProducer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), hgc_names_, LogDebug, and AlCaHLTBitMon_QueryRunRegistry::string.

     : src_(consumes<edm::View<reco::PFRecTrack> >(iConfig.getParameter<edm::InputTag>("src"))),
       trackQuality_((iConfig.existsAs<std::string>("trackQuality"))
                         ? reco::TrackBase::qualityByName(iConfig.getParameter<std::string>("trackQuality"))
                         : reco::TrackBase::highPurity),
       DPtovPtCut_(iConfig.getParameter<std::vector<double> >("DPtOverPtCuts_byTrackAlgo")),
       NHitCut_(iConfig.getParameter<std::vector<unsigned> >("NHitCuts_byTrackAlgo")),
       useIterTracking_(iConfig.getParameter<bool>("useIterativeTracking")) {
   LogDebug("HGCalTrackCollectionProducer")
       << " HGCalTrackCollectionProducer::HGCalTrackCollectionProducer " << std::endl;
 
   const edm::ParameterSet& geoconf = iConfig.getParameterSet("hgcalGeometryNames");
   hgc_names_[0] = geoconf.getParameter<std::string>("HGC_ECAL");
 
   produces<reco::PFRecTrackCollection>("TracksInHGCal");
   produces<reco::PFRecTrackCollection>("TracksNotInHGCal");
 }

Member Function Documentation

void HGCalTrackCollectionProducer::beginLuminosityBlock	(	const edm::LuminosityBlock &	lumi,
		const edm::EventSetup &	es
	)

overrideprivate

Definition at line 95 of file HGCalTrackCollectionProducer.cc.

References funct::abs(), alongMomentum, bField_, constexpr, diskInnerRadius_, diskOuterRadius_, edm::EventSetup::get(), HGCalDDDConstants::getTrForms(), hgc_names_, hgcGeometries_, mps_fire::i, LogDebug, mat_prop_, minusSurface_, plusSurface_, edm::ESHandle< T >::product(), makeMuonMisalignmentScenario::rot, tkGeom_, and DOFs::Z.

                                                                                                                {
   constexpr float m_pion = 0.1396;
   // get dependencies for setting up propagator
   es.get<IdealMagneticFieldRecord>().get(bField_);
   es.get<TrackerDigiGeometryRecord>().get(tkGeom_);
   // get HGC geometries (assume that layers are ordered in Z!)
   for (unsigned i = 0; i < hgcGeometries_.size(); ++i) {
     es.get<IdealGeometryRecord>().get(hgc_names_[i], hgcGeometries_[i]);
   }
 
   // make propagator
   mat_prop_.reset(new PropagatorWithMaterial(alongMomentum, m_pion, bField_.product()));
   // setup HGC layers for track propagation
   Surface::RotationType rot;  //unit rotation matrix
   for (unsigned i = 0; i < hgcGeometries_.size(); ++i) {
     minusSurface_[i].clear();
     plusSurface_[i].clear();
     const HGCalDDDConstants& dddCons = hgcGeometries_[i]->topology().dddConstants();
     std::map<float, float> zrhoCoord;
     std::map<float, float> innerRadiusCoord;
     auto theTrForms = dddCons.getTrForms();
     const auto& firstLayerIt = theTrForms.back();
     float Z(std::abs(firstLayerIt.h3v.z()));
     // use hardcoded radii for now (FIX ME)
     diskInnerRadius_ = 31.5;
     diskOuterRadius_ = 161.0f;
     LogDebug("HGCalTrackCollectionProducer") << "O HAI I'm making a bound disk with Outer R=" << diskOuterRadius_
                                              << " Inner R=" << diskInnerRadius_ << " and Z=" << Z << std::endl;
     minusSurface_[i].push_back(ReferenceCountingPointer<BoundDisk>(new BoundDisk(
         Surface::PositionType(0, 0, -Z), rot, new SimpleDiskBounds(diskInnerRadius_, diskOuterRadius_, -0.001, 0.001))));
     plusSurface_[i].push_back(ReferenceCountingPointer<BoundDisk>(new BoundDisk(
         Surface::PositionType(0, 0, +Z), rot, new SimpleDiskBounds(diskInnerRadius_, diskOuterRadius_, -0.001, 0.001))));
   }
 }

void HGCalTrackCollectionProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 130 of file HGCalTrackCollectionProducer.cc.

References funct::abs(), bField_, DEFINE_FWK_MODULE, diskInnerRadius_, diskOuterRadius_, HLT_2018_cff::distance, DPtovPtCut_, newFWLiteAna::found, edm::Event::getByToken(), TrajectoryStateOnSurface::globalPosition(), PFTrackAlgoTools::goodPtResolution(), mps_fire::i, TrajectoryStateOnSurface::isValid(), LogDebug, mat_prop_, minusSurface_, eostools::move(), NHitCut_, trajectoryStateTransform::outerStateOnSurface(), PV3DBase< T, PVType, FrameType >::perp(), plusSurface_, edm::ESHandle< T >::product(), DiDispStaMuonMonitor_cfi::pt, edm::Event::put(), mathSSE::sqrt(), src_, tkGeom_, HLT_2018_cff::track, trackQuality_, PDWG_EXOHSCP_cff::tracks, useIterTracking_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                      {
   edm::Handle<edm::View<reco::PFRecTrack> > trackHandle;
   evt.getByToken(src_, trackHandle);
   const auto& tracks = *trackHandle;
 
   auto outputInHGCal = std::make_unique<reco::PFRecTrackCollection>();
   auto outputNotInHGCal = std::make_unique<reco::PFRecTrackCollection>();
 
   for (unsigned int i = 0; i < tracks.size(); i++) {
     const auto track = tracks.ptrAt(i);
     bool isGood =
         PFTrackAlgoTools::goodPtResolution(track->trackRef(), DPtovPtCut_, NHitCut_, useIterTracking_, trackQuality_);
     LogDebug("HGCalTrackCollectionProducer") << "HGCalTrackCollectionProducer Track number " << i
                                              << " has a goodPtResolution result of " << isGood << std::endl;
     if (!isGood)
       continue;
     bool found = false;
     const TrajectoryStateOnSurface myTSOS =
         trajectoryStateTransform::outerStateOnSurface(*(track->trackRef()), *(tkGeom_.product()), bField_.product());
     auto detbegin = myTSOS.globalPosition().z() > 0 ? plusSurface_.begin() : minusSurface_.begin();
     auto detend = myTSOS.globalPosition().z() > 0 ? plusSurface_.end() : minusSurface_.end();
     for (auto det = detbegin; det != detend; ++det) {
       LogDebug("HGCalTrackCollectionProducer") << "at HGC detector: " << std::distance(detbegin, det) << std::endl;
       unsigned layer_count = 1;
       for (const auto& layer : *det) {
         LogDebug("HGCalTrackCollectionProducer") << "  at DET layer: " << layer_count++ << std::endl;
         TrajectoryStateOnSurface piStateAtSurface = mat_prop_->propagate(myTSOS, *layer);
         if (piStateAtSurface.isValid()) {
           LogDebug("HGCalTrackCollectionProducer") << "Extrapolation is valid!" << std::endl;
           GlobalPoint pt = piStateAtSurface.globalPosition();
           if (pt.perp() < diskOuterRadius_) {
             if (pt.perp() > diskInnerRadius_) {
               LogDebug("HGCalTrackCollectionProducer")
                   << "(x,y,z,r)=(" << pt.x() << ", " << pt.y() << ", " << pt.z() << ", "
                   << sqrt(pt.x() * pt.x() + pt.y() * pt.y()) << ")" << std::endl;
               if (std::abs(track->trackRef()->eta()) < 1.47)
                 LogDebug("HGCalTrackCollectionProducer") << " ETA IN BARREL REGION: " << track->trackRef()->eta()
                                                          << " (PT: " << track->trackRef()->pt() << ")" << std::endl;
               found = true;
             } else {
               LogDebug("HGCalTrackCollectionProducer")
                   << " but r=" << pt.perp() << " < diskInnerRadius=" << diskInnerRadius_ << " so skipping "
                   << std::endl;
             }
           } else {
             LogDebug("HGCalTrackCollectionProducer")
                 << " but r=" << pt.perp() << " > diskOuterRadius=" << diskOuterRadius_ << " so skipping " << std::endl;
           }
         } else {
           LogDebug("HGCalTrackCollectionProducer") << "Extrapolation is NOT valid!" << std::endl;
         }
       }
     }
     if (found) {
       LogDebug("HGCalTrackCollectionProducer") << " Track going to outputInHGCal pt eta " << track->trackRef()->pt()
                                                << " " << track->trackRef()->eta() << std::endl;
       outputInHGCal->push_back(*track);
     } else {
       outputNotInHGCal->push_back(*track);
     }
   }  // track loop
 
   evt.put(std::move(outputInHGCal), "TracksInHGCal");
   evt.put(std::move(outputNotInHGCal), "TracksNotInHGCal");
 }