TrackerHitAssociator Class Reference

#include <TrackerHitAssociator.h>

Classes
struct	Config

Public Types
typedef std::map < simHitCollectionID, std::vector< PSimHit > >	simhit_collectionMap
typedef std::map< unsigned int, std::vector< PSimHit > >	simhit_map
typedef std::pair < simHitCollectionID, unsigned int >	simhitAddr
typedef std::pair< unsigned int, unsigned int >	simHitCollectionID

Public Member Functions
void	associateCluster (const SiStripCluster *clust, const DetId &detid, std::vector< SimHitIdpr > &simtrackid, std::vector< PSimHit > &simhit) const
std::vector< SimHitIdpr >	associateFastRecHit (const FastTrackerRecHit *rechit) const
std::vector< PSimHit >	associateHit (const TrackingRecHit &thit) const
std::vector< SimHitIdpr >	associateHitId (const TrackingRecHit &thit) const
void	associateHitId (const TrackingRecHit &thit, std::vector< SimHitIdpr > &simhitid, std::vector< simhitAddr > *simhitCFPos=0) const
std::vector< SimHitIdpr >	associateMatchedRecHit (const SiStripMatchedRecHit2D *matchedrechit, std::vector< simhitAddr > *simhitCFPos=0) const
std::vector< PSimHit >	associateMultiRecHit (const SiTrackerMultiRecHit *multirechit) const
std::vector< SimHitIdpr >	associateMultiRecHitId (const SiTrackerMultiRecHit *multirechit, std::vector< simhitAddr > *simhitCFPos=0) const
void	associatePixelRecHit (const SiPixelRecHit *pixelrechit, std::vector< SimHitIdpr > &simhitid, std::vector< simhitAddr > *simhitCFPos=0) const
std::vector< SimHitIdpr >	associateProjectedRecHit (const ProjectedSiStripRecHit2D *projectedrechit, std::vector< simhitAddr > *simhitCFPos=0) const
void	associateSimpleRecHitCluster (const SiStripCluster *clust, const DetId &detid, std::vector< SimHitIdpr > &simtrackid, std::vector< simhitAddr > *simhitCFPos=0) const
template<typename T >
void	associateSiStripRecHit (const T *simplerechit, std::vector< SimHitIdpr > &simtrackid, std::vector< simhitAddr > *simhitCFPos=0) const
	TrackerHitAssociator (const edm::Event &e, const Config &config)
virtual	~TrackerHitAssociator ()

Public Attributes
simhit_collectionMap	SimHitCollMap
simhit_map	SimHitMap

Private Types
typedef std::vector< std::string >	vstring

Private Member Functions
void	makeMaps (const edm::Event &theEvent, const Config &config)

Private Attributes
bool	assocHitbySimTrack_
bool	doPixel_
bool	doStrip_
bool	doTrackAssoc_
edm::Handle< edm::DetSetVector < PixelDigiSimLink > >	pixeldigisimlink
edm::Handle< edm::DetSetVector < StripDigiSimLink > >	stripdigisimlink

Detailed Description

Definition at line 54 of file TrackerHitAssociator.h.

Member Typedef Documentation

typedef std::map<simHitCollectionID, std::vector<PSimHit> > TrackerHitAssociator::simhit_collectionMap

Definition at line 103 of file TrackerHitAssociator.h.

typedef std::map<unsigned int, std::vector<PSimHit> > TrackerHitAssociator::simhit_map

Definition at line 101 of file TrackerHitAssociator.h.

typedef std::pair<simHitCollectionID, unsigned int> TrackerHitAssociator::simhitAddr

Definition at line 74 of file TrackerHitAssociator.h.

typedef std::pair<unsigned int, unsigned int> TrackerHitAssociator::simHitCollectionID

Definition at line 73 of file TrackerHitAssociator.h.

typedef std::vector<std::string> TrackerHitAssociator::vstring

private

Definition at line 107 of file TrackerHitAssociator.h.

Constructor & Destructor Documentation

TrackerHitAssociator::TrackerHitAssociator	(	const edm::Event &	e,
		const Config &	config
	)

Definition at line 83 of file TrackerHitAssociator.cc.

References doPixel_, doStrip_, doTrackAssoc_, edm::Event::getByToken(), makeMaps(), pixeldigisimlink, TrackerHitAssociator::Config::pixelToken_, stripdigisimlink, and TrackerHitAssociator::Config::stripToken_.

                                                                                                      :
  doPixel_(config.doPixel_),
  doStrip_(config.doStrip_),
  doTrackAssoc_(config.doTrackAssoc_),
  assocHitbySimTrack_(config.assocHitbySimTrack_) {
  //if track association there is no need to access the input collections
  if(!doTrackAssoc_) {
    makeMaps(e, config);
  }

  if(doStrip_) e.getByToken(config.stripToken_, stripdigisimlink);
  if(doPixel_) e.getByToken(config.pixelToken_, pixeldigisimlink);
}

virtual TrackerHitAssociator::~TrackerHitAssociator ( )

inlinevirtual

Definition at line 71 of file TrackerHitAssociator.h.

{}

Member Function Documentation

void TrackerHitAssociator::associateCluster	(	const SiStripCluster *	clust,
		const DetId &	detid,
		std::vector< SimHitIdpr > &	simtrackid,
		std::vector< PSimHit > &	simhit
	)		const

Definition at line 390 of file TrackerHitAssociator.cc.

References associateSimpleRecHitCluster(), and SimHitCollMap.

                                                          {
  std::vector<simhitAddr> simhitCFPos;
  associateSimpleRecHitCluster(clust, detid, simtrackid, &simhitCFPos);

  for(auto const& theSimHitAddr : simhitCFPos ) {
    simHitCollectionID theSimHitCollID = theSimHitAddr.first;
    simhit_collectionMap::const_iterator it = SimHitCollMap.find(theSimHitCollID);
    if (it!= SimHitCollMap.end()) {
      unsigned int theSimHitIndex = theSimHitAddr.second;
      if (theSimHitIndex < (it->second).size()) simhit.push_back((it->second)[theSimHitIndex]);
      // const PSimHit& theSimHit = (it->second)[theSimHitIndex];
      // std::cout << "For cluster, simHit detId =  " << theSimHit.detUnitId() << " address = (" << (theSimHitAddr.first).first
      //        << ", " << (theSimHitAddr.first).second << ", " << theSimHitIndex
      //        << "), process = " << theSimHit.processType() << " (" << theSimHit.eventId().bunchCrossing()
      //        << ", " << theSimHit.eventId().event() << ", " << theSimHit.trackId() << ")" << std::endl;
    }
  }
}

std::vector< SimHitIdpr > TrackerHitAssociator::associateFastRecHit

(

const FastTrackerRecHit *

rechit

)

const

Definition at line 632 of file TrackerHitAssociator.cc.

References cmsHarvester::index, FastTrackerRecHit::nSimTrackIds(), FastTrackerRecHit::simTrackEventId(), and FastTrackerRecHit::simTrackId().

Referenced by associateHitId().

{
  vector<SimHitIdpr> simtrackid;
  simtrackid.clear();
  for(size_t index =0, indexEnd = rechit->nSimTrackIds();index<indexEnd;++index){
      SimHitIdpr currentId(rechit->simTrackId(index), EncodedEventId(rechit->simTrackEventId(index)));
      simtrackid.push_back(currentId);
  }
  return simtrackid;
}

std::vector< PSimHit > TrackerHitAssociator::associateHit

(

const TrackingRecHit &

thit

)

const

Check if it's the gluedDet.

Definition at line 175 of file TrackerHitAssociator.cc.

References assocHitbySimTrack_, associateHitId(), associateMultiRecHit(), begin, doTrackAssoc_, end, PSimHit::eventId(), plotBeamSpotDB::first, TrackingRecHit::geographicalId(), DetId::rawId(), mps_fire::result, edm::second(), SimHitCollMap, SimHitMap, and PSimHit::trackId().

Referenced by CkfDebugger::analyseCompatibleMeasurements(), SiPixelRecHitsValid::analyze(), TestHits::analyze(), TestSmoothHits::analyze(), TestTrackHits::analyze(), TestOutliers::analyze(), SiPixelErrorEstimation::analyze(), SiPixelTrackingRecHitsValid::analyze(), CkfDebugger::associated(), associateMultiRecHit(), CkfDebugger::correctTrajectory(), GlobalRecHitsAnalyzer::fillTrk(), GlobalRecHitsProducer::fillTrk(), spr::matchedSimTrack(), CkfDebugger::nextCorrectHits(), SiStripRecHitsValid::rechitanalysis(), SiStripTrackingRecHitsValid::rechitanalysis(), SiStripRecHitsValid::rechitanalysis_matched(), SiStripTrackingRecHitsValid::rechitanalysis_matched(), and CkfDebugger::testSeed().

{

  if (const SiTrackerMultiRecHit * rechit = dynamic_cast<const SiTrackerMultiRecHit *>(&thit)){
    return associateMultiRecHit(rechit);
  }

  //vector with the matched SimHit
  std::vector<PSimHit> result;

  if(doTrackAssoc_) return result;

  // Vectors to contain lists of matched simTracks, simHits
  std::vector<SimHitIdpr> simtrackid;
  std::vector<simhitAddr> simhitCFPos;
  
  //get the Detector type of the rechit
  DetId detid=  thit.geographicalId();
  uint32_t detID = detid.rawId();

  // Get the vectors of simtrackIDs and simHit indices associated with this rechit
  associateHitId(thit, simtrackid, &simhitCFPos);
  // std::cout << "recHit subdet, detID = " << detid.subdetId() << ", " << detID << ", (bnch, evt, trk) = ";
  // for (size_t i=0; i<simtrackid.size(); ++i)
  //   std::cout << ", (" << simtrackid[i].second.bunchCrossing() << ", "
  //          << simtrackid[i].second.event() << ", " << simtrackid[i].first << ")";
  // std::cout << std::endl; 

  // Get the vector of simHits associated with this rechit

  if (!assocHitbySimTrack_ && simhitCFPos.size() > 0) {
    // We use the indices to the simHit collections taken
    //  from the DigiSimLinks and returned in simhitCFPos.
    //  simhitCFPos[i] contains the full address of the ith simhit:
    //   <collection, index> = <<subdet, tofBin>, index>

    //check if the recHit is a SiStripMatchedRecHit2D
    if(dynamic_cast<const SiStripMatchedRecHit2D *>(&thit)) {
      for(auto const& theSimHitAddr : simhitCFPos) {
        simHitCollectionID theSimHitCollID = theSimHitAddr.first;
        simhit_collectionMap::const_iterator it = SimHitCollMap.find(theSimHitCollID);
        if (it!= SimHitCollMap.end()) {
          unsigned int theSimHitIndex = theSimHitAddr.second;
          if (theSimHitIndex < (it->second).size()) {
            const PSimHit& theSimHit = (it->second)[theSimHitIndex];
            // Try to remove ghosts by requiring a match to the simTrack also
            unsigned int simHitid = theSimHit.trackId();
            EncodedEventId simHiteid = theSimHit.eventId();
            for(auto const& id : simtrackid) {
              if(simHitid == id.first && simHiteid == id.second) {
                result.push_back(theSimHit);
              }
            } 
          // std::cout << "by CFpos, simHit detId =  " << theSimHit.detUnitId() << " address = (" << (theSimHitAddr.first).first
      //        << ", " << (theSimHitAddr.first).second << ", " << theSimHitIndex
      //        << "), process = " << theSimHit.processType() << " (" << theSimHit.eventId().bunchCrossing()
      //        << ", " << theSimHit.eventId().event() << ", " << theSimHit.trackId() << ")" << std::endl;
      }
        }
      }
    } else {
      for(auto const& theSimHitAddr : simhitCFPos) {
        simHitCollectionID theSimHitCollID = theSimHitAddr.first;
        simhit_collectionMap::const_iterator it = SimHitCollMap.find(theSimHitCollID);
        if (it!= SimHitCollMap.end()) {
          unsigned int theSimHitIndex = theSimHitAddr.second;
          if (theSimHitIndex < (it->second).size()) {
             result.push_back((it->second)[theSimHitIndex]);
          // std::cout << "by CFpos, simHit detId =  " << theSimHit.detUnitId() << " address = (" << (theSimHitAddr.first).first
      //        << ", " << (theSimHitAddr.first).second << ", " << theSimHitIndex
      //        << "), process = " << theSimHit.processType() << " (" << theSimHit.eventId().bunchCrossing()
      //        << ", " << theSimHit.eventId().event() << ", " << theSimHit.trackId() << ")" << std::endl;
          }
        }
      }
    }
    return result;
  }

  // Get the SimHit from the trackid instead
  std::map<unsigned int, std::vector<PSimHit> >::const_iterator it = SimHitMap.find(detID);
  if (it!= SimHitMap.end()) {
    vector<PSimHit>::const_iterator simHitIter = (it->second).begin();
    vector<PSimHit>::const_iterator simHitIterEnd = (it->second).end();
    for (;simHitIter != simHitIterEnd; ++simHitIter) {
      const PSimHit& ihit = *simHitIter;
      unsigned int simHitid = ihit.trackId();
      EncodedEventId simHiteid = ihit.eventId();
      // std::cout << "by simTk, simHit, process = " << ihit.processType() << " (" << ihit.eventId().bunchCrossing()
      //        << ", " << ihit.eventId().event() << ", " << ihit.trackId() << ")";
      for(auto id : simtrackid) {
        if(simHitid == id.first && simHiteid == id.second) {
//          cout << "Associator ---> ID" << ihit.trackId() << " Simhit x= " << ihit.localPosition().x() 
//                   << " y= " <<  ihit.localPosition().y() << " z= " <<  ihit.localPosition().x() << endl;
      // std::cout << " matches";
      result.push_back(ihit);
      break;
        }
      } 
      // std::cout << std::endl;
    }

  }else{

    std::map<unsigned int, std::vector<PSimHit> >::const_iterator itrphi =
      SimHitMap.find(detID+2);  //iterator to the simhit in the rphi module
    std::map<unsigned int, std::vector<PSimHit> >::const_iterator itster = 
      SimHitMap.find(detID+1);//iterator to the simhit in the stereo module
    if (itrphi!= SimHitMap.end()&&itster!=SimHitMap.end()) {
      std::vector<PSimHit> simHitVector = itrphi->second;
      simHitVector.insert(simHitVector.end(),(itster->second).begin(),(itster->second).end());
      vector<PSimHit>::const_iterator simHitIter = simHitVector.begin();
      vector<PSimHit>::const_iterator simHitIterEnd = simHitVector.end();
      for (;simHitIter != simHitIterEnd; ++simHitIter) {
    const PSimHit& ihit = *simHitIter;
    unsigned int simHitid = ihit.trackId();
    EncodedEventId simHiteid = ihit.eventId();
    for(auto const& id : simtrackid) {
      if(simHitid == id.first && simHiteid == id.second) { 
//    if(simHitid == simtrackid[i].first && simHiteid.bunchCrossing() == simtrackid[i].second.bunchCrossing()) {
        //    cout << "GluedDet Associator ---> ID" << ihit.trackId() << " Simhit x= " << ihit.localPosition().x() 
            //         << " y= " <<  ihit.localPosition().y() << " z= " <<  ihit.localPosition().x() << endl; 
        result.push_back(ihit);
        break;
      }
    }
      }
    }
  }

  return result;
}

std::vector< SimHitIdpr > TrackerHitAssociator::associateHitId

(

const TrackingRecHit &

thit

)

const

Definition at line 309 of file TrackerHitAssociator.cc.

Referenced by TestTrackHits::analyze(), TestOutliers::analyze(), associateHit(), associateMultiRecHitId(), CkfDebugger::assocTrackId(), CkfDebugger::correctMeas(), CkfDebugger::correctTrajectory(), QuickTrackAssociatorByHitsImpl::getAllSimTrackIdentifiers(), TrackAssociatorByHitsImpl::getDoubleCount(), QuickTrackAssociatorByHitsImpl::getDoubleCount(), MuonAssociatorByHitsHelper::getMatchedIds(), and TrackAssociatorByHitsImpl::getMatchedIds().

{
  std::vector< SimHitIdpr > simhitid;
  associateHitId(thit, simhitid);
  return simhitid;
}

void TrackerHitAssociator::associateHitId	(	const TrackingRecHit &	thit,
		std::vector< SimHitIdpr > &	simhitid,
		std::vector< simhitAddr > *	simhitCFPos = 0
	)		const

Definition at line 316 of file TrackerHitAssociator.cc.

References associateFastRecHit(), associateMatchedRecHit(), associateMultiRecHitId(), associatePixelRecHit(), associateProjectedRecHit(), associateSiStripRecHit(), TrackingRecHit::geographicalId(), trackerHitRTTI::isFast(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, and StripSubdetector::TOB.

{

    simtkid.clear();
    
  //get the Detector type of the rechit
    DetId detid=  thit.geographicalId();
    if (const SiTrackerMultiRecHit * rechit = dynamic_cast<const SiTrackerMultiRecHit *>(&thit)){
       simtkid=associateMultiRecHitId(rechit, simhitCFPos);
    }
    
  //cout << "Associator ---> get Detid " << detID << endl;
  //check we are in the strip tracker
    if(detid.subdetId() == StripSubdetector::TIB ||
       detid.subdetId() == StripSubdetector::TOB || 
       detid.subdetId() == StripSubdetector::TID ||
       detid.subdetId() == StripSubdetector::TEC) 
      {
    //check if it is a simple SiStripRecHit2D
    if(const SiStripRecHit2D * rechit = 
       dynamic_cast<const SiStripRecHit2D *>(&thit))
      {   
        associateSiStripRecHit(rechit, simtkid, simhitCFPos);
      }
    //check if it is a SiStripRecHit1D
    else  if(const SiStripRecHit1D * rechit = 
         dynamic_cast<const SiStripRecHit1D *>(&thit))
      {   
        associateSiStripRecHit(rechit, simtkid, simhitCFPos);
      }
    //check if it is a SiStripMatchedRecHit2D
    else  if(const SiStripMatchedRecHit2D * rechit = 
         dynamic_cast<const SiStripMatchedRecHit2D *>(&thit))
      {   
        simtkid = associateMatchedRecHit(rechit, simhitCFPos);
      }
    //check if it is a  ProjectedSiStripRecHit2D
    else if(const ProjectedSiStripRecHit2D * rechit = 
        dynamic_cast<const ProjectedSiStripRecHit2D *>(&thit))
      {   
        simtkid = associateProjectedRecHit(rechit, simhitCFPos);
      }
    else{
      //std::cout << "associate to invalid" << std::endl;
      //throw cms::Exception("Unknown RecHit Type") << "TrackerHitAssociator failed second casting of " << typeid(thit).name() << " type ";
    }
      }
    //check we are in the pixel tracker
    else if( (unsigned int)(detid.subdetId()) == PixelSubdetector::PixelBarrel || 
         (unsigned int)(detid.subdetId()) == PixelSubdetector::PixelEndcap) 
      {
      if(const SiPixelRecHit * rechit = dynamic_cast<const SiPixelRecHit *>(&thit))
    {     
      associatePixelRecHit(rechit, simtkid, simhitCFPos);
    }
      }
    //check if these are GSRecHits (from FastSim)
    if(trackerHitRTTI::isFast(thit))
      {
      simtkid = associateFastRecHit(static_cast<const FastTrackerRecHit *>(&thit));
      }
}

std::vector< SimHitIdpr > TrackerHitAssociator::associateMatchedRecHit	(	const SiStripMatchedRecHit2D *	matchedrechit,
		std::vector< simhitAddr > *	simhitCFPos = 0
	)		const

Definition at line 494 of file TrackerHitAssociator.cc.

References associateSiStripRecHit(), spr::find(), SiStripMatchedRecHit2D::monoHit(), and SiStripMatchedRecHit2D::stereoHit().

Referenced by associateHitId().

{
  std::vector<SimHitIdpr> matched_mono;
  std::vector<SimHitIdpr> matched_st;

  const SiStripRecHit2D mono = matchedrechit->monoHit();
  const SiStripRecHit2D st = matchedrechit->stereoHit();
  //associate the two simple hits separately
  associateSiStripRecHit(&mono, matched_mono, simhitCFPos);
  associateSiStripRecHit(&st, matched_st, simhitCFPos);
  
  //save in a vector all the simtrack-id's that are common to mono and stereo hits
  std::vector<SimHitIdpr> simtrackid;
  if(!(matched_mono.empty() || matched_st.empty())){
    std::vector<SimHitIdpr> idcachev;
    for(auto const& mhit: matched_mono){
      //save only once the ID
      if(find(idcachev.begin(), idcachev.end(), mhit) == idcachev.end()) {
    idcachev.push_back(mhit);
    //save if the stereoID matched the monoID
    if(find(matched_st.begin(), matched_st.end(), mhit) != matched_st.end()) {
      simtrackid.push_back(mhit);
    }
      }
    }
  }
  
  return simtrackid;
}

std::vector< PSimHit > TrackerHitAssociator::associateMultiRecHit

(

const SiTrackerMultiRecHit *

multirechit

)

const

Definition at line 606 of file TrackerHitAssociator.cc.

References associateHit(), i, SiTrackerMultiRecHit::recHits(), findQualityFiles::size, SiTrackerMultiRecHit::weight(), and SiTrackerMultiRecHit::weights().

Referenced by associateHit().

                                                                                                           {
  std::vector<const TrackingRecHit*> componenthits = multirechit->recHits();
  //        std::vector<PSimHit> assimhits;
  int size=multirechit->weights().size(), idmostprobable=0;
  
  for (int i=0; i<size; ++i){
    if(multirechit->weight(i)>multirechit->weight(idmostprobable)) idmostprobable=i;
  }
  
  return associateHit(*componenthits[idmostprobable]);
}

std::vector< SimHitIdpr > TrackerHitAssociator::associateMultiRecHitId	(	const SiTrackerMultiRecHit *	multirechit,
		std::vector< simhitAddr > *	simhitCFPos = 0
	)		const

Definition at line 618 of file TrackerHitAssociator.cc.

References associateHitId(), i, SiTrackerMultiRecHit::recHits(), findQualityFiles::size, SiTrackerMultiRecHit::weight(), and SiTrackerMultiRecHit::weights().

Referenced by associateHitId().

                                                                                                                                                    {
  std::vector<const TrackingRecHit*> componenthits = multirechit->recHits();
  int size=multirechit->weights().size(), idmostprobable=0;
  
  for (int i=0; i<size; ++i){
    if(multirechit->weight(i)>multirechit->weight(idmostprobable)) idmostprobable=i;
  }
  
  std::vector< SimHitIdpr > simhitid;
  associateHitId(*componenthits[idmostprobable], simhitid, simhitCFPos);
  return simhitid;
}

void TrackerHitAssociator::associatePixelRecHit	(	const SiPixelRecHit *	pixelrechit,
		std::vector< SimHitIdpr > &	simhitid,
		std::vector< simhitAddr > *	simhitCFPos = 0
	)		const

Definition at line 537 of file TrackerHitAssociator.cc.

References PixelDigi::channelToPixel(), SiPixelRecHit::cluster(), edm::DetSet< T >::data, spr::find(), TrackingRecHit::geographicalId(), edm::Ref< C, T, F >::isNull(), pixeldigisimlink, DetId::rawId(), and DetId::subdetId().

Referenced by associateHitId().

{
  //
  // Pixel associator
  //
  DetId detid=  pixelrechit->geographicalId();
  uint32_t detID = detid.rawId();

  edm::DetSetVector<PixelDigiSimLink>::const_iterator isearch = pixeldigisimlink->find(detID); 
  if(isearch != pixeldigisimlink->end()) {  //if it is not empty
    edm::DetSet<PixelDigiSimLink> link_detset = (*isearch);
    SiPixelRecHit::ClusterRef const& cluster = pixelrechit->cluster();
    
    //check the reference is valid
    
    if(!(cluster.isNull())){//if the cluster is valid
      
      int minPixelRow = (*cluster).minPixelRow();
      int maxPixelRow = (*cluster).maxPixelRow();
      int minPixelCol = (*cluster).minPixelCol();
      int maxPixelCol = (*cluster).maxPixelCol();    
      //std::cout << "    Cluster minRow " << minPixelRow << " maxRow " << maxPixelRow << std::endl;
      //std::cout << "    Cluster minCol " << minPixelCol << " maxCol " << maxPixelCol << std::endl;
      edm::DetSet<PixelDigiSimLink>::const_iterator linkiter = link_detset.data.begin(), linkEnd = link_detset.data.end();
      int dsl = 0;
      std::vector<SimHitIdpr> idcachev;
      std::vector<simhitAddr> CFposcachev;
      for( ; linkiter != linkEnd; ++linkiter) {
    ++dsl;
    std::pair<int,int> pixel_coord = PixelDigi::channelToPixel(linkiter->channel());
    //std::cout << "    " << dsl << ") Digi link: row " << pixel_coord.first << " col " << pixel_coord.second << std::endl;      
    if(  pixel_coord.first  <= maxPixelRow && 
         pixel_coord.first  >= minPixelRow &&
         pixel_coord.second <= maxPixelCol &&
         pixel_coord.second >= minPixelCol ) {
      //std::cout << "      !-> trackid   " << linkiter->SimTrackId() << endl;
      //std::cout << "          fraction  " << linkiter->fraction()   << endl;
      SimHitIdpr currentId(linkiter->SimTrackId(), linkiter->eventId());
      if(find(idcachev.begin(),idcachev.end(),currentId) == idcachev.end()){
        simtrackid.push_back(currentId);
        idcachev.push_back(currentId);
      }

      if (simhitCFPos != 0) {
      //create a vector that contains all the position (in the MixCollection) of the SimHits that contributed to the RecHit
      //write position only once
        unsigned int currentCFPos = linkiter->CFposition();
        unsigned int tofBin = linkiter->TofBin();
        simHitCollectionID theSimHitCollID = std::make_pair(detid.subdetId(), tofBin);
        simhitAddr currentAddr = std::make_pair(theSimHitCollID, currentCFPos);

        if(find(CFposcachev.begin(), CFposcachev.end(), currentAddr) == CFposcachev.end()) {
          CFposcachev.push_back(currentAddr);
          simhitCFPos->push_back(currentAddr);
        }
      }

    } 
      }
    }
    else{      
      edm::LogError("TrackerHitAssociator")<<"no Pixel cluster reference attached";
      
    }
  }
}

std::vector< SimHitIdpr > TrackerHitAssociator::associateProjectedRecHit	(	const ProjectedSiStripRecHit2D *	projectedrechit,
		std::vector< simhitAddr > *	simhitCFPos = 0
	)		const

Definition at line 525 of file TrackerHitAssociator.cc.

References associateSiStripRecHit(), and ProjectedSiStripRecHit2D::originalHit().

Referenced by associateHitId().

{
  //projectedRecHit is a "matched" rechit with only one component

  std::vector<SimHitIdpr> matched_mono;
 
  const SiStripRecHit2D mono = projectedrechit->originalHit();
  associateSiStripRecHit(&mono, matched_mono, simhitCFPos);
  return matched_mono;
}

void TrackerHitAssociator::associateSimpleRecHitCluster	(	const SiStripCluster *	clust,
		const DetId &	detid,
		std::vector< SimHitIdpr > &	simtrackid,
		std::vector< simhitAddr > *	simhitCFPos = 0
	)		const

Definition at line 412 of file TrackerHitAssociator.cc.

References SiStripCluster::amplitudes(), edm::DetSet< T >::data, spr::find(), plotBeamSpotDB::first, SiStripCluster::firstStrip(), plotBeamSpotDB::last, DetId::rawId(), stripdigisimlink, and DetId::subdetId().

Referenced by associateCluster(), and associateSiStripRecHit().

                                                                      {
  
  uint32_t detID = detid.rawId();
  edm::DetSetVector<StripDigiSimLink>::const_iterator isearch = stripdigisimlink->find(detID); 
  if(isearch != stripdigisimlink->end()) {  //if it is not empty
    edm::DetSet<StripDigiSimLink> link_detset = (*isearch);
    
    if(clust!=0){//the cluster is valid
      int clusiz = clust->amplitudes().size();
      int first  = clust->firstStrip();     
      int last   = first + clusiz;
      
//       std::cout << "CLUSTERSIZE " << clusiz << " first strip = " << first << " last strip = " << last-1 << std::endl;
//       std::cout << " detID = " << detID << " DETSET size = " << link_detset.data.size() << std::endl;
      //use a vector
      std::vector<SimHitIdpr> idcachev;
      std::vector<simhitAddr> CFposcachev;
      int channel;
//      for(auto linkiter : link_detset.data){
      for(edm::DetSet<StripDigiSimLink>::const_iterator linkiter = link_detset.data.begin(), linkerEnd = link_detset.data.end(); linkiter != linkerEnd; ++linkiter){      
        channel = (int)(linkiter->channel()); 
    if( channel >= first  && channel < last ){
      
      //check this digisimlink
//    printf("%s%4d%s%8d%s%3d%s%8.4f\n", "CHANNEL = ", linkiter->channel(), " TrackID = ", linkiter->SimTrackId(),
//       " tofBin = ", linkiter->TofBin(), " fraction = ", linkiter->fraction());
      /*
        std::cout << "CHECKING CHANNEL  = " << linkiter->channel()   << std::endl;
        std::cout << "TrackID  = " << linkiter->SimTrackId()  << std::endl;
        std::cout << "Position = " << linkiter->CFposition()  << std::endl;
        std::cout << " fraction = " << linkiter->fraction() << std::endl;
      */
      
      SimHitIdpr currentId(linkiter->SimTrackId(), linkiter->eventId());
      
      //create a vector with the list of SimTrack ID's of the tracks that contributed to the RecHit
      //write the id only once in the vector
      
      if(find(idcachev.begin(),idcachev.end(),currentId ) == idcachev.end()){
        /*
          std::cout << " Adding track id  = " << currentId.first  
          << " Event id = " << currentId.second.event() 
          << " Bunch Xing = " << currentId.second.bunchCrossing() 
          << std::endl;
        */
        idcachev.push_back(currentId);
        simtrackid.push_back(currentId);
      }
          
      if (simhitCFPos != 0) {
      //create a vector that contains all the position (in the MixCollection) of the SimHits that contributed to the RecHit
      //write position only once
        unsigned int currentCFPos = linkiter->CFposition();
        unsigned int tofBin = linkiter->TofBin();
        simHitCollectionID theSimHitCollID = std::make_pair(detid.subdetId(), tofBin);
        simhitAddr currentAddr = std::make_pair(theSimHitCollID, currentCFPos);

        if(find(CFposcachev.begin(), CFposcachev.end(), currentAddr ) == CFposcachev.end()) {
//            std::cout << "CHECKING CHANNEL  = " << linkiter->channel()   << std::endl;
//        std::cout << "\tTrackID  = " << linkiter->SimTrackId()  << "\tCFPos = " << currentCFPos <<"\ttofBin = " << tofBin << std::endl;
//        simhit_collectionMap::const_iterator it = SimHitCollMap.find(theSimHitCollID);
//        if (it!= SimHitCollMap.end()) {
//      PSimHit theSimHit = it->second[currentCFPos];
//      std::cout << "\tLocal Pos = " << theSimHit.localPosition()
//            << "\tProcess = " << theSimHit.processType() << std::endl;
//        }
          CFposcachev.push_back(currentAddr);
          simhitCFPos->push_back(currentAddr);
        }
      }
    }
      }    
    }
    else {
      edm::LogError("TrackerHitAssociator")<<"no cluster reference attached";
    }
  }
}

template<typename T >

void TrackerHitAssociator::associateSiStripRecHit ( const T *  simplerechit, std::vector< SimHitIdpr > &  simtrackid, std::vector< simhitAddr > *  simhitCFPos = 0  ) const

inline

Definition at line 381 of file TrackerHitAssociator.cc.

References associateSimpleRecHitCluster().

Referenced by associateHitId(), associateMatchedRecHit(), and associateProjectedRecHit().

{
  const SiStripCluster* clust = &(*simplerechit->cluster());
  associateSimpleRecHitCluster(clust, simplerechit->geographicalId(), simtrackid, simhitCFPos);
}

void TrackerHitAssociator::makeMaps ( const edm::Event &  theEvent, const Config &  config  )

private

Definition at line 97 of file TrackerHitAssociator.cc.

References assocHitbySimTrack_, TrackerHitAssociator::Config::cfTokens_, edm::Event::getByToken(), StripDigiSimLink::HighTof, HLT_25ns10e33_v2_cff::labels, edm::Event::labelsForToken(), StripDigiSimLink::LowTof, NULL, edm::Handle< T >::product(), edm::ProductLabels::productInstance, SimHitCollMap, SimHitMap, trackerHits::simHits, and TrackerHitAssociator::Config::simHitTokens_.

Referenced by TrackerHitAssociator().

                                                                                                      {
  // Step A: Get Inputs
  //  The collections are specified via ROUList in the configuration, and can
  //  be either crossing frames (e.g., mix/g4SimHitsTrackerHitsTIBLowTof)
  //  or just PSimHits (e.g., g4SimHits/TrackerHitsTIBLowTof)
  if (assocHitbySimTrack_) {
    for(auto const& cfToken : config.cfTokens_) {
      edm::Handle<CrossingFrame<PSimHit> > cf_simhit;
      //int Nhits = 0;
      if (theEvent.getByToken(cfToken, cf_simhit)) {
        std::unique_ptr<MixCollection<PSimHit> > thisContainerHits(new MixCollection<PSimHit>(cf_simhit.product())); 
        for (auto const& isim : *thisContainerHits) {
          DetId theDet(isim.detUnitId());
          SimHitMap[theDet].push_back(isim);
        //  ++Nhits;
        }
      // std::cout << "simHits from crossing frames; map size = " << SimHitCollMap.size() << ", Hit count = " << Nhits << std::endl;
      }
    }
    for(auto const& simHitToken : config.simHitTokens_) {
      edm::Handle<std::vector<PSimHit> > simHits;
      //int Nhits = 0;
      if(theEvent.getByToken(simHitToken, simHits)) {
        for (auto const& isim : *simHits) {
          DetId theDet(isim.detUnitId());
          SimHitMap[theDet].push_back(isim);
        //++Nhits;
        }
      // std::cout << "simHits from prompt collections; map size = " << SimHitCollMap.size() << ", Hit count = " << Nhits << std::endl;
      }
    }
  } else {
    const char* const highTag = "HighTof";
    unsigned int tofBin;
    edm::EDConsumerBase::Labels labels;
    simHitCollectionID theSimHitCollID;
    for(auto const& cfToken : config.cfTokens_) {
      edm::Handle<CrossingFrame<PSimHit> > cf_simhit;
      //int Nhits = 0;
      if (theEvent.getByToken(cfToken, cf_simhit)) {
        std::unique_ptr<MixCollection<PSimHit> > thisContainerHits(new MixCollection<PSimHit>(cf_simhit.product())); 
        theEvent.labelsForToken(cfToken, labels);
        if(std::strstr(labels.productInstance, highTag) != NULL) {
          tofBin = StripDigiSimLink::HighTof;
        } else {
          tofBin = StripDigiSimLink::LowTof; 
        }    
        for (auto const& isim : *thisContainerHits) {
          DetId theDet(isim.detUnitId());
          theSimHitCollID = std::make_pair(theDet.subdetId(), tofBin);
          SimHitCollMap[theSimHitCollID].push_back(isim);
        //++Nhits;
        }
      // std::cout << "simHits from crossing frames; map size = " << SimHitCollMap.size() << ", Hit count = " << Nhits << std::endl;
      }
    }
    for(auto const& simHitToken : config.simHitTokens_) {
      edm::Handle<std::vector<PSimHit> > simHits;
      //int Nhits = 0;
      if(theEvent.getByToken(simHitToken, simHits)) {
        theEvent.labelsForToken(simHitToken, labels);
        if(std::strstr(labels.productInstance, highTag) != NULL) {
          tofBin = StripDigiSimLink::HighTof;
        } else {
          tofBin = StripDigiSimLink::LowTof; 
        }
        for (auto const& isim : *simHits) {
          DetId theDet(isim.detUnitId());
          theSimHitCollID = std::make_pair(theDet.subdetId(), tofBin);
          SimHitCollMap[theSimHitCollID].push_back(isim);
        //++Nhits;
        }
      // std::cout << "simHits from prompt collections; map size = " << SimHitCollMap.size() << ", Hit count = " << Nhits << std::endl;
      }
    }
  }
}

Member Data Documentation

bool TrackerHitAssociator::assocHitbySimTrack_

private

Definition at line 112 of file TrackerHitAssociator.h.

Referenced by associateHit(), and makeMaps().

bool TrackerHitAssociator::doPixel_

private

Definition at line 112 of file TrackerHitAssociator.h.

Referenced by TrackerHitAssociator().

bool TrackerHitAssociator::doStrip_

private

Definition at line 112 of file TrackerHitAssociator.h.

Referenced by TrackerHitAssociator().

bool TrackerHitAssociator::doTrackAssoc_

private

Definition at line 112 of file TrackerHitAssociator.h.

Referenced by associateHit(), and TrackerHitAssociator().

edm::Handle< edm::DetSetVector<PixelDigiSimLink> > TrackerHitAssociator::pixeldigisimlink

private

Definition at line 111 of file TrackerHitAssociator.h.

Referenced by associatePixelRecHit(), and TrackerHitAssociator().

simhit_collectionMap TrackerHitAssociator::SimHitCollMap

Definition at line 104 of file TrackerHitAssociator.h.

Referenced by associateCluster(), associateHit(), and makeMaps().

simhit_map TrackerHitAssociator::SimHitMap

Definition at line 102 of file TrackerHitAssociator.h.

Referenced by associateHit(), CkfDebugger::hasDelta(), makeMaps(), and CkfDebugger::printSimHits().

edm::Handle< edm::DetSetVector<StripDigiSimLink> > TrackerHitAssociator::stripdigisimlink

private

Definition at line 110 of file TrackerHitAssociator.h.

Referenced by associateSimpleRecHitCluster(), and TrackerHitAssociator().