#include <TrackerCleaner.h>

Inheritance diagram for TrackerCleaner< T >:

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

	~TrackerCleaner () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Types
typedef edmNew::DetSetVector< T >	TrackClusterCollection

Private Member Functions
template<>
bool	match_rechit_type (const TrackingRecHit &murechit)

template<>
bool	match_rechit_type (const TrackingRecHit &murechit)

bool	match_rechit_type (const TrackingRecHit &murechit)

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

Private Attributes
std::map< std::string, edm::EDGetTokenT< TrackClusterCollection > >	inputs_

const edm::EDGetTokenT< edm::View< pat::Muon > >	mu_input_

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

template<typename T>
class TrackerCleaner< T >

Author: Stefan Wayand; Christian Veelken, LLR

Definition at line 40 of file TrackerCleaner.h.

Member Typedef Documentation

template<typename T >

typedef edmNew::DetSetVector<T> TrackerCleaner< T >::TrackClusterCollection

private

Definition at line 50 of file TrackerCleaner.h.

Constructor & Destructor Documentation

template<typename T >

TrackerCleaner< T >::TrackerCleaner ( const edm::ParameterSet & iConfig )

explicit

Definition at line 60 of file TrackerCleaner.h.

References edm::ParameterSet::getParameter(), and TrackerCleaner< T >::inputs_.

                                                                 :
     mu_input_(consumes<edm::View<pat::Muon> >(iConfig.getParameter<edm::InputTag>("MuonCollection")))
   
 {
   std::vector<edm::InputTag> inCollections =  iConfig.getParameter<std::vector<edm::InputTag> >("oldCollection");
   for (auto inCollection : inCollections){
      inputs_[inCollection.instance()] = consumes<TrackClusterCollection >(inCollection);
      produces<TrackClusterCollection>(inCollection.instance());
   }
     
   
 }

template<typename T >

TrackerCleaner< T >::~TrackerCleaner ( )

override

Definition at line 75 of file TrackerCleaner.h.

 {
 // nothing to be done yet...  
 }

Member Function Documentation

template<>

bool TrackerCleaner< SiStripCluster >::match_rechit_type ( const TrackingRecHit & murechit )

private

Definition at line 34 of file TrackerCleaner.cc.

                                                                                     {
   
      const std::type_info &hit_type = typeid(murechit);
      if (hit_type == typeid(SiStripRecHit2D)) return true;
      else if (hit_type == typeid(SiStripRecHit1D)) return true;
      else if (hit_type == typeid(SiStripMatchedRecHit2D)) return true;
      else if (hit_type == typeid(ProjectedSiStripRecHit2D))  return true;
 
      
      
      return false;
    
  }

template<>

bool TrackerCleaner< SiPixelCluster >::match_rechit_type ( const TrackingRecHit & murechit )

private

Definition at line 49 of file TrackerCleaner.cc.

References DEFINE_FWK_MODULE.

                                                                                     {
    
       const std::type_info &hit_type = typeid(murechit);
      if (hit_type == typeid(SiPixelRecHit)) return true;
 
      return false;
    
  }

template<typename T >

bool TrackerCleaner< T >::match_rechit_type ( const TrackingRecHit & murechit )

private

Definition at line 24 of file TrackerCleaner.cc.

Referenced by TrackerCleaner< T >::produce().

                                                                        {
    
      assert(0); // CV: make sure general function never gets called;
              //     always use template specializations
      return false;
    
  }

template<typename T >

void TrackerCleaner< T >::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 82 of file TrackerCleaner.h.

References edm::View< T >::begin(), edm::View< T >::end(), BaseTrackerRecHit::firstClusterRef(), edm::Event::getByToken(), training_settings::idx, TrackerCleaner< T >::inputs_, TrackerCleaner< T >::match_rechit_type(), eostools::move(), TrackerCleaner< T >::mu_input_, extraflags_cff::muons, convertSQLitetoXML_cfg::output, edm::Event::put(), reco::Track::recHitsBegin(), and reco::Track::recHitsEnd().

 {
   
    using namespace edm;
    
    edm::Handle< edm::View<pat::Muon> > muonHandle;
    iEvent.getByToken(mu_input_, muonHandle);
    edm::View<pat::Muon> muons = *muonHandle;
    
   for (auto input_ : inputs_){
    
    edm::Handle<TrackClusterCollection > inputClusters;
    iEvent.getByToken(input_.second, inputClusters);
 
    std::vector<bool> vetodClusters;
 
    vetodClusters.resize(inputClusters->dataSize(), false);
 
    for (edm::View<pat::Muon>::const_iterator iMuon = muons.begin(); iMuon != muons.end(); ++iMuon) {   
     if(!iMuon->isGlobalMuon() ) continue;
     const reco::Track* mutrack = iMuon->globalTrack().get();
   //  reco::Track *mutrack = new reco::Track(*(iMuon->innerTrack() ));
     for (trackingRecHit_iterator hitIt = mutrack->recHitsBegin(); hitIt != mutrack->recHitsEnd(); ++hitIt) {
         const TrackingRecHit &murechit = **hitIt;     
         if(!(murechit).isValid()) continue;
     
     if (match_rechit_type(murechit)){
       auto & thit = reinterpret_cast<BaseTrackerRecHit const&>(murechit);
           auto const & cluster = thit.firstClusterRef();
       vetodClusters[cluster.key()]=true;
     }
     
         
      }  
   }
   std::unique_ptr<TrackClusterCollection > output(new TrackClusterCollection());
 
     int idx = 0;
     for ( typename TrackClusterCollection::const_iterator clustSet = inputClusters->begin(); clustSet != inputClusters->end(); ++clustSet ) { 
       DetId detIdObject( clustSet->detId() );
       typename TrackClusterCollection::FastFiller spc(*output, detIdObject);
       for (typename edmNew::DetSet<T>::const_iterator clustIt = clustSet->begin(); clustIt != clustSet->end(); ++clustIt ) { 
         idx++;  
         if (vetodClusters[idx-1]) continue;
     //if (!vetodClusters[idx-1]) continue; for inverted selction
         spc.push_back(*clustIt);
       }
     }
     iEvent.put(std::move(output),input_.first);
   }
   
 }