Inheritance diagram for HGCalTrackCollectionProducer:

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

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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_

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

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
	)

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 129 of file HGCalTrackCollectionProducer.cc.

References funct::abs(), bField_, debug_, DEFINE_FWK_MODULE, diskInnerRadius_, diskOuterRadius_, SoftLeptonByDistance_cfi::distance, DPtovPtCut_, runEdmFileComparison::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(), EnergyCorrector::pt, edm::Event::put(), mathSSE::sqrt(), src_, tkGeom_, HiIsolationCommonParameters_cff::track, l1t::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");
 }