#include <MuonAssociatorByHits.h>

Public Member Functions
reco::RecoToSimCollection	associateRecoToSim (edm::Handle< edm::View< reco::Track >> &tCH, edm::Handle< TrackingParticleCollection > &tPCH, const edm::Event event, const edm::EventSetup setup) const

reco::RecoToSimCollection	associateRecoToSim (const edm::RefToBaseVector< reco::Track > &, const edm::RefVector< TrackingParticleCollection > &, const edm::Event event=nullptr, const edm::EventSetup setup=nullptr) const
	Association Reco To Sim with Collections. More...

virtual reco::SimToRecoCollection	associateSimToReco (edm::Handle< edm::View< reco::Track >> &tCH, edm::Handle< TrackingParticleCollection > &tPCH, const edm::Event event, const edm::EventSetup setup) const

reco::SimToRecoCollection	associateSimToReco (const edm::RefToBaseVector< reco::Track > &, const edm::RefVector< TrackingParticleCollection > &, const edm::Event event=nullptr, const edm::EventSetup setup=nullptr) const
	Association Sim To Reco with Collections. More...

	MuonAssociatorByHits (const edm::ParameterSet &conf, edm::ConsumesCollector &&iC)

virtual	~MuonAssociatorByHits ()

Private Attributes
CSCHitAssociator::Config	cscHitAssociatorConfig_

std::unique_ptr< muonAssociatorByHitsDiagnostics::InputDumper >	diagnostics_

DTHitAssociator::Config	dtHitAssociatorConfig_

GEMHitAssociator::Config	gemHitAssociatorConfig_

MuonAssociatorByHitsHelper	helper_

RPCHitAssociator::Config	rpcHitAssociatorConfig_

TrackerHitAssociator::Config	trackerHitAssociatorConfig_

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	ttopoToken_

Detailed Description

Definition at line 31 of file MuonAssociatorByHits.h.

Constructor & Destructor Documentation

◆ MuonAssociatorByHits()

MuonAssociatorByHits::MuonAssociatorByHits	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &&	iC
	)

Definition at line 135 of file MuonAssociatorByHits.cc.

References diagnostics_, edm::ParameterSet::getUntrackedParameter(), and eostools::move().

     : helper_(conf),
       trackerHitAssociatorConfig_(conf, std::move(iC)),
       gemHitAssociatorConfig_(conf, iC),
       rpcHitAssociatorConfig_(conf, iC),
       cscHitAssociatorConfig_(conf, iC),
       dtHitAssociatorConfig_(conf, iC),
       ttopoToken_(iC.esConsumes()) {
   // hack for consumes
   if (conf.getUntrackedParameter<bool>("dumpInputCollections")) {
     diagnostics_ = std::make_unique<InputDumper>(conf, std::move(iC));
   }
 }

◆ ~MuonAssociatorByHits()

MuonAssociatorByHits::~MuonAssociatorByHits ( )

virtual

Definition at line 149 of file MuonAssociatorByHits.cc.

149 {}

Member Function Documentation

◆ associateRecoToSim() [1/2]

reco::RecoToSimCollection MuonAssociatorByHits::associateRecoToSim	(	edm::Handle< edm::View< reco::Track >> &	tCH,
		edm::Handle< TrackingParticleCollection > &	tPCH,
		const edm::Event *	event,
		const edm::EventSetup *	setup
	)		const

inline

Definition at line 38 of file MuonAssociatorByHits.h.

References edm::HandleBase::id(), dqmiolumiharvest::j, edm::RefToBaseVector< T >::push_back(), edm::RefVector< C, T, F >::push_back(), and singleTopDQM_cfi::setup.

Referenced by MuonAssociatorEDProducer::produce().

                                                                                  {
     edm::RefToBaseVector<reco::Track> tc;
     for (unsigned int j = 0; j < tCH->size(); j++)
       tc.push_back(tCH->refAt(j));
 
     edm::RefVector<TrackingParticleCollection> tpc(tPCH.id());
     for (unsigned int j = 0; j < tPCH->size(); j++)
       tpc.push_back(edm::Ref<TrackingParticleCollection>(tPCH, j));
 
     return associateRecoToSim(tc, tpc, event, setup);
   }

◆ associateRecoToSim() [2/2]

RecoToSimCollection MuonAssociatorByHits::associateRecoToSim	(	const edm::RefToBaseVector< reco::Track > &	tC,
		const edm::RefVector< TrackingParticleCollection > &	TPCollectionH,
		const edm::Event *	event = `nullptr`,
		const edm::EventSetup *	setup = `nullptr`
	)		const

Association Reco To Sim with Collections.

Definition at line 151 of file MuonAssociatorByHits.cc.

References MuonAssociatorByHitsHelper::associateRecoToSimIndices(), edm::RefToBaseVector< T >::begin(), cscHitAssociatorConfig_, diagnostics_, dtHitAssociatorConfig_, MillePedeFileConverter_cfg::e, edm::RefToBaseVector< T >::end(), gemHitAssociatorConfig_, helper_, mps_fire::resources, rpcHitAssociatorConfig_, singleTopDQM_cfi::setup, trackerHitAssociatorConfig_, and ttopoToken_.

                                       {
   RecoToSimCollection outputCollection(&e->productGetter());
 
   MuonAssociatorByHitsHelper::TrackHitsCollection tH;
   for (auto it = tC.begin(), ed = tC.end(); it != ed; ++it) {
     tH.push_back(std::make_pair((*it)->recHitsBegin(), (*it)->recHitsEnd()));
   }
 
   // Retrieve tracker topology from geometry
   const TrackerTopology *tTopo = &setup->getData(ttopoToken_);
 
   // Tracker hit association
   TrackerHitAssociator trackertruth(*e, trackerHitAssociatorConfig_);
   // CSC hit association
   CSCHitAssociator csctruth(*e, *setup, cscHitAssociatorConfig_);
   // DT hit association
   bool printRtS(true);
   DTHitAssociator dttruth(*e, *setup, dtHitAssociatorConfig_, printRtS);
   // RPC hit association
   RPCHitAssociator rpctruth(*e, rpcHitAssociatorConfig_);
   // GEM hit association
   GEMHitAssociator gemtruth(*e, gemHitAssociatorConfig_);
 
   MuonAssociatorByHitsHelper::Resources resources = {
       tTopo, &trackertruth, &csctruth, &dttruth, &rpctruth, &gemtruth, {}};
 
   if (diagnostics_) {
     resources.diagnostics_ = [this, e](const TrackHitsCollection &hC, const TrackingParticleCollection &pC) {
       diagnostics_->dump(hC, pC, *e);
     };
   }
 
   auto bareAssoc = helper_.associateRecoToSimIndices(tH, TPCollectionH, resources);
   for (auto it = bareAssoc.begin(), ed = bareAssoc.end(); it != ed; ++it) {
     for (auto itma = it->second.begin(), edma = it->second.end(); itma != edma; ++itma) {
       outputCollection.insert(tC[it->first], std::make_pair(TPCollectionH[itma->idx], itma->quality));
     }
   }
 
   outputCollection.post_insert();  // perhaps not even necessary
   return outputCollection;
 }

◆ associateSimToReco() [1/2]

virtual reco::SimToRecoCollection MuonAssociatorByHits::associateSimToReco	(	edm::Handle< edm::View< reco::Track >> &	tCH,
		edm::Handle< TrackingParticleCollection > &	tPCH,
		const edm::Event *	event,
		const edm::EventSetup *	setup
	)		const

inlinevirtual

Definition at line 53 of file MuonAssociatorByHits.h.

References edm::HandleBase::id(), dqmiolumiharvest::j, edm::RefToBaseVector< T >::push_back(), edm::RefVector< C, T, F >::push_back(), and singleTopDQM_cfi::setup.

Referenced by MuonAssociatorEDProducer::produce().

                                                                                          {
     edm::RefToBaseVector<reco::Track> tc;
     for (unsigned int j = 0; j < tCH->size(); j++)
       tc.push_back(tCH->refAt(j));
 
     edm::RefVector<TrackingParticleCollection> tpc(tPCH.id());
     for (unsigned int j = 0; j < tPCH->size(); j++)
       tpc.push_back(edm::Ref<TrackingParticleCollection>(tPCH, j));
 
     return associateSimToReco(tc, tpc, event, setup);
   }

◆ associateSimToReco() [2/2]

SimToRecoCollection MuonAssociatorByHits::associateSimToReco	(	const edm::RefToBaseVector< reco::Track > &	tC,
		const edm::RefVector< TrackingParticleCollection > &	TPCollectionH,
		const edm::Event *	event = `nullptr`,
		const edm::EventSetup *	setup = `nullptr`
	)		const

Association Sim To Reco with Collections.

Definition at line 198 of file MuonAssociatorByHits.cc.

References MuonAssociatorByHitsHelper::associateSimToRecoIndices(), edm::RefToBaseVector< T >::begin(), cscHitAssociatorConfig_, dtHitAssociatorConfig_, MillePedeFileConverter_cfg::e, edm::RefToBaseVector< T >::end(), gemHitAssociatorConfig_, helper_, mps_fire::resources, rpcHitAssociatorConfig_, singleTopDQM_cfi::setup, trackerHitAssociatorConfig_, and ttopoToken_.

                                       {
   SimToRecoCollection outputCollection(&e->productGetter());
   MuonAssociatorByHitsHelper::TrackHitsCollection tH;
   for (auto it = tC.begin(), ed = tC.end(); it != ed; ++it) {
     tH.push_back(std::make_pair((*it)->recHitsBegin(), (*it)->recHitsEnd()));
   }
 
   // Retrieve tracker topology from geometry
   const TrackerTopology *tTopo = &setup->getData(ttopoToken_);
 
   // Tracker hit association
   TrackerHitAssociator trackertruth(*e, trackerHitAssociatorConfig_);
   // CSC hit association
   CSCHitAssociator csctruth(*e, *setup, cscHitAssociatorConfig_);
   // DT hit association
   bool printRtS = false;
   DTHitAssociator dttruth(*e, *setup, dtHitAssociatorConfig_, printRtS);
   // RPC hit association
   RPCHitAssociator rpctruth(*e, rpcHitAssociatorConfig_);
   // GEM hit association
   GEMHitAssociator gemtruth(*e, gemHitAssociatorConfig_);
 
   MuonAssociatorByHitsHelper::Resources resources = {
       tTopo, &trackertruth, &csctruth, &dttruth, &rpctruth, &gemtruth, {}};
 
   auto bareAssoc = helper_.associateSimToRecoIndices(tH, TPCollectionH, resources);
   for (auto it = bareAssoc.begin(), ed = bareAssoc.end(); it != ed; ++it) {
     for (auto itma = it->second.begin(), edma = it->second.end(); itma != edma; ++itma) {
       outputCollection.insert(TPCollectionH[it->first], std::make_pair(tC[itma->idx], itma->quality));
     }
   }
 
   outputCollection.post_insert();  // perhaps not even necessary
   return outputCollection;
 }