#include <MuonDetRecHitMixer.h>

Inheritance diagram for MuonDetRecHitMixer< T1, T2 >:

Classes
struct	TodoListEntry

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

	~MuonDetRecHitMixer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

	EDConsumerBase ()

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Types
typedef std::map< uint32_t, int >	detIdToIntMap

typedef edm::RangeMap< T1, edm::OwnVector< T2 > >	RecHitCollection

Private Member Functions
void	addRecHits (std::map< T1, std::vector< T2 > > &, const RecHitCollection &, bool, const detIdToIntMap &, const detIdToIntMap &, int &)

template<>
uint32_t	getRawDetId (const CSCRecHit2D &recHit)

template<>
uint32_t	getRawDetId (const DTRecHit1DPair &recHit)

template<>
uint32_t	getRawDetId (const RPCRecHit &recHit)

uint32_t	getRawDetId (const T2 &)

void	printHitMapRH (const edm::EventSetup &, const RecHitCollection &)

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

Private Attributes
std::string	moduleLabel_

edm::InputTag	srcHitMapMuMinus_

edm::InputTag	srcHitMapMuPlus_

std::vector< TodoListEntry >	todoList_

int	verbosity_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	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::EDProducer
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

template<typename T1, typename T2>
class MuonDetRecHitMixer< T1, T2 >

Merge collections of hits in muon detectors (CSC, DT and RPC) for original Zmumu event and "embedded" simulated tau decay products

Author: Christian Veelken, LLR

Version

Revision:: 1.2

Id:: MuonDetRecHitMixer.h,v 1.2 2012/12/13 09:52:06 veelken Exp

Definition at line 30 of file MuonDetRecHitMixer.h.

Member Typedef Documentation

template<typename T1 , typename T2 >

typedef std::map<uint32_t, int> MuonDetRecHitMixer< T1, T2 >::detIdToIntMap

private

Definition at line 41 of file MuonDetRecHitMixer.h.

template<typename T1 , typename T2 >

typedef edm::RangeMap<T1, edm::OwnVector<T2> > MuonDetRecHitMixer< T1, T2 >::RecHitCollection

private

Definition at line 39 of file MuonDetRecHitMixer.h.

Constructor & Destructor Documentation

template<typename T1 , typename T2 >

MuonDetRecHitMixer< T1, T2 >::MuonDetRecHitMixer ( const edm::ParameterSet & cfg )

explicit

Definition at line 67 of file MuonDetRecHitMixer.h.

References MuonDetRecHitMixer< T1, T2 >::TodoListEntry::cleanCollection1_, MuonDetRecHitMixer< T1, T2 >::TodoListEntry::cleanCollection2_, Exception, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::InputTag::instance(), MuonDetRecHitMixer< T1, T2 >::TodoListEntry::srcCollection1_, MuonDetRecHitMixer< T1, T2 >::TodoListEntry::srcCollection2_, AlCaHLTBitMon_QueryRunRegistry::string, MuonDetRecHitMixer< T1, T2 >::todoList_, and MuonDetRecHitMixer< T1, T2 >::verbosity_.

   : moduleLabel_(cfg.getParameter<std::string>("@module_label")),        
     srcHitMapMuPlus_(cfg.getParameter<edm::InputTag>("srcHitMapMuPlus")),
     srcHitMapMuMinus_(cfg.getParameter<edm::InputTag>("srcHitMapMuMinus"))
 {
   edm::VParameterSet cfgTodoList = cfg.getParameter<edm::VParameterSet>("todo");
   for ( edm::VParameterSet::const_iterator cfgTodoItem = cfgTodoList.begin();
     cfgTodoItem != cfgTodoList.end(); ++cfgTodoItem ) {
     TodoListEntry todoItem;
     todoItem.srcCollection1_ = cfgTodoItem->getParameter<edm::InputTag>("collection1");
     todoItem.cleanCollection1_ = cfgTodoItem->getParameter<bool>("cleanCollection1");
     todoItem.srcCollection2_ = cfgTodoItem->getParameter<edm::InputTag>("collection2");
     todoItem.cleanCollection2_ = cfgTodoItem->getParameter<bool>("cleanCollection2");
 
     std::string instanceLabel1 = todoItem.srcCollection1_.instance();
     std::string instanceLabel2 = todoItem.srcCollection2_.instance();
     if ( instanceLabel1 != instanceLabel2 ) {
       throw cms::Exception("Configuration") 
     << "Mismatch in Instance labels for collection 1 = " << instanceLabel1 << " and 2 = " << instanceLabel2 << " !!\n";
     }
     
     if ( instanceLabel1 == "" ) produces<RecHitCollection>();
     else produces<RecHitCollection>(instanceLabel1);
     
     todoList_.push_back(todoItem);
   }
 
   if ( todoList_.size() == 0 ) {
     throw cms::Exception("Configuration") 
       << "Empty to-do list !!\n";
   }
 
   verbosity_ = ( cfg.exists("verbosity") ) ?
     cfg.getParameter<int>("verbosity") : 0;
 }

template<typename T1 , typename T2 >

MuonDetRecHitMixer< T1, T2 >::~MuonDetRecHitMixer ( )

Definition at line 104 of file MuonDetRecHitMixer.h.

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

Member Function Documentation

template<typename T1 , typename T2 >

void MuonDetRecHitMixer< T1, T2 >::addRecHits	(	std::map< T1, std::vector< T2 > > &	recHits_output,
		const RecHitCollection &	recHitCollection_input,
		bool	applyCleaning,
		const detIdToIntMap &	hitMapMuPlus,
		const detIdToIntMap &	hitMapMuMinus,
		int &	numHits_cleaned
	)

private

Definition at line 165 of file MuonDetRecHitMixer.h.

References edm::RangeMap< ID, C, P >::begin(), and edm::RangeMap< ID, C, P >::end().

 {
   for ( typename RecHitCollection::const_iterator recHit = recHitCollection_input.begin();
     recHit != recHitCollection_input.end(); ++recHit ) {
     uint32_t rawDetId = getRawDetId(*recHit);
 
     bool isToBeCleaned = false;
     if ( applyCleaning ) {      
       isToBeCleaned |= matchHitMapRH(rawDetId, hitMapMuPlus);
       isToBeCleaned |= matchHitMapRH(rawDetId, hitMapMuMinus);
     }
       
     if ( !isToBeCleaned ) {
       T1 detId(rawDetId);
       recHits_output[detId].push_back(*recHit);
     } else {
       ++numHits_cleaned;
     }
   }
 }

template<>

uint32_t MuonDetRecHitMixer< CSCDetId, CSCRecHit2D >::getRawDetId ( const CSCRecHit2D & recHit )

private

Definition at line 27 of file MuonDetRecHitMixer.cc.

References CSCRecHit2D::cscDetId(), and DetId::rawId().

 {
   return recHit.cscDetId().rawId();
 }

template<>

uint32_t MuonDetRecHitMixer< DTLayerId, DTRecHit1DPair >::getRawDetId ( const DTRecHit1DPair & recHit )

private

Definition at line 33 of file MuonDetRecHitMixer.cc.

References DTRecHit1DPair::geographicalId(), and DetId::rawId().

 {
   return recHit.geographicalId().rawId();
 }

template<>

uint32_t MuonDetRecHitMixer< RPCDetId, RPCRecHit >::getRawDetId ( const RPCRecHit & recHit )

private

Definition at line 39 of file MuonDetRecHitMixer.cc.

References DetId::rawId(), and RPCRecHit::rpcId().

 {
   return recHit.rpcId().rawId();
 }

template<typename T1 , typename T2 >

uint32_t MuonDetRecHitMixer< T1, T2 >::getRawDetId ( const T2 & recHit )

private

Definition at line 19 of file MuonDetRecHitMixer.cc.

References assert().

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

template<typename T1 , typename T2 >

void MuonDetRecHitMixer< T1, T2 >::printHitMapRH	(	const edm::EventSetup &	es,
		const RecHitCollection &	recHits
	)

private

Definition at line 188 of file MuonDetRecHitMixer.h.

References edm::RangeMap< ID, C, P >::begin(), gather_cfg::cout, edm::RangeMap< ID, C, P >::end(), and printMuonDetId().

 {
   std::cout << "detIds:";
   for ( typename RecHitCollection::const_iterator rh = recHits.begin();
     rh != recHits.end(); ++rh ) {
     printMuonDetId(es, getRawDetId(*rh));
   }
 }

template<typename T1 , typename T2 >

void MuonDetRecHitMixer< T1, T2 >::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 110 of file MuonDetRecHitMixer.h.

References gather_cfg::cout, edm::Event::getByLabel(), python.multivaluedict::map(), moduleLabel_(), edm::Event::put(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   edm::Handle<detIdToIntMap> hitMapMuPlus;
   evt.getByLabel(srcHitMapMuPlus_, hitMapMuPlus);
   edm::Handle<detIdToIntMap> hitMapMuMinus;
   evt.getByLabel(srcHitMapMuMinus_, hitMapMuMinus);
 
   for ( typename std::vector<TodoListEntry>::const_iterator todoItem = todoList_.begin();
     todoItem != todoList_.end(); ++todoItem ) {
 
     typedef edm::Handle<RecHitCollection> RecHitCollectionHandle;
     RecHitCollectionHandle recHitCollection1;
     evt.getByLabel(todoItem->srcCollection1_, recHitCollection1);
     RecHitCollectionHandle recHitCollection2;
     evt.getByLabel(todoItem->srcCollection2_, recHitCollection2);
     
     std::map<T1, std::vector<T2> > recHits_output;
     
     int numHits_cleaned = 0;
     addRecHits(recHits_output, *recHitCollection1, todoItem->cleanCollection1_, *hitMapMuPlus, *hitMapMuMinus, numHits_cleaned);
     addRecHits(recHits_output, *recHitCollection2, todoItem->cleanCollection2_, *hitMapMuPlus, *hitMapMuMinus, numHits_cleaned);
 
     std::auto_ptr<RecHitCollection> recHitCollection_output(new RecHitCollection());
     for ( typename std::map<T1, std::vector<T2> >::const_iterator recHit = recHits_output.begin();
       recHit != recHits_output.end(); ++recHit ) {
       recHitCollection_output->put(recHit->first, recHit->second.begin(), recHit->second.end());
     }
     std::string instanceLabel = todoItem->srcCollection1_.instance(); 
     if ( verbosity_ ) {
       std::cout << "<MuonDetRecHitMixer::produce>:" << std::endl;   
       std::cout << " #Hits(input1 = " << todoItem->srcCollection1_.label() << ":" << todoItem->srcCollection1_.instance() << ") = " << recHitCollection1->size() << std::endl;
       if ( verbosity_ >= 2 ) printHitMapRH(es, *recHitCollection1);
       std::cout << " #Hits(input2 = " << todoItem->srcCollection2_.label() << ":" << todoItem->srcCollection2_.instance() << ") = " << recHitCollection2->size() << std::endl;
       if ( verbosity_ >= 2 ) printHitMapRH(es, *recHitCollection2);
       std::cout << " #Hits(output = " << moduleLabel_ << ":" << instanceLabel << ") = " << recHitCollection_output->size() << std::endl;
       if ( verbosity_ >= 2 ) printHitMapRH(es, *recHitCollection_output);
       std::cout << " #Hits = " << numHits_cleaned << " removed during cleaning." << std::endl;   
     }    
     evt.put(recHitCollection_output, instanceLabel);
   }
 }