|
|
|
#include <SimTracker/TrackAssociation/interface/TrackAssociatorByHits.h>
Definition at line 22 of file TrackAssociatorByHits.h.
| TrackAssociatorByHits::TrackAssociatorByHits | ( | const edm::ParameterSet & | conf | ) | [explicit] |
Definition at line 38 of file TrackAssociatorByHits.cc.
00038 : 00039 conf_(conf), 00040 AbsoluteNumberOfHits(conf_.getParameter<bool>("AbsoluteNumberOfHits")), 00041 SimToRecoDenominator(conf_.getParameter<string>("SimToRecoDenominator")), 00042 quality_SimToReco(conf_.getParameter<double>("Quality_SimToReco")), 00043 purity_SimToReco(conf_.getParameter<double>("Purity_SimToReco")), 00044 cut_RecoToSim(conf_.getParameter<double>("Cut_RecoToSim")), 00045 UsePixels(conf_.getParameter<bool>("UsePixels")), 00046 UseGrouped(conf_.getParameter<bool>("UseGrouped")), 00047 UseSplitting(conf_.getParameter<bool>("UseSplitting")), 00048 ThreeHitTracksAreSpecial(conf_.getParameter<bool>("ThreeHitTracksAreSpecial")) 00049 { 00050 }
| TrackAssociatorByHits::~TrackAssociatorByHits | ( | ) |
| RecoToSimCollectionSeed TrackAssociatorByHits::associateRecoToSim | ( | edm::Handle< edm::View< TrajectorySeed > > & | seedCollectionH, | |
| edm::Handle< TrackingParticleCollection > & | TPCollectionH, | |||
| const edm::Event * | event = 0 | |||
| ) | const [virtual] |
Reimplemented from TrackAssociatorBase.
Definition at line 313 of file TrackAssociatorByHits.cc.
References AbsoluteNumberOfHits, edm::View< T >::begin(), conf_, cut_RecoToSim, edm::View< T >::end(), getShared(), edm::AssociationMap< Tag >::insert(), LogTrace, edm::AssociationMap< Tag >::post_insert(), edm::Handle< T >::product(), edm::AssociationMap< Tag >::size(), t, ThreeHitTracksAreSpecial, and TrackerHitAssociator::TrackerHitAssociator().
00315 { 00316 00317 edm::LogVerbatim("TrackAssociator") << "Starting TrackAssociatorByHits::associateRecoToSim - #seeds=" 00318 <<seedCollectionH->size()<<" #TPs="<<TPCollectionH->size(); 00319 int nshared = 0; 00320 float quality=0;//fraction or absolute number of shared hits 00321 std::vector< SimHitIdpr> SimTrackIds; 00322 std::vector< SimHitIdpr> matchedIds; 00323 RecoToSimCollectionSeed outputCollection; 00324 00325 TrackerHitAssociator * associate = new TrackerHitAssociator::TrackerHitAssociator(*e, conf_); 00326 00327 const TrackingParticleCollection tPC = *(TPCollectionH.product()); 00328 00329 const edm::View<TrajectorySeed> sC = *(seedCollectionH.product()); 00330 00331 //get the ID of the recotrack by hits 00332 int tindex=0; 00333 for (edm::View<TrajectorySeed>::const_iterator seed=sC.begin(); seed!=sC.end(); seed++, tindex++) { 00334 matchedIds.clear(); 00335 int ri=0;//valid rechits 00336 int nsimhit = seed->recHits().second-seed->recHits().first; 00337 LogTrace("TrackAssociator") << "\nNEW SEED - seed number " << tindex << " # valid=" << nsimhit; 00338 getMatchedIds<edm::OwnVector<TrackingRecHit>::const_iterator>(matchedIds, SimTrackIds, ri, seed->recHits().first, seed->recHits().second, associate ); 00339 00340 //save id for the track 00341 std::vector<SimHitIdpr> idcachev; 00342 if(!matchedIds.empty()){ 00343 00344 int tpindex =0; 00345 for (TrackingParticleCollection::const_iterator t = tPC.begin(); t != tPC.end(); ++t, ++tpindex) { 00346 LogTrace("TrackAssociator") << "TP number " << tpindex << " pdgId=" << t->pdgId() << " with number of PSimHits: " << nsimhit; 00347 idcachev.clear(); 00348 nshared = getShared(matchedIds, idcachev, t); 00349 00350 if (AbsoluteNumberOfHits) quality = static_cast<double>(nshared); 00351 else if(ri!=0) quality = (static_cast<double>(nshared)/static_cast<double>(ri)); 00352 else quality = 0; 00353 //cut on the fraction 00354 if(quality > cut_RecoToSim && !(ThreeHitTracksAreSpecial && ri==3 && nshared<3) ){ 00355 //if a track has just 3 hits we require that all 3 hits are shared 00356 outputCollection.insert(edm::RefToBase<TrajectorySeed>(seedCollectionH,tindex), 00357 std::make_pair(edm::Ref<TrackingParticleCollection>(TPCollectionH, tpindex),quality)); 00358 LogTrace("TrackAssociator") << "Seed number " << tindex << " with #hits=" << ri 00359 << "associated to TP (pdgId, nb segments, p) = " 00360 << (*t).pdgId() << " " << (*t).g4Tracks().size() 00361 << " " << (*t).momentum() <<" number " << tpindex << " with quality =" << quality; 00362 } else { 00363 LogTrace("TrackAssociator") <<"Seed number " << tindex << " with #hits=" << ri << " NOT associated with quality =" << quality; 00364 } 00365 }//TP loop 00366 } 00367 } 00368 LogTrace("TrackAssociator") << "% of Assoc Seeds=" << ((double)outputCollection.size())/((double)seedCollectionH->size()); 00369 delete associate; 00370 outputCollection.post_insert(); 00371 return outputCollection; 00372 }
| reco::RecoToSimCollection TrackAssociatorByHits::associateRecoToSim | ( | edm::Handle< edm::View< reco::Track > > & | tCH, | |
| edm::Handle< TrackingParticleCollection > & | tPCH, | |||
| const edm::Event * | event = 0 | |||
| ) | const [inline, virtual] |
compare reco to sim the handle of reco::Track and TrackingParticle collections
Reimplemented from TrackAssociatorBase.
Definition at line 40 of file TrackAssociatorByHits.h.
References TrackAssociatorBase::associateRecoToSim(), and event().
00042 { 00043 return TrackAssociatorBase::associateRecoToSim(tCH,tPCH,event); 00044 }
| RecoToSimCollection TrackAssociatorByHits::associateRecoToSim | ( | edm::RefToBaseVector< reco::Track > & | tC, | |
| edm::RefVector< TrackingParticleCollection > & | TPCollectionH, | |||
| const edm::Event * | event = 0 | |||
| ) | const [virtual] |
Association Reco To Sim with Collections.
Implements TrackAssociatorBase.
Definition at line 64 of file TrackAssociatorByHits.cc.
References AbsoluteNumberOfHits, edm::RefToBaseVector< T >::begin(), conf_, cut_RecoToSim, edm::RefToBaseVector< T >::end(), getShared(), edm::AssociationMap< Tag >::insert(), LogTrace, edm::AssociationMap< Tag >::post_insert(), edm::RefVector< C, T, F >::product(), edm::RefToBaseVector< T >::size(), edm::AssociationMap< Tag >::size(), edm::RefVector< C, T, F >::size(), t, ThreeHitTracksAreSpecial, track, and TrackerHitAssociator::TrackerHitAssociator().
00066 { 00067 00068 edm::LogVerbatim("TrackAssociator") << "Starting TrackAssociatorByHits::associateRecoToSim - #tracks="<<tC.size()<<" #TPs="<<TPCollectionH.size(); 00069 int nshared = 0; 00070 float quality=0;//fraction or absolute number of shared hits 00071 std::vector< SimHitIdpr> SimTrackIds; 00072 std::vector< SimHitIdpr> matchedIds; 00073 RecoToSimCollection outputCollection; 00074 00075 TrackerHitAssociator * associate = new TrackerHitAssociator::TrackerHitAssociator(*e, conf_); 00076 00077 TrackingParticleCollection tPC; 00078 if (TPCollectionH.size()!=0) tPC = *const_cast<TrackingParticleCollection*>(TPCollectionH.product()); 00079 00080 //get the ID of the recotrack by hits 00081 int tindex=0; 00082 for (edm::RefToBaseVector<reco::Track>::const_iterator track=tC.begin(); track!=tC.end(); track++, tindex++){ 00083 matchedIds.clear(); 00084 int ri=0;//valid rechits 00085 LogTrace("TrackAssociator") << "\nNEW TRACK - track number " << tindex <<" with pt =" << (*track)->pt() << " # valid=" << (*track)->found(); 00086 getMatchedIds<trackingRecHit_iterator>(matchedIds, SimTrackIds, ri, (*track)->recHitsBegin(), (*track)->recHitsEnd(), associate); 00087 00088 //LogTrace("TrackAssociator") << "MATCHED IDS LIST BEGIN" ; 00089 //for(size_t j=0; j<matchedIds.size(); j++){ 00090 // LogTrace("TrackAssociator") << "matchedIds[j].first=" << matchedIds[j].first; 00091 //} 00092 //LogTrace("TrackAssociator") << "MATCHED IDS LIST END" ; 00093 //LogTrace("TrackAssociator") << "#matched ids=" << matchedIds.size() << " #tps=" << tPC.size(); 00094 00095 //save id for the track 00096 std::vector<SimHitIdpr> idcachev; 00097 if(!matchedIds.empty()){ 00098 00099 int tpindex =0; 00100 for (TrackingParticleCollection::iterator t = tPC.begin(); t != tPC.end(); ++t, ++tpindex) { 00101 int nsimhit = t->trackerPSimHit_end()-t->trackerPSimHit_begin(); 00102 LogTrace("TrackAssociator") << "TP number " << tpindex << " pdgId=" << t->pdgId() << " with number of PSimHits: " << nsimhit; 00103 idcachev.clear(); 00104 nshared = getShared(matchedIds, idcachev, t); 00105 00106 if (AbsoluteNumberOfHits) quality = static_cast<double>(nshared); 00107 else if(ri!=0) quality = (static_cast<double>(nshared)/static_cast<double>(ri)); 00108 else quality = 0; 00109 //cut on the fraction 00110 //float purity = 1.0*nshared/ri; 00111 if(quality > cut_RecoToSim && !(ThreeHitTracksAreSpecial && ri==3 && nshared<3)){ 00112 //if a track has just 3 hits we require that all 3 hits are shared 00113 outputCollection.insert(tC[tindex], 00114 std::make_pair(edm::Ref<TrackingParticleCollection>(TPCollectionH, tpindex), 00115 quality)); 00116 LogTrace("TrackAssociator") << "reco::Track number " << tindex << " with #hits=" << ri <<" pt=" << (*track)->pt() 00117 << " associated to TP (pdgId, nb segments, p) = " 00118 << (*t).pdgId() << " " << (*t).g4Tracks().size() 00119 << " " << (*t).momentum() << " #hits=" << nsimhit 00120 << " TP index=" << tpindex<< " with quality =" << quality; 00121 } else { 00122 LogTrace("TrackAssociator") <<"reco::Track number " << tindex << " with #hits=" << ri 00123 << " NOT associated with quality =" << quality; 00124 } 00125 }//TP loop 00126 } 00127 } 00128 LogTrace("TrackAssociator") << "% of Assoc Tracks=" << ((double)outputCollection.size())/((double)tC.size()); 00129 delete associate; 00130 outputCollection.post_insert(); 00131 return outputCollection; 00132 }
| SimToRecoCollectionSeed TrackAssociatorByHits::associateSimToReco | ( | edm::Handle< edm::View< TrajectorySeed > > & | seedCollectionH, | |
| edm::Handle< TrackingParticleCollection > & | TPCollectionH, | |||
| const edm::Event * | event = 0 | |||
| ) | const [virtual] |
Reimplemented from TrackAssociatorBase.
Definition at line 376 of file TrackAssociatorByHits.cc.
References AbsoluteNumberOfHits, edm::View< T >::begin(), conf_, edm::View< T >::end(), getShared(), edm::AssociationMap< Tag >::insert(), LogTrace, edm::AssociationMap< Tag >::post_insert(), edm::Handle< T >::product(), quality_SimToReco, edm::AssociationMap< Tag >::size(), t, ThreeHitTracksAreSpecial, and TrackerHitAssociator::TrackerHitAssociator().
00378 { 00379 00380 edm::LogVerbatim("TrackAssociator") << "Starting TrackAssociatorByHits::associateSimToReco - #seeds=" 00381 <<seedCollectionH->size()<<" #TPs="<<TPCollectionH->size(); 00382 float quality=0;//fraction or absolute number of shared hits 00383 int nshared = 0; 00384 std::vector< SimHitIdpr> SimTrackIds; 00385 std::vector< SimHitIdpr> matchedIds; 00386 SimToRecoCollectionSeed outputCollection; 00387 00388 TrackerHitAssociator * associate = new TrackerHitAssociator::TrackerHitAssociator(*e, conf_); 00389 00390 TrackingParticleCollection tPC =*const_cast<TrackingParticleCollection*>(TPCollectionH.product()); 00391 00392 const edm::View<TrajectorySeed> sC = *(seedCollectionH.product()); 00393 00394 //get the ID of the recotrack by hits 00395 int tindex=0; 00396 for (edm::View<TrajectorySeed>::const_iterator seed=sC.begin(); seed!=sC.end(); seed++, tindex++) { 00397 int ri=0;//valid rechits 00398 LogTrace("TrackAssociator") << "\nNEW SEED - seed number " << tindex << " # valid=" << seed->recHits().second-seed->recHits().first; 00399 getMatchedIds<edm::OwnVector<TrackingRecHit>::const_iterator>(matchedIds, SimTrackIds, ri, seed->recHits().first, seed->recHits().second, associate ); 00400 00401 //save id for the track 00402 std::vector<SimHitIdpr> idcachev; 00403 if(!matchedIds.empty()){ 00404 int tpindex =0; 00405 for (TrackingParticleCollection::iterator t = tPC.begin(); t != tPC.end(); ++t, ++tpindex) { 00406 idcachev.clear(); 00407 int nsimhit = t->trackerPSimHit().size(); 00408 LogTrace("TrackAssociator") << "TP number " << tpindex << " pdgId=" << t->pdgId() << " with number of PSimHits: " << nsimhit; 00409 nshared = getShared(matchedIds, idcachev, t); 00410 00411 if (AbsoluteNumberOfHits) quality = static_cast<double>(nshared); 00412 else if(ri!=0) quality = ((double) nshared)/((double)ri); 00413 else quality = 0; 00414 //LogTrace("TrackAssociator") << "Final count: nhit(TP) = " << nsimhit 00415 //<< " nshared = " << nshared 00416 //<< " nrechit = " << ri; 00417 if(quality > quality_SimToReco && !(ThreeHitTracksAreSpecial && ri==3 && nshared<3) ){ 00418 outputCollection.insert(edm::Ref<TrackingParticleCollection>(TPCollectionH, tpindex), 00419 std::make_pair(edm::RefToBase<TrajectorySeed>(seedCollectionH,tindex), quality)); 00420 LogTrace("TrackAssociator") << "TrackingParticle number " << tpindex << " with #hits=" << nsimhit 00421 << " associated to seed number " << tindex << " with #hits=" << ri 00422 << " with hit quality =" << quality ; 00423 } 00424 else { 00425 LogTrace("TrackAssociator") << "TrackingParticle number " << tpindex << " with #hits=" << nsimhit << " NOT associated with quality =" << quality; 00426 } 00427 } 00428 } 00429 } 00430 LogTrace("TrackAssociator") << "% of Assoc TPs=" << ((double)outputCollection.size())/((double)TPCollectionH->size()); 00431 delete associate; 00432 outputCollection.post_insert(); 00433 return outputCollection; 00434 }
| reco::SimToRecoCollection TrackAssociatorByHits::associateSimToReco | ( | edm::Handle< edm::View< reco::Track > > & | tCH, | |
| edm::Handle< TrackingParticleCollection > & | tPCH, | |||
| const edm::Event * | event = 0 | |||
| ) | const [inline, virtual] |
compare reco to sim the handle of reco::Track and TrackingParticle collections
Reimplemented from TrackAssociatorBase.
Definition at line 47 of file TrackAssociatorByHits.h.
References TrackAssociatorBase::associateSimToReco(), and event().
00049 { 00050 return TrackAssociatorBase::associateSimToReco(tCH,tPCH,event); 00051 }
| SimToRecoCollection TrackAssociatorByHits::associateSimToReco | ( | edm::RefToBaseVector< reco::Track > & | tC, | |
| edm::RefVector< TrackingParticleCollection > & | TPCollectionH, | |||
| const edm::Event * | event = 0 | |||
| ) | const [virtual] |
Association Sim To Reco with Collections.
Implements TrackAssociatorBase.
Definition at line 136 of file TrackAssociatorByHits.cc.
References AbsoluteNumberOfHits, edm::RefToBaseVector< T >::begin(), conf_, edm::RefToBaseVector< T >::end(), getShared(), edm::AssociationMap< Tag >::insert(), LayerFromDetid(), LogTrace, SiStripDetId::partnerDetId(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, edm::AssociationMap< Tag >::post_insert(), edm::RefVector< C, T, F >::product(), purity_SimToReco, quality_SimToReco, DetId::rawId(), SimToRecoDenominator, edm::RefToBaseVector< T >::size(), edm::AssociationMap< Tag >::size(), edm::RefVector< C, T, F >::size(), DetId::subdetId(), t, SiStripDetId::TEC, ThreeHitTracksAreSpecial, SiStripDetId::TIB, SiStripDetId::TID, SiStripDetId::TOB, track, TrackerHitAssociator::TrackerHitAssociator(), UseGrouped, UsePixels, and UseSplitting.
00138 { 00139 edm::LogVerbatim("TrackAssociator") << "Starting TrackAssociatorByHits::associateSimToReco - #tracks="<<tC.size()<<" #TPs="<<TPCollectionH.size(); 00140 float quality=0;//fraction or absolute number of shared hits 00141 int nshared = 0; 00142 std::vector< SimHitIdpr> SimTrackIds; 00143 std::vector< SimHitIdpr> matchedIds; 00144 SimToRecoCollection outputCollection; 00145 00146 TrackerHitAssociator * associate = new TrackerHitAssociator::TrackerHitAssociator(*e, conf_); 00147 00148 TrackingParticleCollection tPC; 00149 if (TPCollectionH.size()!=0) tPC = *const_cast<TrackingParticleCollection*>(TPCollectionH.product()); 00150 00151 //for (TrackingParticleCollection::const_iterator t = tPC.begin(); t != tPC.end(); ++t) { 00152 // LogTrace("TrackAssociator") << "NEW TP DUMP"; 00153 // for (TrackingParticle::g4t_iterator g4T = t -> g4Track_begin();g4T != t -> g4Track_end(); ++g4T) { 00154 // LogTrace("TrackAssociator") << "(*g4T).trackId()=" <<(*g4T).trackId() ; 00155 // } 00156 //} 00157 00158 //get the ID of the recotrack by hits 00159 int tindex=0; 00160 for (edm::RefToBaseVector<reco::Track>::const_iterator track=tC.begin(); track!=tC.end(); track++, tindex++){ 00161 LogTrace("TrackAssociator") << "\nNEW TRACK - hits of track number " << tindex <<" with pt =" << (*track)->pt() << " # valid=" << (*track)->found(); 00162 int ri=0;//valid rechits 00163 getMatchedIds<trackingRecHit_iterator>(matchedIds, SimTrackIds, ri, (*track)->recHitsBegin(), (*track)->recHitsEnd(), associate); 00164 00165 //save id for the track 00166 std::vector<SimHitIdpr> idcachev; 00167 if(!matchedIds.empty()){ 00168 00169 int tpindex =0; 00170 for (TrackingParticleCollection::iterator t = tPC.begin(); t != tPC.end(); ++t, ++tpindex) { 00171 idcachev.clear(); 00172 int nsimhit = t->trackerPSimHit().size(); 00173 float totsimhit = 0; 00174 std::vector<PSimHit> tphits; 00175 LogTrace("TrackAssociator") << "TP number " << tpindex << " pdgId=" << t->pdgId() << " with number of PSimHits: " << nsimhit; 00176 00177 nshared = getShared(matchedIds, idcachev, t); 00178 00179 //for(std::vector<PSimHit>::const_iterator TPhit = t->trackerPSimHit_begin(); TPhit != t->trackerPSimHit_end(); TPhit++){ 00180 // unsigned int detid = TPhit->detUnitId(); 00181 // DetId detId = DetId(TPhit->detUnitId()); 00182 // LogTrace("TrackAssociator") << " hit trackId= " << TPhit->trackId() << " det ID = " << detid 00183 // << " SUBDET = " << detId.subdetId() << " layer = " << LayerFromDetid(detId); 00184 //} 00185 00186 if (nshared!=0) {//do not waste time recounting when it is not needed!!!! 00187 00188 //count the TP simhit 00189 LogTrace("TrackAssociator") << "recounting of tp hits"; 00190 for(std::vector<PSimHit>::const_iterator TPhit = t->trackerPSimHit_begin(); TPhit != t->trackerPSimHit_end(); TPhit++){ 00191 DetId dId = DetId(TPhit->detUnitId()); 00192 00193 unsigned int subdetId = static_cast<unsigned int>(dId.subdetId()); 00194 if (!UsePixels && (subdetId==PixelSubdetector::PixelBarrel || subdetId==PixelSubdetector::PixelEndcap) ) 00195 continue; 00196 00197 //unsigned int dRawId = dId.rawId(); 00198 SiStripDetId* stripDetId = 0; 00199 if (subdetId==SiStripDetId::TIB||subdetId==SiStripDetId::TOB|| 00200 subdetId==SiStripDetId::TID||subdetId==SiStripDetId::TEC) 00201 stripDetId= new SiStripDetId(dId); 00202 LogTrace("TrackAssociator") << "consider hit SUBDET = " << subdetId 00203 << " layer = " << LayerFromDetid(dId) 00204 << " id = " << dId.rawId(); 00205 bool newhit = true; 00206 for(std::vector<PSimHit>::const_iterator TPhitOK = tphits.begin(); TPhitOK != tphits.end(); TPhitOK++){ 00207 DetId dIdOK = DetId(TPhitOK->detUnitId()); 00208 //unsigned int dRawIdOK = dIdOK.rawId(); 00209 LogTrace("TrackAssociator") << "\t\tcompare with SUBDET = " << dIdOK.subdetId() 00210 << " layer = " << LayerFromDetid(dIdOK) 00211 << " id = " << dIdOK.rawId(); 00212 //no grouped, no splitting 00213 if (!UseGrouped && !UseSplitting) 00214 if (LayerFromDetid(dId)==LayerFromDetid(dIdOK) && 00215 dId.subdetId()==dIdOK.subdetId()) newhit = false; 00216 //no grouped, splitting 00217 if (!UseGrouped && UseSplitting) 00218 if (LayerFromDetid(dId)==LayerFromDetid(dIdOK) && 00219 dId.subdetId()==dIdOK.subdetId() && 00220 (stripDetId==0 || stripDetId->partnerDetId()!=dIdOK.rawId())) 00221 newhit = false; 00222 //grouped, no splitting 00223 if (UseGrouped && !UseSplitting) 00224 if (LayerFromDetid(dId)==LayerFromDetid(dIdOK) && 00225 dId.subdetId()==dIdOK.subdetId() && 00226 stripDetId!=0 && stripDetId->partnerDetId()==dIdOK.rawId()) 00227 newhit = false; 00228 //grouped, splitting 00229 if (UseGrouped && UseSplitting) 00230 newhit = true; 00231 } 00232 if (newhit) { 00233 LogTrace("TrackAssociator") << "\t\tok"; 00234 tphits.push_back(*TPhit); 00235 } 00236 else LogTrace("TrackAssociator") << "\t\tno"; 00237 delete stripDetId; 00238 } 00239 totsimhit = tphits.size(); 00240 } 00241 00242 if (AbsoluteNumberOfHits) quality = static_cast<double>(nshared); 00243 else if(SimToRecoDenominator == "sim" && totsimhit!=0) quality = ((double) nshared)/((double)totsimhit); 00244 else if(SimToRecoDenominator == "reco" && ri!=0) quality = ((double) nshared)/((double)ri); 00245 else quality = 0; 00246 //LogTrace("TrackAssociator") << "Final count: nhit(TP) = " << nsimhit << " re-counted = " << totsimhit 00247 //<< " nshared = " << nshared << " nrechit = " << ri; 00248 00249 float purity = 1.0*nshared/ri; 00250 if (quality>quality_SimToReco && !(ThreeHitTracksAreSpecial && totsimhit==3 && nshared<3) && (AbsoluteNumberOfHits||(purity>purity_SimToReco))) { 00251 //if a track has just 3 hits we require that all 3 hits are shared 00252 outputCollection.insert(edm::Ref<TrackingParticleCollection>(TPCollectionH, tpindex), 00253 std::make_pair(tC[tindex],quality)); 00254 LogTrace("TrackAssociator") << "TrackingParticle number " << tpindex << " with #hits=" << nsimhit 00255 << " re-counted = " << totsimhit << " nshared = " << nshared 00256 << " associated to track number " << tindex << " with pt=" << (*track)->pt() 00257 << " with hit quality =" << quality ; 00258 } else { 00259 LogTrace("TrackAssociator") << "TrackingParticle number " << tpindex << " with #hits=" << nsimhit 00260 << " re-counted = " << totsimhit << " nshared = " << nshared 00261 << " NOT associated with quality =" << quality; 00262 } 00263 } 00264 } 00265 } 00266 LogTrace("TrackAssociator") << "% of Assoc TPs=" << ((double)outputCollection.size())/((double)TPCollectionH.size()); 00267 delete associate; 00268 outputCollection.post_insert(); 00269 return outputCollection; 00270 }
| const TrackingRecHit* TrackAssociatorByHits::getHitPtr | ( | trackingRecHit_iterator | iter | ) | const [inline, private] |
| const TrackingRecHit* TrackAssociatorByHits::getHitPtr | ( | edm::OwnVector< TrackingRecHit >::const_iterator | iter | ) | const [inline, private] |
| void TrackAssociatorByHits::getMatchedIds | ( | std::vector< SimHitIdpr > & | matchedIds, | |
| std::vector< SimHitIdpr > & | SimTrackIds, | |||
| int & | ri, | |||
| iter | begin, | |||
| iter | end, | |||
| TrackerHitAssociator * | associate | |||
| ) | const [inline] |
Definition at line 437 of file TrackAssociatorByHits.cc.
References TrackerHitAssociator::associateHitId(), TrackingRecHit::geographicalId(), getHitPtr(), TrackingRecHit::isValid(), it, iter, j, LogTrace, and DetId::rawId().
00442 { 00443 matchedIds.clear(); 00444 ri=0;//valid rechits 00445 for (iter it = begin; it != end; it++){ 00446 if (getHitPtr(it)->isValid()){ 00447 ri++; 00448 uint32_t t_detID= getHitPtr(it)->geographicalId().rawId(); 00449 SimTrackIds.clear(); 00450 SimTrackIds = associate->associateHitId(*getHitPtr(it)); 00451 //save all the id of matched simtracks 00452 if(!SimTrackIds.empty()){ 00453 for(size_t j=0; j<SimTrackIds.size(); j++){ 00454 LogTrace("TrackAssociator") << " hit # " << ri << " valid=" << getHitPtr(it)->isValid() 00455 << " det id = " << t_detID << " SimId " << SimTrackIds[j].first 00456 << " evt=" << SimTrackIds[j].second.event() 00457 << " bc=" << SimTrackIds[j].second.bunchCrossing(); 00458 matchedIds.push_back(SimTrackIds[j]); 00459 } 00460 } 00462 //****to be used when the crossing frame is in the event and with flag TrackAssociatorByHitsESProducer.associateRecoTracks = false 00463 //std::vector<PSimHit> tmpSimHits = associate->associateHit(*getHitPtr(it)); 00465 //for(size_t j=0; j<tmpSimHits.size(); j++) { 00466 // LogTrace("TrackAssociator") << " hit # " << ri << " valid=" << getHitPtr(it)->isValid() 00467 // << " det id = " << t_detID << " SimId " << SimTrackIds[j].first 00468 // << " evt=" << SimTrackIds[j].second.event() 00469 // << " bc=" << SimTrackIds[j].second.bunchCrossing() 00470 // << " process=" << tmpSimHits[j].processType() 00471 // << " eloss=" << tmpSimHits[j].energyLoss() 00472 // << " part type=" << tmpSimHits[j].particleType() 00473 // << " pabs=" << tmpSimHits[j].pabs() 00474 // << " trId=" << tmpSimHits[j].trackId(); 00475 //} 00477 }else{ 00478 LogTrace("TrackAssociator") <<"\t\t Invalid Hit On "<<getHitPtr(it)->geographicalId().rawId(); 00479 } 00480 }//trackingRecHit loop 00481 }
| int TrackAssociatorByHits::getShared | ( | std::vector< SimHitIdpr > & | matchedIds, | |
| std::vector< SimHitIdpr > & | idcachev, | |||
| TrackingParticleCollection::const_iterator | t | |||
| ) | const |
Definition at line 484 of file TrackAssociatorByHits.cc.
References count, find(), j, and LogTrace.
Referenced by associateRecoToSim(), and associateSimToReco().
00486 { 00487 int nshared = 0; 00488 for(size_t j=0; j<matchedIds.size(); j++){ 00489 //LogTrace("TrackAssociator") << "now matchedId=" << matchedIds[j].first; 00490 if(find(idcachev.begin(), idcachev.end(),matchedIds[j]) == idcachev.end() ){ 00491 //only the first time we see this ID 00492 idcachev.push_back(matchedIds[j]); 00493 00494 for (TrackingParticle::g4t_iterator g4T = t -> g4Track_begin(); g4T != t -> g4Track_end(); ++g4T) { 00495 LogTrace("TrackAssociator") << " TP (ID, Ev, BC) = " << (*g4T).trackId() 00496 << ", " << t->eventId().event() << ", "<< t->eventId().bunchCrossing() 00497 << " Match(ID, Ev, BC) = " << matchedIds[j].first 00498 << ", " << matchedIds[j].second.event() << ", " 00499 << matchedIds[j].second.bunchCrossing() ; 00500 //<< "\t G4 Track Momentum " << (*g4T).momentum() 00501 //<< " \t reco Track Momentum " << track->momentum(); 00502 if((*g4T).trackId() == matchedIds[j].first && t->eventId() == matchedIds[j].second){ 00503 int countedhits = std::count(matchedIds.begin(), matchedIds.end(), matchedIds[j]); 00504 nshared += countedhits; 00505 00506 LogTrace("TrackAssociator") << "hits shared by this segment : " << countedhits; 00507 LogTrace("TrackAssociator") << "hits shared so far : " << nshared; 00508 } 00509 }//g4Tracks loop 00510 } 00511 } 00512 return nshared; 00513 }
Definition at line 272 of file TrackAssociatorByHits.cc.
References align::tpb::layerNumber(), LogTrace, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::rawId(), DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, and StripSubdetector::TOB.
Referenced by associateSimToReco().
00273 { 00274 int layerNumber=0; 00275 unsigned int subdetId = static_cast<unsigned int>(detId.subdetId()); 00276 if ( subdetId == StripSubdetector::TIB) 00277 { 00278 TIBDetId tibid(detId.rawId()); 00279 layerNumber = tibid.layer(); 00280 } 00281 else if ( subdetId == StripSubdetector::TOB ) 00282 { 00283 TOBDetId tobid(detId.rawId()); 00284 layerNumber = tobid.layer(); 00285 } 00286 else if ( subdetId == StripSubdetector::TID) 00287 { 00288 TIDDetId tidid(detId.rawId()); 00289 layerNumber = tidid.wheel(); 00290 } 00291 else if ( subdetId == StripSubdetector::TEC ) 00292 { 00293 TECDetId tecid(detId.rawId()); 00294 layerNumber = tecid.wheel(); 00295 } 00296 else if ( subdetId == PixelSubdetector::PixelBarrel ) 00297 { 00298 PXBDetId pxbid(detId.rawId()); 00299 layerNumber = pxbid.layer(); 00300 } 00301 else if ( subdetId == PixelSubdetector::PixelEndcap ) 00302 { 00303 PXFDetId pxfid(detId.rawId()); 00304 layerNumber = pxfid.disk(); 00305 } 00306 else 00307 LogTrace("TrackAssociator") << "Unknown subdetid: " << subdetId; 00308 00309 return layerNumber; 00310 }
const bool TrackAssociatorByHits::AbsoluteNumberOfHits [private] |
Definition at line 76 of file TrackAssociatorByHits.h.
Referenced by associateRecoToSim(), and associateSimToReco().
const edm::ParameterSet& TrackAssociatorByHits::conf_ [private] |
Definition at line 75 of file TrackAssociatorByHits.h.
Referenced by associateRecoToSim(), and associateSimToReco().
const double TrackAssociatorByHits::cut_RecoToSim [private] |
const double TrackAssociatorByHits::purity_SimToReco [private] |
const double TrackAssociatorByHits::quality_SimToReco [private] |
const std::string TrackAssociatorByHits::SimToRecoDenominator [private] |
const bool TrackAssociatorByHits::ThreeHitTracksAreSpecial [private] |
Definition at line 84 of file TrackAssociatorByHits.h.
Referenced by associateRecoToSim(), and associateSimToReco().
const bool TrackAssociatorByHits::UseGrouped [private] |
const bool TrackAssociatorByHits::UsePixels [private] |
const bool TrackAssociatorByHits::UseSplitting [private] |
1.5.4