#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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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<>
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

template<typename T>
class TrackerCleaner< T >

Author: Stefan Wayand; Christian Veelken, LLR

Definition at line 38 of file TrackerCleaner.h.

Member Typedef Documentation

template<typename T >

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

private

Definition at line 47 of file TrackerCleaner.h.

Constructor & Destructor Documentation

template<typename T >

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

explicit

Definition at line 55 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 (const auto& inCollection : inCollections) {
     inputs_[inCollection.instance()] = consumes<TrackClusterCollection>(inCollection);
     produces<TrackClusterCollection>(inCollection.instance());
   }
 }

template<typename T >

TrackerCleaner< T >::~TrackerCleaner ( )

override

Definition at line 67 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 26 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 41 of file TrackerCleaner.cc.

                                                                                      {
   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 19 of file TrackerCleaner.cc.

References cms::cuda::assert().

                                                                         {
   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 72 of file TrackerCleaner.h.

References runTheMatrix::const, edm::Event::getByToken(), eostools::move(), patZpeak::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);
   }
 }