#include <Calibration/TkAlCaRecoProducers/plugins/CalibrationTrackSelectorFromDetIdList.cc>

Inheritance diagram for CalibrationTrackSelectorFromDetIdList:

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

	~CalibrationTrackSelectorFromDetIdList () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Member Functions
void	beginRun (edm::Run const &run, const edm::EventSetup &) override

TrackCandidate	makeCandidate (const reco::Track &tk, std::vector< std::unique_ptr< TrackingRecHit > >::iterator hitsBegin, std::vector< std::unique_ptr< TrackingRecHit > >::iterator hitsEnd)

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

Private Attributes
std::vector< DetIdSelector >	detidsels_

edm::EDGetTokenT< reco::TrackCollection >	m_label

bool	m_verbose

edm::ESHandle< TrackerGeometry >	theGeometry

edm::ESHandle< MagneticField >	theMagField

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

Description: Selects tracks that have at leaast one valid hit on a given set of Tracker DetIds

Definition at line 40 of file CalibrationTrackSelectorFromDetIdList.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 59 of file CalibrationTrackSelectorFromDetIdList.cc.

References detidsels_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), m_label, m_verbose, and corrVsCorr::selection.

                                                                                                           :  detidsels_() {
 
   std::vector<edm::ParameterSet> selconfigs = iConfig.getParameter<std::vector<edm::ParameterSet> >("selections");
    
   for(std::vector<edm::ParameterSet>::const_iterator selconfig=selconfigs.begin();selconfig!=selconfigs.end();++selconfig) {
     DetIdSelector selection(*selconfig);
     detidsels_.push_back(selection);
   }
 
   m_verbose = iConfig.getUntrackedParameter<bool>("verbose");
   m_label = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("Input"));
   produces<TrackCandidateCollection>();
 }

CalibrationTrackSelectorFromDetIdList::~CalibrationTrackSelectorFromDetIdList ( )

override

Definition at line 74 of file CalibrationTrackSelectorFromDetIdList.cc.

74 {}

Member Function Documentation

void CalibrationTrackSelectorFromDetIdList::beginRun	(	edm::Run const &	run,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 150 of file CalibrationTrackSelectorFromDetIdList.cc.

References DEFINE_FWK_MODULE, TrackerGeometry::detIds(), detidsels_, edm::EventSetup::get(), LogDebug, m_verbose, edm::ESHandle< T >::product(), theGeometry, and theMagField.

                                                                                                    {
   iSetup.get<TrackerDigiGeometryRecord>().get(theGeometry);
   iSetup.get<IdealMagneticFieldRecord>().get(theMagField);
 
   if(m_verbose){
     for (const auto &detidsel : detidsels_){
       auto theDetIds = theGeometry.product()->detIds();
       for(const auto &theDet : theDetIds){
     if(detidsel.isSelected(theDet)) {
       LogDebug("CalibrationTrackSelectorFromDetIdList") << "detid: " << theDet.rawId() << " is taken" << std::endl;
     }
       } 
     }
   }
 
 }

TrackCandidate CalibrationTrackSelectorFromDetIdList::makeCandidate	(	const reco::Track &	tk,
		std::vector< std::unique_ptr< TrackingRecHit > >::iterator	hitsBegin,
		std::vector< std::unique_ptr< TrackingRecHit > >::iterator	hitsEnd
	)

private

Definition at line 123 of file CalibrationTrackSelectorFromDetIdList.cc.

References alongMomentum, anyDirection, Exception, reco::Track::innerDetId(), reco::Track::innerPosition(), trajectoryStateTransform::innerStateOnSurface(), reco::TrackBase::momentum(), reco::Track::outerDetId(), reco::Track::outerPosition(), trajectoryStateTransform::outerStateOnSurface(), trajectoryStateTransform::persistentState(), edm::OwnVector< T, P >::push_back(), edm::OwnVector< T, P >::reserve(), SurveyInfoScenario_cff::seed, reco::Track::seedDirection(), reco::Track::seedRef(), theGeometry, and theMagField.

Referenced by produce().

                                                                                                                                                                                                 {
   
   PropagationDirection   pdir = tk.seedDirection();
   PTrajectoryStateOnDet state;
   if ( pdir == anyDirection ) throw cms::Exception("UnimplementedFeature") << "Cannot work with tracks that have 'anyDirecton' \n";
 
   if ( (pdir == alongMomentum) == (  (tk.outerPosition()-tk.innerPosition()).Dot(tk.momentum()) >= 0 ) ) {
     // use inner state
     TrajectoryStateOnSurface originalTsosIn(trajectoryStateTransform::innerStateOnSurface(tk, *theGeometry, &*theMagField));
     state = trajectoryStateTransform::persistentState( originalTsosIn, DetId(tk.innerDetId()) );
   } else {
     // use outer state
     TrajectoryStateOnSurface originalTsosOut(trajectoryStateTransform::outerStateOnSurface(tk, *theGeometry, &*theMagField));
     state = trajectoryStateTransform::persistentState( originalTsosOut, DetId(tk.outerDetId()) );
   }
   TrajectorySeed seed(state, TrackCandidate::RecHitContainer(), pdir);
   TrackCandidate::RecHitContainer ownHits;
   ownHits.reserve(hitsEnd - hitsBegin);
   for ( ; hitsBegin != hitsEnd; ++hitsBegin) {
     ownHits.push_back( hitsBegin->release() );
   }
 
   TrackCandidate cand(ownHits, seed, state, tk.seedRef());
   
   return cand;
 }

void CalibrationTrackSelectorFromDetIdList::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 81 of file CalibrationTrackSelectorFromDetIdList.cc.

References TrackingRecHit::clone(), detidsels_, TrackingRecHit::geographicalId(), edm::Event::getByToken(), hfClusterShapes_cfi::hits, LogDebug, m_label, makeCandidate(), eostools::move(), convertSQLitetoXML_cfg::output, edm::Event::put(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), and l1t::tracks.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

                                                                                                  {
   using namespace edm;
   using namespace std;
 
   edm::Handle<std::vector<reco::Track> > tracks;
   iEvent.getByToken(m_label, tracks);
   auto output = std::make_unique<TrackCandidateCollection>();
   
   // loop on tracks
   for (std::vector<reco::Track>::const_iterator ittrk = tracks->begin(), edtrk = tracks->end(); ittrk != edtrk; ++ittrk) {
     const reco::Track *trk = &(*ittrk);
 
     std::vector<std::unique_ptr<TrackingRecHit> > hits; 
    
     bool saveTrack(false);
 
     for (trackingRecHit_iterator ith = trk->recHitsBegin(), edh = trk->recHitsEnd(); ith != edh; ++ith) {
       const TrackingRecHit * hit = (*ith); // ith is an iterator on edm::Ref to rechit
       DetId detid = hit->geographicalId();
       
       for (const auto &detidsel : detidsels_){
     if(detidsel.isSelected(detid)) {
       LogDebug("CalibrationTrackSelectorFromDetIdList") << "Selected by selection " << detid;
       saveTrack=true;
       break;
     }
       }
 
       // here there will be the selection
       hits.emplace_back(hit->clone());
 
     }
 
     if(saveTrack){
       output->push_back( makeCandidate ( *ittrk, hits.begin(), hits.end() ) );
     }
 
   }
   iEvent.put(std::move(output));
 }