#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.
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.
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
int	getShared (MapOfMatchedIds &matchedIds, TrackingParticleCollection::const_iterator trpart) const
	MuonAssociatorByHits (const edm::ParameterSet &)
	~MuonAssociatorByHits ()
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 (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::InputTag	simtracksXFTag
const bool	ThreeHitTracksAreSpecial
const bool	UseGrouped
bool	UseMuon
const bool	UsePixels
const bool	UseSplitting
bool	UseTracker

Detailed Description

Definition at line 26 of file MuonAssociatorByHits.h.

Member Typedef Documentation

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

Definition at line 124 of file MuonAssociatorByHits.h.

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

Definition at line 32 of file MuonAssociatorByHits.h.

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

Definition at line 73 of file MuonAssociatorByHits.h.

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

Definition at line 29 of file MuonAssociatorByHits.h.

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

Definition at line 74 of file MuonAssociatorByHits.h.

typedef std::vector<std::pair<trackingRecHit_iterator, trackingRecHit_iterator> > MuonAssociatorByHits::TrackHitsCollection [private]

Definition at line 118 of file MuonAssociatorByHits.h.

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

Definition at line 31 of file MuonAssociatorByHits.h.

Member Enumeration Documentation

enum MuonAssociatorByHits::MuonTrackType

Enumerator:

InnerTk
OuterTk
GlobalTk
Segments

Definition at line 67 of file MuonAssociatorByHits.h.

{ InnerTk, OuterTk, GlobalTk, Segments };

Constructor & Destructor Documentation

MuonAssociatorByHits::MuonAssociatorByHits ( const edm::ParameterSet & conf )

Definition at line 25 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 71 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 1388 of file MuonAssociatorByHits.cc.

References associateRecoToSimIndices(), associateSimToRecoIndices(), edm::RefVector< C, T, F >::begin(), edm::RefToBaseVector< T >::begin(), conf_, edm::RefVector< C, T, F >::end(), edm::RefToBaseVector< T >::end(), TrackerMuonHitExtractor::getMuonHits(), GlobalTk, i, errorMatrix2Lands::indices, TrackerMuonHitExtractor::init(), InnerTk, edm::RefToBase< T >::isNonnull(), n, OuterTk, edm::OwnVector< T, P >::push_back(), edm::RefVector< C, T, F >::push_back(), edm::RefToBaseVector< T >::push_back(), Segments, and edm::OwnVector< T, P >::size().

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 1376 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 76 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 98 of file MuonAssociatorByHits.cc.

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

Referenced by associateMuons(), and associateRecoToSim().

                                                                                                {

  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->trackPSimHit().size()<<" PSimHit"
        <<" (in "<<ITER->matchedHit()<<" 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("g4SimHits",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);
    
    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).trackPSimHit().size()<<" 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 450 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 473 of file MuonAssociatorByHits.cc.

References AbsoluteNumberOfHits_muon, AbsoluteNumberOfHits_track, acceptOneStubMatchings, conf_, CSC(), MuonTruth::cscBadChambers, DetId::det(), EfficiencyCut_muon, EfficiencyCut_track, getMatchedIds(), getShared(), includeZeroHitMuons, CSCBadChambers::isInBadChamber(), LayerFromDetid(), DetId::Muon, NHitCut_muon, NHitCut_track, SiStripDetId::partnerDetId(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, edm::RefVector< C, T, F >::product(), PurityCut_muon, PurityCut_track, DetId::rawId(), edm::RefVector< C, T, F >::size(), python::multivaluedict::sort(), DetId::subdetId(), SiStripDetId::TEC, ThreeHitTracksAreSpecial, SiStripDetId::TIB, SiStripDetId::TID, SiStripDetId::TOB, DetId::Tracker, UseGrouped, UseMuon, UsePixels, UseSplitting, and UseTracker.

Referenced by associateMuons(), and associateSimToReco().

                                                                                                {


  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);
    
    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 

        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 (LayerFromDetid(dId)==LayerFromDetid(dIdOK) &&
                    dId.subdetId()==dIdOK.subdetId()) newhit = false;
              //no grouped, splitting
              if (!UseGrouped && UseSplitting)
                if (LayerFromDetid(dId)==LayerFromDetid(dIdOK) &&
                    dId.subdetId()==dIdOK.subdetId() &&
                    (stripDetId==0 || stripDetId->partnerDetId()!=dIdOK.rawId()))
                  newhit = false;
              //grouped, no splitting
              if (UseGrouped && !UseSplitting)
                if (LayerFromDetid(dId)==LayerFromDetid(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();
        // Handle the case of TrackingParticles that don't have PSimHits inside, e.g. because they were made on RECOSIM only.
        if (trpart->trackPSimHit().empty()) {
            // 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).trackPSimHit().size()<<" 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).trackPSimHit().size()<<" 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 [inline, private]

Definition at line 112 of file MuonAssociatorByHits.h.

{return &*iter;}

const TrackingRecHit* MuonAssociatorByHits::getHitPtr ( trackingRecHit_iterator iter ) const [inline, private]

Definition at line 113 of file MuonAssociatorByHits.h.

{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

Definition at line 922 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(), getHLTprescales::index, TrackingRecHit::isValid(), j, DetId::Muon, DTRecSegment4D::phiSegment(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::rawId(), CSCSegment::recHits(), DTRecSegment2D::recHits(), MuonSubdetId::RPC, RPCDetId, DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, 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;
      
      if (subdet == PixelSubdetector::PixelBarrel) {
        PXBDetId pxbdetid(detid);
        detector_id << pxbdetid;
      }
      else if (subdet == PixelSubdetector::PixelEndcap) {
        PXFDetId pxfdetid(detid);
        detector_id << pxfdetid;
      }
      else if (subdet == StripSubdetector::TIB) {
        TIBDetId tibdetid(detid);
        detector_id << tibdetid;
      }
      else if (subdet == StripSubdetector::TOB) {
        TOBDetId tobdetid(detid);
        detector_id << tobdetid;
      }
      else if (subdet == StripSubdetector::TID) {
        TIDDetId tiddetid(detid);
        detector_id << tiddetid;
      }
      else if (subdet == StripSubdetector::TEC) {
        TECDetId tecdetid(detid);
        detector_id << tecdetid;
      }
      
      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 1318 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 & detId ) const [private]

Definition at line 879 of file MuonAssociatorByHits.cc.

References DetId::det(), PXFDetId::disk(), TIBDetId::layer(), TOBDetId::layer(), PXBDetId::layer(), align::tib::layerNumber(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::rawId(), DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, align::Tracker, TIDDetId::wheel(), and TECDetId::wheel().

Referenced by associateSimToRecoIndices().

{
  int layerNumber=0;
  if (detId.det() != DetId::Tracker) return layerNumber;

  unsigned int subdetId = static_cast<unsigned int>(detId.subdetId()); 
  if ( subdetId == StripSubdetector::TIB) 
    { 
      TIBDetId tibid(detId.rawId()); 
      layerNumber = tibid.layer();
    }
  else if ( subdetId ==  StripSubdetector::TOB )
    { 
      TOBDetId tobid(detId.rawId()); 
      layerNumber = tobid.layer();
    }
  else if ( subdetId ==  StripSubdetector::TID) 
    { 
      TIDDetId tidid(detId.rawId());
      layerNumber = tidid.wheel();
    }
  else if ( subdetId ==  StripSubdetector::TEC )
    { 
      TECDetId tecid(detId.rawId()); 
      layerNumber = tecid.wheel(); 
    }
  else if ( subdetId ==  PixelSubdetector::PixelBarrel ) 
    { 
      PXBDetId pxbid(detId.rawId()); 
      layerNumber = pxbid.layer();  
    }
  else if ( subdetId ==  PixelSubdetector::PixelEndcap ) 
    { 
      PXFDetId pxfid(detId.rawId()); 
      layerNumber = pxfid.disk();  
    }
  else edm::LogWarning("MuonAssociatorByHits") 
    << "*** WARNING in MuonAssociatorByHits::LayerFromDetid: Unknown Tracker subdetector: subdetId = " <<  subdetId;
  
  return layerNumber;
}

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

Definition at line 1352 of file MuonAssociatorByHits.cc.

References j.

                                                                                                      {

  string hitlog;

  if (!SimTrackIds.empty()) {
    hitlog = " matched to SimTrack";
  
    for(size_t j=0; j<SimTrackIds.size(); j++){
      stringstream trackid;  
      trackid<<SimTrackIds[j].first;
      
      stringstream evtid;    
      evtid<<SimTrackIds[j].second.event();
      
      stringstream bunchxid; 
      bunchxid<<SimTrackIds[j].second.bunchCrossing();
      
      hitlog = hitlog+" Id:"+trackid.str()+"/Evt:("+evtid.str()+","+bunchxid.str()+") ";
    }
  } else hitlog = "  *** UNMATCHED ***";

  return hitlog;
}