MuonAssociatorByHits Class Reference

#include <MuonAssociatorByHits.h>

Inheritance diagram for MuonAssociatorByHits:

Classes
struct	IndexMatch
struct	RefToBaseSort

Public Types
typedef boost::ptr_vector < uint_SimHitIdpr_pair >	MapOfMatchedIds
typedef std::map < edm::RefToBase< reco::Muon > , std::vector< std::pair < TrackingParticleRef, double > >, RefToBaseSort >	MuonToSimCollection
enum	MuonTrackType { InnerTk, OuterTk, GlobalTk, Segments }
typedef std::pair< uint32_t, EncodedEventId >	SimHitIdpr
typedef std::map < TrackingParticleRef, std::vector< std::pair < edm::RefToBase< reco::Muon > , double > > >	SimToMuonCollection
typedef std::pair< unsigned int, std::vector< SimHitIdpr > >	uint_SimHitIdpr_pair

Public Member Functions
void	associateMuons (MuonToSimCollection &recoToSim, SimToMuonCollection &simToReco, const edm::RefToBaseVector< reco::Muon > &, MuonTrackType, const edm::RefVector< TrackingParticleCollection > &, const edm::Event *event=0, const edm::EventSetup *setup=0) const
void	associateMuons (MuonToSimCollection &recoToSim, SimToMuonCollection &simToReco, const edm::Handle< edm::View< reco::Muon > > &, MuonTrackType, const edm::Handle< TrackingParticleCollection > &, const edm::Event *event=0, const edm::EventSetup *setup=0) const
reco::RecoToSimCollection	associateRecoToSim (const edm::RefToBaseVector< reco::Track > &, const edm::RefVector< TrackingParticleCollection > &, const edm::Event *event=0, const edm::EventSetup *setup=0) const
	Association Reco To Sim with Collections. More...
reco::SimToRecoCollection	associateSimToReco (const edm::RefToBaseVector< reco::Track > &, const edm::RefVector< TrackingParticleCollection > &, const edm::Event *event=0, const edm::EventSetup *setup=0) const
	Association Sim To Reco with Collections. More...
void	getMatchedIds (MapOfMatchedIds &tracker_matchedIds_valid, MapOfMatchedIds &muon_matchedIds_valid, MapOfMatchedIds &tracker_matchedIds_INVALID, MapOfMatchedIds &muon_matchedIds_INVALID, int &n_tracker_valid, int &n_dt_valid, int &n_csc_valid, int &n_rpc_valid, int &n_tracker_matched_valid, int &n_dt_matched_valid, int &n_csc_matched_valid, int &n_rpc_matched_valid, int &n_tracker_INVALID, int &n_dt_INVALID, int &n_csc_INVALID, int &n_rpc_INVALID, int &n_tracker_matched_INVALID, int &n_dt_matched_INVALID, int &n_csc_matched_INVALID, int &n_rpc_matched_INVALID, trackingRecHit_iterator begin, trackingRecHit_iterator end, TrackerHitAssociator *trackertruth, DTHitAssociator &dttruth, MuonTruth &csctruth, RPCHitAssociator &rpctruth, bool printRts, const TrackerTopology *) const
int	getShared (MapOfMatchedIds &matchedIds, TrackingParticleCollection::const_iterator trpart) const
	MuonAssociatorByHits (const edm::ParameterSet &conf, edm::ConsumesCollector &&iC)
	MuonAssociatorByHits (const edm::ParameterSet &conf)
	~MuonAssociatorByHits ()
Public Member Functions inherited from TrackAssociatorBase
virtual reco::RecoToSimCollection	associateRecoToSim (edm::Handle< edm::View< reco::Track > > &tCH, edm::Handle< TrackingParticleCollection > &tPCH, const edm::Event *event, const edm::EventSetup *setup) const
	compare reco to sim the handle of reco::Track and TrackingParticle collections More...
virtual reco::RecoToSimCollectionSeed	associateRecoToSim (edm::Handle< edm::View< TrajectorySeed > > &, edm::Handle< TrackingParticleCollection > &, const edm::Event *event, const edm::EventSetup *setup) const
virtual reco::RecoToSimCollectionTCandidate	associateRecoToSim (edm::Handle< TrackCandidateCollection > &, edm::Handle< TrackingParticleCollection > &, const edm::Event *event, const edm::EventSetup *setup) const
virtual reco::SimToRecoCollection	associateSimToReco (edm::Handle< edm::View< reco::Track > > &tCH, edm::Handle< TrackingParticleCollection > &tPCH, const edm::Event *event, const edm::EventSetup *setup) const
	compare reco to sim the handle of reco::Track and TrackingParticle collections More...
virtual reco::SimToRecoCollectionSeed	associateSimToReco (edm::Handle< edm::View< TrajectorySeed > > &, edm::Handle< TrackingParticleCollection > &, const edm::Event *event, const edm::EventSetup *setup) const
virtual reco::SimToRecoCollectionTCandidate	associateSimToReco (edm::Handle< TrackCandidateCollection > &, edm::Handle< TrackingParticleCollection > &, const edm::Event *event, const edm::EventSetup *setup) const
	TrackAssociatorBase ()
	Constructor. More...
virtual	~TrackAssociatorBase ()
	Destructor. More...

Private Types
typedef std::map< size_t, std::vector< IndexMatch > >	IndexAssociation
typedef std::vector< std::pair < trackingRecHit_iterator, trackingRecHit_iterator > >	TrackHitsCollection

Private Member Functions
IndexAssociation	associateRecoToSimIndices (const TrackHitsCollection &, const edm::RefVector< TrackingParticleCollection > &, const edm::Event *event=0, const edm::EventSetup *setup=0) const
IndexAssociation	associateSimToRecoIndices (const TrackHitsCollection &, const edm::RefVector< TrackingParticleCollection > &, const edm::Event *event=0, const edm::EventSetup *setup=0) const
const TrackingRecHit *	getHitPtr (edm::OwnVector< TrackingRecHit >::const_iterator iter) const
const TrackingRecHit *	getHitPtr (const trackingRecHit_iterator &iter) const
int	LayerFromDetid (const DetId &) const
std::string	write_matched_simtracks (const std::vector< SimHitIdpr > &) const

Private Attributes
const bool	AbsoluteNumberOfHits_muon
const bool	AbsoluteNumberOfHits_track
const bool	acceptOneStubMatchings
const edm::ParameterSet &	conf_
const bool	crossingframe
const bool	dumpDT
const bool	dumpInputCollections
double	EfficiencyCut_muon
double	EfficiencyCut_track
const bool	includeZeroHitMuons
unsigned int	NHitCut_muon
unsigned int	NHitCut_track
double	PurityCut_muon
double	PurityCut_track
edm::InputTag	simtracksTag
edm::EDGetTokenT < edm::SimTrackContainer >	simtracksToken_
edm::InputTag	simtracksXFTag
edm::EDGetTokenT < CrossingFrame< SimTrack > >	simtracksXFToken_
edm::EDGetTokenT < edm::SimVertexContainer >	simvertsToken_
edm::EDGetTokenT < CrossingFrame< SimVertex > >	simvertsXFToken_
const bool	ThreeHitTracksAreSpecial
const bool	UseGrouped
bool	UseMuon
const bool	UsePixels
const bool	UseSplitting
bool	UseTracker

Detailed Description

Definition at line 30 of file MuonAssociatorByHits.h.

Member Typedef Documentation

typedef std::map<size_t, std::vector<IndexMatch> > MuonAssociatorByHits::IndexAssociation

private

Definition at line 134 of file MuonAssociatorByHits.h.

typedef boost::ptr_vector<uint_SimHitIdpr_pair> MuonAssociatorByHits::MapOfMatchedIds

Definition at line 36 of file MuonAssociatorByHits.h.

typedef std::map<edm::RefToBase<reco::Muon>, std::vector<std::pair<TrackingParticleRef, double> >, RefToBaseSort> MuonAssociatorByHits::MuonToSimCollection

Definition at line 78 of file MuonAssociatorByHits.h.

typedef std::pair<uint32_t, EncodedEventId> MuonAssociatorByHits::SimHitIdpr

Definition at line 33 of file MuonAssociatorByHits.h.

typedef std::map<TrackingParticleRef, std::vector<std::pair<edm::RefToBase<reco::Muon>, double> > > MuonAssociatorByHits::SimToMuonCollection

Definition at line 79 of file MuonAssociatorByHits.h.

typedef std::vector<std::pair<trackingRecHit_iterator, trackingRecHit_iterator> > MuonAssociatorByHits::TrackHitsCollection

private

Definition at line 128 of file MuonAssociatorByHits.h.

typedef std::pair<unsigned int,std::vector<SimHitIdpr> > MuonAssociatorByHits::uint_SimHitIdpr_pair

Definition at line 35 of file MuonAssociatorByHits.h.

Member Enumeration Documentation

enum MuonAssociatorByHits::MuonTrackType

Enumerator
InnerTk
OuterTk
GlobalTk
Segments

Definition at line 72 of file MuonAssociatorByHits.h.

{ InnerTk, OuterTk, GlobalTk, Segments };

Constructor & Destructor Documentation

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

Definition at line 19 of file MuonAssociatorByHits.cc.

References conf_, crossingframe, edm::ParameterSet::dump(), EfficiencyCut_muon, includeZeroHitMuons, NHitCut_muon, PurityCut_muon, simtracksToken_, simtracksXFTag, simtracksXFToken_, simvertsToken_, simvertsXFToken_, UseMuon, and UseTracker.

                                                                                                   :  
  includeZeroHitMuons(conf.getParameter<bool>("includeZeroHitMuons")),
  acceptOneStubMatchings(conf.getParameter<bool>("acceptOneStubMatchings")),
  UseTracker(conf.getParameter<bool>("UseTracker")),
  UseMuon(conf.getParameter<bool>("UseMuon")),
  AbsoluteNumberOfHits_track(conf.getParameter<bool>("AbsoluteNumberOfHits_track")),
  NHitCut_track(conf.getParameter<unsigned int>("NHitCut_track")),
  EfficiencyCut_track(conf.getParameter<double>("EfficiencyCut_track")),
  PurityCut_track(conf.getParameter<double>("PurityCut_track")),
  AbsoluteNumberOfHits_muon(conf.getParameter<bool>("AbsoluteNumberOfHits_muon")),
  NHitCut_muon(conf.getParameter<unsigned int>("NHitCut_muon")),
  EfficiencyCut_muon(conf.getParameter<double>("EfficiencyCut_muon")),
  PurityCut_muon(conf.getParameter<double>("PurityCut_muon")),
  UsePixels(conf.getParameter<bool>("UsePixels")),
  UseGrouped(conf.getParameter<bool>("UseGrouped")),
  UseSplitting(conf.getParameter<bool>("UseSplitting")),
  ThreeHitTracksAreSpecial(conf.getParameter<bool>("ThreeHitTracksAreSpecial")),
  dumpDT(conf.getParameter<bool>("dumpDT")),
  dumpInputCollections(conf.getParameter<bool>("dumpInputCollections")),
  crossingframe(conf.getParameter<bool>("crossingframe")),
  simtracksTag(conf.getParameter<edm::InputTag>("simtracksTag")),
  simtracksXFTag(conf.getParameter<edm::InputTag>("simtracksXFTag")),
  conf_(conf)
{
  edm::LogVerbatim("MuonAssociatorByHits") << "constructing  MuonAssociatorByHits" << conf_.dump();

  // up to the user in the other cases - print a message
  if (UseTracker) edm::LogVerbatim("MuonAssociatorByHits")<<"\n UseTracker = TRUE  : Tracker SimHits and RecHits WILL be counted";
  else edm::LogVerbatim("MuonAssociatorByHits") <<"\n UseTracker = FALSE : Tracker SimHits and RecHits WILL NOT be counted";
  
  // up to the user in the other cases - print a message
  if (UseMuon) edm::LogVerbatim("MuonAssociatorByHits")<<" UseMuon = TRUE  : Muon SimHits and RecHits WILL be counted";
  else edm::LogVerbatim("MuonAssociatorByHits") <<" UseMuon = FALSE : Muon SimHits and RecHits WILL NOT be counted"<<endl;
  
  // check consistency of the configuration when allowing zero-hit muon matching (counting invalid hits)
  if (includeZeroHitMuons) {
    edm::LogVerbatim("MuonAssociatorByHits") 
      <<"\n includeZeroHitMuons = TRUE"
      <<"\n ==> (re)set NHitCut_muon = 0, PurityCut_muon = 0, EfficiencyCut_muon = 0"<<endl;
    NHitCut_muon = 0;
    PurityCut_muon = 0.;
    EfficiencyCut_muon = 0.;
  }

  if (crossingframe) {
    simtracksXFToken_=iC.consumes<CrossingFrame<SimTrack> >(simtracksXFTag);
    simvertsXFToken_=iC.consumes<CrossingFrame<SimVertex> >(simtracksXFTag);
  }
  else{
    simtracksToken_=iC.consumes<edm::SimTrackContainer>(edm::InputTag("g4SimHits"));
    simvertsToken_=iC.consumes<edm::SimVertexContainer>(edm::InputTag("g4SimHits"));
  }

  //hack for consumes
  RPCHitAssociator rpctruth(conf_,std::move(iC));
  TrackerMuonHitExtractor hitExtractor(conf_,std::move(iC)); 
  DTHitAssociator dttruth(conf_,std::move(iC));
  MuonTruth muonTruth(conf_,std::move(iC));
}

MuonAssociatorByHits::MuonAssociatorByHits

(

const edm::ParameterSet &

conf

)

Definition at line 80 of file MuonAssociatorByHits.cc.

References conf_, edm::ParameterSet::dump(), EfficiencyCut_muon, includeZeroHitMuons, NHitCut_muon, PurityCut_muon, UseMuon, and UseTracker.

                                                                       :  
  includeZeroHitMuons(conf.getParameter<bool>("includeZeroHitMuons")),
  acceptOneStubMatchings(conf.getParameter<bool>("acceptOneStubMatchings")),
  UseTracker(conf.getParameter<bool>("UseTracker")),
  UseMuon(conf.getParameter<bool>("UseMuon")),
  AbsoluteNumberOfHits_track(conf.getParameter<bool>("AbsoluteNumberOfHits_track")),
  NHitCut_track(conf.getParameter<unsigned int>("NHitCut_track")),
  EfficiencyCut_track(conf.getParameter<double>("EfficiencyCut_track")),
  PurityCut_track(conf.getParameter<double>("PurityCut_track")),
  AbsoluteNumberOfHits_muon(conf.getParameter<bool>("AbsoluteNumberOfHits_muon")),
  NHitCut_muon(conf.getParameter<unsigned int>("NHitCut_muon")),
  EfficiencyCut_muon(conf.getParameter<double>("EfficiencyCut_muon")),
  PurityCut_muon(conf.getParameter<double>("PurityCut_muon")),
  UsePixels(conf.getParameter<bool>("UsePixels")),
  UseGrouped(conf.getParameter<bool>("UseGrouped")),
  UseSplitting(conf.getParameter<bool>("UseSplitting")),
  ThreeHitTracksAreSpecial(conf.getParameter<bool>("ThreeHitTracksAreSpecial")),
  dumpDT(conf.getParameter<bool>("dumpDT")),
  dumpInputCollections(conf.getParameter<bool>("dumpInputCollections")),
  crossingframe(conf.getParameter<bool>("crossingframe")),
  simtracksTag(conf.getParameter<edm::InputTag>("simtracksTag")),
  simtracksXFTag(conf.getParameter<edm::InputTag>("simtracksXFTag")),
  conf_(conf)
{
  edm::LogVerbatim("MuonAssociatorByHits") << "constructing  MuonAssociatorByHits" << conf_.dump();

  // up to the user in the other cases - print a message
  if (UseTracker) edm::LogVerbatim("MuonAssociatorByHits")<<"\n UseTracker = TRUE  : Tracker SimHits and RecHits WILL be counted";
  else edm::LogVerbatim("MuonAssociatorByHits") <<"\n UseTracker = FALSE : Tracker SimHits and RecHits WILL NOT be counted";
  
  // up to the user in the other cases - print a message
  if (UseMuon) edm::LogVerbatim("MuonAssociatorByHits")<<" UseMuon = TRUE  : Muon SimHits and RecHits WILL be counted";
  else edm::LogVerbatim("MuonAssociatorByHits") <<" UseMuon = FALSE : Muon SimHits and RecHits WILL NOT be counted"<<endl;
  
  // check consistency of the configuration when allowing zero-hit muon matching (counting invalid hits)
  if (includeZeroHitMuons) {
    edm::LogVerbatim("MuonAssociatorByHits") 
      <<"\n includeZeroHitMuons = TRUE"
      <<"\n ==> (re)set NHitCut_muon = 0, PurityCut_muon = 0, EfficiencyCut_muon = 0"<<endl;
    NHitCut_muon = 0;
    PurityCut_muon = 0.;
    EfficiencyCut_muon = 0.;
  }

}

MuonAssociatorByHits::~MuonAssociatorByHits

(

)

Definition at line 128 of file MuonAssociatorByHits.cc.

{
}

Member Function Documentation

void MuonAssociatorByHits::associateMuons	(	MuonToSimCollection &	recoToSim,
		SimToMuonCollection &	simToReco,
		const edm::RefToBaseVector< reco::Muon > &	muons,
		MuonTrackType	trackType,
		const edm::RefVector< TrackingParticleCollection > &	tPC,
		const edm::Event *	event = 0,
		const edm::EventSetup *	setup = 0
	)		const

PART 1: Fill MuonAssociatorByHits::TrackHitsCollection

PART 2: call the association routines

Definition at line 1391 of file MuonAssociatorByHits.cc.

References associateRecoToSimIndices(), associateSimToRecoIndices(), edm::RefToBaseVector< T >::begin(), edm::RefVector< C, T, F >::begin(), conf_, edm::RefToBaseVector< T >::end(), edm::RefVector< C, T, F >::end(), TrackerMuonHitExtractor::getMuonHits(), GlobalTk, reco::Muon::globalTrack(), i, TrackerMuonHitExtractor::init(), InnerTk, edm::Ref< C, T, F >::isNonnull(), reco::Muon::isTrackerMuon(), n, OuterTk, reco::Muon::outerTrack(), edm::RefToBaseVector< T >::push_back(), edm::RefVector< C, T, F >::push_back(), edm::OwnVector< T, P >::push_back(), Segments, edm::OwnVector< T, P >::size(), and reco::Muon::track().

Referenced by RecoMuonValidator::analyze(), associateMuons(), and MuonMCClassifier::produce().

                                                                                 {

    MuonAssociatorByHits::TrackHitsCollection muonHitRefs;
    edm::OwnVector<TrackingRecHit> allTMRecHits;  // this I will fill in only for tracker muon hits from segments
    TrackingRecHitRefVector hitRefVector;              // same as above, plus used to get null iterators for muons without a track
    switch (trackType) {
        case InnerTk: 
            for (edm::RefToBaseVector<reco::Muon>::const_iterator it = muons.begin(), ed = muons.end(); it != ed; ++it) {
                edm::RefToBase<reco::Muon> mur = *it;
                if (mur->track().isNonnull()) { 
                    muonHitRefs.push_back(std::make_pair(mur->track()->recHitsBegin(), mur->track()->recHitsEnd()));
                } else {
                    muonHitRefs.push_back(std::make_pair(hitRefVector.begin(), hitRefVector.end()));
                }
            }
            break;
        case OuterTk: 
            for (edm::RefToBaseVector<reco::Muon>::const_iterator it = muons.begin(), ed = muons.end(); it != ed; ++it) {
                edm::RefToBase<reco::Muon> mur = *it;
                if (mur->outerTrack().isNonnull()) { 
                    muonHitRefs.push_back(std::make_pair(mur->outerTrack()->recHitsBegin(), mur->outerTrack()->recHitsEnd()));
                } else {
                    muonHitRefs.push_back(std::make_pair(hitRefVector.begin(), hitRefVector.end()));
                }
            }
            break;
        case GlobalTk: 
            for (edm::RefToBaseVector<reco::Muon>::const_iterator it = muons.begin(), ed = muons.end(); it != ed; ++it) {
                edm::RefToBase<reco::Muon> mur = *it;
                if (mur->globalTrack().isNonnull()) { 
                    muonHitRefs.push_back(std::make_pair(mur->globalTrack()->recHitsBegin(), mur->globalTrack()->recHitsEnd()));
                } else {
                    muonHitRefs.push_back(std::make_pair(hitRefVector.begin(), hitRefVector.end()));
                }
            }
            break;
        case Segments: {
                TrackerMuonHitExtractor hitExtractor(conf_); 
                hitExtractor.init(*event, *setup);
                // puts hits in the vector, and record indices
                std::vector<std::pair<size_t, size_t> >   muonHitIndices;
                for (edm::RefToBaseVector<reco::Muon>::const_iterator it = muons.begin(), ed = muons.end(); it != ed; ++it) {
                    edm::RefToBase<reco::Muon> mur = *it;
                    std::pair<size_t, size_t> indices(allTMRecHits.size(), allTMRecHits.size());
                    if (mur->isTrackerMuon()) {
                        std::vector<const TrackingRecHit *> hits = hitExtractor.getMuonHits(*mur);
                        for (std::vector<const TrackingRecHit *>::const_iterator ith = hits.begin(), edh = hits.end(); ith != edh; ++ith) {
                            allTMRecHits.push_back(**ith);
                        }
                        indices.second += hits.size();
                    }
                    muonHitIndices.push_back(indices);  
                }
                // puts hits in the ref-vector
                for (size_t i = 0, n = allTMRecHits.size(); i < n; ++i) {
                    hitRefVector.push_back(TrackingRecHitRef(& allTMRecHits, i)); 
                }
                // convert indices into pairs of iterators to references
                typedef std::pair<size_t, size_t> index_pair;
                trackingRecHit_iterator hitRefBegin = hitRefVector.begin();
                for (std::vector<std::pair<size_t, size_t> >::const_iterator idxs = muonHitIndices.begin(), idxend = muonHitIndices.end(); idxs != idxend; ++idxs) {
                    muonHitRefs.push_back(std::make_pair(hitRefBegin+idxs->first, 
                                                         hitRefBegin+idxs->second));
                }
                
            }
            break;
    }

    MuonAssociatorByHits::IndexAssociation recSimColl = associateRecoToSimIndices(muonHitRefs,tPC,event,setup);
    for (MuonAssociatorByHits::IndexAssociation::const_iterator it = recSimColl.begin(), ed = recSimColl.end(); it != ed; ++it) {
        edm::RefToBase<reco::Muon> rec = muons[it->first];
        const std::vector<MuonAssociatorByHits::IndexMatch>  & idxAss = it->second;
        std::vector<std::pair<TrackingParticleRef, double> > & tpAss  = recToSim[rec];
        for (std::vector<MuonAssociatorByHits::IndexMatch>::const_iterator ita = idxAss.begin(), eda = idxAss.end(); ita != eda; ++ita) {
            tpAss.push_back(std::make_pair(tPC[ita->idx], ita->quality));
        }
    }
    MuonAssociatorByHits::IndexAssociation simRecColl = associateSimToRecoIndices(muonHitRefs,tPC,event,setup);
    for (MuonAssociatorByHits::IndexAssociation::const_iterator it = simRecColl.begin(), ed = simRecColl.end(); it != ed; ++it) {
        TrackingParticleRef sim = tPC[it->first];
        const std::vector<MuonAssociatorByHits::IndexMatch>  & idxAss = it->second;
        std::vector<std::pair<edm::RefToBase<reco::Muon>, double> > & recAss = simToRec[sim];
        for (std::vector<MuonAssociatorByHits::IndexMatch>::const_iterator ita = idxAss.begin(), eda = idxAss.end(); ita != eda; ++ita) {
            recAss.push_back(std::make_pair(muons[ita->idx], ita->quality));
        }
    }

}

void MuonAssociatorByHits::associateMuons	(	MuonToSimCollection &	recoToSim,
		SimToMuonCollection &	simToReco,
		const edm::Handle< edm::View< reco::Muon > > &	tCH,
		MuonTrackType	type,
		const edm::Handle< TrackingParticleCollection > &	tPCH,
		const edm::Event *	event = 0,
		const edm::EventSetup *	setup = 0
	)		const

Definition at line 1379 of file MuonAssociatorByHits.cc.

References associateMuons(), event(), edm::HandleBase::id(), j, edm::RefVector< C, T, F >::push_back(), and HcalObjRepresent::setup().

                                                                                                      {

    edm::RefVector<TrackingParticleCollection> tpc(tPCH.id());
    for (unsigned int j=0; j<tPCH->size();j++)
      tpc.push_back(edm::Ref<TrackingParticleCollection>(tPCH,j));
    
    associateMuons(recToSim, simToRec, tCH->refVector(),type,tpc,event,setup);
}

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

virtual

Association Reco To Sim with Collections.

Implements TrackAssociatorBase.

Definition at line 133 of file MuonAssociatorByHits.cc.

References associateRecoToSimIndices(), edm::RefToBaseVector< T >::begin(), edm::RefToBaseVector< T >::end(), edm::AssociationMap< Tag >::insert(), and edm::AssociationMap< Tag >::post_insert().

Referenced by MuonAssociatorEDProducer::produce().

                                                                                {
  RecoToSimCollection  outputCollection;

  TrackHitsCollection tH;
  for (edm::RefToBaseVector<reco::Track>::const_iterator it = tC.begin(), ed = tC.end(); it != ed; ++it) {
    tH.push_back(std::make_pair((*it)->recHitsBegin(), (*it)->recHitsEnd()));
  }
  
  IndexAssociation bareAssoc = associateRecoToSimIndices(tH, TPCollectionH, e, setup);
  for (IndexAssociation::const_iterator it = bareAssoc.begin(), ed = bareAssoc.end(); it != ed; ++it) {
    for (std::vector<IndexMatch>::const_iterator itma = it->second.begin(), edma = it->second.end(); itma != edma; ++itma) {
        outputCollection.insert(tC[it->first], std::make_pair(edm::Ref<TrackingParticleCollection>(TPCollectionH, itma->idx), itma->quality));
    }
  }

  outputCollection.post_insert(); // perhaps not even necessary
  return outputCollection;
}

MuonAssociatorByHits::IndexAssociation MuonAssociatorByHits::associateRecoToSimIndices ( const TrackHitsCollection &  tC, const edm::RefVector< TrackingParticleCollection > &  TPCollectionH, const edm::Event *  event = 0, const edm::EventSetup *  setup = 0  ) const

private

Definition at line 155 of file MuonAssociatorByHits.cc.

References AbsoluteNumberOfHits_muon, AbsoluteNumberOfHits_track, acceptOneStubMatchings, conf_, crossingframe, dumpInputCollections, edm::EventSetup::get(), edm::Event::getByLabel(), getMatchedIds(), getShared(), includeZeroHitMuons, j, gen::k, NHitCut_muon, NHitCut_track, edm::Handle< T >::product(), edm::ESHandle< class >::product(), edm::RefVector< C, T, F >::product(), PurityCut_muon, PurityCut_track, simtracksTag, simtracksXFTag, edm::RefVector< C, T, F >::size(), python.multivaluedict::sort(), AlCaHLTBitMon_QueryRunRegistry::string, and ThreeHitTracksAreSpecial.

Referenced by associateMuons(), and associateRecoToSim().

                                                                            {

  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHand;
  setup->get<IdealGeometryRecord>().get(tTopoHand);
  const TrackerTopology *tTopo=tTopoHand.product();

  int tracker_nshared = 0;
  int muon_nshared = 0;
  int global_nshared = 0;

  double tracker_quality = 0;
  double tracker_quality_cut;
  if (AbsoluteNumberOfHits_track) tracker_quality_cut = static_cast<double>(NHitCut_track); 
  else tracker_quality_cut = PurityCut_track;

  double muon_quality = 0;
  double muon_quality_cut;
  if (AbsoluteNumberOfHits_muon) muon_quality_cut = static_cast<double>(NHitCut_muon); 
  else muon_quality_cut = PurityCut_muon;

  double global_quality = 0;
  
  MapOfMatchedIds tracker_matchedIds_valid, muon_matchedIds_valid;
  MapOfMatchedIds tracker_matchedIds_INVALID, muon_matchedIds_INVALID;

  IndexAssociation     outputCollection;
  bool printRtS(true);

  // Tracker hit association  
  TrackerHitAssociator * trackertruth = new TrackerHitAssociator(*e, conf_);
  // CSC hit association
  MuonTruth csctruth(*e,*setup,conf_);
  // DT hit association
  printRtS = true;
  DTHitAssociator dttruth(*e,*setup,conf_,printRtS);
  // RPC hit association
  RPCHitAssociator rpctruth(*e,*setup,conf_);
  
  TrackingParticleCollection tPC;
  if (TPCollectionH.size()!=0) tPC = *(TPCollectionH.product());

  if (dumpInputCollections) {
    // reco::Track collection
    edm::LogVerbatim("MuonAssociatorByHits")<<"\n"<<"reco::Track collection --- size = "<<tC.size();


    // TrackingParticle collection
    edm::LogVerbatim("MuonAssociatorByHits")<<"\n"<<"TrackingParticle collection --- size = "<<tPC.size();
    int j = 0;
    for(TrackingParticleCollection::const_iterator ITER=tPC.begin(); ITER!=tPC.end(); ITER++, j++) {
      edm::LogVerbatim("MuonAssociatorByHits")
    <<"TrackingParticle "<<j<<", q = "<<ITER->charge()<<", p = "<<ITER->p()
    <<", pT = "<<ITER->pt()<<", eta = "<<ITER->eta()<<", phi = "<<ITER->phi();
      
      edm::LogVerbatim("MuonAssociatorByHits")
    <<"\t pdg code = "<<ITER->pdgId()<<", made of "<<ITER->numberOfHits()<<" PSimHit"
    <<" (in "<<ITER->numberOfTrackerLayers()<<" layers)"
    <<" from "<<ITER->g4Tracks().size()<<" SimTrack:";
      for (TrackingParticle::g4t_iterator g4T=ITER->g4Track_begin(); g4T!=ITER->g4Track_end(); g4T++) {
    edm::LogVerbatim("MuonAssociatorByHits")
      <<"\t\t Id:"<<g4T->trackId()<<"/Evt:("<<g4T->eventId().event()<<","<<g4T->eventId().bunchCrossing()<<")";
      }    
    }

    // SimTrack collection
    edm::Handle<CrossingFrame<SimTrack> > cf_simtracks;
    edm::Handle<edm::SimTrackContainer> simTrackCollection;

    // SimVertex collection
    edm::Handle<CrossingFrame<SimVertex> > cf_simvertices;
    edm::Handle<edm::SimVertexContainer> simVertexCollection;
        
    if (crossingframe) {
      e->getByLabel(simtracksXFTag,cf_simtracks);
      auto_ptr<MixCollection<SimTrack> > SimTk( new MixCollection<SimTrack>(cf_simtracks.product()) );
      edm::LogVerbatim("MuonAssociatorByHits")<<"\n"<<"CrossingFrame<SimTrack> collection with InputTag = "<<simtracksXFTag
                          <<" has size = "<<SimTk->size();
      int k = 0;
      for (MixCollection<SimTrack>::MixItr ITER=SimTk->begin(); ITER!=SimTk->end(); ITER++, k++) {
    edm::LogVerbatim("MuonAssociatorByHits")
      <<"SimTrack "<<k
      <<" - Id:"<<ITER->trackId()<<"/Evt:("<<ITER->eventId().event()<<","<<ITER->eventId().bunchCrossing()<<")"
      <<" pdgId = "<<ITER->type()<<", q = "<<ITER->charge()<<", p = "<<ITER->momentum().P()
      <<", pT = "<<ITER->momentum().Pt()<<", eta = "<<ITER->momentum().Eta()<<", phi = "<<ITER->momentum().Phi()
      <<"\n * "<<*ITER <<endl;
      }
      e->getByLabel(simtracksXFTag,cf_simvertices);
      auto_ptr<MixCollection<SimVertex> > SimVtx( new MixCollection<SimVertex>(cf_simvertices.product()) );
      edm::LogVerbatim("MuonAssociatorByHits")<<"\n"<<"CrossingFrame<SimVertex> collection with InputTag = "<<simtracksXFTag
                          <<" has size = "<<SimVtx->size();
      int kv = 0;
      for (MixCollection<SimVertex>::MixItr VITER=SimVtx->begin(); VITER!=SimVtx->end(); VITER++, kv++){
    edm::LogVerbatim("MuonAssociatorByHits")
      <<"SimVertex "<<kv
      << " : "<< *VITER <<endl;
      }
    }
    else {
      e->getByLabel(simtracksTag,simTrackCollection);
      const edm::SimTrackContainer simTC = *(simTrackCollection.product());
      edm::LogVerbatim("MuonAssociatorByHits")<<"\n"<<"SimTrack collection with InputTag = "<<simtracksTag
                          <<" has size = "<<simTC.size()<<endl;
      int k = 0;
      for(edm::SimTrackContainer::const_iterator ITER=simTC.begin(); ITER!=simTC.end(); ITER++, k++){
    edm::LogVerbatim("MuonAssociatorByHits")
      <<"SimTrack "<<k
      <<" - Id:"<<ITER->trackId()<<"/Evt:("<<ITER->eventId().event()<<","<<ITER->eventId().bunchCrossing()<<")"
      <<" pdgId = "<<ITER->type()<<", q = "<<ITER->charge()<<", p = "<<ITER->momentum().P()
      <<", pT = "<<ITER->momentum().Pt()<<", eta = "<<ITER->momentum().Eta()<<", phi = "<<ITER->momentum().Phi()
      <<"\n * "<<*ITER <<endl;
      }
      e->getByLabel(simtracksTag,simVertexCollection);
      const edm::SimVertexContainer simVC = *(simVertexCollection.product());
      edm::LogVerbatim("MuonAssociatorByHits")<<"\n"<<"SimVertex collection with InputTag = "<<"g4SimHits"
                          <<" has size = "<<simVC.size()<<endl;
      int kv = 0;
      for (edm::SimVertexContainer::const_iterator VITER=simVC.begin(); VITER!=simVC.end(); VITER++, kv++){
    edm::LogVerbatim("MuonAssociatorByHits")
      <<"SimVertex "<<kv
      << " : "<< *VITER <<endl;
      }
    }
  }
  
  int tindex=0;
  for (TrackHitsCollection::const_iterator track=tC.begin(); track!=tC.end(); track++, tindex++) {
    edm::LogVerbatim("MuonAssociatorByHits")
      <<"\n"<<"reco::Track "<<tindex
      <<", number of RecHits = "<< (track->second - track->first) << "\n";
    tracker_matchedIds_valid.clear();
    muon_matchedIds_valid.clear();

    tracker_matchedIds_INVALID.clear();
    muon_matchedIds_INVALID.clear();

    bool this_track_matched = false;
    int n_matching_simhits = 0;

    // all hits = valid +INVALID
    int n_all         = 0;        
    int n_tracker_all = 0;
    int n_dt_all      = 0;     
    int n_csc_all     = 0;    
    int n_rpc_all     = 0;    

    int n_valid         = 0;        
    int n_tracker_valid = 0;
    int n_muon_valid    = 0;   
    int n_dt_valid      = 0;     
    int n_csc_valid     = 0;    
    int n_rpc_valid     = 0;    

    int n_tracker_matched_valid = 0;
    int n_muon_matched_valid    = 0;   
    int n_dt_matched_valid      = 0;     
    int n_csc_matched_valid     = 0;    
    int n_rpc_matched_valid     = 0;    

    int n_INVALID         = 0;        
    int n_tracker_INVALID = 0;
    int n_muon_INVALID    = 0;   
    int n_dt_INVALID      = 0;     
    int n_csc_INVALID     = 0;    
    int n_rpc_INVALID     = 0;    
    
    int n_tracker_matched_INVALID = 0;
    int n_muon_matched_INVALID    = 0;     
    int n_dt_matched_INVALID      = 0;     
    int n_csc_matched_INVALID     = 0;    
    int n_rpc_matched_INVALID     = 0;    
    
    printRtS = true;
    getMatchedIds(tracker_matchedIds_valid, muon_matchedIds_valid,
          tracker_matchedIds_INVALID, muon_matchedIds_INVALID,       
          n_tracker_valid, n_dt_valid, n_csc_valid, n_rpc_valid,
          n_tracker_matched_valid, n_dt_matched_valid, n_csc_matched_valid, n_rpc_matched_valid,
          n_tracker_INVALID, n_dt_INVALID, n_csc_INVALID, n_rpc_INVALID,
          n_tracker_matched_INVALID, n_dt_matched_INVALID, n_csc_matched_INVALID, n_rpc_matched_INVALID,
                  track->first, track->second,
          trackertruth, dttruth, csctruth, rpctruth,
          printRtS,tTopo);
    
    n_matching_simhits = tracker_matchedIds_valid.size() + muon_matchedIds_valid.size() + 
                         tracker_matchedIds_INVALID.size() +muon_matchedIds_INVALID.size(); 

    n_muon_valid   = n_dt_valid + n_csc_valid + n_rpc_valid;
    n_valid        = n_tracker_valid + n_muon_valid;
    n_muon_INVALID = n_dt_INVALID + n_csc_INVALID + n_rpc_INVALID;
    n_INVALID      = n_tracker_INVALID + n_muon_INVALID;

    // all used hits (valid+INVALID), defined by UseTracker, UseMuon
    n_tracker_all = n_tracker_valid + n_tracker_INVALID;
    n_dt_all      = n_dt_valid  + n_dt_INVALID;
    n_csc_all     = n_csc_valid + n_csc_INVALID;
    n_rpc_all     = n_rpc_valid + n_rpc_INVALID;
    n_all         = n_valid + n_INVALID;

    n_muon_matched_valid   = n_dt_matched_valid + n_csc_matched_valid + n_rpc_matched_valid;
    n_muon_matched_INVALID = n_dt_matched_INVALID + n_csc_matched_INVALID + n_rpc_matched_INVALID;

    // selected hits are set initially to valid hits
    int n_tracker_selected_hits = n_tracker_valid;
    int n_muon_selected_hits    = n_muon_valid;
    int n_dt_selected_hits      = n_dt_valid;
    int n_csc_selected_hits     = n_csc_valid;
    int n_rpc_selected_hits     = n_rpc_valid;

    // matched hits are a subsample of the selected hits
    int n_tracker_matched = n_tracker_matched_valid;
    int n_muon_matched    = n_muon_matched_valid;
    int n_dt_matched      = n_dt_matched_valid;
    int n_csc_matched     = n_csc_matched_valid;
    int n_rpc_matched     = n_rpc_matched_valid;

    std::string InvMuonHits, ZeroHitMuon;
    
    if (includeZeroHitMuons && n_muon_valid==0 && n_muon_INVALID!=0) {
      // selected muon hits = INVALID when (useZeroHitMuons == True) and track has no valid muon hits

      InvMuonHits = " ***INVALID MUON HITS***";
      ZeroHitMuon = " ***ZERO-HIT MUON***";

      n_muon_selected_hits = n_muon_INVALID;
      n_dt_selected_hits   = n_dt_INVALID;
      n_csc_selected_hits  = n_csc_INVALID;
      n_rpc_selected_hits  = n_rpc_INVALID;

      n_muon_matched = n_muon_matched_INVALID;
      n_dt_matched   = n_dt_matched_INVALID;
      n_csc_matched  = n_csc_matched_INVALID;
      n_rpc_matched  = n_rpc_matched_INVALID;      
    }

    int n_selected_hits = n_tracker_selected_hits + n_muon_selected_hits;
    int n_matched = n_tracker_matched + n_muon_matched;

    edm::LogVerbatim("MuonAssociatorByHits")
      <<"\n"<<"# TrackingRecHits: "<<(track->second - track->first) 
      <<"\n"<< "# used RecHits     = " << n_all <<" ("<<n_tracker_all<<"/"
      <<n_dt_all<<"/"<<n_csc_all<<"/"<<n_rpc_all<<" in Tracker/DT/CSC/RPC)"<<", obtained from " << n_matching_simhits << " SimHits"
      <<"\n"<< "# selected RecHits = " <<n_selected_hits <<" (" <<n_tracker_selected_hits<<"/"
      <<n_dt_selected_hits<<"/"<<n_csc_selected_hits<<"/"<<n_rpc_selected_hits<<" in Tracker/DT/CSC/RPC)"<<InvMuonHits
      <<"\n"<< "# matched RecHits  = " <<n_matched<<" ("<<n_tracker_matched<<"/"
      <<n_dt_matched<<"/"<<n_csc_matched<<"/"<<n_rpc_matched<<" in Tracker/DT/CSC/RPC)";

    if (n_all>0 && n_matching_simhits == 0)
      edm::LogWarning("MuonAssociatorByHits")
    <<"*** WARNING in MuonAssociatorByHits::associateRecoToSim: no matching PSimHit found for this reco::Track !";

    if (n_matching_simhits != 0) {
      edm::LogVerbatim("MuonAssociatorByHits")
    <<"\n"<< "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"
    <<"\n"<< "reco::Track "<<tindex<<ZeroHitMuon
    <<"\n\t"<< "made of "<<n_selected_hits<<" selected RecHits (tracker:"<<n_tracker_selected_hits<<"/muons:"<<n_muon_selected_hits<<")";

      int tpindex = 0;
      for (TrackingParticleCollection::const_iterator trpart = tPC.begin(); trpart != tPC.end(); ++trpart, ++tpindex) {
    tracker_nshared = getShared(tracker_matchedIds_valid, trpart);
    muon_nshared = getShared(muon_matchedIds_valid, trpart);

    if (includeZeroHitMuons && n_muon_valid==0 && n_muon_INVALID!=0) 
      muon_nshared = getShared(muon_matchedIds_INVALID, trpart);

        global_nshared = tracker_nshared + muon_nshared;

    if (AbsoluteNumberOfHits_track) tracker_quality = static_cast<double>(tracker_nshared);
    else if(n_tracker_selected_hits != 0) tracker_quality = (static_cast<double>(tracker_nshared)/static_cast<double>(n_tracker_selected_hits));
    else tracker_quality = 0;

    if (AbsoluteNumberOfHits_muon) muon_quality = static_cast<double>(muon_nshared);
    else if(n_muon_selected_hits != 0) muon_quality = (static_cast<double>(muon_nshared)/static_cast<double>(n_muon_selected_hits));
    else muon_quality = 0;

    // global_quality used to order the matching TPs
    if (n_selected_hits != 0) {
            if (AbsoluteNumberOfHits_muon && AbsoluteNumberOfHits_track) 
                global_quality = global_nshared;
            else
                global_quality = (static_cast<double>(global_nshared)/static_cast<double>(n_selected_hits));
        } else global_quality = 0;

    bool trackerOk = false;
    if (n_tracker_selected_hits != 0) {
      if (tracker_quality > tracker_quality_cut) trackerOk = true;
      //if a track has just 3 hits in the Tracker we require that all 3 hits are shared
      if (ThreeHitTracksAreSpecial && n_tracker_selected_hits==3 && tracker_nshared<3) trackerOk = false;
    }
    
    bool muonOk = false;
    if (n_muon_selected_hits != 0) {
      if (muon_quality > muon_quality_cut) muonOk = true;
    }
    
    // (matchOk) has to account for different track types (tracker-only, standalone muons, global muons)
    bool matchOk = trackerOk || muonOk;

    // only for global muons: match both tracker and muon stub unless (acceptOneStubMatchings==true)
    if (!acceptOneStubMatchings && n_tracker_selected_hits!=0 && n_muon_selected_hits!=0)
      matchOk = trackerOk && muonOk;
    
    if (matchOk) {

          outputCollection[tindex].push_back(IndexMatch(tpindex, global_quality));
      this_track_matched = true;

      edm::LogVerbatim("MuonAssociatorByHits")
        << "\n\t"<<" **MATCHED** with quality = "<<global_quality<<" (tracker: "<<tracker_quality<<" / muon: "<<muon_quality<<")"
        << "\n\t"<<"   N shared hits = "<<global_nshared<<" (tracker: "<<tracker_nshared<<" / muon: "<<muon_nshared<<")"
        <<"\n"<< "   to: TrackingParticle " <<tpindex<<", q = "<<(*trpart).charge()<<", p = "<<(*trpart).p()
        <<", pT = "<<(*trpart).pt()<<", eta = "<<(*trpart).eta()<<", phi = "<<(*trpart).phi()
        <<"\n\t"<< " pdg code = "<<(*trpart).pdgId()<<", made of "<<(*trpart).numberOfHits()<<" PSimHits"
        <<" from "<<(*trpart).g4Tracks().size()<<" SimTrack:";
      for(TrackingParticle::g4t_iterator g4T=(*trpart).g4Track_begin(); 
          g4T!=(*trpart).g4Track_end(); 
          ++g4T) {
        edm::LogVerbatim("MuonAssociatorByHits")
          <<"\t"<< " Id:"<<(*g4T).trackId()<<"/Evt:("<<(*g4T).eventId().event()<<","<<(*g4T).eventId().bunchCrossing()<<")";        
      }
    }
    else {
      // print something only if this TrackingParticle shares some hits with the current reco::Track
      if (global_nshared != 0) 
        edm::LogVerbatim("MuonAssociatorByHits")
          <<"\n\t"<<" NOT matched to TrackingParticle "<<tpindex
          << " with quality = "<<global_quality<<" (tracker: "<<tracker_quality<<" / muon: "<<muon_quality<<")"
          << "\n"<< "   N shared hits = "<<global_nshared<<" (tracker: "<<tracker_nshared<<" / muon: "<<muon_nshared<<")";
    }
    
      }    //  loop over TrackingParticle

      if (!this_track_matched) {
    edm::LogVerbatim("MuonAssociatorByHits")
      <<"\n"<<" NOT matched to any TrackingParticle";
      }
      
      edm::LogVerbatim("MuonAssociatorByHits")
    <<"%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"<<"\n";
     
    }    //  if(n_matching_simhits != 0)
    
  }    // loop over reco::Track

  if (!tC.size()) 
    edm::LogVerbatim("MuonAssociatorByHits")<<"0 reconstructed tracks (-->> 0 associated !)";

  delete trackertruth;
  for (IndexAssociation::iterator it = outputCollection.begin(), ed = outputCollection.end(); it != ed; ++it) {
    std::sort(it->second.begin(), it->second.end());
  }
  return outputCollection;
}

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

virtual

Association Sim To Reco with Collections.

Implements TrackAssociatorBase.

Definition at line 512 of file MuonAssociatorByHits.cc.

References associateSimToRecoIndices(), edm::RefToBaseVector< T >::begin(), edm::RefToBaseVector< T >::end(), edm::AssociationMap< Tag >::insert(), and edm::AssociationMap< Tag >::post_insert().

Referenced by MuonAssociatorEDProducer::produce().

                                                                            {

  SimToRecoCollection  outputCollection;
  TrackHitsCollection tH;
  for (edm::RefToBaseVector<reco::Track>::const_iterator it = tC.begin(), ed = tC.end(); it != ed; ++it) {
    tH.push_back(std::make_pair((*it)->recHitsBegin(), (*it)->recHitsEnd()));
  }
  
  IndexAssociation bareAssoc = associateSimToRecoIndices(tH, TPCollectionH, e, setup);
  for (IndexAssociation::const_iterator it = bareAssoc.begin(), ed = bareAssoc.end(); it != ed; ++it) {
    for (std::vector<IndexMatch>::const_iterator itma = it->second.begin(), edma = it->second.end(); itma != edma; ++itma) {
        outputCollection.insert(edm::Ref<TrackingParticleCollection>(TPCollectionH, it->first),
                                std::make_pair(tC[itma->idx], itma->quality));
    }
  }

  outputCollection.post_insert(); // perhaps not even necessary
  return outputCollection;
}

MuonAssociatorByHits::IndexAssociation MuonAssociatorByHits::associateSimToRecoIndices ( const TrackHitsCollection &  tC, const edm::RefVector< TrackingParticleCollection > &  TPCollectionH, const edm::Event *  event = 0, const edm::EventSetup *  setup = 0  ) const

private

Definition at line 535 of file MuonAssociatorByHits.cc.

References AbsoluteNumberOfHits_muon, AbsoluteNumberOfHits_track, acceptOneStubMatchings, conf_, EfficiencyCut_muon, EfficiencyCut_track, edm::EventSetup::get(), getMatchedIds(), getShared(), includeZeroHitMuons, NHitCut_muon, NHitCut_track, edm::ESHandle< class >::product(), edm::RefVector< C, T, F >::product(), PurityCut_muon, PurityCut_track, edm::RefVector< C, T, F >::size(), python.multivaluedict::sort(), AlCaHLTBitMon_QueryRunRegistry::string, ThreeHitTracksAreSpecial, UseMuon, and UseTracker.

Referenced by associateMuons(), and associateSimToReco().

                                                                            {


  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHand;
  setup->get<IdealGeometryRecord>().get(tTopoHand);
  const TrackerTopology *tTopo=tTopoHand.product();

  int tracker_nshared = 0;
  int muon_nshared = 0;
  int global_nshared = 0;

  double tracker_quality = 0;
  double tracker_quality_cut;
  if (AbsoluteNumberOfHits_track) tracker_quality_cut = static_cast<double>(NHitCut_track); 
  else tracker_quality_cut = EfficiencyCut_track;
  
  double muon_quality = 0;
  double muon_quality_cut;
  if (AbsoluteNumberOfHits_muon) muon_quality_cut = static_cast<double>(NHitCut_muon); 
  else muon_quality_cut = EfficiencyCut_muon;

  double global_quality = 0;

  double tracker_purity = 0;
  double muon_purity = 0;
  double global_purity = 0;
  
  MapOfMatchedIds tracker_matchedIds_valid, muon_matchedIds_valid;
  MapOfMatchedIds tracker_matchedIds_INVALID, muon_matchedIds_INVALID;

  IndexAssociation  outputCollection;

  bool printRtS(true);

  // Tracker hit association  
  TrackerHitAssociator * trackertruth = new TrackerHitAssociator(*e, conf_);
  // CSC hit association
  MuonTruth csctruth(*e,*setup,conf_);
  // DT hit association
  printRtS = false;
  DTHitAssociator dttruth(*e,*setup,conf_,printRtS);  
  // RPC hit association
  RPCHitAssociator rpctruth(*e,*setup,conf_);
 
  TrackingParticleCollection tPC;
  if (TPCollectionH.size()!=0) tPC = *(TPCollectionH.product());

  bool any_trackingParticle_matched = false;
  
  int tindex=0;
  for (TrackHitsCollection::const_iterator track=tC.begin(); track!=tC.end(); track++, tindex++) {
    if (printRtS) edm::LogVerbatim("MuonAssociatorByHits")
      <<"\n"<<"reco::Track "<<tindex
      <<", number of RecHits = "<< (track->second - track->first) << "\n";
    
    tracker_matchedIds_valid.clear();
    muon_matchedIds_valid.clear();

    tracker_matchedIds_INVALID.clear();
    muon_matchedIds_INVALID.clear();

    int n_matching_simhits = 0;

    // all hits = valid +INVALID
    int n_all         = 0;        
    int n_tracker_all = 0;
    int n_dt_all      = 0;     
    int n_csc_all     = 0;    
    int n_rpc_all     = 0;    

    int n_valid         = 0;        
    int n_tracker_valid = 0;
    int n_muon_valid    = 0;   
    int n_dt_valid      = 0;     
    int n_csc_valid     = 0;    
    int n_rpc_valid     = 0;    

    int n_tracker_matched_valid = 0;
    int n_muon_matched_valid    = 0;   
    int n_dt_matched_valid      = 0;     
    int n_csc_matched_valid     = 0;    
    int n_rpc_matched_valid     = 0;    

    int n_INVALID         = 0;        
    int n_tracker_INVALID = 0;
    int n_muon_INVALID    = 0;   
    int n_dt_INVALID      = 0;     
    int n_csc_INVALID     = 0;    
    int n_rpc_INVALID     = 0;    
    
    int n_tracker_matched_INVALID = 0;
    int n_muon_matched_INVALID    = 0;     
    int n_dt_matched_INVALID      = 0;     
    int n_csc_matched_INVALID     = 0;    
    int n_rpc_matched_INVALID     = 0;    
    
    printRtS = false;
    getMatchedIds(tracker_matchedIds_valid, muon_matchedIds_valid,
          tracker_matchedIds_INVALID, muon_matchedIds_INVALID,       
          n_tracker_valid, n_dt_valid, n_csc_valid, n_rpc_valid,
          n_tracker_matched_valid, n_dt_matched_valid, n_csc_matched_valid, n_rpc_matched_valid,
          n_tracker_INVALID, n_dt_INVALID, n_csc_INVALID, n_rpc_INVALID,
          n_tracker_matched_INVALID, n_dt_matched_INVALID, n_csc_matched_INVALID, n_rpc_matched_INVALID,
                  track->first, track->second,
          trackertruth, dttruth, csctruth, rpctruth,
          printRtS,tTopo);
    
    n_matching_simhits = tracker_matchedIds_valid.size() + muon_matchedIds_valid.size() + 
                         tracker_matchedIds_INVALID.size() +muon_matchedIds_INVALID.size(); 

    n_muon_valid   = n_dt_valid + n_csc_valid + n_rpc_valid;
    n_valid        = n_tracker_valid + n_muon_valid;
    n_muon_INVALID = n_dt_INVALID + n_csc_INVALID + n_rpc_INVALID;
    n_INVALID      = n_tracker_INVALID + n_muon_INVALID;

    // all used hits (valid+INVALID), defined by UseTracker, UseMuon
    n_tracker_all = n_tracker_valid + n_tracker_INVALID;
    n_dt_all      = n_dt_valid  + n_dt_INVALID;
    n_csc_all     = n_csc_valid + n_csc_INVALID;
    n_rpc_all     = n_rpc_valid + n_rpc_INVALID;
    n_all         = n_valid + n_INVALID;

    n_muon_matched_valid   = n_dt_matched_valid + n_csc_matched_valid + n_rpc_matched_valid;
    n_muon_matched_INVALID = n_dt_matched_INVALID + n_csc_matched_INVALID + n_rpc_matched_INVALID;

     // selected hits are set initially to valid hits
    int n_tracker_selected_hits = n_tracker_valid;
    int n_muon_selected_hits    = n_muon_valid;
    int n_dt_selected_hits      = n_dt_valid;
    int n_csc_selected_hits     = n_csc_valid;
    int n_rpc_selected_hits     = n_rpc_valid;

    // matched hits are a subsample of the selected hits
    int n_tracker_matched = n_tracker_matched_valid;
    int n_muon_matched    = n_muon_matched_valid;
    int n_dt_matched      = n_dt_matched_valid;
    int n_csc_matched     = n_csc_matched_valid;
    int n_rpc_matched     = n_rpc_matched_valid;

    std::string InvMuonHits, ZeroHitMuon;

    if (includeZeroHitMuons && n_muon_valid==0 && n_muon_INVALID!=0) {
      // selected muon hits = INVALID when (useZeroHitMuons == True) and track has no valid muon hits
      
      InvMuonHits = " ***INVALID MUON HITS***";
      ZeroHitMuon = " ***ZERO-HIT MUON***";

      n_muon_selected_hits = n_muon_INVALID;
      n_dt_selected_hits   = n_dt_INVALID;
      n_csc_selected_hits  = n_csc_INVALID;
      n_rpc_selected_hits  = n_rpc_INVALID;

      n_muon_matched = n_muon_matched_INVALID;
      n_dt_matched   = n_dt_matched_INVALID;
      n_csc_matched  = n_csc_matched_INVALID;
      n_rpc_matched  = n_rpc_matched_INVALID;
    }

    int n_selected_hits = n_tracker_selected_hits + n_muon_selected_hits;
    int n_matched = n_tracker_matched + n_muon_matched;

    if (printRtS) edm::LogVerbatim("MuonAssociatorByHits")
      <<"\n"<<"# TrackingRecHits: "<<(track->second - track->first) 
      <<"\n"<< "# used RecHits     = " <<n_all    <<" ("<<n_tracker_all<<"/"
      <<n_dt_all<<"/"<<n_csc_all<<"/"<<n_rpc_all<<" in Tracker/DT/CSC/RPC)"<<", obtained from " << n_matching_simhits << " SimHits"
      <<"\n"<< "# selected RecHits = " <<n_selected_hits  <<" (" <<n_tracker_selected_hits<<"/"
      <<n_dt_selected_hits<<"/"<<n_csc_selected_hits<<"/"<<n_rpc_selected_hits<<" in Tracker/DT/CSC/RPC)"<<InvMuonHits
      <<"\n"<< "# matched RecHits = " <<n_matched<<" ("<<n_tracker_matched<<"/"
      <<n_dt_matched<<"/"<<n_csc_matched<<"/"<<n_rpc_matched<<" in Tracker/DT/CSC/RPC)";
    
    if (printRtS && n_all>0 && n_matching_simhits==0)
      edm::LogWarning("MuonAssociatorByHits")
    <<"*** WARNING in MuonAssociatorByHits::associateSimToReco: no matching PSimHit found for this reco::Track !";
    
    if (n_matching_simhits != 0) {
      int tpindex =0;
      for (TrackingParticleCollection::const_iterator trpart = tPC.begin(); trpart != tPC.end(); ++trpart, ++tpindex) {

    //  int n_tracker_simhits = 0;
    int n_tracker_recounted_simhits = 0; 
    int n_muon_simhits = 0; 
    int n_global_simhits = 0; 
    //  std::vector<PSimHit> tphits;

    int n_tracker_selected_simhits = 0;
    int n_muon_selected_simhits = 0; 
    int n_global_selected_simhits = 0; 

    // shared hits are counted over the selected ones
    tracker_nshared = getShared(tracker_matchedIds_valid, trpart);
    muon_nshared = getShared(muon_matchedIds_valid, trpart);

        if (includeZeroHitMuons && n_muon_valid==0 && n_muon_INVALID!=0)
      muon_nshared = getShared(muon_matchedIds_INVALID, trpart);
    
        global_nshared = tracker_nshared + muon_nshared;    
        if (global_nshared == 0) continue; // if this TP shares no hits with the current reco::Track loop over 

    // This does not work with the new TP interface 
    /*
    for(std::vector<PSimHit>::const_iterator TPhit = trpart->pSimHit_begin(); TPhit != trpart->pSimHit_end(); TPhit++) {
          DetId dId = DetId(TPhit->detUnitId());
      DetId::Detector detector = dId.det();
      
      if (detector == DetId::Tracker) {
        n_tracker_simhits++;

        unsigned int subdetId = static_cast<unsigned int>(dId.subdetId());
        if (!UsePixels && (subdetId==PixelSubdetector::PixelBarrel || subdetId==PixelSubdetector::PixelEndcap) )
          continue;

        SiStripDetId* stripDetId = 0;
        if (subdetId==SiStripDetId::TIB||subdetId==SiStripDetId::TOB||
        subdetId==SiStripDetId::TID||subdetId==SiStripDetId::TEC)
          stripDetId= new SiStripDetId(dId);
        
        bool newhit = true;
        for(std::vector<PSimHit>::const_iterator TPhitOK = tphits.begin(); TPhitOK != tphits.end(); TPhitOK++) {
          DetId dIdOK = DetId(TPhitOK->detUnitId());
          //no grouped, no splitting
          if (!UseGrouped && !UseSplitting)
        if (tTopo->layer(dId)==tTopo->layer(dIdOK) &&
            dId.subdetId()==dIdOK.subdetId()) newhit = false;
          //no grouped, splitting
          if (!UseGrouped && UseSplitting)
        if (tTopo->layer(dId)==tTopo->layer(dIdOK) &&
            dId.subdetId()==dIdOK.subdetId() &&
            (stripDetId==0 || stripDetId->partnerDetId()!=dIdOK.rawId()))
          newhit = false;
          //grouped, no splitting
          if (UseGrouped && !UseSplitting)
        if ( tTopo->layer(dId)== tTopo->layer(dIdOK) &&
            dId.subdetId()==dIdOK.subdetId() &&
            stripDetId!=0 && stripDetId->partnerDetId()==dIdOK.rawId())
          newhit = false;
          //grouped, splitting
          if (UseGrouped && UseSplitting)
        newhit = true;
        }
        if (newhit) {
          tphits.push_back(*TPhit);
        }
        delete stripDetId;
      }
      else if (detector == DetId::Muon) {
        n_muon_simhits++;
        
        // discard BAD CSC chambers (ME4/2) from hit counting
        if (dId.subdetId() == MuonSubdetId::CSC) {
              if (csctruth.cscBadChambers->isInBadChamber(CSCDetId(dId))) {
        // edm::LogVerbatim("MuonAssociatorByHits")<<"This PSimHit is in a BAD CSC chamber, CSCDetId = "<<CSCDetId(dId);
        n_muon_simhits--;
          }
        }
        
      }
    }
    */
    //  n_tracker_recounted_simhits = tphits.size();

        // adapt to new TP interface: this gives the total number of hits in tracker
        //   should reproduce the behaviour of UseGrouped=UseSplitting=.true.
    n_tracker_recounted_simhits = trpart->numberOfTrackerHits();
        //   numberOfHits() gives the total number of hits (tracker + muons)
        n_muon_simhits = trpart->numberOfHits() - trpart->numberOfTrackerHits();

        // Handle the case of TrackingParticles that don't have PSimHits inside, e.g. because they were made on RECOSIM only.
        if (trpart->numberOfHits()==0) {
            // FIXME this can be made better, counting the digiSimLinks associated to this TP, but perhaps it's not worth it
            n_tracker_recounted_simhits = tracker_nshared;
            n_muon_simhits = muon_nshared;
        }   
    n_global_simhits = n_tracker_recounted_simhits + n_muon_simhits;

    if (UseMuon) {
      n_muon_selected_simhits = n_muon_simhits;
      n_global_selected_simhits = n_muon_selected_simhits;
    }
    if (UseTracker) {
      n_tracker_selected_simhits = n_tracker_recounted_simhits;
      n_global_selected_simhits += n_tracker_selected_simhits;
    }

    if (AbsoluteNumberOfHits_track) tracker_quality = static_cast<double>(tracker_nshared);
    else if (n_tracker_selected_simhits!=0) 
      tracker_quality = static_cast<double>(tracker_nshared)/static_cast<double>(n_tracker_selected_simhits);
    else tracker_quality = 0;
    
    if (AbsoluteNumberOfHits_muon) muon_quality = static_cast<double>(muon_nshared);
    else if (n_muon_selected_simhits!=0) 
      muon_quality = static_cast<double>(muon_nshared)/static_cast<double>(n_muon_selected_simhits);
    else muon_quality = 0;

    // global_quality used to order the matching tracks
    if (n_global_selected_simhits != 0) {
      if (AbsoluteNumberOfHits_muon && AbsoluteNumberOfHits_track)
        global_quality = global_nshared;
      else 
        global_quality = static_cast<double>(global_nshared)/static_cast<double>(n_global_selected_simhits);
    } 
    else global_quality = 0;

    // global purity
    if (n_selected_hits != 0) {
      if (AbsoluteNumberOfHits_muon && AbsoluteNumberOfHits_track)
        global_purity = global_nshared;
      else
        global_purity = static_cast<double>(global_nshared)/static_cast<double>(n_selected_hits);
    }
    else global_purity = 0;
    
    bool trackerOk = false;
    if (n_tracker_selected_hits != 0) {
      if (tracker_quality > tracker_quality_cut) trackerOk = true;
      
      tracker_purity = static_cast<double>(tracker_nshared)/static_cast<double>(n_tracker_selected_hits);
      if (AbsoluteNumberOfHits_track) tracker_purity = static_cast<double>(tracker_nshared);

      if ((!AbsoluteNumberOfHits_track) && tracker_purity <= PurityCut_track) trackerOk = false;
      
      //if a track has just 3 hits in the Tracker we require that all 3 hits are shared
      if (ThreeHitTracksAreSpecial && n_tracker_selected_hits==3 && tracker_nshared<3) trackerOk = false;
    }
    
    bool muonOk = false;
    if (n_muon_selected_hits != 0) {
      if (muon_quality > muon_quality_cut) muonOk = true;
      
      muon_purity = static_cast<double>(muon_nshared)/static_cast<double>(n_muon_selected_hits);
      if (AbsoluteNumberOfHits_muon) muon_purity = static_cast<double>(muon_nshared);

      if ((!AbsoluteNumberOfHits_muon) &&  muon_purity <= PurityCut_muon) muonOk = false;
    }

    // (matchOk) has to account for different track types (tracker-only, standalone muons, global muons)
    bool matchOk = trackerOk || muonOk;

    // only for global muons: match both tracker and muon stub unless (acceptOneStubMatchings==true)
    if (!acceptOneStubMatchings && n_tracker_selected_hits!=0 && n_muon_selected_hits!=0)
      matchOk = trackerOk && muonOk;
    
    if (matchOk) {
      
          outputCollection[tpindex].push_back(IndexMatch(tindex,global_quality));
      any_trackingParticle_matched = true;
      
      edm::LogVerbatim("MuonAssociatorByHits")
        <<"************************************************************************************************************************"  
        <<"\n"<< "TrackingParticle " << tpindex <<", q = "<<(*trpart).charge()<<", p = "<<(*trpart).p()
        <<", pT = "<<(*trpart).pt()<<", eta = "<<(*trpart).eta()<<", phi = "<<(*trpart).phi()
        <<"\n"<<" pdg code = "<<(*trpart).pdgId()
        <<", made of "<<(*trpart).numberOfHits()<<" PSimHits, recounted "<<n_global_simhits<<" PSimHits"
        <<" (tracker:"<<n_tracker_recounted_simhits<<"/muons:"<<n_muon_simhits<<")"
        <<", from "<<(*trpart).g4Tracks().size()<<" SimTrack:";
      for(TrackingParticle::g4t_iterator g4T=(*trpart).g4Track_begin(); 
          g4T!=(*trpart).g4Track_end(); 
          ++g4T) {
        edm::LogVerbatim("MuonAssociatorByHits")
          <<" Id:"<<(*g4T).trackId()<<"/Evt:("<<(*g4T).eventId().event()<<","<<(*g4T).eventId().bunchCrossing()<<")";
      }
      edm::LogVerbatim("MuonAssociatorByHits")
        <<"\t selected "<<n_global_selected_simhits<<" PSimHits"
        <<" (tracker:"<<n_tracker_selected_simhits<<"/muons:"<<n_muon_selected_simhits<<")"
        << "\n\t **MATCHED** with quality = "<<global_quality<<" (tracker: "<<tracker_quality<<" / muon: "<<muon_quality<<")"
        << "\n\t               and purity = "<<global_purity<<" (tracker: "<<tracker_purity<<" / muon: "<<muon_purity<<")"
        << "\n\t   N shared hits = "<<global_nshared<<" (tracker: "<<tracker_nshared<<" / muon: "<<muon_nshared<<")"
        <<"\n" <<"   to: reco::Track "<<tindex<<ZeroHitMuon
        <<"\n\t"<< " made of "<<n_selected_hits<<" RecHits (tracker:"<<n_tracker_valid<<"/muons:"<<n_muon_selected_hits<<")";
    }
    else {
      // print something only if this TrackingParticle shares some hits with the current reco::Track
      if (global_nshared != 0) {
        if (printRtS) edm::LogVerbatim("MuonAssociatorByHits")
          <<"************************************************************************************************************************"  
          <<"\n"<<"TrackingParticle " << tpindex <<", q = "<<(*trpart).charge()<<", p = "<<(*trpart).p()
          <<", pT = "<<(*trpart).pt()<<", eta = "<<(*trpart).eta()<<", phi = "<<(*trpart).phi()
          <<"\n"<<" pdg code = "<<(*trpart).pdgId()
          <<", made of "<<(*trpart).numberOfHits()<<" PSimHits, recounted "<<n_global_simhits<<" PSimHits"
          <<" (tracker:"<<n_tracker_recounted_simhits<<"/muons:"<<n_muon_simhits<<")"
          <<", from "<<(*trpart).g4Tracks().size()<<" SimTrack:";
        for(TrackingParticle::g4t_iterator g4T=(*trpart).g4Track_begin(); 
        g4T!=(*trpart).g4Track_end(); 
        ++g4T) {
          if (printRtS) edm::LogVerbatim("MuonAssociatorByHits")
        <<" Id:"<<(*g4T).trackId()<<"/Evt:("<<(*g4T).eventId().event()<<","<<(*g4T).eventId().bunchCrossing()<<")";
        }
        if (printRtS) edm::LogVerbatim("MuonAssociatorByHits")
          <<"\t selected "<<n_global_selected_simhits<<" PSimHits"
          <<" (tracker:"<<n_tracker_selected_simhits<<"/muons:"<<n_muon_selected_simhits<<")"
          <<"\n\t NOT matched  to reco::Track "<<tindex<<ZeroHitMuon
          <<" with quality = "<<global_quality<<" (tracker: "<<tracker_quality<<" / muon: "<<muon_quality<<")"
          <<"\n\t and purity = "<<global_purity<<" (tracker: "<<tracker_purity<<" / muon: "<<muon_purity<<")"
          <<"\n\t     N shared hits = "<<global_nshared<<" (tracker: "<<tracker_nshared<<" / muon: "<<muon_nshared<<")";
      }
    }
      }  // loop over TrackingParticle's
    }   // if(n_matching_simhits != 0)
  }    // loop over reco Tracks
  
  if (!any_trackingParticle_matched) {
    edm::LogVerbatim("MuonAssociatorByHits")
      <<"\n"
      <<"************************************************************************************************************************"
      << "\n NO TrackingParticle associated to ANY input reco::Track ! \n"
      <<"************************************************************************************************************************"<<"\n";  
  } else {
    edm::LogVerbatim("MuonAssociatorByHits")
      <<"************************************************************************************************************************"<<"\n";  
  }
  
  delete trackertruth;
  for (IndexAssociation::iterator it = outputCollection.begin(), ed = outputCollection.end(); it != ed; ++it) {
    std::sort(it->second.begin(), it->second.end());
  }
  return outputCollection;
}

const TrackingRecHit* MuonAssociatorByHits::getHitPtr ( edm::OwnVector< TrackingRecHit >::const_iterator  iter ) const

inlineprivate

Definition at line 122 of file MuonAssociatorByHits.h.

References getDQMSummary::iter.

{return &*iter;}

const TrackingRecHit* MuonAssociatorByHits::getHitPtr ( const trackingRecHit_iterator &  iter ) const

inlineprivate

Definition at line 123 of file MuonAssociatorByHits.h.

References getDQMSummary::iter.

{return &**iter;}

void MuonAssociatorByHits::getMatchedIds	(	MapOfMatchedIds &	tracker_matchedIds_valid,
		MapOfMatchedIds &	muon_matchedIds_valid,
		MapOfMatchedIds &	tracker_matchedIds_INVALID,
		MapOfMatchedIds &	muon_matchedIds_INVALID,
		int &	n_tracker_valid,
		int &	n_dt_valid,
		int &	n_csc_valid,
		int &	n_rpc_valid,
		int &	n_tracker_matched_valid,
		int &	n_dt_matched_valid,
		int &	n_csc_matched_valid,
		int &	n_rpc_matched_valid,
		int &	n_tracker_INVALID,
		int &	n_dt_INVALID,
		int &	n_csc_INVALID,
		int &	n_rpc_INVALID,
		int &	n_tracker_matched_INVALID,
		int &	n_dt_matched_INVALID,
		int &	n_csc_matched_INVALID,
		int &	n_rpc_matched_INVALID,
		trackingRecHit_iterator	begin,
		trackingRecHit_iterator	end,
		TrackerHitAssociator *	trackertruth,
		DTHitAssociator &	dttruth,
		MuonTruth &	csctruth,
		RPCHitAssociator &	rpctruth,
		bool	printRts,
		const TrackerTopology *	tTopo
	)		const

Definition at line 957 of file MuonAssociatorByHits.cc.

References MuonTruth::associateCSCHitId(), DTHitAssociator::associateDTHitId(), DTHitAssociator::associateHit(), TrackerHitAssociator::associateHitId(), RPCHitAssociator::associateRecHit(), MuonSubdetId::CSC, CSCDetId, DetId::det(), cond::rpcobgas::detid, MuonSubdetId::DT, end, TrackingRecHit::geographicalId(), DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), cmsHarvester::index, TrackingRecHit::isValid(), j, DetId::Muon, DTRecSegment4D::phiSegment(), TrackerTopology::print(), DetId::rawId(), CSCSegment::recHits(), DTRecSegment2D::recHits(), MuonSubdetId::RPC, RPCDetId, DetId::subdetId(), DetId::Tracker, DTRecHit1D::wireId(), and DTRecSegment4D::zSegment().

Referenced by associateRecoToSimIndices(), and associateSimToRecoIndices().

{
  tracker_matchedIds_valid.clear();
  muon_matchedIds_valid.clear();

  tracker_matchedIds_INVALID.clear();
  muon_matchedIds_INVALID.clear();

  n_tracker_valid = 0;
  n_dt_valid  = 0;
  n_csc_valid = 0;
  n_rpc_valid = 0;

  n_tracker_matched_valid = 0;
  n_dt_matched_valid  = 0;
  n_csc_matched_valid = 0;
  n_rpc_matched_valid = 0;
  
  n_tracker_INVALID = 0;
  n_dt_INVALID  = 0;
  n_csc_INVALID = 0;
  n_rpc_INVALID = 0;
  
  n_tracker_matched_INVALID = 0;
  n_dt_matched_INVALID  = 0;
  n_csc_matched_INVALID = 0;
  n_rpc_matched_INVALID = 0;
  
  std::vector<SimHitIdpr> SimTrackIds;

  // main loop on TrackingRecHits
  int iloop = 0;
  int iH = -1;
  for (trackingRecHit_iterator it = begin;  it != end; it++, iloop++) {
    stringstream hit_index;     
    hit_index<<iloop;

    const TrackingRecHit * hitp = getHitPtr(it);
    DetId geoid = hitp->geographicalId();    

    unsigned int detid = geoid.rawId();    
    stringstream detector_id;
    detector_id<<detid;

    string hitlog = "TrackingRecHit "+hit_index.str();
    string wireidlog;
    std::vector<string> DTSimHits;
    
    DetId::Detector det = geoid.det();
    int subdet = geoid.subdetId();
    
    bool valid_Hit = hitp->isValid();
    
    // Si-Tracker Hits
    if (det == DetId::Tracker && UseTracker) {
      stringstream detector_id;
      detector_id<< tTopo->print(detid);
      
      if (valid_Hit) hitlog = hitlog+" -Tracker - detID = "+detector_id.str();
      else hitlog = hitlog+" *** INVALID ***"+" -Tracker - detID = "+detector_id.str();
      
      iH++;
      SimTrackIds = trackertruth->associateHitId(*hitp);

      if (valid_Hit) {
    n_tracker_valid++;
    
    if(!SimTrackIds.empty()) {
      n_tracker_matched_valid++;
      //tracker_matchedIds_valid[iH] = SimTrackIds;
          tracker_matchedIds_valid.push_back( new uint_SimHitIdpr_pair(iH, SimTrackIds));
    }
      } else {
    n_tracker_INVALID++;

    if(!SimTrackIds.empty()) {
      n_tracker_matched_INVALID++;
      //tracker_matchedIds_INVALID[iH] = SimTrackIds;
          tracker_matchedIds_INVALID.push_back( new uint_SimHitIdpr_pair(iH, SimTrackIds));
    }
      }
    }  
    // Muon detector Hits    
    else if (det == DetId::Muon && UseMuon) {

      // DT Hits      
      if (subdet == MuonSubdetId::DT) {    
    DTWireId dtdetid = DTWireId(detid);
    stringstream dt_detector_id;
    dt_detector_id << dtdetid;
    if (valid_Hit) hitlog = hitlog+" -Muon DT - detID = "+dt_detector_id.str();   
    else hitlog = hitlog+" *** INVALID ***"+" -Muon DT - detID = "+dt_detector_id.str();      
    
    const DTRecHit1D * dtrechit = dynamic_cast<const DTRecHit1D *>(hitp);
    
    // single DT hits
    if (dtrechit) {
      iH++;
      SimTrackIds = dttruth.associateDTHitId(dtrechit);  
      
      if (valid_Hit) {
        n_dt_valid++;

        if (!SimTrackIds.empty()) {
          n_dt_matched_valid++;
          //muon_matchedIds_valid[iH] = SimTrackIds;
              muon_matchedIds_valid.push_back (new uint_SimHitIdpr_pair(iH,SimTrackIds));
        }
      } else {
        n_dt_INVALID++;
        
        if (!SimTrackIds.empty()) {
          n_dt_matched_INVALID++;
          //muon_matchedIds_INVALID[iH] = SimTrackIds;
              muon_matchedIds_INVALID.push_back (new uint_SimHitIdpr_pair(iH,SimTrackIds));
        }
      }

      if (dumpDT) {
        DTWireId wireid = dtrechit->wireId();
        stringstream wid; 
        wid<<wireid;
        std::vector<PSimHit> dtSimHits = dttruth.associateHit(*hitp);
        
        stringstream ndthits;
        ndthits<<dtSimHits.size();
        wireidlog = "\t DTWireId :"+wid.str()+", "+ndthits.str()+" associated PSimHit :";
        
        for (unsigned int j=0; j<dtSimHits.size(); j++) {
          stringstream index;
          index<<j;
          stringstream simhit;
          simhit<<dtSimHits[j];
          string simhitlog = "\t\t PSimHit "+index.str()+": "+simhit.str();
          DTSimHits.push_back(simhitlog);
        }
      }  // if (dumpDT)
    }

    // DT segments  
    else {
      const DTRecSegment4D * dtsegment = dynamic_cast<const DTRecSegment4D *>(hitp);
      
      if (dtsegment) {
        
        std::vector<const TrackingRecHit *> componentHits, phiHits, zHits;
        if (dtsegment->hasPhi()) {
          phiHits = dtsegment->phiSegment()->recHits();
          componentHits.insert(componentHits.end(),phiHits.begin(),phiHits.end());
        }
        if (dtsegment->hasZed()) {
          zHits = dtsegment->zSegment()->recHits();
          componentHits.insert(componentHits.end(),zHits.begin(),zHits.end());
        }
        if (printRtS) edm::LogVerbatim("MuonAssociatorByHits")
          <<"\n\t this TrackingRecHit is a DTRecSegment4D with "
          <<componentHits.size()<<" hits (phi:"<<phiHits.size()<<", z:"<<zHits.size()<<")";
        
        std::vector<SimHitIdpr> i_SimTrackIds;
        int i_compHit = 0;
        for (std::vector<const TrackingRecHit *>::const_iterator ithit =componentHits.begin(); 
         ithit != componentHits.end(); ++ithit) {
          i_compHit++;
          
          const DTRecHit1D * dtrechit1D = dynamic_cast<const DTRecHit1D *>(*ithit);
          
          i_SimTrackIds.clear();
          if (dtrechit1D) {
        iH++;
        i_SimTrackIds = dttruth.associateDTHitId(dtrechit1D);  
        
        if (valid_Hit) {
          // validity check is on the segment, but hits are counted one-by-one
          n_dt_valid++; 

          if (!i_SimTrackIds.empty()) {
            n_dt_matched_valid++;
            //muon_matchedIds_valid[iH] = i_SimTrackIds;
                    muon_matchedIds_valid.push_back (new uint_SimHitIdpr_pair(iH,i_SimTrackIds));
          }
        } else {
          n_dt_INVALID++;
          
          if (!i_SimTrackIds.empty()) {
            n_dt_matched_INVALID++;
            //muon_matchedIds_INVALID[iH] = i_SimTrackIds;
                    muon_matchedIds_INVALID.push_back (new uint_SimHitIdpr_pair(iH,i_SimTrackIds));
                    
          } 
        }
          } else if (printRtS) edm::LogWarning("MuonAssociatorByHits")
        <<"*** WARNING in MuonAssociatorByHits::getMatchedIds, null dynamic_cast of a DT TrackingRecHit !";
          
          unsigned int i_detid = (*ithit)->geographicalId().rawId();
          DTWireId i_dtdetid = DTWireId(i_detid);

          stringstream i_dt_detector_id;
          i_dt_detector_id << i_dtdetid;

          stringstream i_ss;
          i_ss<<"\t\t hit "<<i_compHit<<" -Muon DT - detID = "<<i_dt_detector_id.str();

          string i_hitlog = i_ss.str();
          i_hitlog = i_hitlog + write_matched_simtracks(i_SimTrackIds);
          if (printRtS) edm::LogVerbatim("MuonAssociatorByHits") << i_hitlog;
          
          SimTrackIds.insert(SimTrackIds.end(),i_SimTrackIds.begin(),i_SimTrackIds.end());
        }         
      }  // if (dtsegment)

      else if (printRtS) edm::LogWarning("MuonAssociatorByHits")
        <<"*** WARNING in MuonAssociatorByHits::getMatchedIds, DT TrackingRecHit is neither DTRecHit1D nor DTRecSegment4D ! ";      
    }
      }
      
      // CSC Hits
      else if (subdet == MuonSubdetId::CSC) {
    CSCDetId cscdetid = CSCDetId(detid);
    stringstream csc_detector_id;
    csc_detector_id << cscdetid;
    if (valid_Hit) hitlog = hitlog+" -Muon CSC- detID = "+csc_detector_id.str();      
    else hitlog = hitlog+" *** INVALID ***"+" -Muon CSC- detID = "+csc_detector_id.str();     
    
    const CSCRecHit2D * cscrechit = dynamic_cast<const CSCRecHit2D *>(hitp);
    
    // single CSC hits
    if (cscrechit) {
      iH++;
      SimTrackIds = csctruth.associateCSCHitId(cscrechit);
      
      if (valid_Hit) {
        n_csc_valid++;
        
        if (!SimTrackIds.empty()) {
          n_csc_matched_valid++;
          //muon_matchedIds_valid[iH] = SimTrackIds;
              muon_matchedIds_valid.push_back (new uint_SimHitIdpr_pair(iH,SimTrackIds));
        }
      } else {
        n_csc_INVALID++;
        
        if (!SimTrackIds.empty()) {
          n_csc_matched_INVALID++;
          //muon_matchedIds_INVALID[iH] = SimTrackIds;
              muon_matchedIds_INVALID.push_back (new uint_SimHitIdpr_pair(iH,SimTrackIds));
        }
      }
    }
    
    // CSC segments
    else {
      const CSCSegment * cscsegment = dynamic_cast<const CSCSegment *>(hitp);
      
      if (cscsegment) {
        
        std::vector<const TrackingRecHit *> componentHits = cscsegment->recHits();
        if (printRtS) edm::LogVerbatim("MuonAssociatorByHits")
          <<"\n\t this TrackingRecHit is a CSCSegment with "<<componentHits.size()<<" hits";
        
        std::vector<SimHitIdpr> i_SimTrackIds;
        int i_compHit = 0;
        for (std::vector<const TrackingRecHit *>::const_iterator ithit =componentHits.begin(); 
         ithit != componentHits.end(); ++ithit) {
          i_compHit++;
          
          const CSCRecHit2D * cscrechit2D = dynamic_cast<const CSCRecHit2D *>(*ithit);
          
          i_SimTrackIds.clear();
          if (cscrechit2D) {
        iH++;
        i_SimTrackIds = csctruth.associateCSCHitId(cscrechit2D);

        if (valid_Hit) {
          // validity check is on the segment, but hits are counted one-by-one
          n_csc_valid++;

          if (!i_SimTrackIds.empty()) {
            n_csc_matched_valid++;
            //muon_matchedIds_valid[iH] =  i_SimTrackIds;
                    muon_matchedIds_valid.push_back (new uint_SimHitIdpr_pair(iH,i_SimTrackIds));
          }
        } else {
          n_csc_INVALID++;
          
          if (!i_SimTrackIds.empty()) {
            n_csc_matched_INVALID++;
            //muon_matchedIds_INVALID[iH] =  i_SimTrackIds;
                    muon_matchedIds_INVALID.push_back (new uint_SimHitIdpr_pair(iH,i_SimTrackIds));
          }
        }
          } else if (printRtS) edm::LogWarning("MuonAssociatorByHits")
        <<"*** WARNING in MuonAssociatorByHits::getMatchedIds, null dynamic_cast of a CSC TrackingRecHit !";
          
          unsigned int i_detid = (*ithit)->geographicalId().rawId();
          CSCDetId i_cscdetid = CSCDetId(i_detid);

          stringstream i_csc_detector_id;
          i_csc_detector_id << i_cscdetid;

          stringstream i_ss;
          i_ss<<"\t\t hit "<<i_compHit<<" -Muon CSC- detID = "<<i_csc_detector_id.str();

          string i_hitlog = i_ss.str();
          i_hitlog = i_hitlog + write_matched_simtracks(i_SimTrackIds);
          if (printRtS) edm::LogVerbatim("MuonAssociatorByHits") << i_hitlog;
          
          SimTrackIds.insert(SimTrackIds.end(),i_SimTrackIds.begin(),i_SimTrackIds.end());
        }       
      }  // if (cscsegment)

      else if (printRtS) edm::LogWarning("MuonAssociatorByHits")
        <<"*** WARNING in MuonAssociatorByHits::getMatchedIds, CSC TrackingRecHit is neither CSCRecHit2D nor CSCSegment ! ";
    }
      }
      
      // RPC Hits
      else if (subdet == MuonSubdetId::RPC) {
    RPCDetId rpcdetid = RPCDetId(detid);
    stringstream rpc_detector_id;
    rpc_detector_id << rpcdetid;
    if (valid_Hit) hitlog = hitlog+" -Muon RPC- detID = "+rpc_detector_id.str();      
    else hitlog = hitlog+" *** INVALID ***"+" -Muon RPC- detID = "+rpc_detector_id.str();     
    
    iH++;
    SimTrackIds = rpctruth.associateRecHit(*hitp);
    
    if (valid_Hit) {
      n_rpc_valid++;

      if (!SimTrackIds.empty()) {
        n_rpc_matched_valid++;
        //muon_matchedIds_valid[iH] = SimTrackIds;
            muon_matchedIds_valid.push_back (new uint_SimHitIdpr_pair(iH,SimTrackIds));
            
      }
    } else {
      n_rpc_INVALID++;
      
      if (!SimTrackIds.empty()) {
        n_rpc_matched_INVALID++;
        //muon_matchedIds_INVALID[iH] = SimTrackIds;
            muon_matchedIds_INVALID.push_back (new uint_SimHitIdpr_pair(iH,SimTrackIds));
      }
    }
    
      } else if (printRtS) edm::LogVerbatim("MuonAssociatorByHits")
    <<"TrackingRecHit "<<iloop<<"  *** WARNING *** Unexpected Hit from Detector = "<<det;
    }
    else continue;
    
    hitlog = hitlog + write_matched_simtracks(SimTrackIds);
    
    if (printRtS) edm::LogVerbatim("MuonAssociatorByHits") << hitlog;
    if (printRtS && dumpDT && det==DetId::Muon && subdet==MuonSubdetId::DT) {
      edm::LogVerbatim("MuonAssociatorByHits") <<wireidlog;
      for (unsigned int j=0; j<DTSimHits.size(); j++) {
    edm::LogVerbatim("MuonAssociatorByHits") <<DTSimHits[j];
      }
    }
    
  } //trackingRecHit loop
}

int MuonAssociatorByHits::getShared	(	MapOfMatchedIds &	matchedIds,
		TrackingParticleCollection::const_iterator	trpart
	)		const

Definition at line 1330 of file MuonAssociatorByHits.cc.

References newFWLiteAna::found.

Referenced by associateRecoToSimIndices(), and associateSimToRecoIndices().

                                                                                                                       {
  int nshared = 0;


  // map is indexed over the rechits of the reco::Track (no double-countings allowed)
  for (MapOfMatchedIds::const_iterator iRecH=matchedIds.begin(); iRecH!=matchedIds.end(); ++iRecH) {

    // vector of associated simhits associated to the current rechit
    std::vector<SimHitIdpr> const & SimTrackIds = (*iRecH).second;
    
    bool found = false;
    
    for (std::vector<SimHitIdpr>::const_iterator iSimH=SimTrackIds.begin(); iSimH!=SimTrackIds.end(); ++iSimH) {
      uint32_t simtrackId = iSimH->first;
      EncodedEventId evtId = iSimH->second;
      
      // look for shared hits with the given TrackingParticle (looping over component SimTracks)
      for (TrackingParticle::g4t_iterator simtrack = trpart->g4Track_begin(); simtrack !=  trpart->g4Track_end(); ++simtrack) {
    if (simtrack->trackId() == simtrackId  &&  simtrack->eventId() == evtId) {
      found = true;
      break;
    }
      }
      
      if (found) {
    nshared++;
    break;
      }
    }
  }
  
  return nshared;
}

int MuonAssociatorByHits::LayerFromDetid ( const DetId &  ) const

private

std::string MuonAssociatorByHits::write_matched_simtracks ( const std::vector< SimHitIdpr > &  SimTrackIds ) const

private

Definition at line 1364 of file MuonAssociatorByHits.cc.

References first, j, and edm::second().

                                                                                                      {
  if (SimTrackIds.empty())
    return "  *** UNMATCHED ***";

  string hitlog(" matched to SimTrack");
  
  for(size_t j=0; j<SimTrackIds.size(); j++)
  {
    char buf[64];
    snprintf(buf, 64, " Id:%i/Evt:(%i,%i) ", SimTrackIds[j].first, SimTrackIds[j].second.event(), SimTrackIds[j].second.bunchCrossing());
    hitlog += buf;
  }
  return hitlog;
}

Member Data Documentation

const bool MuonAssociatorByHits::AbsoluteNumberOfHits_muon

private

Definition at line 101 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices(), and associateSimToRecoIndices().

const bool MuonAssociatorByHits::AbsoluteNumberOfHits_track

private

Definition at line 97 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices(), and associateSimToRecoIndices().

const bool MuonAssociatorByHits::acceptOneStubMatchings

private

Definition at line 94 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices(), and associateSimToRecoIndices().

const edm::ParameterSet& MuonAssociatorByHits::conf_

private

Definition at line 114 of file MuonAssociatorByHits.h.

Referenced by associateMuons(), associateRecoToSimIndices(), associateSimToRecoIndices(), and MuonAssociatorByHits().

const bool MuonAssociatorByHits::crossingframe

private

Definition at line 111 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices(), and MuonAssociatorByHits().

const bool MuonAssociatorByHits::dumpDT

private

Definition at line 109 of file MuonAssociatorByHits.h.

const bool MuonAssociatorByHits::dumpInputCollections

private

Definition at line 110 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices().

double MuonAssociatorByHits::EfficiencyCut_muon

private

Definition at line 103 of file MuonAssociatorByHits.h.

Referenced by associateSimToRecoIndices(), and MuonAssociatorByHits().

double MuonAssociatorByHits::EfficiencyCut_track

private

Definition at line 99 of file MuonAssociatorByHits.h.

Referenced by associateSimToRecoIndices().

const bool MuonAssociatorByHits::includeZeroHitMuons

private

Definition at line 93 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices(), associateSimToRecoIndices(), and MuonAssociatorByHits().

unsigned int MuonAssociatorByHits::NHitCut_muon

private

Definition at line 102 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices(), associateSimToRecoIndices(), and MuonAssociatorByHits().

unsigned int MuonAssociatorByHits::NHitCut_track

private

Definition at line 98 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices(), and associateSimToRecoIndices().

double MuonAssociatorByHits::PurityCut_muon

private

Definition at line 104 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices(), associateSimToRecoIndices(), and MuonAssociatorByHits().

double MuonAssociatorByHits::PurityCut_track

private

Definition at line 100 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices(), and associateSimToRecoIndices().

edm::InputTag MuonAssociatorByHits::simtracksTag

private

Definition at line 112 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices().

edm::EDGetTokenT<edm::SimTrackContainer> MuonAssociatorByHits::simtracksToken_

private

Definition at line 118 of file MuonAssociatorByHits.h.

Referenced by MuonAssociatorByHits().

edm::InputTag MuonAssociatorByHits::simtracksXFTag

private

Definition at line 113 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices(), and MuonAssociatorByHits().

edm::EDGetTokenT<CrossingFrame<SimTrack> > MuonAssociatorByHits::simtracksXFToken_

private

Definition at line 116 of file MuonAssociatorByHits.h.

Referenced by MuonAssociatorByHits().

edm::EDGetTokenT<edm::SimVertexContainer> MuonAssociatorByHits::simvertsToken_

private

Definition at line 119 of file MuonAssociatorByHits.h.

Referenced by MuonAssociatorByHits().

edm::EDGetTokenT<CrossingFrame<SimVertex> > MuonAssociatorByHits::simvertsXFToken_

private

Definition at line 117 of file MuonAssociatorByHits.h.

Referenced by MuonAssociatorByHits().

const bool MuonAssociatorByHits::ThreeHitTracksAreSpecial

private

Definition at line 108 of file MuonAssociatorByHits.h.

Referenced by associateRecoToSimIndices(), and associateSimToRecoIndices().

const bool MuonAssociatorByHits::UseGrouped

private

Definition at line 106 of file MuonAssociatorByHits.h.

bool MuonAssociatorByHits::UseMuon

private

Definition at line 96 of file MuonAssociatorByHits.h.

Referenced by associateSimToRecoIndices(), and MuonAssociatorByHits().

const bool MuonAssociatorByHits::UsePixels

private

Definition at line 105 of file MuonAssociatorByHits.h.

const bool MuonAssociatorByHits::UseSplitting

private

Definition at line 107 of file MuonAssociatorByHits.h.

bool MuonAssociatorByHits::UseTracker

private

Definition at line 95 of file MuonAssociatorByHits.h.

Referenced by associateSimToRecoIndices(), and MuonAssociatorByHits().