Inheritance diagram for HGCalTrackCollectionProducer:

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

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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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_

bool	debug_

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 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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 44 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"))),
   debug_(iConfig.getParameter<bool>("debug")),
   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
	)

overrideprivatevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 95 of file HGCalTrackCollectionProducer.cc.

References funct::abs(), alongMomentum, bField_, constexpr, diskInnerRadius_, diskOuterRadius_, edm::EventSetup::get(), HGCalDDDConstants::getTrForms(), hgc_names_, hgcGeometries_, i, LogDebug, mat_prop_, minusSurface_, plusSurface_, edm::ESHandle< class >::product(), makeMuonMisalignmentScenario::rot, tkGeom_, and Gflash::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
	)

overrideprivatevirtual

Implements edm::stream::EDProducerBase.

Definition at line 129 of file HGCalTrackCollectionProducer.cc.

References funct::abs(), bField_, debug_, diskInnerRadius_, diskOuterRadius_, HLT_FULL_cff::distance, DPtovPtCut_, newFWLiteAna::found, edm::Event::getByToken(), TrajectoryStateOnSurface::globalPosition(), PFTrackAlgoTools::goodPtResolution(), i, TrajectoryStateOnSurface::isValid(), LogDebug, mat_prop_, minusSurface_, eostools::move(), NHitCut_, trajectoryStateTransform::outerStateOnSurface(), PV3DBase< T, PVType, FrameType >::perp(), plusSurface_, edm::ESHandle< class >::product(), EnergyCorrector::pt, edm::Event::put(), mathSSE::sqrt(), src_, tkGeom_, testEve_cfg::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_, debug_);
     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");
 }